【AcWing_算法基础课】子矩阵的和

原题链接： AcWing 796. 子矩阵的和

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

---

题目大意：

输入一个矩阵，再输入两个坐标，求以这两个坐标围城的矩阵中数的和

思路：

用a[i][j]存储这个位置的数，用s[i][j]存储前缀和（以a[1][1]和a[i][j]围成的矩阵中数的和）

s[i][j]如何求？

用递归
公式：s[i][j] = s[i-1][j] + s[i][j-1] - s[i-1][j-1] + a[i][j]
可以画图理解一下

以s[4][3] 为例：

• s[3][3] + s[4][2]导致s[3][2]被加了两次，需要减去一次。(下图的黑色部分)
• a[i][j]需要补上（下图红色部分）

如何求a[x1][y1]和a[x2][y2]围成的矩阵中的数之和？

上图是s[2][2]和s[3][5]围城的矩阵。

求出前缀和矩阵后，类似于上面的思想：
ans = s[x2][y2] - s[x1-1][y2] - s[x2][y1-1] + s[x1-1][y1-1]

其中s[x1-1][y1-1]部分减了两次，需要加上

---

C++代码：

#include <iostream>
using namespace std;

const int maxn = 1010;
int a[maxn][maxn], s[maxn][maxn];

int main(){
    int n, m, q;    //n行数 m列数 q询问数
    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 ++ )
            s[i][j] = s[i-1][j] + s[i][j-1] - s[i-1][j-1] + a[i][j];
        
    //q个询问        
    while(q--){
        int x1, y1, x2, y2;
        scanf("%d%d%d%d", &x1, &y1, &x2, &y2);
        
        printf("%d\n", s[x2][y2] - s[x1-1][y2] - s[x2][y1-1] + s[x1-1][y1-1]);
    }
    
    
    return 0;
}

python3代码：

if __name__ == '__main__':
    # n行数 m列数 q询问
    n, m, q = map(int, input().split())
    
    # 原始数组
    a = []
    for i in range(n):
        a.append(list(map(int, input().split())))
    
    # 前缀和数组
    s = [[0 for _ in range(m+1)] for _ in range(n+1)]
    for i in range(n):
        for j in range(m):
            s[i+1][j+1] = s[i][j+1] + s[i+1][j] - s[i][j] + a[i][j]
    
    # q个询问
    while q:
        x1, y1, x2, y2 = map(int,input().split())
        print(s[x2][y2]+s[x1-1][y1-1]-s[x2][y1-1]-s[x1-1][y2])
        q -= 1

复杂度分析：

• 时间复杂度：O(n²)，主要花在维护前缀和数组上
• 空间复杂度：O(n²)，维护一个二维的前缀和数组