如果给定一个图是一个DAG,就可以按拓扑排序给顶点编号,并利用这条递推关系式计算出d。但是,如果图中有圈,就无法依赖这样的顺序进行计算。在这种情况下,记当前顶点i的最短路长度为d[i],并 设初值 d [ s ] = 0 , d [ i [ = I N F ( 足 够 大 的 常 数 ) d[s]=0, d[i[=INF(足够大的常数) d[s]=0,d[i[=INF(足够大的常数),再不断使用这条递推关系式更新d的值,就可以计算出新的d。只要图中不存在负圈,这样的更新操作就是有限的。结束之后的d就是所求的最短距离了。
// 从顶点from指向顶点to的权值为cost的边
struct edge { int from, to, cost; };
edge es[MAX_E]; // 边
int d[MAX_V]; // 最短距离
int V, E; // V 是顶点数,E是边数
// 求解从顶点s出发到所有点的最短距离
void shortest_path(int s)
{
for(int i=0; i < V; i++) d[i] = INF;
d[s] = 0;
while(true)
{
bool update = false;
for(int i=0;i<E;i++)
{
edge e = es[i];
if(d[e.from] != INF && d[e.to] > d[e.from] + e.cost)
{
d[e.to] = d[e.from] + e.cost;
update = true;
}
}
if(!update) break;
}
}
如果在图中不存在从s可达的负圈,那么最短路不会经过同一个顶点两次(也就是说,最多通过|V|-1条边),while(true)
的循环最多执行|V|-1次,因此,复杂度是
O
(
∣
V
∣
×
∣
E
∣
)
O(|V| \times |E|)
O(∣V∣×∣E∣)。反之,如果存在从s可达的负圈,那么在第|V|次循环中也会更新d的值,因此也可以用这个性质来检查负圈。如果一开始对所有的顶点i,都把d[i]初始化为0,那么可以检查出所有的负圈。
// 如果返回true则存在负圈
bool find_negative_loop()
{
memset(d, 0, sizeof(d));
for(int i=0;i<V;i++)
for(int j=0;j<E;j++)
{
edge e = es[j];
if(d[e.to] > d[e.from] + e.cost)
{
d[e.to] = d[e.from] + e.cost;
// 如果第n次仍然更新了,则存在负圈
if (i == V - 1) return true;
}
}
return false;
}
文章出自:《挑战程序设计竞赛第2版》2.5节