差分及其用法

一维差分

在差分当中我们要用到两个数组，原数组a[]和差分数组D[]

D[]的定义D[i]=a[i]-a[i-1],即原数组a[]的相邻元素差，从定义中也可以推出a[i]=D[1]+D[2]+D[3]+……+D[i],也就是说a[]是D[]的前缀和，即”差分是前缀和的逆运算"。

我们可以把每个D[]看成一条直线上的小线段，这些小线段相加就成为了a[]。

如何利用差分数组记录区间修改呢？为什么利用差分数组修改能提升修改的效率呢？

把区间[L,R]内每个元素 a[]+d，只需要把对应的D[]做如下操作

把D[L]加上d：D[L]+=d

把D[R+1]减去d:D[R+1]-=d

因为前缀和：a[x]=D[1]+D[2]+……+D[x],所以有:

1<=x<L , 前缀和a[x]不变；

L<=x<=R , 前缀和a[x]增加了d;

R<x<=N , 前缀和a[x]不变，因为被D[R+1]中减去的d抵消了；

 所以每次操作只需要修改[L,R]的两个端点的D[]值，复杂度是O(1)的。经过这种操作后，原来在a[]上做的复杂度为O(n)的区间修改操作，就变成了在D[]上做得复杂度为O(1)的端点操作，完成区间修改并得到D[]后，最后用D[]计算a[]，复杂度是O(n)的。m次区间修改和1次查询，总复杂度为O(n+m)。

二维差分

一维是线性数组，一个区间[[L,R]有两个端点；二维是矩阵，一个区间由四个端点围成。

a[][]是D[][]的前缀和，即由原点(1,1)和坐标(i,j)围成的矩阵中，所有的D[][]相加等于a[i][j]。

上面的图中，每个格子的面积是一个D[][]，例如阴影格子是D[i][j]，它由4个坐标点定义:（i-1,j)、(i,j)、(i,j-1)、(i-1,j-1) 。坐标点(i,j)的值是a[i][j]，它等于坐标(1,1)和(i,j)所围成的所有格子的总面积，

就等于a[i][j]。因此在二维情况下，差分变成了相邻的a[][]的"面积差",计算公式:D[i][j] =a[i][j]-a[i-1][j]-a[i][j-1]+a[i-1][j-1]。阴影方格表示 D[i][j]D[i][j] 的值，它的面积这样求：大面积 a[i][j] 减去两个小面积 a[i−1][j]、a[i][j−1]，由于两个小面积的公共面积 a[i−1][j−1] 被减了 22次，所以需要加来 1 次。

在一维情况下，做区间修改只需要修改区间两个端点的D[]值。而在二维情况下，一个区间是一个小矩阵，有4个端点，只需要修改这4个端点的D[][]值。例如坐标点(x1,y1)~(x2,y2)定义的区间，对应的四个端点。

D[x1][y1] += d; //二维区间的起点

D[x1][y2+1] -= d; //把x看成常数，y从y1到y2+1

D[x2+1][y1] -= d; //把y看成常数，x从x1到x2+1

D[x2+1][y2+1] += d; //由于前两式把d减了2次，多减了1次，这里加1次回来

 

如上图所示，两个黑色点围成的矩阵是题目给出的区间修改范围。我们只需要改变4个D[][]值，即改变图中4个阴影块的面积即可。

（ps:以上坐标点，说的都是一个面积，例如a[x1][y1]，它代表了左下角那个矩形面积）。

以acwing798 差分矩阵举例

输入一个 nn 行 mm 列的整数矩阵，再输入 qq 个操作，每个操作包含五个整数 x1,y1,x2,y2,cx1,y1,x2,y2,c，其中 (x1,y1)(x1,y1) 和 (x2,y2)(x2,y2) 表示一个子矩阵的左上角坐标和右下角坐标。

每个操作都要将选中的子矩阵中的每个元素的值加上 cc。

请你将进行完所有操作后的矩阵输出。

输入格式

第一行包含整数 n,m,qn,m,q。

接下来 nn 行，每行包含 mm 个整数，表示整数矩阵。

接下来 qq 行，每行包含 55 个整数 x1,y1,x2,y2,cx1,y1,x2,y2,c，表示一个操作。

输出格式

共 nn 行，每行 mm 个整数，表示所有操作进行完毕后的最终矩阵。

数据范围

1≤n,m≤10001≤n,m≤1000,
1≤q≤1000001≤q≤100000,
1≤x1≤x2≤n1≤x1≤x2≤n,
1≤y1≤y2≤m1≤y1≤y2≤m,
−1000≤c≤1000−1000≤c≤1000,
−1000≤矩阵内元素的值≤1000−1000≤矩阵内元素的值≤1000

输入样例：

3 4 3
1 2 2 1
3 2 2 1
1 1 1 1
1 1 2 2 1
1 3 2 3 2
3 1 3 4 1

输出样例：

2 3 4 1
4 3 4 1
2 2 2 2

#include <iostream>
using namespace std;
const int N = 1010;
int a[N][N], b[N][N];
int n, m, q;
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++)
            cin >> 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];//二维前缀和

    for (int i = 1; i <= n; i++)
    {
        for (int j = 1; j <= m; j++)
        {
            printf("%d ", b[i][j]);
        }
        printf("\n");
    }
    return 0;
}