bzoj2253: [2010 Beijing wc]纸箱堆叠

链接

　　http://www.lydsy.com/JudgeOnline/problem.php?id=2253

题解

　　我现在的智商也就….幼儿园？
　　$BIT$还能维护前缀$max$啊涨姿势了，其实只要能搞前缀和的东西好像都能用$BIT$维护吧。
　　这道题目和我以前做的$CDQ$分治不一样，这道题前面的信息对后面是有影响的。我的意思是说，只有$[l,mid]$对$[mid+1,r]$的影响处理完了，才能做$mid+1,r$的。这样的话就是说，处理顺序是先$[l,mid]$再处理区间之间的影响，最后$[mid+1,r]$。
　　离散之后维护一个权值$BIT$，里面维护最大值，像其它三维偏序问题一样用两个指针乱搞，然后就行了..
　　博主知道他写的题解别人都看不懂，但是博主由于太懒不想再改了…

代码

//cdq分治+BIT
#include <cstdio>
#include <algorithm>
#define maxn 200000
#define lowbit(x) (x&-x)
using namespace std;
struct box{int x,y,z,ans;}b[maxn];
int N, tmp[maxn], ndtot, bit[maxn];
inline bool operator<(box a, box b)
{
    if(a.z==b.z)
    {
        if(a.x==b.x)return a.y<b.y;
        return a.x<b.x;
    }
    return a.z<b.z;
}
inline bool cmp(box a, box b){return a.x==b.x?a.y<b.y:a.x<b.x;}
inline void add(int pos, int v)
{for(;pos<=N*3;pos+=lowbit(pos))bit[pos]=max(bit[pos],v);}
inline int max(int pos)
{int ans=0;for(;pos;pos-=lowbit(pos))ans=max(ans,bit[pos]);return ans;}
inline void erase(int pos){for(;pos<=3*N;pos+=lowbit(pos))bit[pos]=0;}
void init()
{
    int i, a, p, t=1;
    scanf("%d%d%d",&a,&p,&N);
    for(i=1;i<=N;i++)
    {
        tmp[3*i-2]=t*a%p, tmp[3*i-1]=tmp[3*i-2]*a%p, tmp[3*i]=tmp[3*i-1]*a%p;
        t=tmp[3*i];
        sort(tmp+3*i-2,tmp+3*i+1);
        b[i].x=tmp[3*i-2], b[i].y=tmp[3*i-1], b[i].z=tmp[3*i];
        b[i].ans=1;
    }
    sort(tmp+1,tmp+3*N+1);
    for(i=1;i<=N;i++)b[i].x=lower_bound(tmp+1,tmp+3*N+1,b[i].x)-tmp;
    for(i=1;i<=N;i++)b[i].y=lower_bound(tmp+1,tmp+3*N+1,b[i].y)-tmp;
    for(i=1;i<=N;i++)b[i].z=lower_bound(tmp+1,tmp+3*N+1,b[i].z)-tmp;
    sort(b+1,b+N+1);
}
void CDQ(int l, int r)
{
    if(l==r)return;
    int mid=(l+r)>>1, i, j, x=1e9;
    for(i=l;i<r;i++)if(b[i].z!=b[i+1].z and abs(i-mid)<abs(x-mid))x=i;
    if(x==1e9)return;
    mid=x;
    CDQ(l,mid), sort(b+l,b+mid+1,cmp), sort(b+mid+1,b+r+1,cmp);
    for(i=l,j=mid+1;j<=r;j++)
    {
        for(;i<=mid and b[i].x<b[j].x;i++)add(b[i].y,b[i].ans);
        b[j].ans=max(b[j].ans,max(b[j].y-1)+1);
    }
    for(j=l;j<i;j++)erase(b[j].y);
    sort(b+mid+1,b+r+1), CDQ(mid+1,r);
}
int main()
{
    int ans=-1, i;
    init();
    CDQ(1,N);
    for(i=1;i<=N;i++)ans=max(ans,b[i].ans);
    printf("%d",ans);
    return 0;
}