单源最短路(SPFA算法)

时间:2023-02-13 09:36:39

原文衔接:http://blog.csdn.net/runninghui/article/details/8895586


解决存在负环的图的单源最短路径,bellman-ford算法是比较经典的一个,但是大家都知道,这个算法的效率并不咋的,因为它只知道要求单源最短路,至多做|v|(j图的结点数)次松弛操作,感觉有点盲目吧,这里介绍一个有西南交通大学段凡丁1994年发明的一个算法即SPFA,很大程度上优化了bellman-ford算法(建议没有学过的,先去了解一下这个算法),算法的时间效率我就不说了,因为我觉得当我们熟悉某个算法之后,分析时间复杂度就没什么问题了,如果盲目的记忆,意义不大。

SPFA算法的精妙之处在于不是盲目的做松弛操作,而是用一个队列保存当前做了松弛操作的结点。只要队列不空,就可以继续从队列里面取点,做松弛操作,想想bellman-ford算法吧,它只知道做|v|次循环就对了。下面讲讲SPFA为什么这样做呢?还是举个例子:

单源最短路(SPFA算法)

当前源点1在队列里面,于是我们取了1结点来做对图进行松弛操作,显然这个时候2,3结点的距离更新了,入了队列,我们假设他们没入队列,即现在队列已经空了,那么还有没有必要继续做松弛操作呢?显然没必要了啊,因为源点1要到其他结点必须经过2或3结点啊。现在懂了吧。

先讲一下SPFA的大致思想

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

判断有无负环:如果某个点进入队列的次数超过V次则存在负环(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算法。个人觉得LLL优化每次要求平均值,不太好,为了简单,我们可以之间用C++STL里面的优先队列来进行SLF优化

下面贴一下用用优先队列来进行SLF优化代码

[cpp] view plain copy print?
  1. #include <iostream>  
  2. #include <cstring>  
  3. #include <queue>  
  4. using namespace std;  
  5. const int INF = 0x3fffffff;  
  6. const int MAX = 100;  
  7. int map[MAX][MAX];  
  8. int dis[MAX];  
  9. bool vis[MAX];  
  10. int num[MAX];//记录每个结点入队的次数  
  11. struct cmp  
  12. {  
  13.      bool operator()(int x,int y)  
  14.      {  
  15.           return x>y;  
  16.      }  
  17. };  
  18. bool SPFA(int s0,int n)  
  19. {  
  20.      priority_queue<int,vector<int>,cmp> q;  
  21.      memset(vis,false,sizeof(vis));  
  22.      memset(num,0,sizeof(num));  
  23.      for(int i=0;i<n;i++)  
  24.           dis[i] = INF;  
  25.      dis[s0] = 0;  
  26.      q.push(s0);  
  27.      vis[s0] = true;  
  28.      num[s0]++;  
  29.      while(!q.empty())  
  30.      {  
  31.           int p = q.top();  
  32.           q.pop();  
  33.           for(int i=0;i<n;i++)  
  34.           {  
  35.                if(dis[p]+map[p][i]<dis[i])  
  36.                {  
  37.                     dis[i] = dis[p]+map[p][i];  
  38.                     if(!vis[i])  
  39.                     {  
  40.                          q.push(i);  
  41.                          num[i]++;  
  42.                          if(num[i]>n)//存在负环  
  43.                          {  
  44.                               return false;  
  45.                          }  
  46.                          vis[i]=true;  
  47.                     }  
  48.                }  
  49.           }  
  50.           vis[p] = false;  
  51.      }  
  52.      return true;  
  53. }