子矩阵的和（二维前缀和）题解【附图、举例、详细】

子矩阵的和（二维前缀和）

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

输入格式
第一行包含三个整数n，m，q。接下来n行，每行包含m个整数，表示整数矩阵。接下来q行，每行包含四个整数x1, y1, x2, y2，表示一组询问。

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

数据范围
1≤n,m≤10001≤n,m≤1000,
1≤q≤2000001≤q≤200000,
1≤x1≤x2≤n1≤x1≤x2≤n,
1≤y1≤y2≤m1≤y1≤y2≤m,
−1000≤矩阵内元素的值≤1000−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 点；

1.首先 初始化二维前缀和数组，也就是知道 s[i][j] 中所有数的和怎么算【通俗点来说，这个 s[i][j] 表示以 i 为长，j 为宽的小矩形，但是这个小矩形是以（1,1）为左上角，以（i，j）为右下角的】 。

2.其次，我们还要知道如何用 s[i][j]中所有数的和 来表示 以（x1,y1）为左上角 ，以 (x2,y2) 为右下角的子矩阵中所有数的和 。

【下面小白为了简约，就直言 s[i][j] 、（…）子矩阵了】。

下面，小白先说一下 求 s[i][j] 的 “背景知识” 。根据二维矩阵（数组）的特征，我们是一行计算完，在累加下一行某一列中的数的（假定有一个数组 a[3][2]，它就表示一个第 3 行、第 2 列的数，我们知道，计算机在进行 “扫描” 时，是先把第一行扫描完，再对下一行进行扫描的，所以这个 a[3][2] 可以表示从第 1 行、第 1 列这个数，也就是 a[1][1] 按照小白之前所说的顺序一直 “扫描” 到了 a[3][2] ，这就是之后我们在求 s[i][j] 时的 “背景知识” ）。

我们选取 n=5，m=5 的矩阵
首先，小白我来讲一下这个 s[i][j] 怎么算（ s[i][j] 对应 1个 i*j 的矩形，这个矩形即我们所要去表示的）【我们以 i=3，j=2 为例】。

所有喷涂颜色的部分表示 s[i][j] 。
s[i-1][j] 表示 浅蓝色 的部分：s[2][2] 。
s[i][j-1] 表示 紫色方框 的部分：s[3][1]。
s[i-1][j-1] 表示 红色方框 的部分：s[2][1]，这个部分在 2 和 3 中分别 +1 次，要 -1 次，避免重复多算。
a[i][j] 表示 最后未算的一个小方块 ：a[3][2] ，我们要把它加上。
最后，根据公式 “s[i][j]=s[i-1][j]+s[i][j-1]-s[i-1][j-1]+a[i][j]” ，我们就可以求得所有的 s[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];

前面，我们已经解决了 s[i][j] 如何和计算的问题，下面我们就可以 表示这个 “ 以（x1,y1）为左上角 ，以 (x2,y2) 为右下角的子矩阵 ” 了。

黄绿色背景 矩形区域表示 s[x2][y2]： 即 s[5][3]
粗体蓝色矩形 区域表示 s[x1-1][y2]： 即 s[2][3]
粗体紫色矩形 区域表示 s[x2][y1-1]： 即 s[5][1]
粗体黑色矩形 区域表示 s[x1-1][y1-1]： 即 s[2][1]，在2和3中共减了2次，多减1次，要+1次

这个我们最终要求的红色矩形即表示为 “s[x2][y2]-s[x2][y1-1]-s[x1-1][y2]+s[x1-1][y1-1]” ，即：
“s[5][3]-s[5][1]-s[2][3]+s[2][1]” 。

求 1 次这个红色子矩形即表示题目中的 1 次询问。

// 询问 
while(q--)
{
	int x1,y1,x2,y2;
	scanf("%d%d%d%d",&x1,&y1,&x2,&y2);
	printf("%d\n",s[x2][y2]-s[x2][y1-1]-s[x1-1][y2]+s[x1-1][y1-1]);
}

好的，最后小白为大家附上这道题目的 AC代码，我们此题的讲解就到此结束了，欢迎大家在下方评论留言哦，我们下次再会 ！

#include <bits/stdc++.h>
using namespace std;
const int N = 1010;

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

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