我的 Floyd-Warshall C++ 实现中的错误

Posted

技术标签:

【中文标题】我的 Floyd-Warshall C++ 实现中的错误【英文标题】:Bug in my Floyd-Warshall C++ implementation 【发布时间】:2011-03-02 22:00:09 【问题描述】:

我的大学有一个任务,已经成功实施了 Dijkstra 和 Bellman-Ford,但我在这方面遇到了麻烦。一切看起来都很好,但它没有给我正确的答案。

代码如下:

void FloydWarshall()

    //Also assume that n is the number of vertices and edgeCost(i,i) = 0

    int path[500][500];

    /* A 2-dimensional matrix. At each step in the algorithm, path[i][j] is the shortest path
       from i to j using intermediate vertices (1..k−1).  Each path[i][j] is initialized to
       edgeCost(i,j) or infinity if there is no edge between i and j.
    */

    for(int i = 0 ; i <= nvertices ; i++)
        for(int j = 0 ; j <= nvertices ; j++)
            path[i][j] = INFINITY;

    for(int j = 0 ; j < narestas ; j++) //narestas = number of edges
    
        path[arestas[j]->v1][arestas[j]->v2] = arestas[j]->peso; //peso = weight of the edge (aresta = edge)
        path[arestas[j]->v2][arestas[j]->v1] = arestas[j]->peso;
    

    for(int i = 0 ; i <= nvertices ; i++) //path(i, i) = 0
        path[i][i] = 0;

    //test print, it's working fine
    //printf("\n\n\nResultado FloydWarshall:\n");
    //for(int i = 1 ; i <= nvertices ; i++)
    //    printf("distancia ao vertice %d:  %d\n", i, path[1][i]);


    // Here's the problem, it messes up, and even a edge who costs 4, and the minimum is 4, it prints 2.

    //for k = 1 to n
    for(int k = 1 ; k <= nvertices ; k++)
       //for i = 1 to n
       for(int i = 1 ; i <= nvertices ; i++)
           //for j := 1 to n
           for(int j = 1 ; j <= nvertices ; j++)
               if(path[i][j] > path[i][k] + path[k][j])
                   path[i][j] = path[i][k] + path[k][j];


    printf("\n\n\nResultado FloydWarshall:\n");
    for(int i = 1 ; i <= nvertices ; i++)
        printf("distancia ao vertice %d:  %d\n", i, path[1][i]);

我正在使用我制作的这个图表示例:

6 7

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

意味着我们有 6 个顶点(1 到 6)和 7 个边(1,2),权重为 4...等等。

如果有人需要更多信息,我愿意提供,只是厌倦了查看此代码并没有发现错误。

【问题讨论】:

【参考方案1】:

另外,你的迭代的开始和结束不是在几个地方相差一个吗?您可能希望它们从 0 运行到 nvertices-1;即for (int i = 0; i &lt; nvertices; i++)

【讨论】:

【参考方案2】:

没关系,我休息一下吃点东西发现了错误。

无穷大被定义为INT_MAX,所以一旦你向它添加一些东西,它就会变成负数。

我只定义了一些大的东西(对于我的问题,比如超过 9000,没有图形路径会占用更多),它工作正常。

但我可以知道你为什么建议尹吗?我没听懂。

谢谢

【讨论】:

【参考方案3】:
 if(path[i][j] > path[i][k] + path[k][j])
  path[i][j] = path[i][k] + path[k][j];

在这里做一些检查。例如如果 path[i][k] 和 path[k][j] 是非无限的,并且 i!=j i!=k 和 k!=j。

【讨论】:

以上是关于我的 Floyd-Warshall C++ 实现中的错误的主要内容,如果未能解决你的问题,请参考以下文章

Floyd-Warshall 找不到问题所在

在 Haskell 中 Floyd-Warshall 的表现——修复空间泄漏

如何在 Floyd-Warshall 算法中输出最短路径?

Floyd-Warshall算法及其并行化实现(基于MPI)

Floyd-Warshall算法及其并行化实现(基于MPI)

Floyd-Warshall:所有最短路径