题目背景
裸体就意味着身体。
题目描述
“第一分钟,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\n",ASK(a-1,b-1)+ASK(c,d)-ASK(a-1,d)-ASK(c,b-1)); 55 } 56 return 0; 57 }