TYVJ1716 上帝造题的七分钟

 

时间: 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 }

 

 

转载于:https://www.cnblogs.com/SilverNebula/p/5994716.html