SPFA算法求单源最短路径
很多时候,给定的图存在负权边,这时类似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]); } } }下载源码: