差分与差分矩阵

差分:通过对差分数组的加减运算,再通过前缀和的方式,从而达到对一个数组或者矩阵数组的部分区域进行加减。

数组的差分:

输入一个长度为 n 的整数序列。

接下来输入 m 个操作,每个操作包含三个整数 l,r,c,表示将序列中 [l,r] 之间的每个数加上 c。

请你输出进行完所有操作后的序列。

输入格式

第一行包含两个整数 n 和 m。

第二行包含 n 个整数,表示整数序列。

接下来 m 行,每行包含三个整数 l,r,c,表示一个操作。

输出格式

共一行,包含 n 个整数,表示最终序列。

数据范围

1≤n,m≤100000,
1≤l≤r≤n,
−1000≤c≤1000,
−1000≤整数序列中元素的值≤1000

输入样例:

6 3
1 2 2 1 2 1
1 3 1
3 5 1
1 6 1

输出样例:

3 4 5 3 4 2

通过对差分数组的加减来对原数组的部分区域进行加减。

 可以求出前缀和数组为: 2  3  3  1  2  1

和在原数组上直接进行暴力查找在相加的结果是相同的。

#include <bits/stdc++.h>

using namespace std;

const int N = 1e5;

int n,m;
int a[N],s[N];

int main ()
{
    int l,r,c;
    scanf("%d%d",&n,&m);
    for(int i = 1; i <= n; i++)
    {
        scanf("%d",&a[i]);//原始的数组
    }
    for(int i=1;i<=n;i++)
    {
        s[i]=a[i]-a[i-1];//差分数组
    }
    while(m--)
    {
        scanf("%d%d%d",&l,&r,&c);
        s[l] += c;
        s[r+1] -= c;
    }
    for(int i = 1; i <= n; i++)
    {
        a[i] = s[i]+a[i-1];
        cout<<a[i]<<" ";
    }
    return 0;
}

 上面的代码还有一种利用同一个函数同时进行差分数组的初始化,以及进行加减,可以看下面差分矩阵的问题。

差分矩阵:和一维数组的差分思路大致的相同。

输入一个 nn 行 mm 列的整数矩阵,再输入 qq 个操作,每个操作包含五个整数 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

 如图的事例,和一维数组的思路基本相同

 然后再进行前缀和的求解,代码如下:

#include <bits/stdc++.h>

using namespace std;

const int N = 1010;

int a[N][N],b[N][N];

void num (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 ()
{
    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 ++)
            num(i, j, i, j, a[i][j]);
            
    while(q --)
    {
        int x1, y1, x2, y2, c;
        cin >> x1>> y1 >> x2 >> y2 >> c;
        num (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[i][j];
    
    for(int i = 1; i <= n; i ++)
    {
        for(int j = 1; j <= m; j ++)
        {
               printf("%d ",b[i][j]);
        }
        printf("\n");
    }
    return 0;
}

 内容还有很多的不足,以后还会不断的补充充实。

 

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值