POJ1195 二维线段树

Mobile phones

 POJ - 1195 

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
——————————————————————————————————————————————————————————————————————————————————————————————————————————————————————————————————————————————
二维线段树，单点修改，区间查询，维护和。
______________________________________________________________________________________________________________________________________________

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
using namespace std;

const int maxn=1028;
struct LIE
{
    int ll,lr,sum;
};
struct HANG
{
    int hl,hr;
    LIE lie[maxn<<2];
}hang[maxn<<2];
int op,n;
void readint(int &x)
{
    char c=getchar();
    int f=1;
    for(;c<'0' || c>'9';c=getchar())if(c=='-')f=-f;
    x=0;
    for(;c<='9'&& c>='0';c=getchar())x=(x<<1)+(x<<3)+c-'0';
    x*=f;
}
void buil(int pre,int cur,int ll,int lr)
{
    hang[pre].lie[cur].ll=ll;hang[pre].lie[cur].lr=lr;
    hang[pre].lie[cur].sum=0;
    if(ll==lr)return ;
    int mid=(ll+lr)>>1;
    buil(pre,cur<<1,ll,mid);
    buil(pre,cur<<1|1,mid+1,lr);
}
void build(int cur,int l,int r,int ll,int rr)
{
    hang[cur].hl=l;hang[cur].hr=r;
    buil(cur,1,ll,rr);
    if(l==r)return ;
    int mid=(l+r)>>1;
    build(cur<<1,l,mid,ll,rr);
    build(cur<<1|1,mid+1,r,ll,rr);
}
void upda(int pre,int cur,int hh,int lh,int dat)
{
    hang[pre].lie[cur].sum+=dat;
    if(hang[pre].lie[cur].ll==hang[pre].lie[cur].lr)return ;
    int mid=(hang[pre].lie[cur].ll+hang[pre].lie[cur].lr)>>1;
    if(lh<=mid) upda(pre,cur<<1,hh,lh,dat);
    else upda(pre,cur<<1|1,hh,lh,dat);
}
void update(int cur,int x,int y,int dat)
{
    upda(cur,1,x,y,dat);
    if(hang[cur].hl==hang[cur].hr)return;
    int mid=(hang[cur].hl+hang[cur].hr)>>1;
    if(x<=mid)update(cur<<1,x,y,dat);
    else update(cur<<1|1,x,y,dat);
}
int qure(int pre,int cur,int yl,int yr)
{
    if(yl<=hang[pre].lie[cur].ll && hang[pre].lie[cur].lr<=yr)return hang[pre].lie[cur].sum;
    int mid=(hang[pre].lie[cur].ll+hang[pre].lie[cur].lr)>>1;
    int sum=0;
    if(yl<=mid)sum+=qure(pre,cur<<1,yl,yr);
    if(mid<yr)sum+=qure(pre,cur<<1|1,yl,yr);
    return sum; 
}
int query(int cur,int xl,int yl,int xr,int yr)
{
    if(xl<=hang[cur].hl && hang[cur].hr<=xr)return qure(cur,1,yl,yr);
    int mid=(hang[cur].hl+hang[cur].hr)>>1;
    int sum=0;
    if(xl<=mid)sum+=query(cur<<1,xl,yl,xr,yr);
    if(mid<xr) sum+=query(cur<<1|1,xl,yl,xr,yr);
    return sum;
}
int main()
{
    readint(op);
    while(op!=3)
    {
        if(!op)
        {
            readint(n);
            build(1,0,n-1,0,n-1);
        }
        else if(op==1)
        {
            int x,y,dat;
            readint(x);readint(y);readint(dat);
            update(1,x,y,dat);
        }
        else if(op==2)
        {
            int xl,xr,yl,yr;
            readint(xl);readint(yl);readint(xr);readint(yr);
            printf("%d
",query(1,xl,yl,xr,yr));
        }
        readint(op);
    }
    return 0;
}