SPFA算法 百科内容来自于: 百度百科

SPFA原理

最短路问题

最短路问题

求单源最短路的SPFA算法的全称是:Shortest Path Faster Algorithm,是西南交通大学段凡丁于1994年发表的。从名字我们就可以看出,这种算法在效率上一定有过人之处。很多时候,给定的图存在负权边,这时类似Dijkstra等算法便没有了用武之地,而Bellman-Ford算法的复杂度又过高,SPFA算法便派上用场了。简洁起见,我们约定加权有向图G不存在负权回路,即最短路径一定存在。如果某个点进入队列的次数超过N次则存在负环(SPFA无法处理带负环的图)。当然,我们可以在执行该算法前做一次拓扑排序,以判断是否存在负权回路,但这不是我们讨论的重点。我们用数组d记录每个结点的最短路径估计值,而且用邻接表来存储图G。我们采取的方法是动态逼近法:设立一个先进先出的队列用来保存待优化的结点,优化时每次取出队首结点u,并且用u点当前的最短路径估计值对离开u点所指向的结点v进行松弛操作,如果v点的最短路径估计值有所调整,且v点不在当前的队列中,就将v点放入队尾。这样不断从队列中取出结点来进行松弛操作,直至队列空为止。定理: 只要最短路径存在,上述SPFA算法必定能求出最小值。
证明:每次将点放入队尾,都是经过松弛操作达到的。换言之,每次的优化将会有某个点v的最短路径估计值d[v]变小。所以算法的执行会使d越来越小。由于我们假定图中不存在负权回路,所以每个结点都有最短路径值。因此,算法不会无限执行下去,随着d值的逐渐变小,直到到达最短路径值时,算法结束,这时的最短路径估计值就是对应结点的最短路径值。
期望的时间复杂度:O(ke), 其中k为所有顶点进队的平均次数,可以证明k一般小于等于2。
对SPFA的一个很直观的理解就是由无权图的BFS转化而来。在无权图中,BFS首先到达的顶点所经历的路径一定是最短路(也就是经过的最少顶点数),所以此时利用数组记录节点访问可以使每个顶点只进队一次,但在带权图中,最先到达的顶点所计算出来的路径不一定是最短路。一个解决方法是放弃数组,此时所需时间自然就是指数级的,所以我们不能放弃数组,而是在处理一个已经在队列中且当前所得的路径比原来更好的顶点时,直接更新最优解。
(“算法编程后实际运算情况表明m一般没有超过2n.事实上顶点入队次数m是一个不容易事先分析出来的数,但它确是一个随图的不同而略有不同的常数.所谓常数,就是与e无关,与n也无关,仅与边的权值分布有关.一旦图确定,权值确定,原点确定,m就是一个确定的常数.所以SPFA算法复杂度为O(e)(证毕)"——SPFA的论文
但事实上这个证明是非常不严谨甚至错误的,事实上在bellman算法的论文中已有这方面的内容,所以国际上一般不承认SPFA算法)
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%。 在实际的应用中SPFA的算法时间效率不是很稳定,为了避免最坏情况的出现,通常使用效率更加稳定的Dijkstra算法。

伪代码

SPFA实际上是Bellman-Ford基础上的队列优化
SPFA(G,w,s)
1. INITIALIZE-SINGLE-SOURCE(G,s)
2. INITIALIZE-QUEUE(Q)
3. ENQUEUE(Q,s)
4. While Not EMPTY(Q)
5. Do u<-DLQUEUE(Q)
6. For 每条边(u,v) in E[G]
7. Do tmp<-d[v]
8. Relax(u,v,w)
9. If (d[v] < tmp) and (v不在Q中)
10. ENQUEUE(Q,v)
一种更容易读懂的伪代码:
ProcedureSPFA;
Begin
initialize-single-source(G,s);
initialize-queue(Q);
enqueue(Q,s);
whilenotempty(Q)dobegin
u:=dequeue(Q);
foreachv∈adj[u]dobegin
tmp:=d[v];
relax(u,v);
if(tmp<>d[v])and(notvinQ)thenenqueue(Q,v);
end;
end;
End;

C++代码

#include<iostream>
#include<deque>
#include<stack>
#include<vector>
usingnamespacestd;

constintMAXN=100;
constintINF=0x7FFFFFFF;

structedge
{
intto;
intweight;
};

vector<edge>adjmap[MAXN];//邻接表
boolin_queue[MAXN];//顶点是否在队列中
intin_sum[MAXN];//顶点入队次数
intdist[MAXN];//源点到各点的最短路径
intpath[MAXN];//存储到达i的前一个顶点
intnodesum;//顶点数
intedgesum;//边数

boolSPFA(intsource)
{
deque<int>dq;
inti,j,x,to;
for(i=1;i<=nodesum;i++)
{
in_sum[i]=0;
in_queue[i]=false;
dist[i]=INF;
path[i]=-1;
}
dq.push_back(source);
in_sum[source]++;
dist[source]=0;
in_queue[source]=true;
//初始化完成

while(!dq.empty())
{
x=dq.front();
dq.pop_front();
in_queue[x]=false;
for(i=0;i<adjmap[x].size();i++)
{
to=adjmap[x][i].to;
if((dist[x]<INF)&&(dist[to]>dist[x]+adjmap[x][i].weight))
{
dist[to]=dist[x]+adjmap[x][i].weight;
path[to]=x;
if(!in_queue[to])
{
in_queue[to]=true;
in_sum[to]++;
if(in_sum[to]==nodesum)
returnfalse;
if(!dq.empty())
{
if(dist[to]>dist[dq.front()])
dq.push_back(to);
else
dq.push_front(to);
}
else
dq.push_back(to);
}
}
}
}
returntrue;
}

voidPrint_Path(intx)
{
stack<int>s;
intw=x;
while(path[w]!=-1)
{
s.push(w);
w=path[w];
}
cout<<"顶点1到顶点"<<x<<"的最短路径长度为:"<<dist[x]<<endl;
cout<<"所经过的路径为:1";
while(!s.empty())
{
cout<<s.top()<<"";
s.pop();
}
cout<<endl;
}

intmain()
{
inti,s,e,w;
edgetemp;
cout<<"输入顶点数和边数:";
cin>>nodesum>>edgesum;
for(i=1;i<=nodesum;i++)
adjmap[i].clear();//清空邻接表
for(i=1;i<=edgesum;i++)
{
cout<<"输入第"<<i<<"条边的起点、终点还有对应的权值:";
cin>>s>>e>>w;
temp.to=e;
temp.weight=w;
adjmap[s].push_back(temp);
}
if(SPFA(1))
{
for(i=2;i<=nodesum;i++)
Print_Path(i);
}
else
cout<<"图中存在负权回路"<<endl;
return0;
}


pascal代码

 const
 maxp=10000;{最大结点数}
 var{变量定义}
 p,c,s,t:longint;{p,结点数;c,边数;s:起点;t:终点}
 a,b:array[1..maxp,0..maxp]oflongint;{a[x,y]存x,y之间边的权;b[x,c]存与x相连的第c个边的另一个结点y}
 d:array[1..maxp]ofinteger;{队列}
 v:array[1..maxp]ofboolean;{是否入队的标记}
 dist:array[1..maxp]oflongint;{到起点的最短路}
 head,tail:longint;{队首/队尾指针}
 procedureinit;
 vari,x,y,z:longint;
 begin
 read(p,c);
 fori:=1tocdo
 begin
 readln(x,y,z);{x,y:一条边的两个结点;z:这条边的权值}
 inc(b[x,0]);b[x,b[x,0]]:=y;a[x,y]:=z;{b[x,0]:以x为一个结点的边的条数}
 inc(b[y,0]);b[y,b[y,0]]:=x;a[y,x]:=z;
 end;
 readln(s,t);{读入起点与终点}
 end;
 procedurespfa(s:longint);{SPFA}
 vari,j,now,sum:longint;
 begin
 fillchar(d,sizeof(d),0);
 fillchar(v,sizeof(v),false);
 forj:=1topdodist[j]:=maxlongint;
 dist[s]:=0;v[s]:=true;d[1]:=s;{队列的初始状态,s为起点}
 head:=1;tail:=1;
 whilehead<=taildo{队列不空}
 begin
 now:=d[head];{取队首元素}
 fori:=1tob[now,0]do
 ifdist[b[now,i]]>dist[now]+a[now,b[now,i]]then
 begin
 dist[b[now,i]]:=dist[now]+a[now,b[now,i]];{修改最短路}
 ifnotv[b[now,i]]then{扩展结点入队}
 begin
 inc(tail);
 d[tail]:=b[now,i];
 v[b[now,i]]:=true;
 end;
 end;
 v[now]:=false;{释放结点,一定要释放掉,因为这节点有可能下次用来松弛其它节点}
 inc(head);{出队}
 end;
 end;
 procedureprint;
 begin
 writeln(dist[t]);
 end;
 begin
 init;
 spfa(s);
 print;
 end.

$firstVoiceSent
- 来自原声例句
小调查
请问您想要如何调整此模块?

感谢您的反馈,我们会尽快进行适当修改!
进来说说原因吧 确定
小调查
请问您想要如何调整此模块?

感谢您的反馈,我们会尽快进行适当修改!
进来说说原因吧 确定