非常漂亮的差分矩阵

#include<iostream>
using namespace std;

const int N = 1010;
int a[N][N],s[N][N];

int main()
{
    int n,m,q;
    cin>>n>>m>>q;
    for(int i=1; i<=n; i++)
        for(int j=1; j<=m; j++){     //下标从1开始
        cin>>s[i][j];
        a[i][j] = s[i][j]-s[i-1][j]-s[i][j-1]+s[i-1][j-1];     //构造差分数组
    }
    while(q--){
        int x1,x2,y1,y2,c;
        cin>>x1>>y1>>x2>>y2>>c;
        a[x1][y1]+=c;
        a[x1][y2+1]-=c;
        a[x2+1][y1]-=c;
        a[x2+1][y2+1]+=c;       //使(x1,y1)到(x2,y2)之间的所有数都加上c
    }
    for(int i=1; i<=n; i++){
        for(int j=1; j<=m; j++){
            a[i][j]+=a[i-1][j]+a[i][j-1]-a[i-1][j-1];
            cout<<a[i][j]<<' ';     //区间和计算(注意数据输入时是用数组s表示，而输出用a表示)
        }
            cout<<endl;
    }


    return 0;
}