算法流程
算法大致流程是用一个队列来进行维护。初始时将源加入队列。每次从队列中取出一个元素,并对所有与他相邻的点进行松弛,若某个相邻的点松弛成功,则将其入队。直到队列为空时算法结束。
这个算法,简单的说就是队列优化的bellman-ford,利用了每个点不会更新次数太多的特点发明的此算法
SPFA——Shortest Path Faster Algorithm,它可以在O(kE)的时间复杂度内求出源点到其他所有点的最短路径,可以处理负边。SPFA的实现甚至比Dijkstra或者Bellman_Ford还要简单:
设Dist代表S到I点的当前最短距离,Fa代表S到I的当前最短路径中I点之前的一个点的编号。开始时Dist全部为+∞,只有Dist[S]=0,Fa全部为0。
维护一个队列,里面存放所有需要进行迭代的点。初始时队列中只有一个点S。用一个布尔数组记录每个点是否处在队列中。
每次迭代,取出队头的点v,依次枚举从v出发的边v->u,设边的长度为len,判断Dist[v]+len是否小于 Dist[u],若小于则改进Dist[u],将Fa[u]记为v,并且由于S到u的最短距离变小了,有可能u可以改进其它的点,所以若u不在队列中,就将它放入队尾。这样一直迭代下去直到队列变空,也就是S到所有的最短距离都确定下来,结束算法。若一个点入队次数超过n,则有负权环。
SPFA 在形式上和宽度优先搜索非常类似,不同的是宽度优先搜索中一个点出了队列就不可能重新进入队列,但是SPFA中一个点可能在出队列之后再次被放入队列,也就是一个点改进过其它的点之后,过了一段时间可能本身被改进,于是再次用来改进其它的点,这样反复迭代下去。设一个点用来作为迭代点对其它点进行改进的平均次数为k,有办法证明对于通常的情况,k在2左右。
由于一般不存在负权环,所以我就不写出现负权的问题了。
本代码采用邻接矩阵的方法储存图,方便理解。
其实我也不知道fa数组是干嘛的。。。。。。。。。
#include <stdio.h>
#include <string.h>
#include <algorithm>
#include <queue>
#define maxn 1000
#define INF 1000000
using namespace std;
int dist[maxn+2],fa[maxn+2],bool_n[maxn+2],use[maxn+2];
int grap[maxn+2][maxn+2];
int n;
queue<int>q;
void clear()
{
memset(fa,0,sizeof(fa));
memset(bool_n,0,sizeof(bool_n));
memset(dist,INF,sizeof(dist));
memset(grap,0,sizeof(grap));
dist[1]=0;
}
void spfa(int x)
{
int m,i;
q.push(x);
use[x]=1;
while(!q.empty())
{
m=q.front();
bool_n[m]=0;
q.pop();
for(i=1;i<=n;i++)
{
if(i!=m && grap[m][i])
{
if(dist[i]>dist[m]+grap[m][i])
{
dist[i]=dist[m]+grap[m][i];
fa[i]=m;
if(!bool_n[i])
q.push(i),bool_n[i]=1;
}
}
}
}
}
int main()
{
int m,x,y,z,i;
while(~scanf("%d%d",&n,&m))
{
clear();
for(i=0;i<m;i++)
{
scanf("%d%d%d",&x,&y,&z);
grap[x][y]=z;
grap[y][x]=z;
}
spfa(1);
while(scanf("%d",&x),x)
printf("%d\n",dist[x]);
}
return 0;
}