bzoj 1176 [Balkan2007]Mokia 【CDQ分治】

W过大，很难在线维护，考虑离线算法

给每个操作加一个时间属性n，显然，对于n=i的询问，对它有影响的修改只在n<i中，所以可以CDQ（因为是按时间序读进来的，所以不用排序了

对于统计矩形和，可以使用二维前缀和的思想，即只需要统计四个点即可

这样就转化成了三维偏序问题，只是询问和修改要分开处理。

（初始值不用管，输出的时候加上就行

 #include<iostream>
 #include<cstdio>
 #include<algorithm>
 using namespace std;
 const int N=;
 int n,s,qu,tot,ans[N],t[N],p[N],cnt,q[N];
 struct qwe
 {
     int o,x,y,v,n;
     qwe(){}
     qwe(int a1,int a2,int a3,int a4,int a5)
     {
         o=a1,x=a2,y=a3,v=a4,n=a5;
     }
 }a[N];
 bool cmp(const qwe &a,const qwe &b)
 {
     return (a.x==b.x&&a.y==b.y&&a.o<b.o)||(a.x==b.x&&a.y<b.y)||(a.x<b.x);
 }
 int read()
 {
     int r=,f=;
     char p=getchar();
     while(p>''||p<'')
     {
         if(p=='-')
             f=-;
         p=getchar();
     }
     while(p>=''&&p<='')
     {
         r=r*+p-;
         p=getchar();
     }
     return r*f;
 }
 int lb(int x)
 {
     return x&(-x);
 }
 void update(int x,int v)
 {
     for(int i=x;i<=n;i+=lb(i))
         t[i]+=v;
 }
 int ques(int x)
 {
     int r=;
     for(int i=x;i>=;i-=lb(i))
         r+=t[i];
     return r;
 }
 void cdq(int l,int r)
 {
     if(l==r)
         return;
     int mid=(l+r)/,tmp=;
     cdq(l,mid);
     cdq(mid+,r);
     sort(a+l,a+mid+,cmp);
     sort(a+mid+,a+r+,cmp);
     int i=l,j=mid+;
     while(j<=r)
     {
         for(;a[i].o==&&i<=mid;i++);
         for(;a[j].o==&&j<=r;j++);
         if(i<=mid&&a[i].x<=a[j].x)
             update(a[i].y,a[i].v),i++,tmp=i-;
         else if(j<=r)
             ans[a[j].n]+=ques(a[j].y),j++;
     }
     for(int t=l;t<=tmp;t++)
         if(a[t].o==)
             update(a[t].y,-a[t].v);
 }
 int main()
 {
     s=read();n=read();
     while()
     {
         qu=read();
         if(qu==)
             break;
         if(qu==)
         {
             int x=read(),y=read(),z=read();
             a[++tot]=qwe(,x,y,z,tot);
         }
         else
         {
             int x1=read(),y1=read(),x2=read(),y2=read();
             p[++cnt]=tot;
             q[cnt]=(x2-x1+)*(y2-y1+)*s;
             a[++tot]=qwe(,x1-,y1-,,tot);
             a[++tot]=qwe(,x2,y2,,tot);
             a[++tot]=qwe(,x1-,y2,,tot);
             a[++tot]=qwe(,x2,y1-,,tot);
         }
     }
     cdq(,tot);
     for(int i=;i<=cnt;i++)
         printf("%d\n",q[i]+ans[p[i]+]+ans[p[i]+]-ans[p[i]+]-ans[p[i]+]);
     return ;
 }