洛谷P4514 上帝造题的七分钟(二维树状数组)

Solution

二维树状数组
(x1,y1) ( x 1 , y 1 ) (x2,y2) ( x 2 , y 2 ) 的和 sum=sum[x2][y2]sum[x11][y2]sum[x2][y11]+sum[x11][y11] s u m = s u m [ x 2 ] [ y 2 ] − s u m [ x 1 − 1 ] [ y 2 ] − s u m [ x 2 ] [ y 1 − 1 ] + s u m [ x 1 − 1 ] [ y 1 − 1 ]

sum[x][y]=i=1xj=1yk=1it=1ja[k][t]=i=1xj=1xa[i][j](xi+1)(yj+1) s u m [ x ] [ y ] = ∑ i = 1 x ∑ j = 1 y ∑ k = 1 i ∑ t = 1 j a [ k ] [ t ] = ∑ i = 1 x ∑ j = 1 x a [ i ] [ j ] ∗ ( x − i + 1 ) ∗ ( y − j + 1 )

sum[x][y]=i=1xj=1yk=1it=1ja[k][t]=i=1xj=1xa[i][j](xi+1)(yj+1)=i=1xj=1xa[i][j](xyxjyi++ij+x+yij+1)=i=1xj=1xa[i][j](xy+x+y+1)a[i][j]i(y+1)a[i][j]j(x+1)+a[i][j]ij s u m [ x ] [ y ] = ∑ i = 1 x ∑ j = 1 y ∑ k = 1 i ∑ t = 1 j a [ k ] [ t ] = ∑ i = 1 x ∑ j = 1 x a [ i ] [ j ] ∗ ( x − i + 1 ) ∗ ( y − j + 1 ) = ∑ i = 1 x ∑ j = 1 x a [ i ] [ j ] ∗ ( x y − x j − y i + + i j + x + y − i − j + 1 ) = ∑ i = 1 x ∑ j = 1 x a [ i ] [ j ] ∗ ( x y + x + y + 1 ) − a [ i ] [ j ] ∗ i ∗ ( y + 1 ) − a [ i ] [ j ] ∗ j ∗ ( x + 1 ) + a [ i ] [ j ] ∗ i ∗ j

因此我们开四个二维树状数组就可以解决了

Code

#include <cstdio>
#define N 2050
int n,m,a,b,c,d,k;
char str[10];
struct BIT{
    int tree[N][N];
    void add(int x,int y,int z){
        for(int i=x;i<=n;i+=i&(-i))
            for(int j=y;j<=m;j+=j&(-j))
                tree[i][j]+=z;
    }
    int query(int x,int y){
        int ans=0;
        for(int i=x;i;i-=i&-i)
            for(int j=y;j;j-=j&-j)
                ans+=tree[i][j];
        return ans;
    }
}A,Ai,Aj,Aij;
inline void add(int x,int y,int z){
    A.add(x,y,z);
    Ai.add(x,y,z*x);
    Aj.add(x,y,z*y);
    Aij.add(x,y,x*y*z);
} 
inline int query(int x,int y){
    return A.query(x,y)*(x*y+x+y+1)-Ai.query(x,y)*(y+1)-Aj.query(x,y)*(x+1)+Aij.query(x,y);
}
int main(){
    scanf("%s%d%d",str,&n,&m);
    while(scanf("%s%d%d%d%d",str,&a,&b,&c,&d)>0){
        if(str[0]=='L'){
            scanf("%d",&k);
            add(a,b,k);
            add(a,d+1,-k);
            add(c+1,b,-k);
            add(c+1,d+1,k);
        }else printf("%d\n",query(c,d)-query(a-1,d)-query(c,b-1)+query(a-1,b-1));
    }
    return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值