C++图论之Bellman-Ford算法和SPFA算法的实现

 更新时间:2022年06月14日 15:55:15   作者:玄澈_  
贝尔曼-福特算法(Bellman-Ford)是由理查德·贝尔曼和莱斯特·福特创立的,求解单源最短路径问题的一种算法。SPFA 算法是 Bellman-Ford算法 的队列优化算法的别称,通常用于求含负权边的单源最短路径。本文将详解两个算法的实现,需要的可以参考一下

给定一张有向图,若对于图中的某一条边(x,y,z),有dist[y]≤dist[x]+z成立,则称该边满足三角形不等式。如果所有边都满足三角形不等式,则dist数组就是所求的最短路。

Bellman-Ford算法

(x,y,z)表示的是一条从 x 出发, 到达 y ,长度为 z 的有向边。

首先介绍基于迭代的Bellman-Ford算法,它的流程如下:

1.扫描所有边(x,y,z),若dist[y]>dist[x]+z, 则用dist[x]+z更新dist[y]

2.重复上述操作,直到没有更新操作发生。

Bellman-Ford算法的时间复杂度是O(nm)

通过Bellman-Ford算法我们可以求解有边数限制的最短路问题。

例题:AcWing 853. 有边数限制的最短路

算法步骤

初始化 dist 数组为正无穷, dist[1] = 0

(外重循环)循环 i 从 1 到 n ,遍历 n 次表示:是不经过超过 i 条边到达终点的最短距离

(内重循环)循环 i 从 1 到 m, 遍历 m 条边,把所有的边都进行松弛操作:

每次取出两点以及以及连接他们的权重 (a,b,w)

用以下公式更新最短距离: dist[b]=min(dist[b],dist[a]+w)

注意点:

需要把dist数组进行一个备份,这样防止每次更新的时候出现串联

由于存在负权边,所以 return -1 的条件是dist[n]>0x3f3f3f/2

代码实现

#include <iostream>
#include <cstring>
using namespace std;
 
const int N = 510, M = 10010;
 
struct Edge
{
    int a, b, w;
}e[M]; // 存下每一条即可
int dist[N];
int back[N]; // 备份数组放置串联
int n, m, k;
 
void bellman_ford()
{
    memset(dist, 0x3f, sizeof dist);
    dist[1] = 0;
    
    for(int i = 0; i < k; i ++ ) // 不超过k条边
    {
        memcpy(back, dist, sizeof back);
        for(int j = 0; j < m; j ++ ) // 遍历所有边
        {
            int a = e[j].a, b = e[j].b, w = e[j].w;
            dist[b] = min(dist[b], back[a] + w);
        }
    }
}
 
int main()
{
    cin >> n >> m >> k;
    for(int i = 0; i < m; i ++ )
    {
        int a, b, w;
        scanf("%d%d%d", &a, &b, &w);
        e[i] = {a, b, w};
    }
    
    bellman_ford();
    if(dist[n] > 0x3f3f3f3f / 2) puts("impossible");
    else cout << dist[n] << endl;
    
    return 0;
}

SPFA算法

SPFA算法在国际上通称为“队列优化的“Bellman-Ford算法”。

SPFA算法的流程如下:

1.建立一个队列,起初队列中只含有起点1

2.取出头结点 x ,扫描它的所有出边(x,y,z),若dist[y]>dist[x]+z,则使dist[y]用dist[x]+z来更新。同时若y不再队列中,则将y入队

在任意时刻,该算法的队列都保持了该拓展的节点。每次入队都相当于完成了一次 dist 数组的更新操作,使其满足三角不等式。一个节点可能会入队、出队多次。最终,图中所有的结点全部收敛到全部满足三角不等式的状态。

这个队列避免了对Bellman-Ford算法中不需要拓展的多余结点的冗余扫描,在随机图上的运行效率O(km)级别,其中 k 是一个很小的常数。

代码实现

SPFA求最短路

#include <cstring>
#include <iostream>
#include <algorithm>
#include <queue>
 
using namespace std;
 
const int N = 1e6 + 10;
 
int n, m;
int h[N], e[N], w[N], ne[N], idx;
int dist[N];
bool st[N];
 
void add(int a, int b, int c)
{
    e[idx] = b, ne[idx] = h[a], w[idx] = c, h[a] = idx ++ ;
}
 
void spfa()
{
    memset(dist, 0x3f, sizeof dist);
    queue<int> q;
    dist[1] = 0;
    st[1] = true;
    q.push(1);
    
    while(q.size())
    {
        int t = q.front();
        q.pop();
        
        st[t] = false;
        
        for(int i = h[t]; ~i; i = ne[i])
        {
            int j = e[i];
            if(dist[j] > dist[t] + w[i])
            {
                dist[j] = dist[t] + w[i];
                if(!st[j])
                {
                    q.push(j);
                    st[j] = true;
                }
            }
        }
    }
}
 
int main()
{
    scanf("%d%d", &n, &m);
 
    memset(h, -1, sizeof h);
    while (m -- )
    {
        int a, b, c;
        scanf("%d%d%d", &a, &b, &c);
        add(a, b, c);
    }
 
    spfa();
    
    if(dist[n] == 0x3f3f3f3f) puts("impossible");
    else printf("%d",dist[n]);
 
    return 0;
}
 

以上就是C++图论之Bellman-Ford算法和SPFA算法的实现的详细内容,更多关于C++ Bellman-Ford SPFA算法的资料请关注脚本之家其它相关文章!

相关文章

  • 详解C++中的常量

    详解C++中的常量

    这篇文章主要介绍了C++中的常量的相关资料,文中示例代码非常详细,帮助大家更好的理解和学习,感兴趣的朋友可以了解下
    2020-07-07
  • C语言sizeof与字符串处理与动态内存分配及main函数参数详解

    C语言sizeof与字符串处理与动态内存分配及main函数参数详解

    这篇文章主要介绍了C语言字符串处理函数、sizeof、动态内存分配函数、main函数参数问题,static在修饰变量的时候,如果是修饰全局变量,则跟全局变量功能一样,通过示例代码给大家介绍的非常详细,需要的朋友可以参考下
    2022-07-07
  • C语言中双向链表和双向循环链表详解

    C语言中双向链表和双向循环链表详解

    这篇文章主要介绍了C语言中双向链表和双向循环链表详解的相关资料,需要的朋友可以参考下
    2017-06-06
  • C语言数据类型和格式说明符基础教程示例

    C语言数据类型和格式说明符基础教程示例

    这篇文章主要为大家介绍了C语言数据类型和格式说明符基础教程示例详解,有需要的朋友可以借鉴参考下,希望能够有所帮助,祝大家多多进步,早日升职加薪
    2023-12-12
  • visual studio 2013中配置opencv图文教程 Opencv2.4.9安装配置教程

    visual studio 2013中配置opencv图文教程 Opencv2.4.9安装配置教程

    这篇文章主要为大家详细介绍了Opencv2.4.9安装教程,以及在visualstudio 2013中opencv的配置步骤,具有一定的参考价值,感兴趣的小伙伴们可以参考一下
    2017-04-04
  • C语言详解如何实现带头双向循环链表

    C语言详解如何实现带头双向循环链表

    带头双向循环链表:结构最复杂,一般用在单独存储数据。实际中使用的链表数据结构,都是带头双向循环链表。另外这个结构虽然结构复杂,但是使用代码实现以后会发现结构会带来很多优势,实现反而简单
    2022-04-04
  • C语言关键字auto与register及static专项详解

    C语言关键字auto与register及static专项详解

    这篇文章主要解释了c语言中什么是数据类型,什么是变量,他们的真正含义是什么。分析了属性关键字auto,register和static的用法
    2022-07-07
  • C语言数据结构与算法之排序总结(二)

    C语言数据结构与算法之排序总结(二)

    这篇文章住要介绍的是选择类排序中的简单、树形和堆排序,归并排序、分配类排序的基数排序,文中的示例代码讲解详细,感兴趣的小伙伴可以了解一下
    2021-12-12
  • C++执行shell命令的多种实现方法

    C++执行shell命令的多种实现方法

    在linux系统下,用C++程序执行shell命令有多种方式,主要介绍了3中方法,具有一定的参考价值,感兴趣的可以了解一下
    2021-11-11
  • Qt中JSON操作的具体使用

    Qt中JSON操作的具体使用

    本文主要介绍了Qt中JSON操作的具体使用,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着小编来一起学习学习吧
    2022-08-08

最新评论