Bellman-Ford算法



贝尔曼-福特算法（Bellman-Ford）是由理查德·贝尔曼（Richard Bellman） 和 莱斯特·福特（英语：Lester Ford） 创立的，求解单源最短路径问题的一种算法。有时候这种算法也被称为 Moore-Bellman-Ford 算法，因为 Edward F. Moore 也为这个算法的发展做出了贡献。它的原理是对图进行V-1次松弛操作，得到所有可能的最短路径。其优于迪科斯彻算法的方面是边的权值可以为负数、实现简单，缺点是时间复杂度过高，高达。但算法可以进行若干种优化，提高了效率。

贝尔曼-福特算法与迪科斯彻算法类似，都以松弛操作为基础，即估计的最短路径值渐渐地被更加准确的值替代，直至得到最优解。在两个算法中，计算时每个边之间的估计距离值都比真实值大，并且被新找到路径的最小长度替代。 然而，迪科斯彻算法以贪心法选取未被处理的具有最小权值的节点，然后对其的出边进行松弛操作；而贝尔曼-福特算法简单地对所有边进行松弛操作，共|V | − 1次，其中 |V |是图的点的数量。在重复地计算中，已计算得到正确的距离的边的数量不断增加，直到所有边都计算得到了正确的路径。这样的策略使得贝尔曼-福特算法比迪科斯彻算法适用于更多种类的输入。

松弛操作伪代码：

Bellman-Ford算法伪代码：

关于实现

//========================================================
// Bellman-Ford算法实现
//========================================================
#include <stdio.h>
#include <stdlib.h>

#define INFINITE 0xff
#define VERTEXS 5
#define EDGES 10

typedef struct edge {
	char u;
	char v;
	int weight;
} EDGE;

// s t x y z
EDGE edges[EDGES] = {{0, 1, 6}, {0, 3, 7}, {1, 2, 5}, {1, 3, 8},
					{1, 4, -4}, {2, 1, -2}, {3, 2, -3}, {3, 4, 9}, 
					{4, 0, 2}, {4, 2, 7}};

int dis[VERTEXS]; // 最短路径估计distance
int prev[VERTEXS]; // v的前驱结点

bool BellmanFord(int s)
{
	// Step 1: initialize graph
	for (int i = 0; i < VERTEXS; i++) {// 对dis，prev进行初始化
		dis[i] = INFINITE;
		prev[i] = 0;
	}

	dis[s] = 0; 
	// Step 2: relax edges repeatedly
	for (int i = 0; i < VERTEXS; i++) {   // 对每条边进行|V|-1次松弛操作
		for (int j = 0; j < EDGES; j++) { // 对每条边进行松弛操作
			if (dis[edges[j].v] > dis[edges[j].u] + edges[j].weight) {
				dis[edges[j].v] = dis[edges[j].u] + edges[j].weight;
				prev[edges[j].v] = edges[j].u;
			}
		}
	}
	
	// Step 3: check for negative-weight cycles
	for (int j = 0; j < EDGES; j++) { // 判断是否有负权回路
		if (dis[edges[j].v] > dis[edges[j].u] + edges[j].weight) 
			return false;
	}

	return true;
}

void PrintPath(int s)
{
	while (s != prev[s]) {
		printf("%c-->", s + 'a');
		s = prev[s];
	}
	if (s == prev[s])
		printf("%c\n", s + 'a');
}

int main()
{
	bool result = BellmanFord(0);
	if (!result) {
		printf("图中有权值为负的环路!\n");
		return -1;
	}

	for (int i = 0; i < VERTEXS; i++) 
		printf("d[%d] = %d, p[%d] = %d\n", i, dis[i], i, prev[i]);

	for (int i = 1; i < VERTEXS; i++)
		PrintPath(i);

	system("pause");
	return 0;
}