BZOJ2253 2010 Beijing wc 纸箱堆叠 CDQ分治

这题之前度娘上没有CDQ分治做法，gerwYY出来以后写了一个。不过要sort3遍，常数很大。

gerw说可以类似划分树的思想优化复杂度，但是蒟蒻目前不会划分树（会了主席树就懒得去弄了）。

嗯 将memset改成手动clear会快很多。

还有就是第一维相同的情况，划分为两个不存在第一维相同的两个区间即可。

#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<cstring>
#include<iostream>
using namespace std;
const int Maxn= + ;
int C[Maxn*],ans=,A,P,n=,h[Maxn],tot=;
struct Node{
    int a,b,c,ans;
    void init(int x,int y,int z){
        if(x>y)swap(x,y);
        if(x>z)swap(x,z);
        if(y>z)swap(y,z);
        a=x,b=y,h[++tot]=c=z;
        ans=;
    }
}p[Maxn],q[Maxn];
inline void Add(int x,const int&y){
    for(;<x&&x<=n;x+=x&-x)C[x]=max(C[x],y);
}
inline void HashIt(int&x){
    x=lower_bound(h+,h+n+,x)-h;
}
inline int Query(int x){
    int ret=;
    for(;<x&&x<=n;x-=x&-x)ret=max(ret,C[x]);
    return ret;
}
inline void Clear(int x) {
    for(;x<=n;x+=x&-x)C[x]=;
}
inline bool cmpa(const Node&x,const Node&y){
    if(x.a!=y.a)return x.a<y.a;
    if(x.b!=y.b)return x.b<y.b;
    return x.c<y.c;
}
inline bool cmpb(const Node&x,const Node&y){
    if(x.b!=y.b)return x.b<y.b;
    return x.c<y.c;
}
void init(){
    scanf("%d%d%d",&A,&P,&n);
    for(int a,b,c,t=,i=;i<=n;i++){
        a=(t=(long long)A*t%P);
        b=(t=(long long)A*t%P);
        c=(t=(long long)A*t%P);
        p[i].init(a,b,c);
    }
    sort(p+,p+n+,cmpa);
    sort(h+,h+n+);
    for(int i=;i<=n;i++)HashIt(p[i].c);
}

void CDQ(int l,int r){
    if(l==r)return;
    int mid=-,L=(l+r)>>,R=L+;
    while(l<=L || R<=r){
        if(l<=L && p[L].a!=p[L+].a){mid=L;break;}
        if(R<=r && p[R].a!=p[R-].a){mid=R-;break;}
        L--,R++;
    }
    if(mid==-)return;
    CDQ(l,mid);
    sort(p+l,p+mid+,cmpb);
    sort(p+mid+,p+r+,cmpb);
    int i=l;
    for(int j=mid+;j<=r;j++){
        for(;i<=mid && p[i].b<p[j].b;i++)Add(p[i].c,p[i].ans);
        p[j].ans=max(p[j].ans,Query(p[j].c-)+);
    }
    for(int j=l;j<=i;j++)Clear(p[j].c);
    sort(p+mid+,p+r+,cmpa);
    CDQ(mid+,r);
}
int main(){
    init();
    CDQ(,n);
    for(int i=;i<=n;i++)
        if(p[i].ans>ans)ans=p[i].ans;
    printf("%d\n",ans);
    return ;
}