很多时候,给定的图存在负权边,这时类似Dijkstra等算法便没有了用武之地,而Bellman-Ford算法的复杂度又过高,SPFA算法便派上用场了。简洁起见,我们约定有向加权图G不存在负权回路,即最短路径一定存在。当然,我们可以在执行该算法前做一次拓扑排序,以判断是否存在负权回路。
我们用数组dist记录每个结点的最短路径估计值,而且用邻接表(或邻接矩阵)来存储图G。采取的方法是动态逼近法:
设立一个先进先出的队列用来保存待优化的结点,优化时每次取出队首结点u,并且用u点当前的最短路径估计值对离开u点所指向的结点v进行松弛操作(看程序注释),如果v点的最短路径估计值有所调整,且v点不在当前的队列中,就将v点放入队尾。这样不断从队列中取出结点来进行松弛操作,直至队列空为止。
定理: 只要最短路径存在,上述SPFA算法必定能求出最小值。
证明:每次将点放入队尾,都是经过松弛操作达到的。换言之,每次的优化将会有某个点v的最短路径估计值dist[v]变小。所以算法的执行会使dist越来越小。由于我们假定图中不存在负权回路, 所以每个结点都有最短路径值。因此,算法不会无限执行下去,随着dist值的逐渐变小,直到到达最短路径值时,算法结束,这时的最短路径估计值就是对应结点的最短路径值。(证毕)
判断有无负环:如果某个点进入队列的次数超过N(图的顶点数)次则存在负环(SPFA无法处理带负环的图)
在一幅图中,我们仅仅知道结点一个节点到另一个节点的最短路径长度是多少,有时候意义不大。这附图如果是地图的模型的话,在算出最短路径长度后,我们总要说明“怎么走”才算真正解决了问题。如何在计算过程中记录下来最短路径是怎么走的,并在最后将它输出呢?
Path[] 数组,Path[i]表示从S到i的最短路径中,结点i之前的结点的编号。注意,是“之前”,不是“之后”。
我们只需要在借助结点u对结点v进行松弛的同时,标记下Path[v] = u,记录的工作就完成了,可递归输出。
// 该注意的是有些点可能重复入队,所以出队的点也要重新置未标记
例(HDU1874)
某省自从实行了很多年的畅通工程计划后,终于修建了很多路。不过路多了也不好,每次要从一个城镇到另一个城镇时,都有许多种道路方案可以选择,而某些方案要比另一些方案行走的距离要短很多。这让行人很困扰。
现在,已知起点和终点,请你计算出要从起点到终点,最短需要行走多少距离。
Input
本题目包含多组数据,请处理到文件结束。
每组数据第一行包含两个正整数N和M(0<N<200,0<M<1000),分别代表现有城镇的数目和已修建的道路的数目。城镇分别以0~N-1编号。
接下来是M行道路信息。每一行有三个整数A,B,X(0<=A,B<N,A!=B,0<X<10000),表示城镇A和城镇B之间有一条长度为X的双向道路。
再接下一行有两个整数S,T(0<=S,T<N),分别代表起点和终点。
Output
对于每组数据,请在一行里输出最短需要行走的距离。如果不存在从S到T的路线,就输出-1.
Sample Input
3 3
0 1 1
0 2 3
1 2 1
0 2
3 1
0 1 1
1 2
Sample Output
2
-1
import java.util.Scanner;
import java.util.Queue;
import java.util.LinkedList;
public class Main{
private int map[][]; //邻接矩阵
private int dist[]; // s到终点的最短路径
//private int[] Path;// 记录前驱点 。若Path[i]=k,表示从s到i的最短路径中i的前一个点是k
static final int INF=10000;
private int n;//顶点个数
private int s;//起点
private int e;//终点
public Main(int n,int s,int e,int[][] map){
this.n=n;
this.s=s;
this.e=e;
this.map=map;
dist = new int[n];
}
public int[] getDist(){
return dist;
}
private void spfa(){
int x;
Queue<Integer> q=new LinkedList<Integer>() ;
boolean visited[]=new boolean[n];
for(int i=0;i<n;i++)
dist [i]=INF;
dist[s]=0;
q.offer(s);
visited[s]=true;
while(!q.isEmpty()){
x=q.poll();
visited[x]=false; // 置出队的点未标记
for(int i=0;i<n;i++)
if(dist[x]+map[x][i]<dist[i]) //这里就是所谓的松弛操作了
{
dist[i]=map[x][i]+dist[x]; //更新路径
//Path[i]=x;
if(!visited[i]) // 未被访问
{
q.offer(i);
visited[i]=true;
}
}
}
}
public static void main(String[] args) {
Scanner sc = new Scanner(System.in);
int n,m,s,e,u,v,w;
int[][] map;
while (sc.hasNext()) {
n = sc.nextInt();//城镇数目即图的顶点数
m = sc.nextInt();//边数
map = new int[n][n];
for (int i = 0; i < n; i++)
for (int j = 0; j < n; j++)
map[i][j] = INF;
while (m-- > 0) {
u = sc.nextInt();
v = sc.nextInt();
w = sc.nextInt();
if (map[u][v] > w || map[v][u] > w) {//注意是双向的
map[u][v] = w;
map[v][u] = w;
}
}
s = sc.nextInt();
e = sc.nextInt();
Main ma=new Main(n,s,e,map);
ma.spfa();
int[] dist=ma.getDist();
if (dist[e] == INF)
System.out.println(-1);
else System.out.println(dist[e]);
}
}
}
下载源码: