BZOJ3132: 上帝造题的七分钟

3132: 上帝造题的七分钟

Time Limit: 20 Sec  Memory Limit: 128 MB
Submit: 483  Solved: 222
[Submit][Status]

Description

“第一分钟，X说，要有矩阵，于是便有了一个里面写满了0的n×m矩阵。

第二分钟，L说，要能修改，于是便有了将左上角为(a,b)，右下角为(c,d)的一个矩形区域内的全部数字加上一个值的操作。

第三分钟，k说，要能查询，于是便有了求给定矩形区域内的全部数字和的操作。

第四分钟，彩虹喵说，要基于二叉树的数据结构，于是便有了数据范围。

第五分钟，和雪说，要有耐心，于是便有了时间限制。

第六分钟，吃钢琴男说，要省点事，于是便有了保证运算过程中及最终结果均不超过32位有符号整数类型的表示范围的限制。

第七分钟，这道题终于造完了，然而，造题的神牛们再也不想写这道题的程序了。”

——《上帝造裸题的七分钟》

所以这个神圣的任务就交给你了。

Input

输入数据的第一行为X n m，代表矩阵大小为n×m。

从输入数据的第二行开始到文件尾的每一行会出现以下两种操作：

　　L a b c d delta —— 代表将(a,b),(c,d)为顶点的矩形区域内的所有数字加上delta。

　　k a b c d　　 —— 代表求(a,b),(c,d)为顶点的矩形区域内所有数字的和。

请注意，k为小写。

Output

针对每个k操作，在单独的一行输出答案。

Sample Input

X 4 4
L 1 1 3 3 2
L 2 2 4 4 1
k 2 2 3 3

Sample Output

12

HINT

对于100%的数据，1 ≤ n ≤ 2048, 1 ≤ m ≤ 2048, 1 ≤ abs(delta) ≤ 500,操作不超过200000个,保证运算过程中及最终结果均不超过32位带符号整数类型的表示范围。

Source

tyvj

题解：
线段树做不了吧，MLE。。。
然后树状数组：
对于矩阵A，A[i][j]表示[i,j]-[n,m]的增量。那么子矩阵[1,1]-[x,y]的总和为：

代码：

 #include<cstdio>

 #include<cstdlib>

 #include<cmath>

 #include<cstring>

 #include<algorithm>

 #include<iostream>

 #include<vector>

 #include<map>

 #include<set>

 #include<queue>

 #include<string>

 #define inf 1000000000

 #define maxn 2050

 #define maxm 500+100

 #define eps 1e-10

 #define ll long long

 #define pa pair<int,int>

 #define for0(i,n) for(int i=0;i<=(n);i++)

 #define for1(i,n) for(int i=1;i<=(n);i++)

 #define for2(i,x,y) for(int i=(x);i<=(y);i++)

 #define for3(i,x,y) for(int i=(x);i>=(y);i--)

 #define mod 1000000007

 using namespace std;

 inline int read()

 {

     int x=,f=;char ch=getchar();

     while(ch<''||ch>''){if(ch=='-')f=-;ch=getchar();}

     while(ch>=''&&ch<=''){x=*x+ch-'';ch=getchar();}

     return x*f;

 }
 int n,m,s[][maxn][maxn];
 inline void add(int z,int x,int y,int t)
 {
     for(int i=x;i<=n;i+=i&(-i))
         for(int j=y;j<=m;j+=j&(-j))
             s[z][i][j]+=t;
 }
 inline int sum(int z,int x,int y)
 {
     int t=;
     for(int i=x;i;i-=i&(-i))
         for(int j=y;j;j-=j&(-j))
             t+=s[z][i][j];
     return t;
 }
 inline int getsum(int x,int y)
 {
     return (x+)*(y+)*sum(,x,y)-(y+)*sum(,x,y)-(x+)*sum(,x,y)+sum(,x,y);
 }

 int main()

 {

     freopen("input.txt","r",stdin);

     freopen("output.txt","w",stdout);

     n=read();m=read();char ch[];
     while(scanf("%s",ch)!=EOF)
     {
         if(ch[]=='L')
         {
             int x1=read(),y1=read(),x2=read(),y2=read(),z=read();
             add(,x1,y1,z);      add(,x2+,y2+,z);              add(,x2+,y1,-z);          add(,x1,y2+,-z);
             add(,x1,y1,z*x1);   add(,x2+,y2+,z*(x2+));       add(,x2+,y1,-z*(x2+));   add(,x1,y2+,-z*x1);
             add(,x1,y1,z*y1);   add(,x2+,y2+,z*(y2+));       add(,x2+,y1,-z*y1);       add(,x1,y2+,-z*(y2+));
             add(,x1,y1,z*x1*y1);add(,x2+,y2+,z*(x2+)*(y2+));add(,x2+,y1,-z*(x2+)*y1);add(,x1,y2+,-z*x1*(y2+));
         }
         else
         {
             int x1=read(),y1=read(),x2=read(),y2=read();
             printf("%d\n",getsum(x1-,y1-)+getsum(x2,y2)-getsum(x2,y1-)-getsum(x1-,y2));
         }
     }

     return ;

 }