2 求含有负权边的图的单源最短路径——Bellman Ford算法与SPFA算法综合分析(下).docVIP

求含有负权边的图的单源最短路径——Bellman Ford算法与SPFA算法综合分析（下）二、SPFA算法 ??? 前面已经说到，Bellman-Ford存在大量的废运算，即使优化过仍是这样。SPFA作为Bellman-Ford的一个改进，尽可能地减少了废运算，从而提高了效率。 ???? Bellman-Ford是枚举每条路径从而算出结果的，也就是一个点是否需要松弛依赖于其连接的点。SPFA的思想是用一个队列存放所有需要松弛的点，初始时队列中只有源点s，进行计算时，若队列为空则算法结束，否则将队首顶点取出后对其相邻的点基于取出地点进行松弛。如果某邻点松弛成功（改变了d中对应的值）并且那个点不在队列中，则将其入队。在进行以上操作时需要记录每个点进入队列的次数。如果有点进入队列超过|V|次则说明图存在负权回路，退出算法并返回false。 算法核心代码如下：（针对有向图，Pascal语言） type point=^node; node=record d,l:longint; next:point; end; var edge:array[1..maxn] of point; closest,d:array[1..maxn] of longint; function spfa(s:longint):boolean; const inf=100000000; var i,j,k,l,r,d1:longint; p:point; q:array[0..maxn] of longint; f:array[1..maxn] of boolean; t:array[1..maxn] of longint; begin l:=1;r:=1;q[1]:=s; fillchar(f,sizeof(f),false); f[s]:=true; fillchar(t,sizeof(t),0);t[s]:=1; for i:=1 to n do begin d[i]:=inf; closest[i]:=-1; end; d[s]:=0;d1:=0; while l=r do begin i:=q[l mod maxn];inc(l); while d[i]*(r-l+1)d1 do begin inc(r); q[r mod maxn]:=i; i:=q[l mod maxn];inc(l); end; p:=edge[i]; d1:=d1-d[i]; while pnil do begin k:=p^.d;j:=p^.l; if (d[k]d[i]+j) then begin d[k]:=d[i]+j; closest[k]:=i; if not f[k] then begin f[k]:=true; inc(t[k]); if t[k]n then exit(false); if d[k]d[i] then begin inc(r); q[r mod maxn]:=k; end else begin dec(l); q[l mod maxn]:=k; end; d1:=d1+d[k]; end; end; p:=p^.next; end; f[i]:=false; end; exit(true); end; 为了方便计算连到某个点的所有边，此处的edge数组用邻接表表示。edge[i]存放着连着i的所有边的链表。 队列中元素的上限为n，所以用了滚动数组来把p的长度压为n。 SPFA的复杂度很难计算，一般公认为O（kE），k为常数。 三、Bellman-Ford、SPFA、朴素Dijkstra效率实测 ?????? 评测环境：WindowsXP，FreePascal2.40，Pentium（R） Dual-Core CPU T4300@2.10GHz，2G内存 ?????  可以看出，SPFA是优于优化后的Bellman-Ford的，所以尽可能地用SPFA吧，记得滚动那个队列数