bzoj 4553: [Tjoi2016&Heoi2016]序列 （CDQ分治+DP+树状数组）

题目描述

传送门

题目大意：有一个长为n的序列，m个变换，每次变换把x位置的权值变成v.求一个子序列满足在任意变换中都单调不降（任意变换不能同时进行），问子序列的最长长度。

题解

先考虑暴力DP。如果点j 可以用来更新i的答案。那么所有包含j的变换以及j初值的最大值一定小于等于i的初值。所有包含i的变换以及i的初值的最小值一定大于等于j的初值。也就是max[j]<=a[i],a[j]<=min[i],j<=i。
那么这实际上就是对于j（j,max[j],a[j]） i（i,a[i],min[i]）求三维偏序。我们对于位置分治，然后将第二维排序，第三维用树状数组维护即可。

代码

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#define N 300003
using namespace std;
int n,m,cnt,tr[N],b[N],f[N];
struct data{
    int x,y,v,mx,mn,id,pd;
}a[N];
int cmp(data a,data b)
{
    return a.x<b.x||a.x==b.x&&a.id<b.id;
}
int cmp1(data a,data b)
{
    return a.id<b.id;
}
int lowbit(int x)
{
    return x&(-x);
}
void change(int x,int v)
{
    for (int i=x;i<=cnt;i+=lowbit(i))
     tr[i]=max(tr[i],v);
}
int query(int x)
{
    int ans=0;
    if (!x) return ans;
    for (int i=x;i>=1;i-=lowbit(i))
     ans=max(ans,tr[i]);
    return ans;
}
void clear(int x)
{
    for (int i=x;i<=cnt;i+=lowbit(i)) tr[i]=0;
}
void divide(int l,int r)
{
    if (l==r) return;
    int mid=(l+r)/2;
    divide(l,mid);
    for (int i=l;i<=mid;i++) a[i].x=a[i].mx,a[i].y=a[i].v,a[i].pd=1;
    for (int i=mid+1;i<=r;i++) a[i].x=a[i].v,a[i].y=a[i].mn,a[i].pd=0;
    sort(a+l,a+r+1,cmp);
    for (int i=l;i<=r;i++) 
     if (a[i].pd) change(a[i].y,f[a[i].id]);
     else f[a[i].id]=max(f[a[i].id],query(a[i].y)+1);
    for (int i=l;i<=r;i++) 
     if (a[i].pd) clear(a[i].y);
    sort(a+l,a+r+1,cmp1);
    divide(mid+1,r);
}
int main()
{
    freopen("a.in","r",stdin);
//  freopen("my.out","w",stdout);
    scanf("%d%d",&n,&m);
    for (int i=1;i<=n;i++){
     scanf("%d",&a[i].v);
     b[++cnt]=a[i].v;
     a[i].mx=a[i].v; a[i].mn=a[i].v; a[i].id=i;
    }
    for (int i=1;i<=m;i++) {
        int x,v; scanf("%d%d",&x,&v);
        a[x].mx=max(a[x].mx,v);
        a[x].mn=min(a[x].mn,v);
        b[++cnt]=v;
    }
    sort(b+1,b+cnt+1);
    cnt=unique(b+1,b+cnt+1)-b-1;
    for (int i=1;i<=n;i++) {
        a[i].v=lower_bound(b+1,b+cnt+1,a[i].v)-b;
        a[i].mx=lower_bound(b+1,b+cnt+1,a[i].mx)-b;
        a[i].mn=lower_bound(b+1,b+cnt+1,a[i].mn)-b;
    }
    for (int i=1;i<=n;i++) f[i]=1;
    divide(1,n);
    int ans=0;
    for (int i=1;i<=n;i++) ans=max(ans,f[i]);
    printf("%d\n",ans);
}