【ybtoj 4.2 树状数组课堂过关】F. 矩阵修改矩阵查询

题目链接：矩阵修改矩阵查询

分析

又是矩阵修改，参照D题，可以使用二维差分。
后面维护出现次数中出现的四个情况即可。

上代码

#include<iostream>
#include<cstdio>
#include<algorithm>
using namespace std;
typedef long long ll;

int n,m;
ll c[3001][3001],c1[3001][3001],c2[3001][3001],c3[3001][3001];

int lowbit(int x)
{
	return x&(-x);
}

void update(int x,int y,int k)
{
	for(int i=x;i<=n;i+=lowbit(i))
	{
		for(int j=y;j<=m;j+=lowbit(j))
		{
			c[i][j]+=k;
			c1[i][j]+=x*k;
			c2[i][j]+=y*k;
			c3[i][j]+=x*y*k;
		}
	}
}

ll getsum(int x,int y)
{
	ll sum=0;
	for(int i=x;i>0;i-=lowbit(i))
	{
		for(int j=y;j>0;j-=lowbit(j))
		{
			sum+=c[i][j]*(x+1)*(y+1)
			-(y+1)*c1[i][j]
			-(x+1)*c2[i][j]
			+c3[i][j];
		}
	}
	return sum;
}

int main()
{
    cin>>n>>m;
    int t;
	while(scanf("%d",&t)!=EOF)
	{
		int x1,x2,y1,y2;
		scanf("%d%d%d%d",&x1,&y1,&x2,&y2);
		if(t==1)
		{
			int k;
			scanf("%d",&k);
			update(x1,y1,k);
			update(x1,y2+1,-k);
			update(x2+1,y1,-k);
			update(x2+1,y2+1,k);
		}
		else
		{
			cout<<getsum(x2,y2)-getsum(x2,y1-1)-getsum(x1-1,y2)+getsum(x1-1,y1-1)<<endl;
		}
	} 
    return 0;
}