SPFA
求单源最短路的SPFA算法的全称是:Shortest Path Faster Algorithm
从名字我们就可以看出,这种算法在效率上一定有过人之处。很多时候,给定的图存在负权边,这时类似Dijkstra等算法便没有了用武之地,而Bellman-Ford算法的复杂度又过高,SPFA算法便派上用场了。
【伪代码】
伪代码 :
function Dijkstra(Graph, source):
dist[source] ← 0 // Distance from source to source
for each vertex v in Graph: // Initialization
if v ≠ source: // Where v has not yet been removed from Q (unvisited
dist[v] ← infinity // Unknown distance function from source to v
end if
add v to Q // All nodes initially in Q (unvisited nodes)
end for
while Q is not empty:
u ← vertex in Q with min dist[u] // Source node in first case
remove u from Q
for each neighbor v of u: // where v is still in Q.
alt ← dist[u] + length(u, v) //松弛操作
if alt < dist[v]: // A shorter path to v has been found
dist[v] ← alt;
end if
end for
end while
return dist[]
这段伪代码可以明白,这里采取的思想是 不断的将s的邻接点加入队列,取出不断的进行松弛操作(替换节点的最短距离),直到队列为空 。
如果一个结点被加入队列超过n-1次,那么显然图中有负环。
这里给出一个简单的例题 HDU2544
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=2544
最短路
Time Limit: 5000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total Submission(s): 42397 Accepted Submission(s): 18551
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条商店到赛场的路线。
输入保证至少存在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
【题目代码】
//SPFA shortest path faster algorithm
#include <cstdio>
#include <vector>
#include <queue>
#define INF 0xfffffff
using namespace std;
struct node{
int v,len;
node(int v=0,int len=0):v(v),len(len){}
};
vector<node>G[110];
int minDist[110];
int inqueue[110];
int n,m;
void init(){
for(int i=0;i<110;i++){
minDist[i]=INF;
inqueue[i]=0;
G[i].clear();
}
}
int Kijkstra(){
inqueue[1]=true;
queue<int >Q;
minDist[1]=0;
Q.push(1);
while(!Q.empty()){
int vex = Q.front();
Q.pop();
inqueue[vex]=0;
for(int i=0;i<G[vex].size();i++){
int v = G[vex][i].v;
if(G[vex][i].len+minDist[vex]<minDist[v]){
minDist[v]=G[vex][i].len+minDist[vex];
if(!inqueue[v])
{
inqueue[v]=1;
Q.push(v);
}
}
}
}
return minDist[n];
}
int main(){
while(scanf("%d%d",&n,&m)!=EOF&&(n||m)){
init();
int a,b,len;
for(int i=0;i<m;i++){
scanf("%d%d%d",&a,&b,&len);
G[a].push_back(node(b,len));
G[b].push_back(node(a,len));
}
int ans=Kijkstra();
printf("%d\n",ans);
}
return 0;
}