cdq分治入门--BZOJ1176: [Balkan2007]Mokia

对w*w，w<=2000000的矩形，一开始全是0（或一开始全是s），n<=170000个操作，每次操作：矩阵内某点加上一个数，查某一个子矩阵的和，保证修改数<=160000，询问数<=10000。

这还是一个比较明显的三维偏序：时间维，以及x和y。由于现在时间维是一个Ti<Tj，而x和y是要查x1<=x<=x2，y1<=y<=y2，查一个范围答案在归并排序直接查不方便，所以一个询问拆4个，就变成普通的三维偏序了。

 #include<stdio.h>
 #include<string.h>
 #include<stdlib.h>
 #include<algorithm>
 //#include<iostream>
 using namespace std;

 int n,m,s;
 #define maxn 300011
 #define maxm 2000011
 struct Point
 {
     int x,y,v;bool type;
     //type=1:Q type=0:A
 }a[maxn];
 int ans[maxn];

 struct BIT
 {
     int a[maxm],n;
     void clear(int m) {n=m;memset(a,,sizeof(a));}
     void add(int x,int v) {for (;x<=n;x+=x&-x) a[x]+=v;}
     int query(int x) {int ans=;for (;x;x-=x&-x) ans+=a[x];return ans;}
 }t;

 int ord[maxn],tmpord[maxn];
 void solve(int L,int R)
 {
     if (L==R) {ord[L]=L;return;}
     const int mid=(L+R)>>;
     solve(L,mid);
     solve(mid+,R);
     int i=L,j=mid+,k=L;
     while (i<=mid && j<=R)
     {
         if (a[ord[i]].x<=a[ord[j]].x)
         {
             if (a[ord[i]].type==)
                 t.add(a[ord[i]].y,a[ord[i]].v);
             tmpord[k++]=ord[i++];
         }
         else
         {
             if (a[ord[j]].type)
                 ans[ord[j]]+=t.query(a[ord[j]].y);
             tmpord[k++]=ord[j++];
         }
     }
     for (;j<=R;j++)
     {
         if (a[ord[j]].type) ans[ord[j]]+=t.query(a[ord[j]].y);
         tmpord[k++]=ord[j];
     }
     for (int ii=L;ii<i;ii++) if (!a[ord[ii]].type) t.add(a[ord[ii]].y,-a[ord[ii]].v);
     for (;i<=mid;i++) tmpord[k++]=ord[i];
     for (int x=L;x<=R;x++) ord[x]=tmpord[x];
 }

 int main()
 {
     scanf("%d%d",&s,&m);
     t.clear(m);
     n=;
     int op,x,y,z,w;
     while (scanf("%d",&op) && op!=)
     {
         if (op-)
         {
             scanf("%d%d%d%d",&x,&y,&z,&w);
             a[++n]=(Point){x-,y-,,};
             a[++n]=(Point){x-,w,-,};
             a[++n]=(Point){z,y-,-,};
             a[++n]=(Point){z,w,,};
         }
         else
         {
             scanf("%d%d%d",&x,&y,&z);
             a[++n]=(Point){x,y,z,};
         }
     }
     solve(,n);

     for (int i=;i<=n;i++) if (a[i].type)
     {
         int Ans=;
         for (int to=i+;i<to;i++) Ans+=ans[i]*a[i].v;
         i--;
         printf("%d\n",Ans);
     }
     return ;
 }

这题面的s没有用，怎么题面也不改。。