题目背景
裸体就意味着身体。
题目描述
“第一分钟，X说，要有矩阵，于是便有了一个里面写满了 000 的 n×mn×mn×m 矩阵。
第二分钟，L说，要能修改，于是便有了将左上角为 (a,b)(a,b)(a,b) ，右下角为 (c,d)(c,d)(c,d) 的一个矩形区域内的全部数字加上一个值的操作。
第三分钟，k说，要能查询，于是便有了求给定矩形区域内的全部数字和的操作。
第四分钟，彩虹喵说，要基于二叉树的数据结构，于是便有了数据范围。
第五分钟，和雪说，要有耐心，于是便有了时间限制。
第六分钟，吃钢琴男说，要省点事，于是便有了保证运算过程中及最终结果均不超过32位有符号整数类型的表示范围的限制。
第七分钟，这道题终于造完了，然而，造题的神牛们再也不想写这道题的程序了。”
——《上帝造裸题的七分钟》
所以这个神圣的任务就交给你了。
输入输出格式
输入格式：
输入数据的第一行为X n m，代表矩阵大小为 n×mn×mn×m 。
从输入数据的第二行开始到文件尾的每一行会出现以下两种操作：
L a b c d delta—— 代表将 (a,b),(c,d)(a,b),(c,d)(a,b),(c,d) 为顶点的矩形区域内的所有数字加上delta。k a b c d—— 代表求 (a,b),(c,d)(a,b),(c,d)(a,b),(c,d) 为顶点的矩形区域内所有数字的和。
请注意， kkk 为小写。
输出格式：
针对每个k操作，在单独的一行输出答案。
输入输出样例
说明
对于10%的数据， 1≤n≤16,1≤m≤161 ≤ n ≤ 16, 1 ≤ m ≤ 161≤n≤16,1≤m≤16 , 操作不超过200个.
对于60%的数据， 1≤n≤512,1≤m≤5121 ≤ n ≤ 512, 1 ≤ m ≤ 5121≤n≤512,1≤m≤512 .
对于100%的数据， 1≤n≤2048,1≤m≤2048,−500≤delta≤5001 ≤ n ≤ 2048, 1 ≤ m ≤ 2048, -500 ≤ delta ≤ 5001≤n≤2048,1≤m≤2048,−500≤delta≤500 ,操作不超过200000个,保证运算过程中及最终结果均不超过32位带符号整数类型的表示范围。
【思路分析】
一道描述很清奇的题目，而且居然是黑题！！！？？？机房有人选这道题做难题选讲，结果题目一给就基本会做了（至少会口胡）。如果上一道是很水的LCT的话，那么这道题就是更更水的二维区间修改，显然对于二维区间修改我们可以想到的算法有：二维树状数组，CDQ分治，KD—Tree。（然而感觉这辈子都不会写KD—Tree了）
鉴于这道题的数据范围很小，我们就选择码量较小的二维树状数组（可能不会少很多，而且应该是最慢的）。
不会二维树状数组的读者可以看这篇博客，感觉讲得蛮不错。
【代码实现】
1 #include<cstdio> 2 #include<algorithm> 3 #include<cctype> 4 #define lowbit(x) x&(-x) 5 using namespace std; 6 const int maxn=2050; 7 int t[maxn][maxn][4],n,m,a,b,c,d,w; 8 inline int read() 9 { 10 char c; 11 int f,x=0; 12 while(!isdigit(c=getchar())&&c!='-'); c=='-'?(f=-1,x=0):(f=1,x=c-'0'); 13 while(isdigit(c=getchar())) x=x*10+c-'0'; 14 return x*f; 15 } 16 void add(int x,int y,int pp,int ww) 17 { 18 for(int i=x;i<=n;i+=lowbit(i)) 19 for(int j=y;j<=m;j+=lowbit(j)) 20 t[i][j][pp]+=ww; 21 } 22 int ask(int x,int y,int pp) 23 { 24 int ans=0; 25 for(int i=x;i;i-=lowbit(i)) 26 for(int j=y;j;j-=lowbit(j)) 27 ans+=t[i][j][pp]; 28 return ans; 29 } 30 void ADD(int x,int y,int ww) 31 { 32 add(x,y,0,ww); 33 add(x,y,1,ww*x); 34 add(x,y,2,ww*y); 35 add(x,y,3,ww*x*y); 36 } 37 int ASK(int x,int y) 38 { 39 int ans=ask(x,y,0)*(x*y+x+y+1)-ask(x,y,1)*(y+1)-ask(x,y,2)*(x+1)+ask(x,y,3); 40 return ans; 41 } 42 int main() 43 { 44 char ch[3]; 45 scanf("%s%d%d",ch,&n,&m); 46 while(~scanf("%s",ch)) 47 { 48 a=read(),b=read(),c=read(),d=read(); 49 if(ch[0]=='L') 50 { 51 w=read(); 52 ADD(a,b,w),ADD(a,d+1,-w),ADD(c+1,b,-w),ADD(c+1,d+1,w); 53 } 54 if(ch[0]=='k') printf("%d ",ASK(a-1,b-1)+ASK(c,d)-ASK(a-1,d)-ASK(c,b-1)); 55 } 56 return 0; 57 }