[bzoj1176]Mokia[CDQ分治]

啃了一天论文，发现CDQ分治的原理其实很简单，大概就是这样的一类分治：将左右区间按一定规律排序后分开处理，递归到底时直接计算答案，对于一个区间，按照第二关键字split成两个区间，先处理左区间，之后因为整个区间是有序的，就可以根据左区间来推算右区间的答案，最后递归处理右区间即可。拿此题做比方，先把全区间按照x坐标排序，然后自左向右用前一半（按时间排序）的修改来更新后一半的查询，之后将整个区间按照时间分成两个部分，递归处理。归纳起来就是split->left->push->right->merge->return。

 #include <iostream>
 #include <algorithm>
 #include <cstdio>
 #include <cstdlib>
 #include <cstring>
 #include <ctime>
 #include <cmath>

 using namespace std;

 struct node { int    v,x,y,op,d,pos; }a[],temp[];

 int    s,w,cnt;
 int    L[],Ans[];

 bool    operator<(const node temp1,const node temp2)
 {
     if(temp1.x!=temp2.x)return temp1.x<temp2.x;
     if(temp1.y!=temp2.y)return temp1.y<temp2.y;
     return temp1.op<temp2.op;
 }

 void    Add(const int x,const int d) { for(int i=x;i<=w;i+=i&-i)L[i]+=d; return ; }
 int    Query(const int x) { int t=; for(int i=x;i;i-=i&-i) t+=L[i]; return t; }

 void    CDQ(const int l,const int r)
 {
     if(l==r)return ;

     int    i,mid=l+((r-l)>>);

     for(i=l;i<=r;++i)
     {
         if(a[i].v<=mid && a[i].op==)Add(a[i].y,a[i].d);
         if(a[i].v>mid  && a[i].op==)Ans[a[i].pos]+=Query(a[i].y)*a[i].d;
     }
     for(i=l;i<=r;++i)
         if(a[i].v<=mid && a[i].op==)Add(a[i].y,-a[i].d);

     int l1=l,l2=mid+;
     for(i=l;i<=r;++i)
     {
         if(a[i].v<=mid)temp[l1++]=a[i];
         else    temp[l2++]=a[i];
     }
     for(i=l;i<=r;++i)a[i]=temp[i];
     CDQ(l,mid),CDQ(mid+,r);
     return ;
 }

 int main()
 {
     int    i,op,x,y,z,xa,xb,ya,yb;

     scanf("%d%d",&s,&w);
     while()
     {
         scanf("%d",&op);
         if(op==)
         {
             scanf("%d%d%d",&x,&y,&z);
             a[++cnt]=(node){cnt,x,y,,z,};
         }
         else if(op==)
         {
             scanf("%d%d%d%d",&xa,&ya,&xb,&yb);
             a[++cnt]=(node){cnt,xb  ,yb  ,,+,++Ans[]};
             a[++cnt]=(node){cnt,xb  ,ya-,,-,Ans[]};
             a[++cnt]=(node){cnt,xa-,yb  ,,-,Ans[]};
             a[++cnt]=(node){cnt,xa-,ya-,,+,Ans[]};
         }
         else    break;
     }

     sort(a+,a+cnt+);

     CDQ(,cnt);

     for(i=;i<=Ans[];++i)
         printf("%d\n",Ans[i]);

     return ;
 }