原题链接: 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(n2),主要花在维护前缀和数组上
- 空间复杂度:O(n2),维护一个二维的前缀和数组