【BZOJ3132】上帝造题的七分钟

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位带符号整数类型的表示范围。

Source

tyvj

合理差分,拯救世界
刷空水题,目空一切

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define MAXN 2050
#define lowbit(x)   (x&-x)
#define GET (ch>='0'&&ch<='9')
using namespace std;
char ch[2];
int n,m;
int x1,y1,x2,y2,w;
void in(int &x)
{
    char ch=getchar();x=0;int flag=1;
    while (!GET)    flag=ch=='-'?-1:1,ch=getchar();
    while (GET) x=x*10+ch-'0',ch=getchar();x*=flag;
}
struct BIT
{
    int c[MAXN][MAXN];
    inline void add(int x,int y,int delta)
    {
        for (int i=x;i<=n;i+=lowbit(i)) for (int j=y;j<=m;j+=lowbit(j)) c[i][j]+=delta;
    }
    inline int query(int x,int y)
    {
        int ret=0;
        for (int i=x;i;i-=lowbit(i))    for (int j=y;j;j-=lowbit(j))    ret+=c[i][j];
        return ret;
    }
}a,b,c,d;
inline int sum(int x,int y) {   return (x+1)*(y+1)*a.query(x,y)-(y+1)*b.query(x,y)-(x+1)*c.query(x,y)+d.query(x,y); }
inline void modify(int x1,int y1,int x2,int y2,int delta)
{
    a.add(x1,y1,delta);a.add(x2+1,y1,-delta);a.add(x1,y2+1,-delta);a.add(x2+1,y2+1,delta);
    b.add(x1,y1,x1*delta);b.add(x2+1,y1,-(x2+1)*delta);b.add(x1,y2+1,-x1*delta);b.add(x2+1,y2+1,(x2+1)*delta);
    c.add(x1,y1,y1*delta);c.add(x2+1,y1,-y1*delta);c.add(x1,y2+1,-(y2+1)*delta);c.add(x2+1,y2+1,(y2+1)*delta);
    d.add(x1,y1,x1*y1*delta);d.add(x2+1,y1,-(x2+1)*y1*delta);d.add(x1,y2+1,-x1*(y2+1)*delta);d.add(x2+1,y2+1,(x2+1)*(y2+1)*delta);
}
int main()
{
    scanf("%s",ch);
    while (ch[0]!='X')  scanf("%s",ch);
    in(n);in(m);
    while (~scanf("%s",ch))
    {
        in(x1);in(y1);in(x2);in(y2);
        if (ch[0]=='L') in(w),modify(x1,y1,x2,y2,w);
        if (ch[0]=='k') printf("%d\n",sum(x2,y2)-sum(x1-1,y2)-sum(x2,y1-1)+sum(x1-1,y1-1));
    }
}