树状数组模板

• 单点更新，区间查询

#define lowbit(x) (x)&(-x)
const int maxn=10005;
int sum[10005];
void add(int x,int val)
{
    while(x<maxn)
    {
        sum[x]+=val;
        x+=lowbit(x);
    }
}
int query(int x)
{
    int ans=0;
    while(x)
    {
        ans+=sum[x];
        x-=lowbit(x);
    }
    return ans;
}

• 区间更新,单点查询

#define lowbit(x) (x)&(-x)
const int maxn=10005;
int sum[10005];
void add(int x,int val)
{
    while(x<maxn)
    {
        sum[x]+=val;
        x+=lowbit(x);
    }
}

void update(int x,int y,int val)
{
    add(x,val);
    add(y+1,-val);
}

int query(int x)
{
    int ans=0;
    while(x)
    {
        ans+=sum[x];
        x-=lowbit(x);
    }
    return ans;
}

• 二维树状数组的区间修改单点查询

#define lowbit(x) (x)&(-x)
typedef long long ll;
ll sum[2505][2505];
int n,m,q;
void update(int x,int y,ll val)
{
    for(int i=x;i<=n;i+=lowbit(i))
        for(int j=y;j<=m;j+=lowbit(j))
        sum[i][j]+=val;
}
void update(int x1,int y1,int x2,int y2,ll val)
{
    update(x1,y1,val);
    update(x2+1,y1,-val);
    update(x1,y2+1,-val);
    update(x2+1,y2+1,val);
}
ll query(int x1,int y1)
{
    ll ans=0;
    for(int i=x1;i>0;i-=lowbit(i))
        for(int j=y1;j>0;j-=lowbit(j))
        ans+=sum[i][j];
    return ans;
}