bzoj 2683: 简单题

Description

你有一个N*N的棋盘，每个格子内有一个整数，初始时的时候全部为0，现在需要维护两种操作：

 

命令	参数限制	内容
1 x y A	1<=x,y<=N，A是正整数	将格子x,y里的数字加上A
2 x1 y1 x2 y2	1<=x1<= x2<=N 1<=y1<= y2<=N	输出x1 y1 x2 y2这个矩形内的数字和
3	无	终止程序

Input

输入文件第一行一个正整数N。

接下来每行一个操作。

 

Output

对于每个2操作，输出一个对应的答案。

 

Sample Input

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

 

数据规模和约定

测试点分布如下：

1，2 1<=N<=100,操作数不超过100个，内存限制512M。 

3，4 1<=N<=1000,操作数不超过100000个，内存限制512M。 

5，6 1<=N<=200000,操作数不超过100000个，内存限制512M。 

7-10 1<=N<=500000,操作数不超过200000个，内存限制20M。

对于100%的数据，操作1中的A不超过2000。

同1176

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),然后把左边扫两遍把右边的询问计算两次就可以了。

然后再递归处理

【似乎稍微超了5M的内存限制。。】

#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[500001];
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()
{
     scanf("%d",&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]);
     return 0;
}