BZOJ3173: [Tjoi2013]最长上升子序列

1~n<=1e5依次插入序列中某一个位置，求每次插入后的最长上升子序列。

因为新插入的数对前面插入的答案没影响，所以只要能想方设法构造出最终序列即可。

方法一：平衡树！。。。。

方法二：在树状数组上倍增求第K大

然而LIS写残了。。记得是lowerbound不是upper。。

然后最后输出记得和前面小的答案取max因为新插入后最长上升子序列并不一定涉及新插入的数。。

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cstdlib>
#include<math.h>
//#include<queue>
//#include<iostream>
using namespace std;

int n;
#define maxn 100011
int a[maxn];
struct BIT
{
    int a[maxn];
    BIT() {memset(a,0,sizeof(a));}
    void add(int x,int v) {for (;x<=n;x+=x&-x) a[x]+=v;}
    int find(int x)
    {
        int now=0,ans=0;
        for (int i=19;i>=0;i--)
        {
            now+=(1<<i);
            if (now>=n || ans+a[now]>=x) now-=(1<<i);
            else ans+=a[now];
        }
        return now+1;
    }
}t;
int num[maxn],g[maxn],last[maxn];
int main()
{
    scanf("%d",&n);
    for (int i=1;i<=n;i++) scanf("%d",&a[i]),a[i]++,t.a[i]++,t.a[i+(i&-i)]+=t.a[i];
    for (int i=n;i;i--)
    {
        int tmp=t.find(a[i]);
        num[tmp]=i;
        t.add(tmp,-1);
    }
    int ans=0;
    for (int i=1;i<=n;i++)
    {
        int now=lower_bound(g,g+ans+1,num[i])-g;
        if (now<=ans) g[now]=min(g[now],num[i]);
        else ans++,g[ans]=num[i];
        last[num[i]]=now;
    }
    for (int i=1;i<=n;i++) last[i]=max(last[i-1],last[i]),printf("%d
",last[i]);
    return 0;
}