算法描述》关于SPFA和Dijkstra算法的两三事

  本来我是想把这两个算法分开写描述的,但是SPFA其实就是Dijkstra的稀疏图优化,所以其实代码差不多,所以就放在一起写了。

  因为SPFA是Dijkstra的优化,所以我想来讲讲Dijkstra。

  什么是Dijkstra

  Dijkstra是一种求单源最短路的基础算法,时间复杂度在不加堆优化的情况下是o(n^2)的,加了堆优化就能简化到o(nlogn),而且算法稳定性很强(从这点上来说比SPFA好多了,具体怎么好下面再讲),基础思路如下:

  首先,把所有点到源的距离设为最大,然后把源加入队列,接着每次对队列首做这种操作:用邻接表找到从这个点能链到的所有所有点,接着要按大小顺序找当前距离和当前队首最近的点(这个按大小顺序就是可以用堆优化的地方),做松弛操作,松弛成功,将这个点入队,否则按大小顺序找下一条边,然后重复做松弛操作,直到当前点所有能链到的边松弛过,弹出当前点,重复以上操作,直到队列里没有其他元素。结束程序,输出单源到某点距离。

  加了堆优化之后,有如此6的复杂度的Dijkstra,一般来说就是稳定单源最短路的最佳程序,但是在全图有负环的时候,Dijkstra算法就显得很无力了,在这时,我们就要用SPFA这种神奇的乱搞算法了。

  什么是SPFA

  SPFA算法全称是Shortest Path Faster Algorithm算法,人家都叫“Faster”了,那肯定有过人之处,基础思路如下:

  首先,还是把所有点到源的距离设为最大,把源加入队列,接着每次对队首做这种操作:用邻接表找到从这个点能链到的所有所有点,然后对每条边做松弛操作,只要松弛成功,我们就判断当前松弛成功的点是否在队列,如果不在,就入队,否则就只是做松弛,然后弹出当前点,对下一点再次做以上操作,直到队列为空。

  当然,和Dijkstra相比,SPFA看起来相当的乱搞,但是,对于绝大多数(稀疏)图来说,加了SLF和LLL两种优化的SPFA肯定比Dijkstra快得多(SLF和LLL就是两种决策优化,前者对于最值优化,后者是均值优化)。

  顺带附上一道codevs1021的SPFA题目与代码

  麦克找了个新女朋友,玛丽卡对他非常恼火并伺机报复。

    因为她和他们不住在同一个城市,因此她开始准备她的长途旅行。

    在这个国家中每两个城市之间最多只有一条路相通,并且我们知道从一个城市到另一个城市路上所需花费的时间。

    麦克在车中无意中听到有一条路正在维修,并且那儿正堵车,但没听清楚到底是哪一条路。无论哪一条路正在维修,从玛丽卡所在的城市都能到达麦克所在的城市。

    玛丽卡将只从不堵车的路上通过,并且她将按最短路线行车。麦克希望知道在最糟糕的情况下玛丽卡到达他所在的城市需要多长时间,这样他就能保证他的女朋友离开该城市足够远。

编写程序,帮助麦克找出玛丽卡按最短路线通过不堵车道路到达他所在城市所需的最长时间(用分钟表示)。

输入描述 Input Description

第一行有两个用空格隔开的数N和M,分别表示城市的数量以及城市间道路的数量。1≤N≤1000,1≤M≤N*(N-1)/2。城市用数字1至N标识,麦克在城市1中,玛丽卡在城市N中。

接下来的M行中每行包含三个用空格隔开的数A,B和V。其中1≤A,B≤N,1≤V≤1000。这些数字表示在A和城市B中间有一条双行道,并且在V分钟内是就能通过。

 

输出描述 Output Description

   输出文件的第一行中写出用分钟表示的最长时间,在这段时间中,无论哪条路在堵车,玛丽卡应该能够到达麦克处,如果少于这个时间的话,则必定存在一条路,该条路一旦堵车,玛丽卡就不能够赶到麦克处。

样例输入 Sample Input

5 7

1 2 8

1 4 10

2 3 9

2 4 10

2 5 1

3 4 7

3 5 10

样例输出 Sample Output

27

  题解如下

  这道题关键在于最优解可能会无法连通,所以我们求出最优解之后,要枚举最优解中哪条路会断,取其中最坏情况(题目所需),得答案就好,因为算次优解一定不能使数据再包括最优解的某条边,不然就和最优解一样了,所以以上方法是正确的。

  贴出代码

 1 #include<stdio.h>
 2 struct shit{
 3     int aim;
 4     int get;
 5     int next;
 6     int lon;
 7 }e[3010000];
 8 int max(int x,int y)
 9 {
10     return x>y?x:y;
11 }
12 int pre[3010000];
13 int d[1010000],a,b,n,m,num,star,ass,quq[60100000],point,head[1010000],ans;
14 bool f[10101000];
15 void fuck(int x,int y,int s)//建边
16 {
17     e[++point].aim=x;
18     e[point].lon=s;
19     e[point].get=y;
20     e[point].next=head[y];
21     head[y]=point;
22 }
23 int main()
24 {
25     scanf("%d%d",&n,&m);
26     for(int i=1;i<=m;i++)
27     {
28         scanf("%d%d%d",&a,&b,&num);
29         fuck(a,b,num);
30         fuck(b,a,num);
31     }
32     quq[++star]=1;
33     ass=1;
34     f[1]=true;
35     for(int i=2;i<=n;i++)
36         d[i]=214748364;
37     while(star<=ass)
38     {
39         for(int k=head[quq[star]];k!=0;k=e[k].next)
40             if(d[quq[star]]+e[k].lon<d[e[k].aim])
41             {
42                 d[e[k].aim]=d[quq[star]]+e[k].lon;
43                 pre[e[k].aim]=k;//第一次的SPFA要记录路径(点的位置)
44                 if(!f[e[k].aim])
45                 {
46                     quq[++ass]=e[k].aim;
47                     f[e[k].aim]=true;
48                 }
49             }
50         f[quq[star++]]=false;
51     }
52     ans=max(d[n],ans);
53     int q=n;
54     while(q!=0)
55     {
56     a=e[pre[q]].lon;
57     e[pre[q]].lon=214748364;//直接设置当前边数值极大(即设置为不连通)
58     if(pre[q]%2==1)//根据读入边的方法来双向删边
59         e[pre[q]+1].lon=214748364;
60     else
61         e[pre[q]-1].lon=214748364;
62     star=1;
63     for(int i=1;i<=ass;i++)
64     quq[i]=0;
65     quq[star]=1;
66     ass=1;
67     f[1]=true;
68     for(int i=2;i<=n;i++)
69     {
70         d[i]=214748364;
71         f[i]=false;
72     }
73     while(star<=ass)
74     {
75         for(int k=head[quq[star]];k!=0;k=e[k].next)
76             if(d[quq[star]]+e[k].lon<d[e[k].aim])
77             {
78                 d[e[k].aim]=d[quq[star]]+e[k].lon;
79                 if(!f[e[k].aim])
80                 {
81                     quq[++ass]=e[k].aim;
82                     f[e[k].aim]=true;
83                 }
84             }
85         f[quq[star++]]=false;
86     }
87     if(d[n]!=214748364)ans=max(d[n],ans);//如果数值就为我们设置的最大,即切去当前边的图无法连通
88     e[pre[q]].lon=a;
89     if(pre[q]%2==1)//复原切去的边
90         e[pre[q]+1].lon=a;
91     else
92         e[pre[q]-1].lon=a;
93         q=e[pre[q]].get;
94     }
95     printf("%d",ans);
96     return 0;
97 }

  PS:没有加LLL和SLF的主要原因是我还不会打→_→

时间: 2024-10-26 06:02:14

算法描述》关于SPFA和Dijkstra算法的两三事的相关文章

算法导论-第24章 Dijkstra算法

Dikstra算法解决的是有向图上单源最短路径问题(无向图可以看成有相反的两条有向边),且要求边的权重都是非负值. 算法导论用了很多引理,性质来证明Dijstra算法的正确性,这里不说了,也表达不明白,只说我理解的过程. 有一个图G( V,E) ,选定一个源点s,维护一个集合Q=V-s,  Q中点有一个d值表示此时从s到该点的已知距离,s.d=0 :初始化都为正无穷,表明不可达.然后对s点所连接的点(设为点集M)进行松弛操作,就是设点m属于M, m.d > s.d+ w(s,m) 则更新 m.d

python数据结构与算法——图的最短路径(Dijkstra算法)

1 # Dijkstra算法——通过边实现松弛 2 # 指定一个点到其他各顶点的路径——单源最短路径 3 4 # 初始化图参数 5 G = {1:{1:0, 2:1, 3:12}, 6 2:{2:0, 3:9, 4:3}, 7 3:{3:0, 5:5}, 8 4:{3:4, 4:0, 5:13, 6:15}, 9 5:{5:0, 6:4}, 10 6:{6:0}} 11 12 13 # 每次找到离源点最近的一个顶点,然后以该顶点为重心进行扩展 14 # 最终的到源点到其余所有点的最短路径 15

[C++] 多源最短路径(带权有向图):【Floyd算法(动态规划法)】 VS nX Dijkstra算法(贪心算法)

1 Floyd算法 1.1 Code /** * 弗洛伊德(Floyd)算法: 1 解决问题: 多源最短路径问题 求每一对顶点之间的最短路径 Background: 带权有向图 2 算法思想: 动态规划(DP, Dynamic Programming) 3 时间复杂度: O(n^3) */ #include<stdio.h> #include<iostream> using namespace std; // 1 定义图模型(邻接矩阵表示法)的基本存储结构体 # define Ma

【算法】祭奠spfa 最短路算法dijspfa

题目链接 本题解来源 其他链接 卡spfa的数据组 题解堆优化的dijkstra 题解spfa讲解 来自以上题解的图片来自常暗踏阴 使用前向星链表存图 直接用队列优化spfa struct cmp { bool operator()(int a,int b) { return dist[a]>dist[b]; } }; priority_queue<int,vector<int>,cmp> q;void dijspfa() { q.push(s); memset(inq,0,

图论(三) (一)最短路径算法 2.Dijkstra算法

Dijkstra 算法解决的是带权重的有向图上单源最短路径问题,该算法要求所有边的权重都为非负值.该算法的时间复杂度是O(N2),相比于处理无负权的图时,比Bellmad-Ford算法效率更高. 算法描述: 首先引用<算法导论>中的一段比较官方的话,如果可以看懂,那下一部分就可以跳过了: "Dijkstra算法在运行过程中维持的关键信息是一组结点集合S.从源结点s到该集合中每个结点之间的最短路径已经被找到.算法重复从结点集 V - S 中算则最短路径估计的最小的结点 u ,将 u 加

最短路 Dijkstra算法

Dijksitra算法求最短路仅仅适用于不存在右边是负权的情况(Bellman-Ford算法没有这一个限制).主要特点是从起点为中心向外层层扩展,直到扩展到终点为止. 最短路的最优子结构性质 即一个最短路路径中经过的所有点这条路均是其最短路.(反证法易证) Dijkstra基本思路: ①找到最短距离已经确定的顶点,从它出发更新相邻顶点的最短距离 ②此后不需要再关心1中的"最短距离已经确定的顶点" 在最开始的时候,只有起点的最短距离是确定的.而在尚未使用过的顶点中,距离d[i]最小的顶点

Dijkstra算法(三)之 Java详解

前面分别通过C和C++实现了迪杰斯特拉算法,本文介绍迪杰斯特拉算法的Java实现. 目录 1. 迪杰斯特拉算法介绍 2. 迪杰斯特拉算法图解 3. 迪杰斯特拉算法的代码说明 4. 迪杰斯特拉算法的源码 转载请注明出处:http://www.cnblogs.com/skywang12345/ 更多内容:数据结构与算法系列 目录 迪杰斯特拉算法介绍 迪杰斯特拉(Dijkstra)算法是典型最短路径算法,用于计算一个节点到其他节点的最短路径. 它的主要特点是以起始点为中心向外层层扩展(广度优先搜索思想

【裸单源最短路:Dijkstra算法两种版本】hdu 1874 畅通工程续

Source : hdu 1874 畅通工程续 http://acm.hdu.edu.cn/showproblem.php?pid=1874 Problem Description 某省自从实行了很多年的畅通工程计划后,终于修建了很多路.不过路多了也不好,每次要从一个城镇到另一个城镇时,都有许多种道路方案可以选择,而某些方案要比另一些方案行走的距离要短很多.这让行人很困扰. 现在,已知起点和终点,请你计算出要从起点到终点,最短需要行走多少距离. Input 本题目包含多组数据,请处理到文件结束.

双向广搜的DIJKSTRA算法--简易的北京地铁导航实现

本学期的课程设计,实现最短路的算法,于是采用了DIJKSTRA算法,并用双向广搜优化了. 实现了简易的北京地铁导航.于是把代码分享出来. (核心代码是find_min(),Dijkstra()部分) 转载或者用到里面的代码请注明博主姓名以及出处! (注:只输入了图片里的地铁站信息,所用到的文件最下面有下载,因为这些文件是我和同学一条一条的录入的,所以如果你用到请务必注明这些文件的出处) 代码: /**************************************************