bzoj2253 [2010 Beijing wc]纸箱堆叠（cdq分治+树状数组）

bzoj2253 [2010 Beijing wc]纸箱堆叠

原题地址：http://www.lydsy.com/JudgeOnline/problem.php?id=2253

题意：
P 工厂是一个生产纸箱的工厂。纸箱生产线在人工输入三个参数 n p a , , 之后，
即可自动化生产三边边长为

(a mod P,a^2 mod p,a^3 mod P)
(a^4 mod p,a^5 mod p,a^6 mod P)
….
(a^(3n-2) mod p,a^(3n-1) mod p,a^(3n) mod p)

的n个纸箱。在运输这些纸箱时，为了节约空间，必须将它们嵌套堆叠起来。
一个纸箱可以嵌套堆叠进另一个纸箱当且仅当它的最短边、次短边和最长边
长度分别严格小于另一个纸箱的最短边、次短边和最长边长度。这里不考虑
任何旋转后在对角线方向的嵌套堆叠。
你的任务是找出这n个纸箱中数量最多的一个子集，使得它们两两之间都可
嵌套堆叠起来。

数据范围
2<=P<=2000000000
1<=a<=p-1
a^k mod p<>0
ap<=2000000000
1<=N<=50000

题解：
三位偏序的最长上升子序列。

先处理左然后中然后右。

本来想写道水题结果调这么久。

注意：
是最短边、次短边和最长边长度分别严格小于另一个纸箱的最短边、次短边和最长边长度，就是说内部不是按a^(3n-2),a^(3n-1),a^(3n)的顺序，而是要按大小排一遍。
要注意一些全局变量的值在递归回来后改变了记得清零。

代码：

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cstring>
#include<map>
using namespace std;
const int N=50005;
map<pair<int,pair<int,int> >,int> mp;
struct node
{
    int x,y,z,val;
    node(){}
    node(int x,int y,int z,int val):x(x),y(y),z(z),val(val){}
}a[N],b[N];
int q[N],cnt=0,c[N];
bool cmpx(const node &A,const node &B) {return A.x<B.x;}
bool cmpz(const node &A,const node &B) {return A.z<B.z;}
int n,base,mod,aa[3*N];
bool cmpy(const int &A,const int &B) {return a[A].y<a[B].y;}
void add(int x,int d) {for(int i=x;i<=n;i=i+(i&(-i))) c[i]=max(c[i],d);}
void res(int x) {for(int i=x;i<=n;i=i+(i&(-i))) c[i]=0;}
int query(int x) {int ret=0; for(int i=x;i;i=i-(i&(-i))) ret=max(ret,c[i]); return ret;}
void solve(int lf,int rg)
{   
    if(lf==rg) return;
    int mid=(lf+rg)>>1; 
    solve(lf,mid); 
    int tmp=lf; cnt=0;
    for(int i=mid+1;i<=rg;i++) q[++cnt]=i;  
    sort(q+1,q+cnt+1,cmpy);
    for(int i=1;i<=cnt;i++)
    {
        while(tmp<=mid&&a[tmp].y<a[q[i]].y) {add(a[tmp].z,a[tmp].val); tmp++;}
        a[q[i]].val=max(a[q[i]].val,query(a[q[i]].z)+1);
    }
    tmp--; for(int i=tmp;i>=lf;i--) res(a[i].z);
    solve(mid+1,rg);
    int p1=lf; int p2=mid+1; cnt=0;
    while(p1<=mid&&p2<=rg)
    {
        if(a[p1].y<a[p2].y) {b[++cnt]=a[p1]; p1++;}
        else {b[++cnt]=a[p2]; p2++;}
    }
    while(p1<=mid) {b[++cnt]=a[p1]; p1++;}
    while(p2<=rg) {b[++cnt]=a[p2]; p2++;}
    for(int i=1;i<=cnt;i++) a[lf+i-1]=b[i];
}
void init(int i)
{
    if(a[i].x>a[i].z) swap(a[i].x,a[i].z);
    if(a[i].x>a[i].y) swap(a[i].x,a[i].y);
    if(a[i].y>a[i].z) swap(a[i].y,a[i].z);
}
int main()
{
    scanf("%d%d%d",&base,&mod,&n);
    aa[0]=1; int cnt=0;
    for(int i=1;i<=3*n;i++)
    {
        aa[i]=(1LL*aa[i-1]*base)%mod;
        if(i%3==0)
        {
            if(!mp[make_pair(aa[i-2],make_pair(aa[i-1],aa[i]))]) a[++cnt]=node(aa[i-2],aa[i-1],aa[i],1);
            mp[make_pair(aa[i-2],make_pair(aa[i-1],aa[i]))]=1;
        }
    } n=cnt;
    for(int i=1;i<=n;i++) init(i);
    sort(a+1,a+n+1,cmpz); for(int i=1;i<=n;i++) a[i].z=i;
    sort(a+1,a+n+1,cmpx);
    solve(1,n);
    int ret=0; for(int i=1;i<=n;i++) ret=max(ret,a[i].val);
    printf("%d\n",ret);
    return 0;
}