弗洛伊德(Floyd)算法(C/C++)

2024-08-28 18:12
文章标签 算法 c++ floyd 弗洛伊德

本文主要是介绍弗洛伊德(Floyd)算法(C/C++),希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

弗洛伊德算法(Floyd's algorithm),又称为弗洛伊德-沃尔什算法(Floyd-Warshall algorithm),是一种用于在加权图中找到所有顶点对之间最短路径的算法。这个算法适用于有向图和无向图,并且可以处理负权重边,但不能处理负权重循环。

弗洛伊德算法(Floyd-Warshall Algorithm)是一种用于计算图中所有顶点对之间最短路径的动态规划算法。本文将详细介绍弗洛伊德算法的原理,并提供一个C++实现的示例,以帮助读者理解算法的工作原理和编程技巧。

算法原理

弗洛伊德算法的核心思想是通过逐步寻找并更新所有顶点对之间的最短路径来解决问题。算法使用一个距离矩阵来存储顶点之间的距离,并在每一步中考虑通过一个新的中间顶点来更新这些距离。跟上一篇Dijkstra算法一样的原理,也是通过中转点去更新最短距离。不过Floyd算法处理的是多源的最短路问题


算法步骤

  1. 初始化一个距离矩阵,其中dist[i][j]表示顶点i到顶点j的直接距离。如果ij不直接相连,则dist[i][j]为无穷大。
  2. 对于每个顶点k,作为中间顶点,更新dist[i][j]min(dist[i][j], dist[i][k] + dist[k][j])

Floyd是经典三重for循环,所以它的时间复杂度为o(n^3),n是图中顶点的数量。第一层遍历中转点,第二层遍历起点,第三层遍历终点,对于图中点的数量多的情况,Floyd算法的时间复杂度是很高的。

图解算法:

下面我们将以4个点的图进行讲解,图的连边为有向边和无向边的结合。以邻接矩阵的方式进行存储,如果大家喜欢用邻接表存储,也可以使用邻接表,下面介绍两个矩阵,矩阵A表示(i,j)i->j的最短距离,初始化为inf。矩阵B表示i->j路径由i到j的中转点,也就是路径上除去起点的第一个点,初始化为-1。

初始:

按照图中的点距离给其赋值,A矩阵i->i距离都为0,inf为无法到达。B矩阵初始为-1。

第一步:

我们选取一个点(按照顺序选取)把它作为中转点,看看以它为中转点,所能到达的点中有没有产生更小的距离,如果产生了,则更新A矩阵的距离,更新B矩阵的中转点。我们先选取1号点,那么位于1号点的行跟列的值都是不可能变化的,还有就是自己到自己的点也是不会变化的永远是0,图中黄颜色标记的是此步不会改变的点,其他的可能会变。在更新距离的时候我们可以不看图就能更新矩阵,例如下图中2号点到3号点本来为10,我们可以连一个矩阵,以1号点画的两条蓝线为两条边,红色线为剩余2边,我们既然把1号点当作中转点,路径必然为2-1-3,此时距离就是副对角线的顶点值相加2+6=8<10,那么通过1号点绕路的方式距离更短。类似的还有3->2号点,6+2=8<inf。3->4号点,10+6=16<inf。4->3号点,10+6=16<inf。顺便把B矩阵更新完。

更新完后(红色标记为变化的值): 

 第二步:

此时把2号结点作为中转结点,看一看能够更新哪一个最短路径,还是跟上一步一样直接看图更新就可以。如下图,4->1号点,2+4=6<10。1->4号点,2+4=6<10。3->4号点,8+4=12<16。4->3号点,8+4=12<16。对于一些不能更新的值,例如1->3号点,2+8=10>6,这样的则不能更新。

对于B矩阵,要注意3->4跟4->3的路径是相反的,更新是则不能直接修改为2,对于3->4号点第一个中转点还是1号点。更新完后(红色标记为变化的值): 

第三步:

把3号点作为中转结点,跟前几步一样,继续寻找最短距离。经过更新我们发现3号点作为中转点不能更新任意一个距离,所以A、B矩阵不需要更新。在图中,经过验证我们发现3号点中转距离反而变大,所以不更新。

第四步:

把4号点作为中转点,继续更新最短距离。我们发现跟3号点一样,不能更新任何距离,在A矩阵中除了黄色的点之外,所能连起来的矩形,主对角线顶点值相加都比当前值要大。在图中也可以验证,所以不给予更新。

这样我们就更新完所有点,把所有点都当作中转点更新完一遍,这样就完成了Floyd算法,更新时每次按照顺序把点当作中转点,遍历寻找路径的起点,再遍历寻找终点,算法时间复杂度为o(n^3)。

视频讲解可以看一下B站这位UP主的讲解,点击直达


算法实现:

以下是弗洛伊德算法的C++实现示例:

#include <iostream>
#include <vector>
#include <limits>
using namespace std;// 定义图的顶点数
const int N = 100;
// 定义无穷大的初始距离
const int INF = numeric_limits<int>::max();// 弗洛伊德算法的实现
void floydWarshall(vector<vector<int>>& dist) {int n = dist.size();// 遍历所有顶点作为中间顶点for (int k = 0; k < n; k++) {// 遍历所有顶点作为起点for (int i = 0; i < n; i++) {// 遍历所有顶点作为终点for (int j = 0; j < n; j++) {// 如果通过顶点k可以找到更短的路径,则更新dist[i][j]if (dist[i][k] != INF && dist[k][j] != INF && dist[i][k] + dist[k][j] < dist[i][j]) {dist[i][j] = dist[i][k] + dist[k][j];}}}}
}int main() {int n; // 顶点的数量cin >> n;vector<vector<int>> dist(n, vector<int>(n, INF)); // 初始化距离矩阵// 读取邻接矩阵for (int i = 0; i < n; i++) {dist[i][i] = 0; // 自己到自己的距离是0for (int j = i; j < n; j++) {int w;cin >> w;dist[i][j] = w;dist[j][i] = w; // 如果是无向图,需要设置对称的权重}}// 执行弗洛伊德算法floydWarshall(dist);// 打印所有顶点对之间的最短路径for (int i = 0; i < n; i++) {for (int j = 0; j < n; j++) {if (dist[i][j] == INF) {cout << "INF" << " ";} else {cout << dist[i][j] << " ";}}cout << endl;}return 0;
}

Floyd与Dijkstra算法比较 

迪杰斯特拉算法(Dijkstra's algorithm)和弗洛伊德算法(Floyd-Warshall algorithm)都是图论中用于计算图中最短路径的著名算法。它们在某些方面有相似之处,但在设计和应用上存在显著差异,下面我们将对这两种算法的相同跟不同进行解释。

相同点:
  1. 目的两者都旨在解决最短路径问题。
  2. 适用性:它们都可以用于加权图中的最短路径计算,无论是正权还是负权(只有弗洛伊德算法)。
不同点:
  1. 问题范围:
    • 迪杰斯特拉算法:主要用于单元路径的最短路问题,即从单一源点到所有其他顶点的最短路径。
    • 弗洛伊德算法:解决的是所有顶点对之间的最短路径问题,即计算图中每一对顶点之间的最短路径。
  2. 时间复杂度:
    • 迪杰斯特拉算法:具有较高的效率,时间复杂度为O(V^2)(使用朴素实现)或O((V+E) log V)(使用优先队列优化)。(V顶点E条边)
    • 弗洛伊德算法:时间复杂度为O(V^3),因为它需要计算所有顶点对的最短路径。
  3. 实现方式:
    • 迪杰斯特拉算法:通常使用贪心策略,从一个顶点开始,逐步扩展到邻接顶点,直到找到所有顶点的最短路径。
    • 弗洛伊德算法:使用动态规划,通过三层循环迭代地改进路径长度,直到达到最优解。
  4. 对负权边的处理:
    • 迪杰斯特拉算法:不能处理负权边,因为负权边会破坏算法的贪心选择性质。
    • 弗洛伊德算法:可以处理负权边,但图中不能有负权环,否则最短路径问题没有解。
  5. 初始化:
    • 迪杰斯特拉算法:从源点到其他所有顶点的距离初始化为无穷大,源点到自身的距离为0。
    • 弗洛伊德算法:所有顶点到自身的距离初始化为0,其他顶点间的距离初始化为边的权重或无穷大(如果无直接连接)。

本篇详解Floyd算法,如果想看Dijkstra算法的话,可以看博主上一篇博客,针对于Dijkstra算法的详解:迪杰斯特拉(Dijkstra)算法(C/C++)-CSDN博客

执笔至此,感触彼多,全文将至,落笔为终,感谢大家的支持。  

这篇关于弗洛伊德(Floyd)算法(C/C++)的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/1115513

相关文章

C++中unordered_set哈希集合的实现

《C++中unordered_set哈希集合的实现》std::unordered_set是C++标准库中的无序关联容器,基于哈希表实现,具有元素唯一性和无序性特点,本文就来详细的介绍一下unorder... 目录一、概述二、头文件与命名空间三、常用方法与示例1. 构造与析构2. 迭代器与遍历3. 容量相关4

C++中悬垂引用(Dangling Reference) 的实现

《C++中悬垂引用(DanglingReference)的实现》C++中的悬垂引用指引用绑定的对象被销毁后引用仍存在的情况,会导致访问无效内存,下面就来详细的介绍一下产生的原因以及如何避免,感兴趣... 目录悬垂引用的产生原因1. 引用绑定到局部变量,变量超出作用域后销毁2. 引用绑定到动态分配的对象,对象

深入理解Mysql OnlineDDL的算法

《深入理解MysqlOnlineDDL的算法》本文主要介绍了讲解MysqlOnlineDDL的算法,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着小... 目录一、Online DDL 是什么?二、Online DDL 的三种主要算法2.1COPY(复制法)

C++读写word文档(.docx)DuckX库的使用详解

《C++读写word文档(.docx)DuckX库的使用详解》DuckX是C++库,用于创建/编辑.docx文件,支持读取文档、添加段落/片段、编辑表格,解决中文乱码需更改编码方案,进阶功能含文本替换... 目录一、基本用法1. 读取文档3. 添加段落4. 添加片段3. 编辑表格二、进阶用法1. 文本替换2

C++中处理文本数据char与string的终极对比指南

《C++中处理文本数据char与string的终极对比指南》在C++编程中char和string是两种用于处理字符数据的类型,但它们在使用方式和功能上有显著的不同,:本文主要介绍C++中处理文本数... 目录1. 基本定义与本质2. 内存管理3. 操作与功能4. 性能特点5. 使用场景6. 相互转换核心区别

C++右移运算符的一个小坑及解决

《C++右移运算符的一个小坑及解决》文章指出右移运算符处理负数时左侧补1导致死循环,与除法行为不同,强调需注意补码机制以正确统计二进制1的个数... 目录我遇到了这么一个www.chinasem.cn函数由此可以看到也很好理解总结我遇到了这么一个函数template<typename T>unsigned

C++统计函数执行时间的最佳实践

《C++统计函数执行时间的最佳实践》在软件开发过程中,性能分析是优化程序的重要环节,了解函数的执行时间分布对于识别性能瓶颈至关重要,本文将分享一个C++函数执行时间统计工具,希望对大家有所帮助... 目录前言工具特性核心设计1. 数据结构设计2. 单例模式管理器3. RAII自动计时使用方法基本用法高级用法

深入解析C++ 中std::map内存管理

《深入解析C++中std::map内存管理》文章详解C++std::map内存管理,指出clear()仅删除元素可能不释放底层内存,建议用swap()与空map交换以彻底释放,针对指针类型需手动de... 目录1️、基本清空std::map2️、使用 swap 彻底释放内存3️、map 中存储指针类型的对象

C++ STL-string类底层实现过程

《C++STL-string类底层实现过程》本文实现了一个简易的string类,涵盖动态数组存储、深拷贝机制、迭代器支持、容量调整、字符串修改、运算符重载等功能,模拟标准string核心特性,重点强... 目录实现框架一、默认成员函数1.默认构造函数2.构造函数3.拷贝构造函数(重点)4.赋值运算符重载函数

C++ vector越界问题的完整解决方案

《C++vector越界问题的完整解决方案》在C++开发中,std::vector作为最常用的动态数组容器,其便捷性与性能优势使其成为处理可变长度数据的首选,然而,数组越界访问始终是威胁程序稳定性的... 目录引言一、vector越界的底层原理与危害1.1 越界访问的本质原因1.2 越界访问的实际危害二、基