CDQ分治——BZOJ4553 [Tjoi2016&Heoi2016]序列

传送门
和带三维的关键字最长不降挺像
所以我们可以先搞出所有数的原值，最小值和最大值
然后可以CDQ求解
我们首先l~mid按照最大值排序，mid+1~r按照原值排序（每个排列最多只能改动一个数）
然后计算前面的对后面答案的影响即可
跟LIS差不多，只不过最后一维我们用一个树状数组来维护
最后清空一下（和二维矩阵差不多）
对就是这样，答案记录一下最大值就好了
然后关于树状数组和排序的时候写得很丑。。。见谅

#include<bits/stdc++.h>using namespace std;struct ppap{int x,y,v,w;}a[200001],c[200001];int f[200001],g[200001],n,m,ans=0,maxx;inline bool cmp(ppap a,ppap b){return a.y<b.y;}inline bool cmpp(ppap a,ppap b){return a.v<b.v;}inline int lowbit(int x){return x&-x;}inline void ad(int x){for(;x<=maxx;x+=lowbit(x))f[x]=0;}inline void add(int x,int v){for(;x<=maxx;x+=lowbit(x))f[x]=max(f[x],v);}inline int sum(int x){int ans=0;for(;x;x-=lowbit(x))ans=max(ans,f[x]);return ans;}void cdq(int l,int r){    if(l==r){        g[l]=max(g[l],1);ans=max(ans,g[l]);return;    }    int i,cc=l,mid=l+r>>1;    cdq(l,mid);    for(i=l;i<=r;++i)c[i]=a[i];    sort(c+l,c+mid+1,cmp);sort(c+mid+1,c+r+1,cmpp);    for(i=mid+1;i<=r;++i){        for(;c[cc].y<=c[i].v&&cc<=mid;++cc)add(c[cc].v,g[c[cc].w]);        g[c[i].w]=max(g[c[i].w],sum(c[i].x)+1);    }    for(i=l;i<cc;++i)ad(c[i].v);    cdq(mid+1,r);}int main(){    int i,x,y;scanf("%d%d",&n,&m);    for(i=1;i<=n;++i)a[i].w=i,scanf("%d",&a[i].v),a[i].x=a[i].y=a[i].v,maxx=max(maxx,a[i].v);    for(i=1;i<=m;++i){        scanf("%d%d",&x,&y);        a[x].x=min(a[x].x,y);        a[x].y=max(a[x].y,y);        maxx=max(maxx,a[x].y);    }    cdq(1,n);    printf("%d",ans);    return 0;}