首页
学习
活动
专区
圈层
工具
发布

【简单】差分矩阵

输入一个 nm 列的整数矩阵,再输入 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 整数,表示所有操作进行完毕后的最终矩阵。

数据范围

{\rm{1}} \le {\rm{q}} \le {\rm{100000}}

输入样例

代码语言:javascript
复制
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

输出样例

代码语言:javascript
复制
2 3 4 1
4 3 4 1
2 2 2 2

题解

(二维差分)O(1)

二维前缀和(基础扩展) O(1)

求前缀和,以 S_{i,j} 为例(a_{i,j} 为当前位置的值):

S_{i,j}=S_{i,j-1}+S_{i-1,j}-S_{i-1,j-1}+a_{i,j}

一个左上角和右下角坐标为 (x1,y1),(x2,y2) 子矩阵的和可认为:

S_{x_2,y_2}-S_{x_2,y_1-1}-S_{x_1-1,y_2}+S_{x_1-1,y_2-1}

S_{x_2,y_2} 减去两个矩形再加上两个矩形重叠部分多减去的一次。

二维差分(即前缀和的逆运算)O(1)

构造 b 使得 ab 数组的前缀和,即 ba 的差分:

a_{i,j}=b_{1,1}+b_{1,2}+\ldots +b_{2,1}+b_{2,2}+\ldots+b_{i,j}

具体到此题,要使得 a 中间的子矩阵全部加上 c ,即是让其差分 b_{x_1,y_1} 加上 c ,此时,该坐标之后的矩阵(b 的前缀和子矩阵)全部加上 c ,也就多加了一个倒“L”型的区域,将该区域减去即可完成题目操作。详见 insert() 函数。

C++ 代码

代码语言:javascript
复制
#include <iostream>
#include <algorithm>

using namespace std;

const int N = 1010;
int n, m, q;
int a[N][N], b[N][N];

void insert(int x1, int y1, int x2, int y2, int c)
{
    b[x1][y1] += c;
    b[x2 + 1][y1] -= c;
    b[x1][y2 + 1] -= c;
    b[x2 + 1][y2 + 1] += c;
}

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++)
            insert(i, j, i, j, a[i][j]);//将读入的矩阵构造差分更新到b中
    while(q--)
    {
        int x1, y1, x2, y2, c;
        cin >> x1 >> y1 >> x2 >> y2 >> c;
        insert(x1, y1, x2, y2, c);
    }
    for(int i = 1; i <= n; i++)
        for(int j = 1; j <= m; j++)
            b[i][j] += b[i - 1][j] + b[i][j - 1] - b[i - 1][j - 1];//求二维差分矩阵b的前缀和
    for(int i = 1; i <= n; i++)
    {
        for(int j = 1; j <= m; j++)
            printf("%d ", b[i][j]);
        cout << endl;
    }
    return 0;
}

本作品采用知识共享署名-非商业性使用-相同方式共享 4.0 国际许可协议进行许可。

下一篇
举报
领券