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\n",ques(1,1,n));
}
return 0;
}