POJ 1195 二维树状数组

Mobile phones

Description

Suppose that the fourth generation mobile phone base stations in the Tampere area operate as follows. The area is divided into squares. The squares form an S * S matrix with the rows and columns numbered from 0 to S-1. Each square contains a base station. The number of active mobile phones inside a square can change because a phone is moved from a square to another or a phone is switched on or off. At times, each base station reports the change in the number of active phones to the main base station along with the row and the column of the matrix.

Write a program, which receives these reports and answers queries
about the current total number of active mobile phones in any
rectangle-shaped area.

Input

The
input is read from standard input as integers and the answers to the
queries are written to standard output as integers. The input is encoded
as follows. Each input comes on a separate line, and consists of one
instruction integer and a number of parameter integers according to the
following table.

The values will always be in range, so there is no need to check
them. In particular, if A is negative, it can be assumed that it will
not reduce the square value below zero. The indexing starts at 0, e.g.
for a table of size 4 * 4, we have 0 <= X <= 3 and 0 <= Y <=
3.

Table size: 1 * 1 <= S * S <= 1024 * 1024

Cell value V at any time: 0 <= V <= 32767

Update amount: -32768 <= A <= 32767

No of instructions in input: 3 <= U <= 60002

Maximum number of phones in the whole table: M= 2^30

Output

Your
program should not answer anything to lines with an instruction other
than 2. If the instruction is 2, then your program is expected to answer
the query by writing the answer as a single line containing a single
integer to standard output.

Sample Input

0 4
1 1 2 3
2 0 0 2 2
1 1 1 2
1 1 2 -1
2 1 1 2 3
3

Sample Output

3
4

Source

IOI 2001

题意：

在一个二维平面内的所有点，共有三种操作，0 S表示把S*S区域内的点的值清零；X，Y，A表示将点（X,Y）的值加A；L,B,R,T表示求（L,B）~（R,T）之间的所有的点的值之和

代码：

 1 //树状数组模板
 2 #include<iostream>
 3 #include<cstdio>
 4 #include<cstring>
 5 using namespace std;
 6 int flag,X,Y,A,L,B,R,T,S;
 7 int squ[1050][1050];
 8 int lowbit(int x)
 9 {
10     return x&(-x);
11 }
12 void add(int idx,int idy,int val)
13 {
14     for(int i=idx;i<=1024;i+=lowbit(i))
15     {
16         for(int j=idy;j<=1024;j+=lowbit(j))
17         {
18             squ[i][j]+=val;
19         }
20     }
21 }
22 int sum(int idx,int idy)
23 {
24     int s=0;
25     for(int i=idx;i>0;i-=lowbit(i))
26     {
27         for(int j=idy;j>0;j-=lowbit(j))
28         {
29             s+=squ[i][j];
30         }
31     }
32     return s;
33 }
34 int ans(int x1,int y1,int x2,int y2)
35 {
36     return (sum(x2,y2)-sum(x1-1,y2)-sum(x2,y1-1)+sum(x1-1,y1-1));
37 }
38 int main()
39 {
40     while(scanf("%d",&flag)&&flag!=3)
41     {
42         if(flag==0)
43         {
44             scanf("%d",&S);
45             memset(squ,0,sizeof(squ));
46         }
47         else if(flag==1)
48         {
49             scanf("%d%d%d",&X,&Y,&A);
50             if(A>=0)
51             add(X+1,Y+1,A);
52             else
53             {
54                 int num=ans(X+1,Y+1,X+1,Y+1);
55                 add(X+1,Y+1,num+A>0?A:-num);
56             }
57         }
58         else if(flag==2)
59         {
60             scanf("%d%d%d%d",&L,&B,&R,&T);
61             printf("%d\n",ans(L+1,B+1,R+1,T+1));
62         }
63     }
64     return 0;
65 }