hdu1892树状数组

按理来说应该会有二维的树状数组的算法吧。。可是我不知道。。于是就用一维的方式在累加了。。时间复杂度大概在o（m*n*log（n））。。直接暴力是o（m*n^2）

#include<iostream>
using namespace std;

const int MAXN=1003;
int data[MAXN][MAXN];
void swap(int &a,int &b)
{
    int t;
    t=a;
    a=b;
    b=t;
}
int lowbit(int idx)
{
    return idx & (-idx);
}
void add(int x,int y,int inc,int n)
{
    for (int i = x; i <= n; i += lowbit(i))
    {
        data[y][i] += inc;
    }
    return ;
}
int sum_of(int pos,int y) 
{
    int sum = 0;
    for (int i = pos; i > 0; i -= lowbit(i)) 
    {
        sum += data[y][i];
    }
    return sum;
}
void init()
{
    int i,j;
    for(i=1;i<=MAXN-1;i++)
        for(j=1;j<=MAXN-1;j++)
        {
            data[i][j]=lowbit(j);
            //cout<<data[i][j]<<"  ";
        }
            
}
int main()
{
    int t;
    cin>>t;
    int i=0;
    for(i=1;i<=t;i++)
    {
        cout<<"Case "<<i<<":"<<endl;
        init();
        int q;
        cin>>q;
        while(q--)
        {
            char c;
            cin>>c;
            switch(c)
            {
                int x1,y1,x2,y2,n1;
                case 'S':
                    {
                        int suum=0;
                        int j;
                        cin>>x1>>y1>>x2>>y2;
                        x1++;
                        x2++;
                        y1++;
                        y2++;
                        if(y1>y2)
                        {
                            swap(x1,x2);
                            swap(y1,y2);
                        }
                        
                        for(j=y1;j<=y2;j++)
                        {
                            if(x1>x2)
                            {
                                swap(x1,x2);
                            }
                            suum=suum+sum_of(x2,j)-sum_of(x1-1,j);
                        }
                        cout<<suum<<endl;
                    }
                    break;
                case 'A':
                    cin>>x1>>y1>>n1;
                    x1++;
                    y1++;
                    add(x1,y1,n1,MAXN);
                    break;
                case 'D':
                {
                    cin>>x1>>y1>>n1;
                    x1++;
                    y1++;
                    int numxy=sum_of(x1,y1)-sum_of(x1-1,y1);
                    if(numxy-n1<0)
                    {
                        add(x1,y1,-numxy,MAXN);
                    }else
                    {
                        add(x1,y1,-n1,MAXN);
                    }
                    break;
                }
                case 'M':
                {
                    cin>>x1>>y1>>x2>>y2>>n1;
                    x1++;
                    y1++;
                    x2++;
                    y2++;
                    int numxy=sum_of(x1,y1)-sum_of(x1-1,y1);
                    if(numxy-n1<0)
                    {
                        add(x1,y1,-numxy,MAXN);
                        add(x2,y2,numxy,MAXN);
                    }else
                    {
                        add(x1,y1,-n1,MAXN);
                        add(x2,y2,n1,MAXN);
                    }                    
                    break;
                }
            }
        } 
    }
    return 0;
}