hdu 1892 二维树状数组

#include<cstdio>
#include<algorithm>
using namespace std;
#define maxn 1111
int c[maxn][maxn],n;
int lb(int x){return x&(-x);}
void update(int x,int y,int w)
{
    for(int i=x;i<=n;i+=lb(i))
        for(int j=y;j<=n;j+=lb(j))
            c[i][j]+=w;
}
int sum(int x,int y)
{
    int ans=0;
    for(int i=x;i>0;i-=lb(i))
        for(int j=y;j>0;j-=lb(j))
            ans+=c[i][j];
    return ans;
}
int main()
{
    int t,q;
    char ss[5];
    int x1,y1,x2,y2,n1;
    while(~scanf("%d",&t))
        for(int z=1;z<=t;++z)
        {
            printf("Case %d:\n",z);
            
            n=1005;
            memset(c,0,sizeof c);
            scanf("%d",&q);
            for(int i=1;i<=n;++i)
                for(int j=1;j<=n;++j)
//                     update(i,j,1);
//                    c[i][j]=lb(i)*lb(j);
            for(int i=0;i<q;++i)
            {
                scanf("%s",ss);
                if(ss[0]=='A')
                {
                    scanf("%d%d%d",&x1,&y1,&n1);
                    x1++,y1++;
                    update(x1,y1,n1);
                }
                else if(ss[0]=='D')
                {
                    scanf("%d%d%d",&x1,&y1,&n1);
                    x1++,y1++;
                    int ans=min(n1,sum(x1,y1)+sum(x1-1,y1-1)-sum(x1,y1-1)-sum(x1-1,y1));
                    update(x1,y1,-ans);
                    
                }
                else if(ss[0]=='S')
                {
                    scanf("%d%d%d%d",&x1,&y1,&x2,&y2);
                    x1++;x2++;y1++;y2++;
                    if(x1>x2)swap(x1,x2);
                    if(y1>y2)swap(y1,y2);
                    printf("%d\n",sum(x2,y2)+sum(x1-1,y1-1)-sum(x1-1,y2)-sum(x2,y1-1));
                }
                else if(ss[0]=='M')
                {
                    scanf("%d%d%d%d%d",&x1,&y1,&x2,&y2,&n1);
                    x1++,x2++,y1++,y2++;
                    int ans=min(sum(x1,y1)+sum(x1-1,y1-1)-sum(x1-1,y1)-sum(x1,y1-1),n1);
                    update(x1,y1,-ans);
                    update(x2,y2,ans);
                }
                
            }
            
        }
    return 0;
    
}

纯裸的二维 关键点反而在于初始化 怎么来表示树状组的每个位置都被add 1 了呢

直接c[i][j]=lowbit(i)*lowbit(j);

这里的lowbit可以看成是每个维度上分别有几组

因此在相加了之后就会变成两者的乘积

当然也可以暴力add 1 不过这样做的话效率比较低而已