一文带你深入了解算法笔记中的前缀与差分(附源码)

简介: 一文带你深入了解算法笔记中的前缀与差分(附源码)

文章目录

一、前缀和

A、一维前缀和

1、什么是一维前缀和

2、一维前缀和的作用

3、习题:Acwing 795. 前缀和

输入格式

输出格式

数据范围

输入样例:

输出样例:

4、代码详解

B、二维前缀和(矩阵和)

1、二维前缀和推导

2、习题:Acwing 796. 子矩阵的和

输入格式

输出格式

数据范围

输入样例:

输出样例:

3、代码详解

二、差分

A、一维差分

1、什么是差分

2、如何构建差分数组

3、差分数组有什么作用

4、练习 Acwing 797. 差分

输入格式

输出格式

数据范围

输入样例:

输出样例:

5、代码详解

B、二维差分(差分矩阵)

3、习题:Acwing 798. 差分矩阵

输入格式

输出格式

数据范围

输入样例:

输出样例:

4、代码详解


一、前缀和

A、一维前缀和

1、什么是一维前缀和

原数组: a[1], a[2], a[3], a[4], a[5], …, a[n]

前缀和 S[i]为数组的前 i项和

前缀和: S[i] = a[1] + a[2] + a[3] + … + a[i]


注意:前缀和的下标一定要从 1开始, 避免进行下标的转换,即定义S[0]等于0

2、一维前缀和的作用

快速求出元素组中某段区间的和


例如要求出l-r区间的和,则只需要执行sum[r]-sum[l-1]


原理如下:


sum[r] = a[1] + a[2] + a[3] + a[l-1] + a[l] + a[l+1] ...... a[r]; sum[l - 1] = a[1] + a[2] + a[3] + a[l - 1]; sum[r] - sum[l - 1] = a[l] + a[l + 1]+......+ a[r];


并且时间复杂度为O(1)

3、习题:Acwing 795. 前缀和

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

接下来再输入 m个询问,每个询问输入一对 l,r。

对于每个询问,输出原序列中从第 l个数到第 r个数的和

输入格式

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


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


接下来 m 行,每行包含两个整数 l和 r,表示一个询问的区间范围。


输出格式

共 m行,每行输出一个询问的结果。


数据范围

1≤l≤r≤n

1≤n,m≤100000

−1000≤数列中元素的值≤1000


输入样例:5 3

               2 1 3 6 4

              1 2

              1 3

              2 4

输出样例:
3
6
10
4、代码详解

#include<iostream>

using namespace std;

const int N=1e5+10;

int a[N],sum[N];  //原数组与求和数组

int main()

{

   int n,m;

   cin>>n>>m;

   for(int i=1;i<=n;i++)

   {

       scanf("%d",&a[i]);

       sum[i]=sum[i-1]+a[i];   //构建前缀和

   }

   while(m--)

   {

       int L,R;

       cin>>L>>R;

       cout<<sum[R]-sum[L-1]<<endl;

   }

   return 0;

}


B、二维前缀和(矩阵和)

1、二维前缀和推导

微信截图_20230425230646.png

接下来回归问题去求以(x1,y1)为左上角和以(x2,y2)为右下角的矩阵的元素的和。

如图

微信截图_20230425230816.png

不难看出绿色的面积=s[x2][y2]+s[x2][y1-1]+s[x1-1][y2]-s[x1-1][y1-1]


且时间复杂度也为O(1)


2、习题:Acwing 796. 子矩阵的和

输入一个 n行 m列的整数矩阵,再输入 q个询问,每个询问包含四个整数 x1,y1,x2,y2,表示一个子矩阵的左上角坐标和右下角坐标。对于每个询问输出子矩阵中所有数的和。


输入格式

第一行包含三个整数 n,m,q。


接下来 n行,每行包含 m个整数,表示整数矩阵。


接下来 q行,每行包含四个整数 x1,y1,x2,y2,表示一组询问。


输出格式

共 q 行,每行输出一个询问的结果。


数据范围

1≤n,m≤1000

1≤q≤200000

1≤x1≤x2≤n

1≤y1≤y2≤m

−1000≤矩阵内元素的值≤1000


输入样例:

3 4 3
1 7 2 4
3 6 2 8
2 1 2 3
1 1 2 2
2 1 3 4
1 3 3 4

输出样例:

17

27

21

3、代码详解

#include<iostream>

using namespace std;

const int N=1010;

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

int main()

{

   int n,m,q;

   cin>>n>>m>>q;

   for(int i=1;i<=n;i++)

       for(int j=1;j<=m;j++){

           scanf("%d",&a[i][j]);

           //构建前缀和

           s[i][j]=s[i-1][j]+s[i][j-1]-s[i-1][j-1]+a[i][j];

       }

   while(q--)

   {

       int x1,y1,x2,y2;

       scanf("%d%d%d%d",&x1,&y1,&x2,&y2);

       //求并输出矩形前缀和

       cout<<s[x2][y2]-s[x1-1][y2]-s[x2][y1-1]+s[x1-1][y1-1]<<endl;

   }

   return 0;

 

}

————————————————


二、差分

A、一维差分

1、什么是差分

类似于数学中的求导和积分,差分可以看成前缀和的逆运算。


这里举一个例子:


先给一个原数组a:a[1]、a[2]、a[3]、a[4]、a[5]...a[n]


然后这里再构建一个数组b:b[1]、b[2]、b[3]...b[n]


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


即a是b的前缀和,反过来b就是a的差分


2、如何构建差分数组

这里我们采取最暴力的做法,直接构造,并且需要注意下标,所以我们这里还是令a[0]=0

注意:下面的构造中,a为前缀和数组,b为差分数组

a[0]=0;

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

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

......

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

3、差分数组有什么作用

我们先来看一个问题

给定区间[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的一项就会改变a中有b的所有项

但是在有边界之外的也多加上一个C,如图所示

微信截图_20230425231409.png

所以这里我们需要执行两步操作

b[l]+=c;

b[r+1]-=c;

因此我们得出一维差分结论:给a数组中的[ l, r]区间中的每一个数都加上c,只需对差分数组b做 b[l] + = c, b[r+1] - = c。时间复杂度为O(1), 大大提高了效率。


4、练习 Acwing 797. 差分

输入一个长度为 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

5、代码详解

注意使用scanf进行输入,效率比较高

#include<iostream>

using namespace std;

const int N = 1e5 + 10;

int a[N], b[N];

int main()

{

   int n, m;

   scanf("%d%d", &n, &m);

   for (int i = 1; i <= n; i++)

   {

       scanf("%d", &a[i]);

       b[i] = a[i] - a[i - 1];      //构建差分数组

   }

   int l, r, c;

   while (m--)

   {

       scanf("%d%d%d", &l, &r, &c);

       b[l] += c;       //将序列中[l, r]之间的每个数都加上c

       b[r + 1] -= c;

   }

   for (int i = 1; i <= n; i++)

   {

       a[i] = b[i] + a[i - 1];    //前缀和运算

       printf("%d ", a[i]);

   }

   return 0;

}

B、二维差分(差分矩阵)

类似于前面的二维前缀和,这里引用大佬林小鹿的图解

微信截图_20230425231807.png3、习题:Acwing 798. 差分矩阵

输入一个 n行 m列的整数矩阵,再输入 q 个操作,每个操作包含五个整数 x1,y1,x2,y2,c,其中 (x1,y1)和 (x2,y2) 表示一个子矩阵的左上角坐标和右下角坐标。每个操作都要将选中的子矩阵中的每个元素的值加上 c。请你将进行完所有操作后的矩阵输出。

输入格式

第一行包含整数 n,m,q。

接下来 n行,每行包含 m个整数,表示整数矩阵。

接下来 q行,每行包含 55 个整数 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

————————————————

4、代码详解

#include<iostream>

using namespace std;

const int N=1010;

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

int main()

{

   int n,m,q;

   cin>>n>>m>>q;

   for(int i=1;i<=n;i++)

       for(int j=1;j<=m;j++)

           scanf("%d",&s[i][j]);

   //构建差分矩阵

   for(int i=1;i<=n;i++)

       for(int j=1;j<=m;j++)

           a[i][j]=s[i][j]-s[i][j-1]-s[i-1][j]+s[i-1][j-1];  

   while(q--)

   {

       int x1,y1,x2,y2,c;

       scanf("%d%d%d%d%d",&x1,&y1,&x2,&y2,&c);

       a[x1][y1]+=c;

       a[x1][y2+1]-=c;

       a[x2+1][y1]-=c;

       a[x2+1][y2+1]+=c;

   }

   //求和

   for(int i=1;i<=n;i++)

       for(int j=1;j<=m;j++)

          s[i][j] = a[i][j] + s[i - 1][j] + s[i][j - 1] - s[i - 1][j - 1];

   for (int i = 1; i <= n; i ++ )

   {

       for (int j = 1; j <= m; j ++ )  

       printf("%d ", s[i][j]);

       cout << endl;

   }

   return 0;

}

版权声明:本文为CSDN博主「热爱编程的小K」的原创文章,遵循CC 4.0 BY-SA版权协议,转载请附上原文出处链接及本声明。

原文链接:https://blog.csdn.net/qq_72157449/article/details/130029246

相关文章
|
4天前
|
机器学习/深度学习 前端开发 算法
婚恋交友系统平台 相亲交友平台系统 婚恋交友系统APP 婚恋系统源码 婚恋交友平台开发流程 婚恋交友系统架构设计 婚恋交友系统前端/后端开发 婚恋交友系统匹配推荐算法优化
婚恋交友系统平台通过线上互动帮助单身男女找到合适伴侣,提供用户注册、个人资料填写、匹配推荐、实时聊天、社区互动等功能。开发流程包括需求分析、技术选型、系统架构设计、功能实现、测试优化和上线运维。匹配推荐算法优化是核心,通过用户行为数据分析和机器学习提高匹配准确性。
22 3
|
1月前
|
搜索推荐 算法 C语言
【排序算法】八大排序(上)(c语言实现)(附源码)
本文介绍了四种常见的排序算法:冒泡排序、选择排序、插入排序和希尔排序。通过具体的代码实现和测试数据,详细解释了每种算法的工作原理和性能特点。冒泡排序通过不断交换相邻元素来排序,选择排序通过选择最小元素进行交换,插入排序通过逐步插入元素到已排序部分,而希尔排序则是插入排序的改进版,通过预排序使数据更接近有序,从而提高效率。文章最后总结了这四种算法的空间和时间复杂度,以及它们的稳定性。
94 8
|
1月前
|
搜索推荐 算法 C语言
【排序算法】八大排序(下)(c语言实现)(附源码)
本文继续学习并实现了八大排序算法中的后四种:堆排序、快速排序、归并排序和计数排序。详细介绍了每种排序算法的原理、步骤和代码实现,并通过测试数据展示了它们的性能表现。堆排序利用堆的特性进行排序,快速排序通过递归和多种划分方法实现高效排序,归并排序通过分治法将问题分解后再合并,计数排序则通过统计每个元素的出现次数实现非比较排序。最后,文章还对比了这些排序算法在处理一百万个整形数据时的运行时间,帮助读者了解不同算法的优劣。
105 7
|
2月前
|
算法 API 计算机视觉
人脸识别笔记(一):通过yuface调包(参数量54K更快更小更准的算法) 来实现人脸识别
本文介绍了YuNet系列人脸检测算法的优化和使用,包括YuNet-s和YuNet-n,以及通过yuface库和onnx在不同场景下实现人脸检测的方法。
73 1
|
2月前
|
JSON 算法 数据可视化
测试专项笔记(一): 通过算法能力接口返回的检测结果完成相关指标的计算(目标检测)
这篇文章是关于如何通过算法接口返回的目标检测结果来计算性能指标的笔记。它涵盖了任务描述、指标分析(包括TP、FP、FN、TN、精准率和召回率),接口处理,数据集处理,以及如何使用实用工具进行文件操作和数据可视化。文章还提供了一些Python代码示例,用于处理图像文件、转换数据格式以及计算目标检测的性能指标。
74 0
测试专项笔记(一): 通过算法能力接口返回的检测结果完成相关指标的计算(目标检测)
|
2月前
|
算法
❤️算法笔记❤️-(每日一刷-160、相交链表)
❤️算法笔记❤️-(每日一刷-160、相交链表)
18 1
|
2月前
|
数据可视化 搜索推荐 Python
Leecode 刷题笔记之可视化六大排序算法:冒泡、快速、归并、插入、选择、桶排序
这篇文章是关于LeetCode刷题笔记,主要介绍了六大排序算法(冒泡、快速、归并、插入、选择、桶排序)的Python实现及其可视化过程。
21 0
|
2月前
|
算法
❤️算法笔记❤️-(每日一刷-83、删除排序链表中的重复项)
❤️算法笔记❤️-(每日一刷-83、删除排序链表中的重复项)
34 0
|
2月前
|
算法
❤️算法笔记❤️-(每日一刷-26、删除有序数组的重复项)
❤️算法笔记❤️-(每日一刷-26、删除有序数组的重复项)
27 0
|
16天前
|
算法
基于WOA算法的SVDD参数寻优matlab仿真
该程序利用鲸鱼优化算法(WOA)对支持向量数据描述(SVDD)模型的参数进行优化,以提高数据分类的准确性。通过MATLAB2022A实现,展示了不同信噪比(SNR)下模型的分类误差。WOA通过模拟鲸鱼捕食行为,动态调整SVDD参数,如惩罚因子C和核函数参数γ,以寻找最优参数组合,增强模型的鲁棒性和泛化能力。
下一篇
DataWorks