SPFA可以解决权值有负值的图的单源最短路径,但不能存在从源点可达的权为负的回路,能够检测负圈。判断负环:如果某个顶点入队次数超过V次,则说明有负环。
算法期望的时间复杂度为O(k*V),k为所有定点进队的平均次数(在稀疏图中k一般小于2,SPFA算法稳定性较差,在稠密图中SPFA算法时间复杂度会退化)。
SPFA算法是对Bellman-Ford算法的一个优化,使用队列减少不必要的计算。
实现方法:建立一个队列,初始时队列里只有起始点,在建立一个表格记录起始点到所有点的最短路径(该表格的初始值要赋为极大值,该点到他本身的路径赋为0)。然后执行松弛操作,用队列里有的点去刷新起始点到所有点的最短路,如果刷新成功且被刷新点不在队列中则把该点加入到队列最后。重复执行直到队列为空。
SPFA算法有两个优化算法 SLF 和 LLL:
SLF:Small Label First 策略,设要加入的节点是j,队首元素为i,若dist(j)<dist(i),则将j插入队首, 否则插入队尾。
LLL:Large Label Last 策略,设队首元素为i,队列中所有dist值的平均值为x,若dist(i)>x则将i插入 到队尾,查找下
一元素,直到找到某一i使得dist(i)<=x,则将i出对进行松弛操作。
引用网上资料,SLF 可使速度提高 15 ~ 20%;SLF + LLL 可提高约 50%。
代码:
#include <stdio.h>
#include <queue>
#define INF 0xfffffff
#define N 220
using namespace std;
int CityNum, RoadNum;
int map[N][N], dis[N];
bool vis[N]; // 标记已经入队的元素
void Init(){ // 初始化
for(int i = 0; i < CityNum; i++){
for(int j = 0; j < CityNum; j ++){
map[i][j] = INF;
}
dis[i] = INF;
vis[i] = 0;
}
}
void SPFA(int choose){
queue<int>q;
q.push(choose);
dis[choose] = 0;
while(!q.empty()){
choose = q.front();
q.pop();
vis[choose] = 0;
for(int i = 0; i < CityNum; i ++){ // 遍历所有的城市结点
if(dis[i] > dis[choose] + map[choose][i]){
dis[i] = dis[choose] + map[choose][i];
if(!vis[i]){
q.push(i);
vis[i] = 1; // 入队元素标记
}
}
}
}
}
int main()
{
int start, end, len;
while(scanf("%d%d", &CityNum, &RoadNum) != EOF){
Init(); // 初始化
for(int i = 0; i < RoadNum; i ++){ // 建无向图
scanf("%d%d%d", &start, &end, &len);
if(len < map[start][end]){
map[start][end] = map[end][start] = len;
}
}
scanf("%d%d", &start, &end);
SPFA(start);
if(dis[end] == INF){
printf("-1\n");
}
else{
printf("%d\n", dis[end]);
}
}
return 0;
}
从源点开始,依次将路径变短的那个点入队列,再从这个点开始查询是否有相邻的点路径会变短,直到所有点都是最短路径为止。正因为这样 ,如果存在负圈的话将无限进队列,所以SPFA无法处理有负圈的