C++ 一维差分

C++一维差分算法详解

一维差分

C++中的差分是一种非常实用的算法技巧,主要用于高效处理区间更新操作。

1. 差分的基本概念

差分是前缀和的逆运算,核心思想是:

  • 对原数组的区间操作可以转化为对差分数组的单点操作
  • 将O(n)的区间更新优化为O(1)的单点更新

(1)差分数组:

首先给定一个原数组a:a[1], a[2], a[3], a[n];

然后我们构造一个数组b : b[1], b[2], b[3], b[i];

使得 a[i] = b[1] + b[2] + b[3] + , + b[i]

也就是说,a数组是b数组的前缀和数组,反过来我们把b数组叫做a数组的差分数组。换句话说,每一个a[i]都是b数组中从头开始的一段区间和。

(2)考虑如何构造差分b数组?

最为直接的方法
如下:
a[0 ]= 0;
b[1] = a[1] - a[0];
b[2] = a[2] - a[1];
b[3] = a [3] - a[2];

b[n] = a[n] - a[n - 1];

(3)知道了差分数组有什么用呢?

在这里插入图片描述
//diff差分数组 arr 原数组,也是前缀和数组

话说有这么一个问题:(对差分写法的解释说明)
给定区间[l, r ],让我们把a数组中的[l, r] 区间中的每一个数都加上c,即 a[l] + c , a[l + 1] + c , a[l + 2] + c , a[r] + c;
暴力做法是for循环l到r区间,时间复杂度O(n),如果我们需要对原数组执行m次这样的操作,时间复杂度就会变成O(n * m)。有没有更高效的做法吗? 考虑差分做法,(差分数组派上用场了)。

始终要记得,a数组是b数组的前缀和数组,比如对b数组的b[i]的修改,会影响到a数组中从a[i]及往后的每一个数。
首先让差分b数组中的 b[l] + c ,通过前缀和运算,a数组变成 a[l] + c ,a[l + 1] + c, a[n] + c;

然后我们打个补丁,b[r + 1] - c, 通过前缀和运算,a数组变成 a[r + 1] - c,a[r + 2] - c,a[n] - c;

(4)为什么只需要改变起点和终点??

直观理解:水波效应

想象你在一个平静的湖面(全0数组)上:

初始: [0, 0, 0, 0, 0, 0]
场景1:在位置1扔石头(产生波纹)
diff[1] += 3;// 效果:从位置1开始,后面所有元素都+3
恢复后: [0, 3, 3, 3, 3, 3]

场景2:在位置4设置"抵消点"
diff[1] += 3;
diff[4] -= 3; // 抵消点!// 恢复后: [0, 3, 3, 3, 0, 0]
终点后的抵消让增加效果只在 [1,3] 区间内生效!

逐步演示
我们要实现:[1,3] 区间加2
步骤1:起点加2
在这里插入图片描述

步骤2:终点+1处减2(关键!)
在这里插入图片描述

完美!

二维差分是前缀和运算的逆过程,它可用于高效处理矩阵区域的修改操作。以下从原理和实现方法两方面详细介绍: ### 原理 差分算法是前缀和运算的逆过程,它通过计算原序列中相邻元素的差值,得到一个差分序列。二维差分数组 `diff` 可表示相邻元素之间的差值关系,其构建规则为:`diff[i][j] = matrix[i][j] - matrix[i-1][j] - matrix[i][j-1] + matrix[i-1][j-1]`。通过对差分序列进行逐个累加,可以得到原序列的前缀和序列。当对矩阵中以 `(x1, y1)` 为左上角,`(x2, y2)` 为右下角的子矩阵中的数值都加 `c` 时,只需对差分数组进行如下操作:`b[x1][y1]+=c`,`b[x1][y2+1]-=c`,`b[x2+1][y1]-=c`,`b[x2+1][y2+1]+=c`。最后通过累加差分数组,就能得到更新后的原矩阵 [^1][^2][^3]。 ### 实现方法 以下是一个 C++ 实现二维差分的示例代码: ```cpp #include <iostream> using namespace std; int n, m, q; const int N = 1010; int a[N][N], b[N][N]; // 插入函数,用于修改差分数组 void insert(int x1, int y1, int x2, int y2, int c) { b[x1][y1] += c; b[x1][y2 + 1] -= c; b[x2 + 1][y1] -= 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]); // 进行操作 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]; // 输出更新后的矩阵 for (int i = 1; i <= n; i++) { for (int j = 1; j <= m; j++) printf("%d ", b[i][j]); puts(""); } return 0; } ``` 在上述代码中,`insert` 函数用于修改差分数组,`main` 函数完成输入、初始化差分数组、进行操作、求前缀和以及输出结果的操作 [^5]。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值