bzoj 3173: [Tjoi2013]最长上升子序列【dp+线段树】

我也不知道为什么把题看成以插入点为结尾的最长生生子序列……还WA了好几次
先把这个序列最后的样子求出来，具体就是倒着做，用线段树维护点数，最开始所有点都是1，然后线段树上二分找到当前数的位置，把这个点标为0（相当于对于这之前的序列这个点是不存在的），把每个数的位置记为p[i]
然后用另一颗线段树维护每个位置上的LIS，根据时间序，每次插入数的时候求一下以他结尾的LIS然后放进线段树上对应的p[i]（因为按照数从小到大所以直接查这个数位置之前的即可），然后再取全部点的max即可

#include<iostream>
#include<cstdio>
using namespace std;
const int N=100005;
int n,m,a[N],p[N],f[N];
struct xds
{
	int l,r,s,mx;
}t[N<<2];
int read()
{
	int r=0,f=1;
	char p=getchar();
	while(p>'9'||p<'0')
	{
		if(p=='-')
			f=-1;
		p=getchar();
	}
	while(p>='0'&&p<='9')
	{
		r=r*10+p-48;
		p=getchar();
	}
	return r*f;
}
void jian(int ro,int l,int r)
{
	t[ro].l=l,t[ro].r=r,t[ro].s=1;
	if(l==r)
		return;
	int mid=(l+r)>>1;
	jian(ro<<1,l,mid);
	jian(ro<<1|1,mid+1,r);
	t[ro].s=t[ro<<1].s+t[ro<<1|1].s;
}
int ef(int ro,int k)
{
	if(t[ro].l==t[ro].r)
		return t[ro].l;
	if(t[ro<<1].s>=k)
		return ef(ro<<1,k);
	else
		return ef(ro<<1|1,k-t[ro<<1].s);
}
void gai(int ro,int p)
{
	if(t[ro].l==t[ro].r)
	{
		t[ro].s=0;
		return;
	}
	int mid=(t[ro].l+t[ro].r)>>1;
	if(p<=mid)
		gai(ro<<1,p);
	else
		gai(ro<<1|1,p);
	t[ro].s=t[ro<<1].s+t[ro<<1|1].s;
}
void build(int ro,int l,int r)
{
	t[ro].l=l,t[ro].r=r,t[ro].mx=0;
	if(l==r)
		return;
	int mid=(l+r)>>1;
	build(ro<<1,l,mid);
	build(ro<<1|1,mid+1,r);
}
int ques(int ro,int l,int r)
{
	if(t[ro].l==l&&t[ro].r==r)
		return t[ro].mx;
	int mid=(t[ro].l+t[ro].r)>>1;
	if(r<=mid)
		return ques(ro<<1,l,r);
	else if(l>mid)
		return ques(ro<<1|1,l,r);
	else
		return max(ques(ro<<1,l,mid),ques(ro<<1|1,mid+1,r));
}
void update(int ro,int p,int v)
{
	if(t[ro].l==t[ro].r)
	{
		t[ro].mx=v;
		return;
	}
	int mid=(t[ro].l+t[ro].r)>>1;
	if(p<=mid)
		update(ro<<1,p,v);
	else
		update(ro<<1|1,p,v);
	t[ro].mx=max(t[ro<<1].mx,t[ro<<1|1].mx);
}
int main()
{
	n=read();
	for(int i=1;i<=n;i++)
		a[i]=read()+1;
	jian(1,1,n);
	for(int i=n;i>=1;i--)
	{
		p[i]=ef(1,a[i]);
		gai(1,p[i]);
	}
	build(1,1,n);
	for(int i=1;i<=n;i++)
	{
		int nw=ques(1,1,p[i])+1;
		update(1,p[i],nw);
		printf("%d
",ques(1,1,n));
	}
	return 0;
}