COGS1532. [IOI2001]移动电话

1532. [IOI2001]移动电话

★☆ 输入文件：mobilephones.in 输出文件：mobilephones.out 简单对比
时间限制：5 s 内存限制：256 MB

【题目描述】

假设Tampere地区的4G移动通信基站以如下方式运行。整个地区被划分成若干正方形格子。这些格子构成一个S*S的矩阵，它们的行，列编号都是从0到S-1.每一个格子中都有一个基站。每个格子中激活的手机数量可能改变，因为一部手机可能从一个格子移动到另一个格子，打开或者关闭。有时，某一座基站会向总站报告自己的行列坐标，以及该格中激活手机数目的变化。

【输入格式】

输入指令编码如下。

每个指令占一行，包含一个指令码和一些参数，见下表。

指令码	参数	意义
0	S	初始化一个S*S的全零矩阵。这个指令只会出现一次，即第一条指令。
1	X Y A	将格子(X,Y)中激活的手机数量增加A。A有可能是正数或负数。
2	L B R T	询问当前所有坐标(X,Y)满足：L<=X<=R,B<=Y<=T的格子中激活的手机数量之和。
3		结束程序。这个指令只会出现一次，即最后一条指令。

【输出格式】

你的程序不应该对指令2外的所有指令进行回答。对于每个指令2，你的程序需要输出一行一个正整数，即该指令的答案。

【样例输入】

0 4

1 1 2 3

2 0 0 2 2

1 1 1 2

1 1 2 -1

2 1 1 2 3

【样例

【提示】

矩阵大小：1<=S<=1024

任意时刻，每个格子中的激活手机数量V：0<=V<=32767

格子中激活手机数量的变化值：-32768<=A<=32767

输入的指令数目：3<=U<=60002

整个矩阵中的最大手机数量：M=2^30

【来源】

IOI 2001

二维BIT维护区间和

#include<cstdio>
const int N=3005;
int opt,n,c[N][N];
int read(){
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
int lowbit(int x){
    return x&-x;
}
void updata(int x,int y,int val){
    for(int i=x;i<=n;i+=lowbit(i)){
        for(int j=y;j<=n;j+=lowbit(j)){
            c[i][j]+=val;
        } 
    }
}
int query(int x,int y){
    int res=0;
    for(int i=x;i>=1;i-=lowbit(i)){
        for(int j=y;j>=1;j-=lowbit(j)){
            res+=c[i][j];
        } 
    }
    return res;
}
int main(){
    freopen("mobilephones.in","r",stdin);
    freopen("mobilephones.out","w",stdout);
    opt=read();n=read();
    for(int x1,x2,y1,y2,val;(opt=read())!=3;){
        if(opt==1){
            x1=read();y1=read();val=read();x1++;y1++;
            updata(x1,y1,val);
        }
        else{
            x1=read();y1=read();x1++;y1++;
            x2=read();y2=read();x2++;y2++;
            printf("%d
",query(x2,y2)-query(x2,y1-1)-query(x1-1,y2)+query(x1-1,y1-1));
        }
    }
    return 0;
}