BZOJ3173 TJOI2013最长上升子序列（splay）

　　容易发现如果求出最后的序列，只要算一下LIS就好了。序列用平衡树随便搞一下，这里种一棵splay。

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstdlib>
#include<cstring>
#include<algorithm>
using namespace std;
int read()
{
    int x=,f=;char c=getchar();
    while (c<''||c>'') {if (c=='-') f=-;c=getchar();}
    while (c>=''&&c<='') x=(x<<)+(x<<)+(c^),c=getchar();
    return x*f;
}
#define N 100010
int n,tree[N],cnt=;
struct data
{
    int x,f;
    bool operator <(const data&a) const
    {
        return x<a.x;
    }
}a[N];
namespace Splay
{
    struct data{int fa,ch[],s;}tree[N];
    int root=;
    int whichson(int k){return tree[tree[k].fa].ch[]==k;}
    void up(int k){tree[k].s=tree[tree[k].ch[]].s+tree[tree[k].ch[]].s+;}
    void move(int k)
    {
        int fa=tree[k].fa,gf=tree[fa].fa,p=whichson(k);
        if (fa) tree[gf].ch[whichson(fa)]=k;tree[k].fa=gf;
        tree[fa].ch[p]=tree[k].ch[!p],tree[tree[k].ch[!p]].fa=fa;
        tree[k].ch[!p]=fa,tree[fa].fa=k;
        up(fa),up(k);
    }
    void splay(int k)
    {
        while (tree[k].fa)
        {
            int fa=tree[k].fa;
            if (tree[fa].fa)
                if (whichson(fa)^whichson(k)) move(k);
                else move(fa);
            move(k);
        }
        root=k;
    }
    void ins(int &k,int s,int x,int from)
    {
        if (!k) {k=x,tree[k].s=,tree[k].fa=from;return;}
        if (tree[tree[k].ch[]].s>=s) ins(tree[k].ch[],s,x,k);
        else ins(tree[k].ch[],s-tree[tree[k].ch[]].s-,x,k);
    }
    void dfs(int k)
    {
        if (tree[k].ch[]) dfs(tree[k].ch[]);
        a[++cnt].x=k;
        if (tree[k].ch[]) dfs(tree[k].ch[]);
    }
}
int query(int k){int s=;while (k) s=max(s,tree[k]),k-=k&-k;return s;}
void ins(int k,int x){while (k<=n) tree[k]=max(tree[k],x),k+=k&-k;}
int main()
{
#ifndef ONLINE_JUDGE
    freopen("bzoj3173.in","r",stdin);
    freopen("bzoj3173.out","w",stdout);
    const char LL[]="%I64d\n";
#else
    const char LL[]="%lld\n";
#endif
    n=read();
    for (int i=;i<=n;i++)
    {
        int x=read();
        Splay::ins(Splay::root,x,i,);Splay::splay(i);
    }
    Splay::dfs(Splay::root);
    for (int i=;i<=n;i++)
    a[i].f=query(a[i].x)+,ins(a[i].x,a[i].f);
    sort(a+,a+n+);
    for (int i=;i<=n;i++) a[i].f=max(a[i].f,a[i-].f),printf("%d\n",a[i].f);
    return ;
}