zoukankan      html  css  js  c++  java
  • Bellman-Ford算法

    参考资料:

    最短路径算法—Bellman-Ford(贝尔曼-福特)算法分析与实现(C/C++)
    Bellman-Ford算法详讲

    算法推导建议参考第一篇文章,讲的很详细。

    解析(摘自参考资料):

    Dijkstra算法是处理单源最短路径的有效算法,但它局限于边的权值非负的情况,若图中出现权值为负的边,Dijkstra算法就会失效,求出的最短路径就可能是错的。

    这时候,就需要使用其他的算法来求解最短路径,Bellman-Ford算法就是其中最常用的一个。该算法由美国数学家理查德•贝尔曼(Richard Bellman, 动态规划的提出者)和小莱斯特•福特(Lester Ford)发明。

    适用条件&范围:

    单源最短路径(从源点s到其它所有顶点v);

    有向图&无向图(无向图可以看作(u,v),(v,u)同属于边集E的有向图);

    边权可正可负(如有负权回路输出错误提示);

    差分约束系统;

    Bellman-Ford算法的流程如下:

    给定图G(V, E)(其中V、E分别为图G的顶点集与边集),源点s,数组Distant[i]记录从源点s到顶点i的路径长度,初始化数组Distant[n]为无穷大, Distant[s]为0;

    以下操作循环执行至多n-1次,n为顶点数
    对于每一条边e(u, v),如果Distant[u] + w(u, v) < Distant[v],则另Distant[v] = Distant[u]+w(u, v)。w(u, v)为边e(u,v)的权值;
    若上述操作没有对Distant进行更新,说明最短路径已经查找完毕,或者部分点不可达,跳出循环。否则执行下次循环;

    为了检测图中是否存在负环路,即权值之和小于0的环路。对于每一条边e(u, v),如果存在Distant[u] + w(u, v) < Distant[v]的边,则图中存在负环路,即是说改图无法求出单源最短路径。否则数组Distant[n]中记录的就是源点s到各顶点的最短路径长度。

    可知,Bellman-Ford算法寻找单源最短路径的时间复杂度为O(V*E).

    Bellman-Ford算法可以大致分为三个部分

    • 第一,初始化所有点。每一个点保存一个值,表示从原点到达这个点的距离,将原点的值设为0,其它的点的值设为无穷大(表示不可达)。

    • 第二,进行循环,循环下标为从1到n-1(n等于图中点的个数)。在循环内部,遍历所有的边,进行松弛计算。

    • 第三,遍历途中所有的边(edge(u,v)),判断是否存在这样情况:
      d(v) > d (u) + w(u,v)
      则返回false,表示途中存在从源点可达的权为负的回路。

    代码:

    //
    //  main.cpp
    //  Bellman-Ford
    //
    //  Created by wasdns on 16/11/20.
    //  Copyright © 2016年 wasdns. All rights reserved.
    //
    
    #include <iostream>
    #include <cstdio>
    #include <algorithm>
    #include <cstdlib>
    using namespace std;
    
    #define maxn 100000001;                 //无穷大
    
    struct Edge                             //存储边的结构体
    {
        int u;
        int v;
        int weight;
    };
    
    Edge edge[10001];                       //存储边的数组
    
    int Fordlen[10001];                     //源点到各个点的最短路径长度
    
    /*
     用于初始化Fordlen数组的函数
        注意:源点到自身最短路径长度为0。
     */
    
    void Initial(int n, int source)
    {
        for (int i = 1; i <= n; i++)
        {
            Fordlen[i] = maxn;
        }
        
        Fordlen[source] = 0;
    }
    
    /*
     用于输入边,初始化图的函数
     由两部分组成:
            (1)调用Initial函数初始化;
            (2)输入边,同时维护Fordlen数组
     */
    
    void CreatGraph(int n, int edgenum, int source)
    {
        Initial(n, source);
        
        for (int i = 1; i <= edgenum; i++)
        {
            cin >> edge[i].u >> edge[i].v >> edge[i].weight;
            
            if (edge[i].u == source) {
                Fordlen[edge[i].v] = edge[i].weight;
            }
            
            if (edge[i].v == source) {
                Fordlen[edge[i].u] = edge[i].weight;
            }
        }
    }
    
    /*
     松弛计算函数
     */
    
    void relax_cal(int u, int v, int weight)
    {
        if (Fordlen[u] + weight < Fordlen[v]) {
            Fordlen[v] = Fordlen[u] + weight;
        }
    }
    
    /*
     Bellman-Ford算法主体:
        (1)初始化、建图;
        (2)松弛计算,维护Fordlen数组;
        (3)判断是否有负权环,有的话返回false。
     */
    
    bool Alg_Bellman_Ford(int n, int edgenum, int source)
    {
        CreatGraph(n, edgenum, source); //(1)初始化
        
        for (int i = 1; i <= n-1; i++)  //(2)循环最多执行n-1次,每一次遍历所有边
        {
            for (int j = 1; j <= edgenum; j++)
            {
                relax_cal(edge[j].u, edge[j].v, edge[j].weight);
            }
        }
        
        bool flag = true;               //(3)判断是否有负权边
        
        for (int i = 1; i <= edgenum; i++)
        {
            if (Fordlen[edge[i].u] + edge[i].weight < Fordlen[edge[i].v]) {
                
                flag = false;
                
                break;
            }
        }
        
        return flag;
    }
    
    /*
     输出函数:
        (1)如果没有负权环,输出Fordlen数组
        (2)有负权环,输出"The graph has negative circle."
     */
    
    void Print(bool flag, int n)
    {
        if (flag)
        {
            for (int i = 1; i <= n; i++)
            {
                cout << Fordlen[i] << endl;
            }
        }
        
        else cout << "The graph has negative circle." << endl;
    }
    
    int main()
    {
        int n, edgenum, source;
        
        cin >> n >> edgenum >> source;
        
        bool flag = Alg_Bellman_Ford(n, edgenum, source);
        
        Print(flag, n);
        
        return 0;
    }
    
    

    测试:

    一:

    4 6 1
    1 2 20
    1 3 5
    4 1 -200
    2 4 4
    4 2 4
    3 4 2
    

    结果:

    The graph has negative circle.
    

    二:

    4 6 1
    1 2 2
    1 3 5
    4 1 10
    2 4 4
    4 2 4
    3 4 2
    

    结果:

    0
    2
    5
    6
    

    2016/11/20

  • 相关阅读:
    Ubuntu下deb包的安装方法
    IP地址的分类——a,b,c 类是怎样划分的
    GDB十分钟教程
    【Github教程】史上最全github用法:github入门到精通
    java实现各种数据统计图(柱形图,饼图,折线图)
    MyEclipse 10官网下载地址以及破解方法
    二分查找
    游起来吧!超妹!
    [Python]linux自己定义Python脚本命令
    python无私有成员变量
  • 原文地址:https://www.cnblogs.com/qq952693358/p/6084058.html
Copyright © 2011-2022 走看看