理解前缀和与差分：高效计算数组前缀和与差分技巧

介绍： 在计算机科学和算法中，前缀和（PrefixSum）是一种常见而强大的技巧。它能够高效地计算数组或序列中某个位置之前的所有元素的总和。本文将介绍前缀和的概念、计算方法和应用场景，帮助读者理解并掌握这一重要的计算技巧。

一、什么是前缀和与差分呢？

前缀和是指数组或序列中每个位置之前的元素总和。给定一个长度为n的数组A，前缀和P的第i个元素表示A中前i个元素的总和。通常，P[0]设置为0，表示前0个元素的总和为0。通过计算前缀和，我们可以快速获得任意位置之前的元素总和。
而差分可以看成前缀和的逆运算。
合理的使用前缀和与差分，可以将某些复杂的问题简单化。

二、计算前缀和的方法

有多种方法可以计算前缀和，其中最简单直观的方法是使用循环遍历数组，累加每个元素到前一个元素的和中。然而，这种方法的时间复杂度为O(n)，并不是最优解。下面介绍一种更高效的方法。

三、动态规划求前缀和

1、一维前缀和

使用动态规划的思想，我们可以利用一个辅助数组来存储前缀和。假设我们有一个长度为n的数组A和一个长度为n+1的辅助数组P。以下是计算一维前缀和的步骤：

模板：
1、初始化P[0] = 0。
2、对于i从1到n，计算P[i] = P[i-1] + A[i-1]。

通过以上步骤，我们可以在一次遍历中计算出所有前缀和，而时间复杂度为O(n)。这种方法非常高效，特别适用于需要频繁计算前缀和的场景。

以下是一个示例代码，实现对一维数组的前缀和计算：

#include<iostream>

using namespace std;

int n , sum[1<<20];//1右移20位大致等于 1,048,576

int main()
{
	cin >> n;

	for(int i = 1 ; i <= n ; i++)
	{
		cin >> arr[i];//读入数组 

		sum[i] = sim[i-1] + a[i];//进行操作 
	}

	for(int i = 1 ; i <= n ; i++)
	{
		cout<<sum[i];//输出前缀和 
	}

	return 0;
}

2、二维前缀和

在已经掌握一维前缀和的基础上，我们接下来可以思考在二维的平面上的问题了。在一维前缀和的问题中，我们讨论的是某一个区间内的元素之和，那如果放在一个二维平面，我们要的就是某个矩阵的和。

模板：
1、 s[i][j] = s[i - 1][j] + s[i][j - 1] - s[i - 1][j - 1] + a[i][j];
2、s[x2][y2] - s[x1 - 1][y2] - s[x2][y1 - 1] + s[x1 - 1][y1 - 1])

图解：

如果要算绿色矩形的面积，就需要用 整体的面积 减去 黑加蓝色 和 黑加橙色的面积，这样的话黑色的面积就会被剪掉两次，再加上一个黑色的面积解释绿色的面积

以下是一个示例代码，实现对二维数组的前缀和计算：

#include <iostream>

using namespace std;

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

int main () {
    cin >> n >> m >> q;
    
    for(int i = 1; i <= n; i ++) {
        for(int j = 1; j <= m; j ++) {
            cin >> a[i][j];
            s[i][j] = s[i - 1][j] + s[i][j - 1] - s[i - 1][j - 1] + a[i][j];
        }
    }
    
    while (q -- ) {
        int x1, x2, y1, y2;
        cin >> x1 >> y1 >> x2 >> y2;    
        cout << s[x2][y2] -  s[x1 - 1][y2] -  s[x2][y1 - 1] + s[x1 - 1][y1 - 1] << endl;
    }
    return 0;
}

四、应用场景

前缀和在许多算法和数据结构问题中都有广泛应用。以下是一些常见的应用场景：

1、子数组和的计算：通过前缀和，可以快速计算任意子数组的和，从而解决一系列相关问题，如最大子数组和、最小子数组和等。

2、区间和的查询：如果需要频繁查询某个区间的和，可以利用前缀和提前计算出所有区间的和，并存储在辅助数组中，以实现快速查询。

3、数组元素更新：当数组中的元素需要频繁更新时，通过前缀和可以减少更新的时间复杂度，提高算法效率。

例题1 一维前缀和

原题链接：https://www.acwing.com/problem/content/797/

题面： 输入一个长度为 n的整数序列。 接下来再输入 m个询问，每个询问输入一对 l,r。 对于每个询问，输出原序列中从第 l个数到第
r个数的和。

输入格式 第一行包含两个整数 n和 m。 第二行包含 n个整数，表示整数数列。 接下来 m行，每行包含两个整数 l 和
r，表示一个询问的区间范围。

输出格式 共 m行，每行输出一个询问的结果。

数据范围 1≤l≤r≤n, 1≤n,m≤100000, −1000≤数列中元素的值≤1000

输入样例：

5 3 2 1 3 6 4 1 2 1 3 2 4

输出样例：

3 6 10

代码：

#include<bits/stdc++.h>

using namespace std;

const int N = 100050;

int n,l,r,m;
int a[N],b[N];

int main()
{
	cin>>n>>m;
	for(int i=1;i<=n;i++)
	{
		cin>>a[i];
		b[i] = b[i-1] + a[i];
	}
	while(m--)
	{
		cin>>l>>r;
		cout<<b[r] - b[l - 1]<<endl; 
	}
	return 0;
}

例题2 二维前缀和

原题链接：https://www.acwing.com/problem/content/798/

输入一个 n 行 m 列的整数矩阵，再输入 q 个询问，每个询问包含四个整数 x1,y1,x2,y2
，表示一个子矩阵的左上角坐标和右下角坐标。

对于每个询问输出子矩阵中所有数的和。
输入格式

第一行包含三个整数 n，m，q

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

接下来 q 行，每行包含四个整数 x1,y1,x2,y2，表示一组询问。
输出格式

共 q 行，每行输出一个询问的结果。
数据范围

1≤n,m≤1000 , 1≤q≤200000, 1≤x1≤x2≤n, 1≤y1≤y2≤m, −1000≤矩阵内元素的值≤1000

输入样例：

3 4 3
1 7 2 4
3 6 2 8
2 1 2 3
1 1 2 2
2 1 3 4
1 3 3 4

输出样例：

17 27 21

代码：

#include <iostream>
using namespace std;
const int N = 1010;
int n,m,q;
int a[N][N],sum[N][N];
int get (int x1,int y1,int x2,int y2) {
    return sum[x2][y2] - sum[x1 - 1][y2] - sum[x2][y1 - 1] + sum[x1 - 1][y1 - 1];//结合图来理解
}
int main () {
    cin >> n >> m >> q;
    for (int i = 1;i <= n;i++) {
        for (int j = 1;j <= m;j++) {
            cin >> a[i][j];
            sum[i][j] = sum[i][j - 1] + sum[i - 1][j] - sum[i - 1][j - 1] + a[i][j];//结合图来理解
        }
    }
    while (q--) {
        int x1,y1,x2,y2;
        cin >> x1 >> y1 >> x2 >> y2;
        cout << get (x1,y1,x2,y2) << endl;
    }
    return 0;
}

五、差分

1、什么是差分？

差分就是将数列中的每一项分别与前一项数做差。

一维差分

创建一数组b，使得数组a为数组b的前缀和，数组b为数组a的差分

构造方法：b[i] = a[i] - a[i - 1]

此处使用了一个虚拟的构造方式(在数组一个位置加上一个数，那么在它的下一个位置减去这一数)

应用：对于a数组的任意区间[l, r]，令其加上一个数，而不改变其它值

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

差分操作和前缀和一样数组下标都从1开始。

二维差分

理解二维前缀和的思路来去理解二维差分（图见上上上边）

构造方法：b[i][j] += c, b[i + 1][j] -= c, b[i][j + 1] -= c, b[i + 1][j + 1]+= c

2、例题

原题链接：https://www.acwing.com/problem/content/description/800/

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

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

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

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

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

接下来 q 行，每行包含 5 个整数 x1,y1,x2,y2,c，表示一个操作。
输出格式

共 n 行，每行 m个整数，表示所有操作进行完毕后的最终矩阵。
数据范围

1≤n,m≤1000 ,
1≤q≤100000,
1≤x1≤x2≤n,
1≤y1≤y2≤m,
−1000≤c≤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>
#include <algorithm>

using namespace std;

const int N = 1010;

int n, m, q;
int a[N][N], s[N][N];

int main()
{
    cin >> n >> m >> q;
    for (int i = 1; i <= n; i ++ )
        for (int j = 1; j <= m; j ++ )
            scanf("%d", &s[i][j]);

    for (int i = 1; i <= n; i ++ )
        for (int j = 1; j <= m; j ++ )
            a[i][j] = s[i][j] - s[i - 1][j] - s[i][j - 1] + s[i - 1][j - 1];

    while (q -- )
    {
        int x1, y1, x2, y2, c;
        scanf("%d%d%d%d%d", &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;
    }

    for (int i = 1; i <= n; i ++ )
        for (int j = 1; j <= m; j ++ )
            s[i][j] = a[i][j] + s[i - 1][j] + s[i][j - 1] - s[i - 1][j - 1];

    for (int i = 1; i <= n; i ++ )
    {
        for (int j = 1; j <= m; j ++ ) printf("%d ", s[i][j]);
        cout << endl;
    }

    return 0;
}

六、结论

前缀和是一种强大而实用的计算技巧，能够高效地计算数组或序列中某个位置之前的所有元素的总和。通过掌握前缀和与差分的概念、计算方法和应用场景，我们可以在算法和数据处理中发挥其优势，提高代码的效率和性能。

希望本文能帮助读者更好地理解和应用前缀和与差分，为解决实际问题提供有力的工具和思路。

如果有错误请欢迎指正~~，感谢！！！