bzoj 4553: [Tjoi2016&Heoi2016]序列 cdq分治+树状数组

题意

给出一个长度为n的序列，有m个时刻，每个时刻序列的某一个数都会改变，然后马上又变回来。求一个最长的序列，使得其在任意时刻都不下降。
n,m<=100000

分析

我们设原序列为v[i]，mn[i]表示位置i在所有时刻中的最小值，mx[i]表示最大值。
那么若一个j可以放在i的后面，那么必须要满足mx[j]<=v[i]且v[j]<=mn[i]
那么就变成了一个二维数点问题。
考虑cdq分治
二分第一维，设其为mid，然后把第一维不大于mid的点都放到左边，然后递归处理。
然后用两个指针+树状数组统计左边的点对右边的点的贡献，然后递归右边即可。

代码

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#define N 200005
using namespace std;

int n,v[N],mn[N],c[N],f[N],a[N],a1,m,mx[N],tot,x[N],y[N];
struct query{int x,y,id,tim,op;}q[N*2],tmp[N*2];

void ins(int x,int y)
{
    while (x<=a1)
    {
        c[x]=max(c[x],y);
        x+=x&(-x);
    }
}

int find(int x)
{
    int ans=0;
    while (x)
    {
        ans=max(ans,c[x]);
        x-=x&(-x);
    }
    return ans;
}

void reset(int x)
{
    while (x<=a1)
    {
        c[x]=0;
        x+=x&(-x);
    }
}

bool cmptim(query a,query b)
{
    return a.tim<b.tim;
}

void cdq(int l,int r,int L,int R)
{
    if (l>=r) return;
    if (L==R)
    {
        for (int i=l;i<=r;i++)
            if (q[i].op==1) ins(q[i].y,f[q[i].id]);
            else f[q[i].id]=max(f[q[i].id],find(q[i].y)+1);
        for (int i=l;i<=r;i++)
            if (q[i].op==1) reset(q[i].y);
        return;
    }
    int s=0,mid=(L+R)/2;
    for (int i=l;i<=r;i++) if (q[i].x<=mid) s++;
    int u=l,v=l+s;
    for (int i=l;i<=r;i++)
        if (q[i].x<=mid) tmp[u++]=q[i];
        else tmp[v++]=q[i];
    for (int i=l;i<=r;i++) q[i]=tmp[i];
    cdq(l,u-1,L,mid);
    sort(q+l,q+u,cmptim);
    int p1=l,p2=l+s;
    while (p1<u&&p2<v)
        if (q[p1].tim<q[p2].tim)
        {
            if (q[p1].op==1) ins(q[p1].y,f[q[p1].id]);
            p1++;
        }
        else
        {
            if (q[p2].op==2) f[q[p2].id]=max(f[q[p2].id],find(q[p2].y)+1);
            p2++;
        }
    for (int i=p2;i<v;i++)
        if (q[i].op==2) f[q[i].id]=max(f[q[i].id],find(q[i].y)+1);
    for (int i=l;i<u;i++) if (q[i].op==1) reset(q[i].y);
    cdq(u,r,mid+1,R);
}

int main()
{
    scanf("%d%d",&n,&m);
    for (int i=1;i<=n;i++) scanf("%d",&v[i]),a[++a1]=v[i];
    for (int i=1;i<=m;i++) scanf("%d%d",&x[i],&y[i]),a[++a1]=y[i];
    sort(a+1,a+a1+1);
    a1=unique(a+1,a+a1+1)-a-1;
    for (int i=1;i<=n;i++) mn[i]=mx[i]=v[i]=lower_bound(a+1,a+a1+1,v[i])-a;
    for (int i=1;i<=m;i++)
    {
        y[i]=lower_bound(a+1,a+a1+1,y[i])-a;
        mn[x[i]]=min(mn[x[i]],y[i]);mx[x[i]]=max(mx[x[i]],y[i]);
    }
    for (int i=1;i<=n;i++)
    {
        q[++tot].op=2;q[tot].x=v[i];q[tot].y=mn[i];q[tot].id=i;q[tot].tim=tot;
        q[++tot].op=1;q[tot].x=mx[i];q[tot].y=v[i];q[tot].id=i;q[tot].tim=tot;
        f[i]=1;
    }
    cdq(1,tot,1,a1);
    int ans=0;
    for (int i=1;i<=n;i++) ans=max(ans,f[i]);
    printf("%d",ans);
    return 0;
}