SPFA 算法

SPFA(Shortest Path Faster Algorithm)是Bellman-Ford算法的一种队列实现，减少了不必要的冗余计算。

算法大致流程是用一个队列来进行维护。 初始时将源加入队列。 每次从队列中取出一个元素，并对所有与他相邻的点进行松弛，若某个相邻的点松弛成功，则将其入队。 直到队列为空时算法结束。

这个算法，简单的说就是队列优化的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，则有负权环。

以HDU2544为例说明：

【BestCoder Round #2 来了！】7月27号19:00~21:00（赛前30分钟停止注册比赛）

最短路 Time Limit: 5000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) Total Submission(s): 30849    Accepted Submission(s): 13290 Problem Description 在每年的校赛里，所有进入决赛的同学都会获得一件很漂亮的t-shirt。但是每当我们的工作人员把上百件的衣服从商店运回到赛场的时候，却是非常累的！所以现在他们想要寻找最短的从商店到赛场的路线，你可以帮助他们吗？   Input 输入包括多组数据。每组数据第一行是两个整数N、M（N<=100，M<=10000），N表示成都的大街上有几个路口，标号为1的路口是商店所在地，标号为N的路口是赛场所在地，M则表示在成都有几条路。N=M=0表示输入结束。接下来M行，每行包括3个整数A，B，C（1<=A,B<=N,1<=C<=1000）,表示在路口A与路口B之间有一条路，我们的工作人员需要C分钟的时间走过这条路。 输入保证至少存在1条商店到赛场的路线。   Output 对于每组输入，输出一行，表示工作人员从商店走到赛场的最短时间   Sample Input 2 1 1 2 3 3 3 1 2 5 2 3 5 3 1 2 0 0   Sample Output 3 2

#include<stdio.h>
#include<queue>
#include<string.h>
#define inf 999999999
using namespace std;
struct node     //邻接表存变的起始，权值，和和连接起始位置相同的前一条边
{
	int u,v,w,next;
}stu[20000];
int t;
int head[110],use[110],dis[110];
void add(int u,int v,int w)    
{
	stu[t].v=v;
	stu[t].w=w;
	stu[t].u=u;
	stu[t].next=head[u];      <span style="font-family: Arial, Helvetica, sans-serif;">//每个起始点第一次出现时所对应的head[u]都为0，表示此为以u为起的第一条边；</span>
	head[u]=t++;              //head[u]=t++,t实际上是边的序号，这样处理后每个stu[t].next指的是上一条已u为起点的边；
}
void spfa(int x)
{
	int i;
	queue<int>q;
	memset(use,0,sizeof(use));
    q.push(x);
	dis[x]=0;
	use[x]=1;
	while(!q.empty())
	{
		int u=q.front();
		q.pop();
		use[u]=0;
		for(i=head[u];i!=0;i=stu[i].next)     从head[u]找到最后一条以u为起点的边，stu[i].next找前一条边
		{
			int v=stu[i].v;
			if(dis[v]>dis[u]+stu[i].w)   //边值被更新后才入队
			{
				dis[v]=dis[u]+stu[i].w;
				if(use[v]==0)
				{
					use[v]=1;
					q.push(v);
				}
			}
		}
	}
}
int main()
{
	int n,m,i,j,k;
	while(scanf("%d%d%",&n,&m),(m+n))
	{
		memset(head,0,sizeof(head));
		memset(use,0,sizeof(use));
		for(i=1;i<=n;i++)
			dis[i]=inf;
		t=1;
		while(m--)
		{
			int a,b,c;
			scanf("%d%d%d",&a,&b,&c);
			add(a,b,c);
			add(b,a,c);

		}
		spfa(1);
		printf("%d\n",dis[n]);
	}
	return 0;
}