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;
}