算法基础集训（第13天，共106天）------＞带你彻底弄懂【差分】之【差分矩阵】：许多高校的的比赛甚至面试题都有考过手写代码题目+分析+万年模板应有尽有

本文链接：https://blog.csdn.net/mary19831/article/details/126404806

一：二维差分

从一维扩展到二维，将某个子矩阵的数值都加上c，并且输出最后的二维矩阵，也同样需要构造差分二维矩阵，就可以达到O（1）的时间复杂度

二：定义

我们定义：a【】数组是b【】数组的前缀和，也就是a【i】【j】的值是以b【i】【j】为右下角端点的矩阵的的和

三：题目描述

四：核心函数

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

对b数组执行以上操作等价于：

for( int i=x1; i<= x2; i ++)
for( int j=y1;j <= y2; j ++)
a【i】【j】+=c;

我们通过下图来理解这个插入过程

b[x1][ y1 ] +=c ; 对应下方图1 ,让整个a数组中蓝色矩形面积的元素都加上了c。
b[x1,][y2+1]-=c ; 对应下方图2 ,让整个a数组中绿色矩形面积的元素都减去c，使其内元素不发生改变。
b[x2+1][y1]- =c ; 对应下方图3 ,让整个a数组中紫色矩形面积的元素都减去c，使其内元素不发生改变。
b[x2+1][y2+1]+=c; 对应下方图4,,让整个a数组中红色矩形面积的元素都加上c，红色内的相当于被减了两次，再加上一次c，才能保持不被影响。

五：构造思路

1.首先是输入二维数组a【】

2.构造差分数组b【】；（这里有一个思维难点，我们不需要搞清楚b数组的具体构造过程，我们只需要假设a数组为空但是其实已经存在，我们每次在【i】【j】的位置插a【i】【j】那么由于a数组是b数组的前缀和，b数组会自行构造完成）

3.在输入所需要操作的子矩阵范围和数字c

4.对数组b进行前缀和

5.输出

注：我们容易发现，差分矩阵是前缀和矩阵的逆过程；

a数组在此题目中没有什么意义，单纯是为了构造初次差分矩阵b而存在的

六：万年无误代码模板

#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[x2 + 1][y1] -= c;
    b[x1][y2 + 1] -= c;
    b[x2 + 1][y2 + 1] += c;
}

int main()
{
    scanf("%d%d%d", &n, &m, &q);

    for (int i = 1; i <= n; i ++ )
        for (int j = 1; j <= m; j ++ )
            scanf("%d", &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 - 1][j] + b[i][j - 1] - b[i - 1][j - 1];//二维前缀和
            printf("%d ", b[i][j]);
            if(j == m)
            puts("");
        }

    return 0;
}