bzoj3132 上帝造题的七分钟(差分+二维树状数组)
Time Limit: 20 Sec
Memory Limit: 128 MBDescription
“第一分钟,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位带符号整数类型的表示范围。
题目地址:bzoj3132 上帝造题的七分钟
题目大意: 题目很简洁了:)
题解:
dalao写的文章真的好 RabbitHu
刷历年noip ing
自己代码写太好了过来贴一发自恋中
AC代码
#include <cstdio>
using namespace std;
const int N=2050;
int n,m,x0,y0,x1,y1,delta;
int A[N][N],B[N][N],C[N][N],D[N][N];
char ch;
void change(int x,int y,int val){
x++,y++;
for(int i=x;i<=n;i+=i&(-i))
for(int j=y;j<=m;j+=j&(-j)){
A[i][j]+=val;
B[i][j]+=val*i;
C[i][j]+=val*j;
D[i][j]+=val*i*j;
}
}
int query(int x,int y){
x++,y++;
int res=0;
for(int i=x;i>=1;i-=i&(-i))
for(int j=y;j>=1;j-=j&(-j))
res+=A[i][j]*(x+1)*(y+1)-B[i][j]*(y+1)-C[i][j]*(x+1)+D[i][j];
return res;
}
int main(){
scanf("%c%d%d\n",&ch,&n,&m);
while(scanf("%c",&ch)!=EOF){
if(ch=='L'){
scanf("%d%d%d%d%d\n",&x0,&y0,&x1,&y1,&delta);
change(x0,y0,delta);
change(x0,y1+1,-delta);
change(x1+1,y0,-delta);
change(x1+1,y1+1,delta);
}else{
scanf("%d%d%d%d\n",&x0,&y0,&x1,&y1);
printf("%d\n",query(x1,y1)-query(x1,y0-1)-query(x0-1,y1)+query(x0-1,y0-1));
}
}
return 0;
}
作者:skl_win
出处:https://www.cnblogs.com/shaokele/
本文版权归作者和博客园共有,欢迎转载,但未经作者同意必须保留此段声明,且在文章页面明显位置给出原文连接,否则保留追究法律责任的权利。
原文地址:https://www.cnblogs.com/shaokele/p/9910997.html