bzoj 1176: [Balkan2007]Mokia

Description

维护一个W*W的矩阵，初始值均为S.每次操作可以增加某格子的权值,或询问某子矩阵的总权值.修改操作数M<=160000,询问数Q<=10000,W<=2000000.

Input

第一行两个整数,S,W;其中S为矩阵初始值;W为矩阵大小

接下来每行为一下三种输入之一(不包含引号):

"1 x y a"

"2 x1 y1 x2 y2"

"3"

输入1:你需要把(x,y)(第x行第y列)的格子权值增加a

输入2:你需要求出以左上角为(x1,y1),右下角为(x2,y2)的矩阵内所有格子的权值和,并输出

输入3:表示输入结束

Output

对于每个输入2,输出一行,即输入2的答案

Sample Input

0 4
1 2 3 3
2 1 1 3 3
1 2 2 2
2 2 2 3 4
3

Sample Output

3
5

HINT

保证答案不会超过int范围

CDQ分治。。我们读入所有操作后按照操作时间排序【其实就是读入顺序】。这样可以保证右边不对左边产生影响

然后solve(l,r)分成solve(l,mid)和(mid+1,r)来做，把左边的修改按照x排序。然后用树状数组记录y轴上的总和。

每个询问(x1,y1) (x2,y2)变成(1,y1)(x1-1,y2)和(1,y1)(x2,y2),然后把左边扫两遍把右边的询问计算两次就可以了。

然后再递归处理

<pre name="code" class="cpp">#include<cstdio>
#include<algorithm>
using namespace std;
struct add
{
	 int s;
     int x,y,a;
     int x1,y1,x2,y2;
     int p;
}a[200001],b[200001],c[200001];
int ans[200001];
inline bool cmp1(add x,add y)
{
     if(x.x<y.x)
          return true;
     //if(x.x==y.x&&x.y<y.y)
       //   return true;
     return false;
}
inline bool cmp2(add x,add y)
{
     if(x.x1<y.x1)
          return true;
     return false;
}
inline bool cmp3(add x,add y)
{
     if(x.x2<y.x2)
          return true;
     return false;
}
int tr[2000001];
int w;
inline int lowbit(int x)
{
     return x&(-x);
}
inline void add(int x,int xx)
{
     int i;
     for(i=x;i<=w;i+=lowbit(i))
          tr[i]+=xx;
}
inline int ask(int x)
{
     int i;
     int ans=0;
     for(i=x;i>=1;i-=lowbit(i))
          ans+=tr[i];
     return ans;
}
inline void solve(int l,int r)
{
     if(l!=r)
     {
     	  int mid=(l+r)/2;
          solve(l,mid);
          solve(mid+1,r);
          int i;
          int p=0;
          for(i=l;i<=mid;i++)
          {
               if(a[i].s==1)
               {
                    p++;
                    b[p]=a[i];
               }
          }
          sort(b+1,b+1+p,cmp1);
          int p1=0;
          for(i=mid+1;i<=r;i++)
          {
               if(a[i].s==2)
               {
                    p1++;
                    c[p1]=a[i];
               }
          }
          sort(c+1,c+1+p1,cmp2);
          int i1=1,i2=1;
          while(i1<=p&&i2<=p1)
          {
               if(b[i1].x<=c[i2].x1-1)
               {
               	    add(b[i1].y,b[i1].a);
                    i1++;
               }
               else
               {
                    ans[c[i2].p]-=(ask(c[i2].y2)-ask(c[i2].y1-1));
                    i2++;
               }
          }
          if(i2<=p1)
               for(i=i2;i<=p1;i++)
                    ans[c[i].p]-=(ask(c[i].y2)-ask(c[i].y1-1));
          for(i=1;i<=i1-1;i++)
               add(b[i].y,-b[i].a);
          sort(c+1,c+1+p1,cmp3);
          i1=1;
		  i2=1;
          while(i1<=p&&i2<=p1)
          {
               if(b[i1].x<=c[i2].x2)
               {
               	    add(b[i1].y,b[i1].a);
                    i1++;
               }
               else
               {
                    ans[c[i2].p]+=(ask(c[i2].y2)-ask(c[i2].y1-1));
                    i2++;
               }
          }
          if(i2<=p1)
               for(i=i2;i<=p1;i++)
                    ans[c[i].p]+=(ask(c[i].y2)-ask(c[i].y1-1));
          for(i=1;i<=i1-1;i++)
               add(b[i].y,-b[i].a);
     }
}
int main()
{
     int s;
     scanf("%d%d",&s,&w);
     int x;
     int p=0;
     scanf("%d",&x);
     while(x!=3)
     {
          if(x==1)
          {
               p++;
               a[p].s=1;
               scanf("%d%d%d",&a[p].x,&a[p].y,&a[p].a);
          }
          else
          {
               p++;
               a[p].s=2;
               scanf("%d%d%d%d",&a[p].x1,&a[p].y1,&a[p].x2,&a[p].y2);
          }
          a[p].p=p;
          scanf("%d",&x);
     }
    // sort(a+1,a+1+p1,cmp1);
     solve(1,p);
     int i;
     for(i=1;i<=p;i++)
          if(a[i].s==2)
               printf("%d\n",ans[i]+(a[i].x2-a[i].x1+1)*(a[i].y2-a[i].y1+1)*s);
     return 0;
}