HDU 3564 （树状数组+LIS）

题意：
一个序列一开始为空，将1到n依次按照给定的位置插入序列，求每次插入后的LIS。

思路：
分成两个问题，
1、按照给定位置插入元素
2、插入之后求LIS
第一个问题，插入元素可以用插空的方法，倒序插入位置，用二分+树状数组的方式，根据当前位置前面的空位数，二分判断，确认插入的位置在最终序列中对应的位置。
第二个问题，因为插入的元素顺序是由小到大的，当一个元素插入时，插入元素即是当前序列中最大的元素，所以后面插入的元素不影响当前的LIS。
dp时，每次使用当前元素时当前元素为最后一个元素的子序列即为当前的最长上升子序列。

代码：

#include <iostream>
#include <iomanip>
#include <algorithm>
#include <cstring>
#include <cctype>
#include <cstdlib>
#include <cstdio>
#include <cmath>
#include <ctime>
#include <map>
#include <list>
#include <set>
#include <stack>
#include <queue>
#include <string>
#include <sstream>
#define pb push_back
#define X first
#define Y second
#define ALL(x) x.begin(),x.end()
#define INS(x) inserter(x,x.begin())
#define pii pair<int,int>
#define qclear(a) while(!a.empty())a.pop();
#define lowbit(x) (x&-x)
#define sd(n) scanf("%d",&n)
#define sdd(n,m) scanf("%d%d",&n,&m)
#define sddd(n,m,k) scanf("%d%d%d",&n,&m,&k)
#define mst(a,b) memset(a,b,sizeof(a))
#define cout3(x,y,z) cout<<x<<" "<<y<<" "<<z<<endl
#define cout2(x,y) cout<<x<<" "<<y<<endl
#define cout1(x) cout<<x<<endl
#define IOS std::ios::sync_with_stdio(false)
#define SRAND srand((unsigned int)(time(0)))
typedef long long ll;
typedef unsigned long long ull;
typedef unsigned int uint;
using namespace std;
const double PI=acos(-1.0);
const int INF=0x3f3f3f3f;
const ll INFF=0x3f3f3f3f3f3f3f3f;
const ll mod=998244353;
const double eps=1e-5;
const int maxn=100005;
const int maxm=20005;

int n;
int bt[maxn<<1];
int sum(int x){
    int ret=0;
    while(x>0){
        ret+=bt[x];
        x-=lowbit(x);
    }
    return ret;
}
void add(int x,int d){
    while(x<=n){
        bt[x]+=d;
        x+=lowbit(x);
    }
}
int k[maxn];
int a[maxn];
int dp[maxn];
int ans[maxn];
void init(){
    for(int i=1;i<=n;i++){
        add(i,1);
    }
    mst(dp,0);
}
void solve() {
    int t;
    sd(t);
    for(int kase=1;kase<=t;kase++){
        sd(n);
        init();
        for(int i=1;i<=n;i++){
            sd(k[i]);k[i]++;
        }
        for(int i=n;i>=1;i--){
            int now=k[i];
            int l,r;
            l=1;r=n;
            while(l!=r){
                int mid=(l+r)>>1;
                if(sum(mid)<now)l=mid+1;
                else r=mid;
            }
            a[l]=i;
            add(l,-1);
        }
        int len=0;
        printf("Case #%d:\n",kase);
        for(int i=1;i<=n;i++){
            if(a[i]>dp[len]){
                dp[++len]=a[i];
                ans[a[i]]=len;
            }else{
                int pos=lower_bound(dp,dp+len,a[i])-dp;
                dp[pos]=a[i];
                ans[a[i]]=pos;
            }
        }
        for(int i=1;i<=n;i++){
            ans[i]=max(ans[i],ans[i-1]);
            printf("%d\n",ans[i]);
        }
        printf("\n");
    }
    return ;
}
int main() {
#ifdef LOCAL
    freopen("in.txt","r",stdin);
//    freopen("out.txt","w",stdout);
#else
    //    freopen("","r",stdin);
    //    freopen("","w",stdout);
#endif
    solve();
    return 0;
}