假设b[N][N]是我们的差分数组:
二维矩阵的差分主要运用于:在二维数组内的某一个区间内加上或者减去一个数。例如要在以x1,y1,x2,y2为对角线的的一个小矩阵内加上一个数c的差分公式为:
b[x1][y1]+=c;
b[x2+1][y1]-=c;
b[x1][y2+1]-=c;
b[x2+1][y2+1]+=c;
然后再对该数组进行求前缀和的操作即可完成在二维数组内部的一个小矩阵内加上一个数的操作。
求二维前缀和的操作为:
b[i][j]=b[i][j-1]+b[i-1][j]-b[i-1][j-1]+b[i][j];
下面是具体的题目例子:
输入一个 n 行 m 列的整数矩阵,再输入 w 个操作,每个操作包含五个整数 x1,y1,x2,y2,c,其中 (x1,y1)(x1,y1) 和 (x2,y2)(x2,y2) 表示一个子矩阵的左上角坐标和右下角坐标。
每个操作都要将选中的子矩阵中的每个元素的值加上 c。
请你将进行完所有操作后的矩阵输出。
输入格式
第一行包含整数 n,m,q。
接下来 n 行,每行包含 m 个整数,表示整数矩阵。
接下来 w 行,每行包含 5 个整数 x1,y1,x2,y2,c,表示一个操作。
输出格式
共 n 行,每行 m 个整数,表示所有操作进行完毕后的最终矩阵。
数据范围
1≤n,m≤1000,
1≤q≤100000,
1≤x1≤x2≤n,
1≤y1≤y2≤m,
−1000≤c≤1000,
−1000≤矩阵内元素的值≤1000
代码:
#include <iostream>
#include <cstring>
#include <algorithm>
using namespace std;
const int N=1010;
int a[N][N];
int b[N][N];
int n,m,q;
int main()
{
cin>>n>>m>>q;
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
scanf("%d",&a[i][j]);
while(q--)
{
int x1,y1,x2,y2,c;
scanf("%d%d%d%d%d",&x1,&y1,&x2,&y2,&c);
b[x1][y1]+=c;
b[x2+1][y1]-=c;
b[x1][y2+1]-=c;
b[x2+1][y2+1]+=c;
}
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
b[i][j]=b[i][j-1]+b[i-1][j]-b[i-1][j-1]+b[i][j];
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
a[i][j]+=b[i][j];
for(int i=1;i<=n;i++)
{
for(int j=1;j<=m;j++) printf("%d ",a[i][j]);
puts(" ");
}
return 0;