图论------贝尔曼-福德(Bellman-Ford)算法

算法概述:

Bellman-Ford算法核心代码如下

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

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

if(dicv\[j]> dicu\[j] + wj]

dicv\[j] = dicu\[j] + wj;

首先我们要了解一个点就是我们这次不再使用邻接矩阵来存储图的信息,而是定义三个一维数组来存储图的信息。

首先定义 un 来存储边的起点,vn来存储边的终点,wn来存储边的长(权重)。

例如存储如下元素

4 4
1 2 3
2 4 7
3 4 5
3 1 3


和Dijkstra算法一样,这段代码也是求单源最短路径,所以我们同样也要定义一个dicn数组来存储所有点到1的最短距离。注意,这里因为使用的并非邻接矩阵,所以初始化dic时默认除了dic1 = 0以外,暂时全部看作不连通,也就是假设以以上数据为例,dic内容要初始化为

【0 ,inf, inf, inf】,其中inf看作无穷大,表示不通。

那么我们再来看核心代码的最后两句。

if(dicv\[j]> dicu\[j] + wj]

dicv\[j] = dicu\[j] + wj;

dicv\[j]的值是1 到 vj这个点的值, dicu\[j] 是 1 到 uj 这个点的值 , wj 是 uj 到 vj 的值,意思是如果如果 1 通过 1 -> uj -> vj 这条路比 1 -> vj 值小的话,更新dicv\[j]的值。也就是通过uj ->vj 这条边进行松弛来优化路径。

初始化示意图如下:

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

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

if(dicv\[j]> dicu\[j] + wj]

dicv\[j] = dicu\[j] + wj;

当第一轮循环时:

当先通过 j = 1 通过第一条边进行优化时, if(dic 2 > dic 1 + 3],发现dic2 需优化,然后

重置 dic2 = dic1+3 = 3。

当 j = 2 时 通过第二条边进行优化时, if(dic 4 > dic 2 + 7],发现dic4 可以优化,然后

重置 dic 4 = dic 2 + 7 = 3+7 = 10。

当 j = 3时,通过第三条边进行优化时,if(dic 4 > dic 3 + 5],发现dic4无法优化,因为

dic 4 = 10 < dic 3 + 5 =inf + 5。

**当 j = 4 时 通过第四条边进行优化时, if(dic 1 > dic 3 + 3],**发现dic1无法优化,因为

dic 1 = 0 < dic 3 + 3 =inf + 3。

最后经过一轮松弛后得到 dic = 0, 3, inf ,10;


当 i= 2 进行第二轮松弛时,

当先通过 j = 1 通过第一条边进行优化时, if(dic 2 > dic 1 + 3],发现dic2 无需优化,因为****dic2 = dic1+3 = 3。

当 j = 2 时 通过第二条边进行优化时, if(dic 4 > dic 2 + 7],发现dic4 无需优化

当 j = 3时,通过第三条边进行优化时,if(dic 4 > dic 3 + 5],发现dic4无需优化

**当 j = 4 时 通过第四条边进行优化时, if(dic 1 > dic 3 + 3],**发现dic1无需优化

这样我们就完成了这组数据,因为这是构建的有向图,所以1 才无法到达 3.

既然很清晰流程了 我再丢给大家一组数据:

5 7
1 2 8
1 3 1
1 4 2
3 4 2
2 4 3
3 5 3
4 5 3

现在请找到一张草稿纸来试着执行刚才所讲的流程。

具体代码:

#include<stdio.h>
int main(void)
{
int u10, v10, w10;
int n, m, inf = 99999999;
int dic10 = { 0 };

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

for (int i = 1; i <= m; i++)
scanf_s("%d%d%d", &ui, &vi, &wi);****//存储边的信息

for (int i = 2; i <= n; i++)
dici = inf;****//初始化1 的单源最短路径值。

//核心代码

for (int i = 1; i <= n - 1; i++)
for (int j = 1; j <= m; j++)
if (dicv\[j] > dicu\[j] + wj)
dicv\[j] = dicu\[j] + wj;

for (int i = 1; i <= n; i++)
printf("%d ", dici);****//输出结果

}

注意事项:

最外层代码循环n-1次的原因是,n 个顶点,1到任何顶点的最多经过n-1个顶点。每一轮松弛循环如果没有完成单源最短路径的最终结果,都至少会有一次松弛。

细心的同学可能会发现我们刚才使用的数据:

5 7
1 2 8
1 3 1
1 4 2
3 4 2
2 4 3
3 5 3
4 5 3

是上一篇文章的内容,在上一篇中结果是:0 5 1 2 4 。

而在这里答案是:0 8 1 2 4。

那么为什么会造成这种影响呢,因为我们使用的这段代码是处理有向图的,而在上一篇文章中是处理无向图的,那么有没有什么办法来改进这段代码使代码能够处理无向图呢?答案当然是可以的,而且十分简单。

聪明的同学注意到数据

4 4
1 2 3
2 4 7
3 4 5
3 1 3

在第二轮就能发现没有优化证明已经得出结果,每必要进行n-1轮循环,那么大家思考下,怎么提前退出循环呢?

课后作业:

改进这段代码使代码能够处理无向图。

如何证明达到最优解提前退出循环。

明天带来答案。

相关推荐
hetao17338372 小时前
联盟场-提高组-R7 & ZYZSC-S-Round 11
c++·算法
啵啵啵鱼2 小时前
第二章 线性表
java·数据结构·算法
虾..2 小时前
计算机图形学:几何表示
算法
无敌贵点大王2 小时前
嵌入式软件面试经(自用)
c语言·stm32·单片机·学习·mongodb·面试·职场和发展
局部有解2 小时前
Hot100 - 三数之和
数据结构·c++·算法
辉艺棉麻小张2 小时前
用 JSON 建模面料规格:成分、克重、支数与公定回潮率字段设计
数据结构
学逆向的3 小时前
C 数据溢出与精度丢失
c语言·开发语言·网络安全
longlongzihan3 小时前
C++ 常用数据结构速查:vector、stack、queue、list、树、图、哈希表等声明与核心操作
数据结构·c++·list
思茂信息3 小时前
CST技术案例:加载正弦曲线的平面螺旋天线仿真研究
算法·平面·emc·cst·电磁仿真·仿真辐射
嫂子开门我是_我哥3 小时前
敏感性分析
算法·机器学习·论文笔记