bzoj3132 上帝造题的七分钟（差分+二维树状数组）

---

bzoj3132 上帝造题的七分钟（差分+二维树状数组）##

　　Time Limit: 20 Sec
　　Memory Limit: 128 MB

Description###

　　 “第一分钟，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){
    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*x;
            C[i][j]+=val*y;
            D[i][j]+=val*x*y;
        }
}
int query(int x,int y){
    if(x==0 || y==0)return 0;
    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(){
    freopen("1.in","r",stdin);
    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/
　　本文版权归作者和博客园共有，欢迎转载，但未经作者同意必须保留此段声明，且在文章页面明显位置给出原文连接，否则保留追究法律责任的权利。