差分

差分矩阵

原题链接

https://www.acwing.com/problem/content/800/

初始化的公式

差分跟前缀和的关系，差分是前缀和的逆运算

前缀和是将原本所有的数列都当成是一堆差分，然后对其进行求前缀和

差分是将原本所有的数列每一个数都当成一个前缀和，然后求他们的差分

而他们所使用的公式是一样的，只是所求的变量是不一样的，前缀和是根据原数列的数求前缀和，而现在是根据前缀和求原数列的数

差分矩阵初始化公式

由上图的前缀和公式可得，s[x][y] = (s[x][y-1] - s[x-1][y-1]) + s[x-1][y-1] + (s[x-1][y] - s[x-1][y-1]) + a[x][y]

数组s表示各个坐标的前缀和，而差分中原数列的每个坐标的值都表示差分的前缀和，所以，可以将原数列的数当成前缀和公式中的s数组，而a数组的数表示前缀和的数

所以a[x][y] = s[x][y] -(s[x][y-1] - s[x-1][y-1]) - s[x-1][y-1] - (s[x-1][y] - s[x-1][y-1])

从运算的式子上面来看，就是互逆运算。

差分数列初始化公式

应用的公式

快速对子矩阵进行增加或者减去相同的数

差分矩阵公式

给以(x1, y1)为左上角，(x2, y2)为右下角的子矩阵中的所有元素加上c：
b[x1][y1] += c;
b[x1][y2+1] -= c;
b[x2+1][y1] -= c;
b[x2+1][y2+1] += c;

数列差分公式

给区间[l, r]中的每个数加上c：

B[l] += c, B[r + 1] -= c

​

代码

#include<iostream>
using namespace std;

const int N = 1010;
int n,m,q;
int a[N][N],b[N][N];

void insert(int x1,int y1,int x2,int y2,int c)
{
    b[x1][y1] += c;
    b[x1][y2+1] -= c;
    b[x2+1][y1] -= c;
    b[x2+1][y2+1] += c;
}

int main()
{
    cin>>n>>m>>q;
    for(int i = 1; i <= n; i++)
        for(int j = 1; j <= m; j++)
            scanf("%d",&a[i][j]);
    //求差分矩阵
    //对当前位置加上a[i][j] 然后再对其后面的数的减去a[i][j] 相当于只对当前位置增加
    for(int i = 1; i <= n; i++)
        for(int j = 1; j <= m; j++)
            insert(i,j,i,j,a[i][j]);
    //对区间里面进行操作
    while(q--)
    {
        int x1,y1,x2,y2,c;
        cin>>x1>>y1>>x2>>y2>>c;
        insert(x1,y1,x2,y2,c);
    }
    //求差分数组的前缀和
    for(int i = 1; i <= n; i++)
        for(int j = 1; j <= m; j++)
            b[i][j] += b[i][j-1] + b[i-1][j] - b[i-1][j-1];
            
    //输出答案
    for(int i = 1; i <= n; i++)
    {
        for(int j = 1; j <= m; j++)cout<<b[i][j]<<' ';
        cout<<endl;
    }
    return 0;
}

        for(int j = 1; j <= m; j++)cout<<b[i][j]<<' ';
        cout<<endl;
    }
    return 0;
}