AcWing796. 子矩阵的和
案列:
输入一个 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
思路:
在了解思路之前,我们应该先明白,什么是二维数组的前缀和,即二维数组的前n项和所指的是哪几项的和?
![](https://img-blog.csdnimg.cn/img_convert/c5b1784d33a56a3ecda2045b0f5591ee.jpeg)
如图中所示,二维数组的前缀和就是图中方形所覆盖的an的和,包括边界!!!!
那么我们如何利用递推公式写出前缀和呢?
如下图所示:
![](https://img-blog.csdnimg.cn/img_convert/fd9eae71fdd9f69aaf68efbdd7435c06.jpeg)
那么在理解了二维数组的前缀和的概念后,我们看下面这个问题:
我们如何计算这个紫色方形范围内的an的面积呢?
![](https://img-blog.csdnimg.cn/img_convert/528160ceebf86bfcd1aa8d63aa77f589.jpeg)
其中的s紫即可表示从[x1,y1]-[x2,y2]内子矩阵的和
图中所示的求法类似于我们高中所学的概率内容中的容斥原理。
即我们先减去两部分,然后再加上重复减去的部分。但是我们要时刻注意边界问题,图中的式子之所以减一,就是因为紫色方形的边界也要算到前缀和中。
代码:
#include<iostream>
using namespace std;
const int N=1010;
int a[N][N];
int S[N][N];
int main()
{
int 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];//求出(1,1)到(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]);//求出(x1,y1)到(x2,y2)的矩阵内元素和,容斥原理,注意·边界
}
return 0;
}