前往小程序,Get更优阅读体验!
立即前往
首页
学习
活动
专区
工具
TVP
发布
社区首页 >专栏 >【简单】差分矩阵

【简单】差分矩阵

作者头像
字节星球Henry
发布2021-08-09 16:49:56
1.1K0
发布2021-08-09 16:49:56
举报

输入一个 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 国际许可协议进行许可。

本文参与 腾讯云自媒体分享计划,分享自作者个人站点/博客。
如有侵权请联系 cloudcommunity@tencent.com 删除

本文分享自 作者个人站点/博客 前往查看

如有侵权,请联系 cloudcommunity@tencent.com 删除。

本文参与 腾讯云自媒体分享计划  ,欢迎热爱写作的你一起参与!

评论
登录后参与评论
0 条评论
热度
最新
推荐阅读
目录
  • 输入格式
  • 输出格式
  • 数据范围
  • 输入样例
  • 输出样例
  • 题解
  • C++ 代码
领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档