bzoj3173最长上升子序列

Posted 2020-11-18 yoyoball

Portal --> bzoj3173

Solution

　　感觉自己需要智力康复qwq

　　首先题目给的这个序列肯定是一个(1-n)的排列，并且插入的顺序是从小到大

　　仔细思考一下会发现如果知道了最终的序列，问题就比较好解决了，这里提供一种用线段树的做法：

　　如果知道了最终的序列，记数字(i)在该序列中的位置为(loc[i])，那么我们按照(i)从小到大的顺序，查询结尾在([1,loc[i]))的这段位置中的最长上升子序列的最大值(mx)，并将(mx+1)作为以(loc[i])位置为结尾的答案，插入到线段树中(loc[i])对应的节点里，复杂度是(O(nlogn))

　　然后现在的问题是怎么求最终的序列

　　这个可以用平衡树来写，不过其实也可以用线段树来写

　　考虑反过来确定每一个数在最终序列中的位置，因为是反过来考虑的，所以一开始的时候每一个位置都有一个数，然后我们根据读入的插入位置，按照(n-1)的顺序，找到当前这个数的位置，然后将它删掉（也就是对应的线段树节点的(sum-1)）

?　　具体一点就是比如当前考虑到第(i)个数，读入这个数应该要插入在(a[i])的位置后面，也就是应该在当前这个序列的第(a[i]+1)个位置，那么找到这个位置，然后把这个位置删掉，这样就可以得到还没有插入这个数之前的序列的位置集合了，这部分的复杂度也是(O(nlogn))的

?　　然后就十分愉快地做完啦

　　

　　代码大概长这个样子

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
const int MAXN=100010,SEG=MAXN*4;
namespace Seg{/*{{{*/
    int ch[SEG][2],sum[SEG],mx[SEG];
    int n,tot;
    void pushup(int x){
        sum[x]=sum[ch[x][0]]+sum[ch[x][1]];
        mx[x]=max(mx[ch[x][0]],mx[ch[x][1]]);
    }
    void _build(int x,int l,int r){
        sum[x]=0; mx[x]=0;
        if (l==r){sum[x]=0; return;}
        int mid=l+r>>1;
        ch[x][0]=++tot; _build(ch[x][0],l,mid);
        ch[x][1]=++tot; _build(ch[x][1],mid+1,r);
        pushup(x);
    }
    void build(int _n){n=_n; tot=1; _build(1,1,n);}
    void _update(int x,int d,int lx,int rx,int delta){
        if (lx==rx) {sum[x]+=delta;mx[x]+=delta;return;}
        int mid=lx+rx>>1;
        if (d<=mid) _update(ch[x][0],d,lx,mid,delta);
        else _update(ch[x][1],d,mid+1,rx,delta);
        pushup(x);
    }
    void update(int d,int delta){_update(1,d,1,n,delta);}
    int _query_mx(int x,int l,int r,int lx,int rx){
        if (l<=lx&&rx<=r) return mx[x];
        int mid=lx+rx>>1,ret=0;
        if (l<=mid) ret=max(ret,_query_mx(ch[x][0],l,r,lx,mid));
        if (r>mid) ret=max(ret,_query_mx(ch[x][1],l,r,mid+1,rx));
        return ret;
    }
    int query(int l,int r){return _query_mx(1,l,r,1,n);}
    int _get_loc(int x,int lx,int rx,int k){
        if (lx==rx) return lx;
        int mid=lx+rx>>1;
        if (sum[ch[x][0]]>=k) return _get_loc(ch[x][0],lx,mid,k);
        else return _get_loc(ch[x][1],mid+1,rx,k-sum[ch[x][0]]);
    }
    int get_loc(int k){return _get_loc(1,1,n,k);}
};/*}}}*/
int loc[MAXN],a[MAXN],b[MAXN];
int n,m,ans;

int main(){
#ifndef ONLINE_JUDGE
    freopen("a.in","r",stdin);
#endif
    scanf("%d",&n);
    for (int i=1;i<=n;++i) scanf("%d",b+i),++b[i];
    Seg::build(n);
    for (int i=1;i<=n;++i) Seg::update(i,1);
    for (int i=n;i>=1;--i){
        loc[i]=Seg::get_loc(b[i]);
        Seg::update(loc[i],-1);
    }
    for (int i=1;i<=n;++i) a[loc[i]]=i;

    Seg::build(n);
    ans=0;  
    int tmp;
    for (int i=1;i<=n;++i){
        if (loc[i]>1)
            tmp=Seg::query(1,loc[i]-1);
        else
            tmp=0;
        Seg::update(loc[i],tmp+1);
        ans=max(ans,tmp+1);
        printf("%d
",ans);
    }
}