ET算法【前缀和与差分】

已于 2023-11-05 22:24:04 修改
阅读量84 收藏 1
点赞数 2
分类专栏： ET基础算法 文章标签： 算法
于 2023-10-17 23:02:27 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_39530733/article/details/133895822
版权

最近比较忙其实根本不知道在忙啥，好久没写代码了，今天在熄灯前20min水一篇博客估计要鸽到明天咯，近期刷洛谷题单有些刷不动了，感觉是因为那些题需要会的东西有点多，看题解也越来越看不懂了，还得找个课来学学新算法，等a协的培训又有些太慢了马上就要打新生赛了哇。这不就来了嘛，挑了挑还是这个便宜，内容好像大概有一半已经学过了？不懂，但我觉得还是选这个比较好，嗯，我觉得ok就ok，难得的freedom。

一、序

前缀和与差分，主要运用于多次求大区间范围内的和或一起调整元素大小，基础分为一维与二维

diff     --前缀和-->   a   --前缀和-->   prefix                             diff     <--差分--   a   <--差分--   prefix

注意a[ 0 ]不用，值设为0，从a[ 1 ]开始

二、一维

1.一维前缀和

P1049 - 【模板】前缀和 - ETOJ (eriktse.com)

数列总学过吧，an总会输入吧，Sn总会算吧。啊没错，前缀和就是Sn。

核心式：prefix[ i ]=prefix[ i-1 ]+a[ i ]

完整代码如下：

#include<bits/stdc++.h>
using namespace std;
using ll = long long;
#define qio ios::sync_with_stdio(false), cin.tie(0), cout.tie(0);

const int N = 1e5 + 5;
ll a[N], prefix[N];       //考虑数据范围，Sn应该会超int的
int main() {
	qio
	int n;
	cin >> n;
	for (int i = 1; i <= n; i++)cin >> a[i];
	for (int i = 1; i <= n; i++)prefix[i] = prefix[i - 1] + a[i];
	int q, l, r;
	cin >> q;
	for (int i = 1; i <= q; i++) {
		cin >> l >> r;
		cout << prefix[r] - prefix[l-1] << "\n";  //注意l-1，第一次就漏了
	}
}

2.一维差分

P1050 - 【模板】差分 - ETOJ (eriktse.com)

差分就是前后两个数的差值，高中好像也经常用吧（？），差分数组的求和就是a[ i ]的值

核心式：diff[ i ] = a[ i ] - a[ i - 1 ]

手敲图解：

       啊这，还没敲就熄灯咯，阿巴阿巴，真鸽了啊，还不知道啥时候有时间接着写呢，刚开学一个月就要期中了，就离谱，咱也不知道什么时候有时间接着写，咕咕咕    接着写

    a0        a1        a2        a3        a4        a5        a6

                d1        d1        d1        d1        d1        d1

                            d2        d2        d2        d2        d2

                                        d3        d3        d3        d3

                                                    d4        d4        d4

                                                                d5        d5

                                                                            d6

想得到a的值，只需要把a下面的diff全都加起来就成（根据diff的定义），并且可以快速的区间修改，例如要给a2到a4这一区间加上x，只需让d2+x，d5-x即可，修改完后要使用前缀和重新累加得到新的a[ i ]。

完整代码如下：

#include<bits/stdc++.h>
using namespace std;
using ll = long long;
#define qio ios::sync_with_stdio(false), cin.tie(0), cout.tie(0);
const int N = 1e5 + 5;
ll a[N], prefix[N], diff[N];
int main() {
	qio
	int n;
	cin >> n;
	for (int i = 1; i <= n; i++) cin >> a[i];
	for (int i = 1; i <= n; i++) diff[i] = a[i] - a[i - 1];
	int m;
	cin >> m;
	int l, r, x;
	for (int i = 1; i <= m; i++) {
		cin >> l >> r >> x;
		diff[l] += x;
		diff[r + 1] -= x;
	}
	for (int i = 1; i <= n; i++) a[i] = a[i - 1] + diff[i];
	for (int i = 1; i <= n; i++) prefix[i] = prefix[i - 1] + a[i];	
	int q;
	cin >> q;
	for (int i = 1; i <= q; i++) {
		cin >> l >> r;
		cout << prefix[r] - prefix[l - 1] << "\n";
	}
	return 0;
}

三、二维

二维就稍微复杂一点了，我得想想怎么画出抽象而又美丽的图了~

1.一维前缀和

P1060 - 【模板】二维前缀和 - ETOJ (eriktse.com)

计算prefix核心式：

for (int i = 1; i <= n; i++) {
		for (int j = 1; j <= m; j++) {
			prefix[i][j] = prefix[i - 1][j] + prefix[i][j - 1] - prefix[i - 1][j - 1] + a[i][j];
		}
	}

核心式：prefix[x2][y2] - prefix[x2][y1 - 1] - prefix[x1 - 1][y2] + prefix[x1 - 1][y1 - 1]

完整代码：

#include<bits/stdc++.h>
using namespace std;
using ll = long long;
#define qio ios::sync_with_stdio(false), cin.tie(0), cout.tie(0);
const int N = 1e3 + 5;
ll a[N][N], prefix[N][N], diff[N][N];
int main() {
	qio
	int n, m, q;
	cin >> n >> m >> q;
	for (int i = 1; i <= n; i++) {
		for (int j = 1; j <= m; j++) {
			cin >> a[i][j];
		}
	}
	for (int i = 1; i <= n; i++) {
		for (int j = 1; j <= m; j++) {
			prefix[i][j] = prefix[i - 1][j] + prefix[i][j - 1] - prefix[i - 1][j - 1] + a[i][j];
		}
	}
	for (int i = 1; i <= q; i++) {
		int x1, y1, x2, y2, ans;
		cin >> x1 >> y1 >> x2 >> y2;
		ans = prefix[x2][y2] - prefix[x2][y1 - 1] - prefix[x1 - 1][y2] + prefix[x1 - 1][y1 - 1];
		cout << ans << "\n";
	}
}

接着咕咕咕，写的好慢怎么办啊，突然想起劳动通识课还没看呢，一天天的也不知道都干了啥就过完了，阿巴阿巴

2.二维差分

P1061 - 【模板】二维差分 - ETOJ (eriktse.com)

计算diff核心式：

for (int i = 1; i <= n; i++) {
		for (int j = 1; j <= m; j++) {
			diff[i][j] = a[i - 1][j - 1] - a[i][j - 1] - a[i - 1][j] + a[i][j];
/*或者		diff[i][j]+=a[i][j];
			diff[i+1][j]-=a[i][j];
			diff[i][j+1]-=a[i][j];
			diff[i+1][j+1]+=a[i][j];		*/
		}
	}

核心式：

diff[x1][y1] += c;
diff[x1][y2 + 1] -= c;
diff[x2 + 1][y1] -= c;
diff[x2 + 1][y2 + 1] += c;

完整代码：

#include<bits/stdc++.h>
using namespace std;
using ll = long long;
#define qio ios::sync_with_stdio(false), cin.tie(0), cout.tie(0);

const int N = 1e3 + 5;
ll a[N][N], prefix[N][N], diff[N][N];
int main() {
	qio
	int n, m, q;
	cin >> n >> m >> q;
	for (int i = 1; i <= n; i++) {
		for (int j = 1; j <= m; j++) {
			cin >> a[i][j];
		}
	}
	for (int i = 1; i <= n; i++) {
		for (int j = 1; j <= m; j++) {
			diff[i][j] = a[i - 1][j - 1] - a[i][j - 1] - a[i - 1][j] + a[i][j];
/*			diff[i][j]+=a[i][j];
			diff[i+1][j]-=a[i][j];
			diff[i][j+1]-=a[i][j];
			diff[i+1][j+1]+=a[i][j];		*/
		}
	}
	for (int i = 1; i <= q; i++) {
		int x1, y1, x2, y2, c;
		cin >> x1 >> y1 >> x2 >> y2 >> c;
		diff[x1][y1] += c;
		diff[x1][y2 + 1] -= c;
		diff[x2 + 1][y1] -= c;
		diff[x2 + 1][y2 + 1] += c;
	}
	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] + diff[i][j];
			cout << a[i][j] << " ";
		}
		cout << "\n";
	}
}

四、应用

1.鼠鼠我鸭

P1034 - 鼠鼠我鸭 - ETOJ (eriktse.com)

最大字段和的求法，以及三目运算符的使用

#include<bits/stdc++.h>
using namespace std;
using ll = long long;
#define qio ios::sync_with_stdio(false), cin.tie(0), cout.tie(0);
void solve();

const int N = 1e5 + 5;
ll a[N], w[N];
int main() {
    qio
	int p;
	cin >> p;
	while (p--) {
		solve();
	}
}

void solve() {
	int n;
	cin >> n;
	for (int i = 0; i < n; i++)cin >> a[i];
	for (int i = 0; i < n; i++)cin >> w[i];
	//修改前的基本量
	ll ess = 0;
	for (int i = 0; i < n; i++) {
		ess += a[i] * w[i];
	}
	//ma为修改后区间的字段和，fix为所有修改后字段和中最大的一段
	ll fix = 0, ma = 0;
	for (int i = 0; i < n; i++) {
		ma = max(0ll, ma + (a[i] ? -1 : 1) * w[i]);
		fix = max(ma, fix);
	}
	cout << ess + fix << "\n";
}

2.卡牌游戏

P6625 [省选联考 2020 B 卷] 卡牌游戏 - 洛谷 | 计算机科学教育新生态 (luogu.com.cn)

思维题，合并卡牌并没有改变其排序以及前缀和

#include<bits/stdc++.h>
using namespace std;
#define qio ios::sync_with_stdio(0), cin.tie(0),cout.tie(0);
typedef long long ll;
const int N = 2e5 + 10;
ll a[N];
ll prefix[N];
int main() {
	qio
	int n;
	cin >> n;
	for (int i = 1; i <= n; i++) cin >> a[i];
	for (int i = 1; i <= n; i++) prefix[i] = prefix[i - 1] + a[i];
	ll sum = 0;
	for (int i = 1; i <= n; i++) {
		if (i != 1 && prefix[i] >= 0) sum += prefix[i];
	}
	cout << sum << '\n';
}

ov0_
关注
  • 2
    点赞
  • 踩
  • 1
    收藏
    觉得还不错? 一键收藏
  • 3
    评论
专栏目录
STL 自定义sort 前缀和差分练习
FourTO_的博客
05-18 406
公开课第一期 文章目录公开课第一期[A HDU1280 前m大的数](http://acm.hdu.edu.cn/showproblem.php?pid=1280)[B HDU1872 稳定排序](http://acm.hdu.edu.cn/showproblem.php?pid=1872)[C HDU1234 开门人和关门人](http://acm.hdu.edu.cn/showproblem.php?pid=1234)[D HDU1862 EXCEL排序](http://acm.hdu.edu.cn/s
ET-Campus算法竞赛_1_ 前缀和与差分
Neutronium_Zero的博客
01-21 453
有q个操作，每个操作格式为：x1​,y1​,x2​,y2​,c，其中(x1,y1)、(x2,y2)分别表示一个子矩阵的左上角和右下角的坐标，每个操作将对应的子矩阵的每个元素加上c。有q次询问，每次询问格式为：x1​,y1​,x2​,y2​，表示一个子矩阵的左上角和右下角的坐标。给定一个长度为n的数组a，有q次询问，每次询问数组a在区间 [ l , r ] 的和。在修改完成后有q次询问，每次询问数组在区间 [ l , r ] 的和。有m次修改，每次修改让数组a在区间 [ l , r ]增加x。
基础算法篇——前缀和与差分
Huangjiazhen711的博客
11-07 364
我们首先定义一个长度为n的数组，然后我们希望求这个数组的部分长度的总和复杂度为O（n）我们将数组的第i个值定义为ai我们将数组的前n个值的和定义为Sn其实就是类似于我们数学上的基本算法// 由于我们需要S[r] - S[l-1]完成计算// 那么当我们的l=0时，我们需要S[r]-S[-1],这明显是不可行的，但是如果我们将整体往前移动一位// 我们直接让数组从1开始，让S数组也从1开始，并将S[0]=0，这样我们在计算[1,k]之间的数时就可以直接使用S[r]-S[l-1]了。
【基础算法】前缀和与差分
qq_46162246的博客
08-08 231
实现 1.开两个数组s[],a[]，a[]用于存储输入数据，s[]用于存储a[]中前i项和 eg：s[i] = a[1]+a[2]+...+a[i] 2.将两个数组赋值，a[]输入赋值，s[]利用公式s[i] = s[i-1] + a[i]赋值，对于二维数组s[i][j] = s[i-1][j] + s[i][j-1] - s[i-1][j-1] + a[i][j]（这个很好理解，画图即可） 3.求子区间的和，一维公式res = s[endIndex] - s[startIndex - 1]。二维公式res
LADE:使用领导和伴随种群的增强的差分进化算法
01-24
来源：Li Y, Wang S, Yang H, et al. Enhancing differential evolution algorithm using leader-adjoint populations[J]. Information Sciences, 2023, 622: 235-268. 内容：CEC2017测试集, LADE. 注释：本算法为...
自适应差分进化算法代码
04-22
该代码是DE（差分进化）算法的改进，是把参数编码到个体中同时进行进化的参数自适应差分算法。 具体算法细节，可以参考文章： Brest J.Greiner S,Boskovic B.et al. Self-adapting control parameters in ...
西门子ET200SP基座单元的区别与分类以及注意事项.docx
12-26
接口模块是ET200SP系统的核心部分，负责与主控制器（如SIMATIC S7-1500 PLC）之间的通信。根据通信协议的不同，接口模块可分为PN接口模块（Profinet）和DP接口模块（DeviceNet）。PN接口模块支持以太网通信，提供...
基于粒子群和差分进化算法的含分布式电源配电网故障区段定位.pdf
08-10
"基于粒子群和差分进化算法的含分布式电源配电网故障区段定位" 一、分布式电源对配电网故障区段定位的影响 随着分布式电源的引入，传统的故障区段定位方法不再适用。分布式电源的引入将导致配电网中的电流和电压的...
S7-300与ET200SP通讯共18页.pdf.zip
11-21
使用SIMATIC Step 7软件进行配置和编程，可以实现S7-300与ET200SP之间的通讯。在编程时，需要设置CPU与ET200SP之间的通讯参数，如站地址、波特率等，并定义数据交换的逻辑。 5. **诊断与故障排除** 通过PLC的诊断...
图解大顶堆的构建、排序过程
sagima_sdu的博客
08-12 497
开始，
排序算法之桶排序
rd_w_csdn的博客
08-12 818
桶排序（bucket sort）是分治策略的一个典型应用。它通过设置一些具有大小顺序的桶，每个桶对应一个数据范围，将数据平均分配到各个桶中；然后，在每个桶内部分别执行排序；最终按照桶的顺序将所有数据合并。
备战秋招60天算法挑战，Day17
溜达虎的博客
08-13 365
给定两个字符串 和 ，编写一个函数来判断 是否是 的字母异位词。注意： 若 和 中每个字符出现的次数都相同，则称 和 互为字母异位词。举个例子： 视频题解 有效的字母异位词 思路来源 思路来源 思路解析 方法一 排序 分别对和进行排序，然后比较排序后的字符串是否相等。 java代码 python代码 复杂度分析 时间复杂度： 快排的时间复杂度是O(nlogn)，其中是字符串的长度。空间复杂度： 快排的空间复杂度是O(logn)。假设， 。 java代码 python代码 复杂度分析 时
算法：魔法字典
Garson的博客
08-12 548
设计一个使用单词列表进行初始化的数据结构，单词列表中的单词。如果给出一个单词，请判定能否只将这个单词中字母换成另一个字母，使得所形成的新单词存在于你构建的字典中。实现dictionarydictionarysearchWordtruefalse// 返回 False// 将第二个 'h' 替换为 'e' 可以匹配 "hello" ，所以返回 True// 返回 False// 返回 FalsedictionarybuildDictsearch100search。
【秋招笔试】8.11大疆秋招(第二套)-测开岗
清隆学长的博客
08-12 661
> 🍒 本专栏已收集 70+套笔试题，`笔试真题` 会在第一时间跟新 > > 🍹 近期会进行一波价格调整，没订阅的朋友抓紧冲！！！
数据结构——排序（4）：归并排序+计数排序
2301_81888634的博客
08-16 611
至此我们学习了顺序表、链表、栈、队列、二叉树、排序，初阶数据结构完结~撒花！我们C++见！！
代码随想录算法训练营第十七天
b1ankwall的博客
08-16 573
参加代码随想录第43期算法训练营day17|654.最大二叉树，617.合并二叉树，700.二叉搜索树中的搜索，98.验证二叉搜索树
【Python】生成二维迷宫的算法
最新发布
HarryXYC的博客
08-16 165
基于Pygame和PIL的迷宫生成器
算法之递归和迭代
qq_28576837的博客
08-16 131
斐波那契数列（Fibonacci sequence），又称黄金分割数列，这个数列最早是由印度数学家提出来的。如果是用迭代，那么就是用循环从前往后计算，前两个数是固定的，先根据前两个数可以计算出第三个数，然后把第三个数和第二个数再赋值给之前的变量，再相加即可。使用递归时，如果要计算第100个数，那么要先直到第99和第98个，要知道第99个和第98个，就要分别进行递归，直到最后一层，再返回递归回来进行计算。这个序列的前几项是这样的：0,1,1,2,3,5,8,13,21,34,⋯。自然，用递归函数最容易实现。
et200sp底座为什么分白色和蓝色
01-27
ET200SP底座分为白色和蓝色是为了区分不同功能或版本的ET200SP设备。 首先，白色底座通常用于普通的数字输入输出模块或传感器模块。这些模块通常用于控制系统中的信号输入和输出，用于接收和发送各种信号。白色底座...
评论 3
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值