51nod1210

题解：

二维树状数组，再矩阵推一下

代码：

#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
const int N=;
int n,m,q;
LL T[][N][N];
int TN,TM;
void add(LL (*T)[N],int x,int y,LL v)
{
    if (!x||!y)return;
    for (int i=x;i<=TN;i+=i&-i)
     for (int j=y;j<=TM;j+=j&-j)T[i][j]+=v;
}
LL sum(LL (*T)[N],int x,int y)
{
    if (!x||!y)return 0LL;
    LL v=;
    for (int i=x;i;i-=i&-i)
     for (int j=y;j;j-=j&-j)v+=T[i][j];
    return v;
}
void add(int x,int y,LL v)
{
    add(T[],x,y,v);
    add(T[],x,y,v*(-x));
    add(T[],x,y,v*(-y));
    add(T[],x,y,v*(-x)*(-y));
}
LL sum(int x, int y)
{
    return (x+)*(y+)*sum(T[],x,y)+(y+)*sum(T[],x,y
    )+(x+)*sum(T[],x,y)+sum(T[],x,y);
}
void add(int x1,int y1,int x2,int y2,LL v)
{
    add(x1,y1,v);
    add(x2+,y1,-v);
    add(x1,y2+,-v);
    add(x2+,y2+,v);
}
LL sum(int x1, int y1, int x2, int y2)
{
    return sum(x2,y2)-sum(x1-,y2)-sum(x2,y1-)+sum(x1-,y1-);
}
int main()
{
    scanf("%d%d%d",&n,&m,&q);
    TN=n;TM=m;
    for (int i=;i<=m;i++)
     for (int j=;j<=n;j++)
      {
        int x;
        scanf("%d",&x);
        add(j,i,j,i,x);
      }
    for (int i=;i<q;i++)
     {
        int x1,y1,x2,y2,z,a;
        scanf("%d%d%d%d%d%d",&x1,&y1,&x2,&y2,&z,&a);
        LL s = sum(x1, y1, x2, y2);
        if (s<(LL)z*(x2-x1+)*(y2-y1+))add(x1,y1,x2,y2,a);
     }
    for (int i=;i<=m;i++)
     for (int j=;j<=n;j++)
      printf("%lld%c",sum(j,i,j,i),j==n?'\n':' ');
}