前缀和算法

序言

 在这篇文章中，我们将认识一维前缀和以及二维前缀和，前缀和主要用于数据结构（如数组）的区间求和问题。

1. 一维前缀和

1. 题目描述

输入一个长度为 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

---

2. 暴力解法

#include <iostream>
#include <vector>
using namespace std;


const int N = 100010;
int nums[N];

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

遍历每一个查询区间，然后输出结果，显而易见时间超时了！😭

2. 构造前缀和数组

各位请问还记得高中的单调数列求和吗？

Sn = a1 + a2 + ... + a3 + an-1 + an

在这里，如果我想要得到某一段特定区间的求和的话，就比如第 3 个 到 第 5 个：

a3 + a4 + a5 = S5 - S2
S5 = a1 + a2 + a3 + a4 + a5
S2 = a1 + a2

所以我们可以得出任意一段区间的和可以表示为：

al + ... + ar = Sr - S(l - 1)

我们浪费一个空间，下标从 1 开始较直观，所以我们可以表示为：

#include <iostream>
#include <vector>
using namespace std;


const int N = 100010;
int a[N];
int s[N];

int main(){
    int m, n;
    cin >> m >> n;
    
    // 边记录数据边构造前缀和数组
    for(int i = 1; i <= m; ++i){
        cin >> a[i];
        s[i] = s[i - 1] + a[i];
    }
    
    while(n--){
        int l, r;
        cin >> l >> r;
        
        cout << s[r] - s[l - 1] << endl;
    }
    
    
    return 0;
}

甚至我们还可以优化一下，其实在这里我们需要只是前缀和数组，另一个数组可以舍去：

const int N = 100010;
// 构造一个数组就够了
int s[N];

int main(){
    int m, n;
    cin >> m >> n;
    
    // 边记录数据边构造前缀和数组
    int temp; // 记录当前的元素
    for(int i = 1; i <= m; ++i){
        cin >> temp;
        s[i] = s[i - 1] + temp;
    }
....

---

2. 二维前缀和

1. 题目描述

输入一个 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

2. 构造前缀和数组

是的，这道题我们也需要构造前缀和数组，但是怎么构造呢？
我们设 S[i][j] 代表从 S[0][0] - S[i][j] 构造的矩阵，就比如我们表示一下 S[2][2]

那怎么随着元素的增加，来更新前缀和数组（这里代表矩阵的面积），公式我直接给大家了，大家可以照着图推导一下：

S[i][j] = S[i - 1][j] + S[i][j - 1] - S[i - 1][j - 1] + a[i][j]; // 更新规则

在一维前缀和中，我们获取区间的值是使用两个下标对应的前缀和直接相减，但在二维中我们需要求的是指定矩形的面积，这里就会稍显的复杂一下，比如：

在这里我们需要求得以 x1 = 2, y1 = 2 为左上角，x2 = 3, y2 = 2 为右下角的矩形的面积，这也是有固定方法求的，但是稍显复杂，公式为：

area = S[x2][y2] - S[x1 - 1][y2] - S[x2][y1 - 1] + S[x1 - 1][y1 - 1]

大家也是可以通过带入值来验证一下。

所以说整道题的思路就出来啦，先求前缀和数组（这里是每个原点到该坐标的矩形的面积），再根据公式求具体的两个坐标之间的矩形，代码得出是:

#include <iostream>
using namespace std;

const int N = 1010;
int arr[N][N];

int main(){
    int n, m, q;
    cin >> n >> m >> q;
    
    // 初始化前缀和数组
    int temp;
    for(int i = 1; i <= n; i++){
        for(int j = 1; j <= m; j++){
            cin >> temp;
            arr[i][j] = arr[i - 1][j] + arr[i][j - 1] - arr[i - 1][j - 1] + temp; // 更新规则
        }
    }
    
    while(q--){
        int x1, y1, x2, y2;
        cin >> x1 >> y1 >> x2 >> y2;
        
        cout << arr[x2][y2] - arr[x1 - 1][y2] - arr[x2][y1 - 1] + arr[x1 - 1][y1 - 1] << endl;
    }
    
    return 0;
}