时间: 1000ms / 空间: 131072KiB / Java类名: Main
背景
裸体就意味着身体。
描述
“第一分钟,X说,要有矩阵,于是便有了一个里面写满了0的n×m矩阵。
第二分钟,L说,要能修改,于是便有了将左上角为(a,b),右下角为(c,d)的一个矩形区域内的全部数字加上一个值的操作。
第三分钟,k说,要能查询,于是便有了求给定矩形区域内的全部数字和的操作。
第四分钟,彩虹喵说,要基于二叉树的数据结构,于是便有了数据范围。
第五分钟,和雪说,要有耐心,于是便有了时间限制。
第六分钟,吃钢琴男说,要省点事,于是便有了保证运算过程中及最终结果均不超过32位有符号整数类型的表示范围的限制。
第七分钟,这道题终于造完了,然而,造题的神牛们再也不想写这道题的程序了。”
——《上帝造裸题的七分钟》
所以这个神圣的任务就交给你了。
输入格式
输入数据的第一行为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为小写。
输出格式
针对每个k操作,在单独的一行输出答案。
测试样例1
输入
X 4 4
L 1 1 3 3 2
L 2 2 4 4 1
k 2 2 3 3
输出
12
备注
对于10%的数据,1 ≤ n ≤ 16, 1 ≤ m ≤ 16, 操作不超过200个.
对于60%的数据,1 ≤ n ≤ 512, 1 ≤ m ≤ 512.
对于100%的数据,1 ≤ n ≤ 2048, 1 ≤ m ≤ 2048, 1 ≤ delta ≤ 500,操作不超过200000个,保证运算过程中及最终结果均不超过32位带符号整数类型的表示范围。
by XLk
树状数组区间修改区间查询,模板题。
然而并不能看懂,姑且死记模板233
1 /*by SilverN*/ 2 #include<algorithm> 3 #include<iostream> 4 #include<cstring> 5 #include<cstdio> 6 #include<cmath> 7 #include<vector> 8 using namespace std; 9 const int mxn=2560; 10 int read(){ 11 int x=0,f=1;char ch=getchar(); 12 while(ch<‘0‘ || ch>‘9‘){if(ch==‘-‘)f=-1;ch=getchar();} 13 while(ch>=‘0‘ && ch<=‘9‘){x=x*10+ch-‘0‘;ch=getchar();} 14 return x*f; 15 } 16 int n,m; 17 int lmn; 18 int t[mxn][mxn][5]; 19 void add(int x,int y,int v,int k){ 20 while(x<=n){ 21 int tmp=y; 22 while(tmp<=m){ t[x][tmp][k]+=v; tmp+=tmp&-tmp; } 23 x+=x&-x; 24 } 25 return; 26 } 27 int smm(int x,int y,int k){ 28 int res=0; 29 while(x){ 30 int tmp=y; 31 while(tmp){ res+=t[x][tmp][k]; tmp-=tmp&-tmp; } 32 x-=x&-x; 33 } 34 return res; 35 } 36 int clc(int x,int y){ 37 return smm(x,y,1)*x*y-smm(x,y,2)*x-smm(x,y,3)*y+smm(x,y,4); 38 } 39 int main(){ 40 getchar(); 41 n=read();m=read(); 42 char op[2]; 43 int a,b,c,d,w;//x1,y1,x2,y2 44 while(scanf("%s",op)!=EOF){ 45 a=read();b=read();c=read();d=read(); 46 if(b>d)swap(b,d); 47 if(a>c)swap(a,c); 48 if(op[0]==‘L‘){ 49 w=read(); 50 add(a,b,w,1);add(a,d+1,-w,1); 51 add(c+1,b,-w,1);add(c+1,d+1,w,1); 52 // 53 add(a,b,w*(b-1),2);add(a,d+1,-w*d,2); 54 add(c+1,b,-w*(b-1),2);add(c+1,d+1,w*d,2); 55 // 56 add(a,b,w*(a-1),3);add(a,d+1,-w*(a-1),3); 57 add(c+1,b,-w*c,3);add(c+1,d+1,w*c,3); 58 // 59 add(a,b,w*(a-1)*(b-1),4);add(a,d+1,-w*(a-1)*d,4); 60 add(c+1,b,-w*c*(b-1),4);add(c+1,d+1,w*c*d,4); 61 // 62 } 63 else{ 64 int ans=clc(c,d)+clc(a-1,b-1)-clc(a-1,d)-clc(c,b-1); 65 printf("%d\n",ans); 66 } 67 } 68 return 0; 69 }