输入一个 n 行 m 列的整数矩阵,再输入 q 个操作,每个操作包含五个整数 x1,y1,x2,y2,c,其中 (x1,y1) 和 (x2,y2) 表示一个子矩阵的左上角坐标和右下角坐标。
每个操作都要将选中的子矩阵中的每个元素的值加上 c。
请你将进行完所有操作后的矩阵输出。
输入格式
第一行包含整数 n,m,q。
接下来 n 行,每行包含 m 个整数,表示整数矩阵。
接下来 q 行,每行包含 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
输入样例:
3 4 3
1 2 2 1
3 2 2 1
1 1 1 1
1 1 2 2 1
1 3 2 3 2
3 1 3 4 1
输出样例:
2 3 4 1
4 3 4 1
2 2 2 2
思路
本体可以看作是二位的差分,和一维的差分类似,通过反求出差分矩阵;然后通过对矩阵进行O(1)的操作就可以实现对原矩阵的的+c。
原理如下所示:
具体核心代码,这里需要对着图去理解
int x1, y1, x2, y2, c;
while (q --)
{
cin >> x1 >> y1 >> x2 >> y2 >> c;
B[x1][y1] += c;
B[x2+1][y2+1] += c;
B[x1][y2+1] -= c;
B[x2+1][y1] -= c;
}
代码示例
#include <iostream>
using namespace std;
const int N = 1010;
int n , m, q;
int A[N][N], B[N][N];
int main()
{
cin >> n >> m >> q;
for (int i = 1; i <= n; i ++)
{
for (int j = 1; j <= m; j ++)
{
cin >> A[i][j];
}
}
// 求差分数组B
for (int i = 1; i <= n; i ++)
{
for (int j = 1; j <= m; j ++)
{
B[i][j] = A[i][j] + A[i-1][j-1] - A[i][j-1] - A[i-1][j];
}
}
int x1, y1, x2, y2, c;
while (q --)
{
cin >> x1 >> y1 >> x2 >> y2 >> c;
B[x1][y1] += c;
B[x2+1][y2+1] += c;
B[x1][y2+1] -= c;
B[x2+1][y1] -= c;
}
// 处理完再次求A矩阵
for (int i = 1; i <= n; i ++)
{
for (int j = 1; j <= m; j ++)
{
A[i][j] = A[i-1][j] + A[i][j-1] - A[i-1][j-1] + B[i][j];
cout << A[i][j] << ' ';
}
cout << endl;
}
return 0;
}