bzoj4553 [Tjoi2016&Heoi2016]序列

　　首先记录每个点原来的数字V、最小的变化L,最大的变化R，在子序列中两个点i,j(i<j)如果相邻，应满足A[I]<=L[J] && R[I]<=A[J],然后就维护一下这个东西，有点类似二维最长不降子序列，可以用cdq分治做，复杂度O(nlognlogn)。

　　代码

#include<cstdio>
#include<vector>
#include<algorithm>
#define lb(x) (x&-x)
#define pb push_back
using namespace std;
const int N = 301010;
int l[N],r[N],s[N],L[N],R[N],n,m,i,a,b,v[N],ans[N],Ans,c[N];
struct g{
    int l,r,id;
}d[N];
void cc(int x,int w)
{
    while (x<=100000)
    {
        if (w) c[x]=max(c[x],w);else c[x]=0;
        x+=lb(x);
    }
}
int get(int x)
{
    int ans=0;
    while (x)
    {
        ans=max(ans,c[x]);
        x-=lb(x);
    }
    return ans;
}
bool cmp(g a,g b)
{
    if (a.l==b.l)
        return a.id<b.id;
    return a.l<b.l;
}
void solve(int l,int r)
{
    if (l==r)
    {
        ans[l]=max(ans[l],1);
        return;
    }
    int m=(l+r)>>1;
    solve(l,m);
    for (i=l;i<=r;i++)
    {
        if (i<=m)
        {
            d[i].l=v[i];
            d[i].r=R[i];    
        }
        else
        {
            d[i].l=L[i];
            d[i].r=v[i];
        }
        d[i].id=i;
    }
    sort(d+l,d+1+r,cmp);
    for (i=l;i<=r;i++)
    if (d[i].id<=m)
        cc(d[i].r,ans[d[i].id]); 
    else
        ans[d[i].id]=max(ans[d[i].id],get(d[i].r)+1);
    for (i=l;i<=r;i++)
    if (d[i].id<=m)
        cc(d[i].r,0);
    solve(m+1,r);
}
int main()
{
    scanf("%d%d",&n,&m);
    for (i=1;i<=n;i++) 
    {
        scanf("%d",&v[i]);
        R[i]=L[i]=v[i];
    }
    for (i=1;i<=m;i++)
    {
        scanf("%d%d",&a,&b);
        R[a]=max(R[a],b);
        L[a]=min(L[a],b);
    }
    solve(1,n); 
    for (i=1;i<=n;i++)
    if (ans[i]>Ans) Ans=ans[i];
    printf("%d
",Ans);
}