hdu 1892 See you 二维树状数组，很好的一道题

1. 
   
2. 题目给定4种操作:  S x1 y1 x2 y2 询问以(x1 , y1) - (x2 , y2)为对角线的矩形的面积，但是这个对角线不一定是正对角线。A x1 y1 n 把点(x1 , y1)加上n。D x1 y1 n点(x1 , y1)减去n如果不足n就全部删除即可。M x1 y1 x2 y2 n 把点(x1 , y1)点值中扣除n加到(x2 , y2)，如果不过n则把(x1 , y1)值全部加到(x2 , y2)
   
3. #include<cstdio>  
4. #include<cstring>  
5. #include<iostream>  
6. #include<algorithm>  
7. using namespace std;  
8.   
9. const int MAXN = 1e3+10;  
10.   
11. int n;  
12. int num[MAXN][MAXN];  
13. int treeNum[MAXN][MAXN];  
14.   
15. int lowbit(int x){  
16.     return x&(-x);  
17. }  
18.   
19. int getSum(int x , int y){  
20.     int sum = 0;  
21.     for(int i = x ; i > 0 ; i -= lowbit(i))   
22.         for(int j = y ; j > 0 ; j -= lowbit(j))   
23.             sum += treeNum[i][j];  
24.     return sum;  
25. }  
26.   
27. void add(int x , int y , int val){  
28.     for(int i = x ; i < MAXN ; i += lowbit(i))   
29.         for(int j = y ; j < MAXN ; j += lowbit(j))   
30.             treeNum[i][j] += val;  
31. }  
32.   
33. void init(){  
34.     memset(treeNum , 0 , sizeof(treeNum));  
35.     for(int i = 1 ; i < MAXN ; i++){  
36.        for(int j = 1 ; j < MAXN ; j++){  
37.            num[i][j] = 1;  
38.            treeNum[i][j] = lowbit(i)*lowbit(j);  
39.        }  
40.     }  
41. }  
42.   
43. void solve(){  
44.     init();  
45.     char c;  
46.     int x , y , val;  
47.     int x1 , y1 , x2 , y2;  
48.     int x3 , y3 , x4 , y4;  
49.     while(n--){  
50.         scanf("%c" , &c);   
51.         if(c == 'S'){  
52.             scanf("%d%d" , &x1 , &y1);   
53.             scanf("%d%d%*c" , &x2 , &y2);   
54.             x1++ , y1++ , x2++ , y2++;  
55.   
56.             x3 = min(x1 , x2) , y3 = min(y1 , y2);  
57.             x4 = max(x1 , x2) , y4 = max(y1 , y2);  
58.             int ans = getSum(x4 , y4);   
59.             ans -= getSum(x3-1 , y4);  
60.             ans -= getSum(x4 , y3-1);  
61.             ans += getSum(x3-1 , y3-1);  
62.             printf("%d\n" , ans);  
63.         }  
64.         else if(c == 'A'){  
65.             scanf("%d%d%d%*c" , &x , &y , &val);   
66.             x++ , y++;  
67.             num[x][y] += val;  
68.             add(x , y , val);  
69.         }  
70.         else if(c == 'D'){  
71.             scanf("%d%d%d%*c" , &x , &y , &val);   
72.             x++ , y++;  
73.             val = min(num[x][y] , val);  
74.             num[x][y] -= val;  
75.             add(x , y , -val);  
76.         }  
77.         else{  
78.             scanf("%d%d" , &x1 , &y1);   
79.             scanf("%d%d%d%*c" , &x2 , &y2 , &val);   
80.             x1++ , y1++ , x2++ , y2++;  
81.             val = min(num[x1][y1] , val);  
82.             num[x1][y1] -= val;  
83.             num[x2][y2] += val;  
84.             add(x1 , y1 , -val);  
85.             add(x2 , y2 , val);  
86.         }  
87.     }  
88. }  
89.   
90. int main(){  
91.     int cas;  
92.     int Case = 1;  
93.     scanf("%d" , &cas);  
94.     while(cas--){  
95.         scanf("%d%*c" , &n);       
96.         printf("Case %d:\n" , Case++);  
97.         solve();  
98.     }  
99.     return 0;  
100. }