CDQ分治——BZOJ1176 [Balkan2007]Mokia

传送门
纪念我第一个CDQ分治代码（当然调了很久）还有这个很像某手机品牌的题目名称
二位数点但是差值太大树套树MLE+TLE
但是离线很重要，这个时候CDQ分治派上用场了（%%%陈丹琦dalao发明了这个东西）
我们先把所有操作按照x坐标先y坐标后排序（双关键字保平安）
然后区间l,r取mid，计算l~mid的修改对mid+1~r的查询的影响
我的做法是先把编号小于等于mid的换到左边去，然后直接找右边的查询计算左边贡献，最后向下递归
求贡献的话我们用一个树状数组维护一下合法的值（最后别忘记清空，几乎很多CDQ的题都要套上一两个数据结构吧。。。）
这是CDQ的主要思想，这题求区间矩阵和我们可以用矩阵前缀和的方式来实现
也就是说询问操作一分四然后分别操作即可
其实这个向下递归有点奇怪，我把它放前面就会出问题（照理来说应该不会有问题把。。。）希望大神解答>_<
（其实本蒟蒻对于CDQ的理解不深慢慢长进）
贴一下我拙劣的代码

#include<bits/stdc++.h>
using namespace std;
struct ppap{
    int b,j,x,y,v,k;
}c[1000001],d[1000001];
int f[2000001],ans[1000001],S,n,cnt,q;
inline bool cmp(ppap a,ppap b){return a.x==b.x?a.y<b.y:a.x<b.x;}
inline int lowbit(int x){return x&-x;}
inline void add(int x,int v){for(int i=x;i<=n;i+=lowbit(i))f[i]+=v;}
inline int sum(int x){int ans=0;for(int i=x;i;i-=lowbit(i))ans+=f[i];return ans;}
void cdq(int l,int r){
    if(l==r)return;
    int mid=l+r>>1,dd1=l-1,dd2=mid;
    for(int i=l;i<=r;i++)if(c[i].b<=mid)d[++dd1]=c[i];
    else d[++dd2]=c[i];
    for(int i=l;i<=r;i++)c[i]=d[i];
    int cc=l;
    for(int i=mid+1;i<=r;i++)if(c[i].j==2){
        for(;cc<=mid&&c[cc].x<=c[i].x;cc++)if(c[cc].j==1)add(c[cc].y,c[cc].v);
        ans[c[i].k]+=sum(c[i].y)*c[i].v;
    }
    for(int i=l;i<cc;i++)if(c[i].j==1)add(c[i].y,-c[i].v);
    cdq(l,mid);cdq(mid+1,r);
}
int main()
{
    scanf("%d%d",&S,&n);int k;
    while(1){
        scanf("%d",&k);if(k==3)break;
        if(k==1){
            int x,y,z;scanf("%d%d%d",&x,&y,&z);
            c[++cnt].b=cnt;c[cnt].j=1;c[cnt].x=x;c[cnt].y=y;c[cnt].v=z;c[cnt].k=0;
        }else{
            q++;int x,y,xx,yy;scanf("%d%d%d%d",&x,&y,&xx,&yy);
        //  ans[q]=S*(xx-x)*(yy-y);  
        //还有这题的S并没有用。。。上面的预处理也就不用啦
            c[++cnt].b=cnt;c[cnt].j=2;c[cnt].x=xx;c[cnt].y=yy;c[cnt].v=1;c[cnt].k=q;
            c[++cnt].b=cnt;c[cnt].j=2;c[cnt].x=x-1;c[cnt].y=y-1;c[cnt].v=1;c[cnt].k=q;
            c[++cnt].b=cnt;c[cnt].j=2;c[cnt].x=x-1;c[cnt].y=yy;c[cnt].v=-1;c[cnt].k=q;
            c[++cnt].b=cnt;c[cnt].j=2;c[cnt].x=xx;c[cnt].y=y-1;c[cnt].v=-1;c[cnt].k=q;
        }
    }
    sort(c+1,c+cnt+1,cmp);
    cdq(1,cnt);
    for(int i=1;i<=q;i++)printf("%d\n",ans[i]);
    return 0;
}

呵呵，我是不会告诉你这题是双倍经验的。。。QAQ
双倍经验传送门（BZOJ2683非权限号的就别想啦。。。）