c++ spfa+洛谷2296 寻找道路

Hi 大家好今天是2020/8/26蒟蒻来给大家讲spfa

注意：加粗的可以不用看

一.算法简介
SPFA(Shortest Path Faster Algorithm)算法是求单源最短路径的一种算法，它是Bellman-ford（贝尔曼福特）的队列优化，它是一种十分高效的最短路算法。很多时候，给定的图存在负权边，这时类似Dijkstra等算法便没有了用武之地，而Bellman-Ford算法的复杂度又过高，SPFA算法便派上用场了。SPFA的复杂度大约是O(kE),k是每个点的平均进队次数(一般的，k是一个常数，在稀疏图中小于2)。但是，SPFA算法稳定性较差，在稠密图中SPFA算法时间复杂度会退化。
二、原理及证明
SPFA算法的全称是:Shortest Path Faster Algorithm，是西南交通大学段凡丁于 1994 年发表的论文中的名字。不过，段凡丁的证明是错误的，且在 Bellman-Ford 算法提出后不久(1957 年)已有队列优化内容，所以国际上不承认 SPFA 算法是段凡丁提出的。为了避免最坏情况的出现，在正权图上应使用效率更高的Dijkstra算法。若给定的图存在负权边，类似Dijkstra算法等算法便没有了用武之地，SPFA算法便派上用场了。简洁起见，我们约定加权有向图G不存在负权回路，即最短路径一定存在。用数组d记录每个结点的最短路径估计值，而且用邻接表来存储图G。我们采取的方法是动态逼近法:设立一个先进先出的队列用来保存待优化的结点，优化时每次取出队首结点u，并且用u点当前的最短路径估计值对离开u点所指向的结点v进行松弛操作，如果v点的最短路径估计值有所调整，且v点不在当前的队列中，就将v点放入队尾。这样不断从队列中取出结点来进行松弛操作，直至队列空为止。
定理:只要最短路径存在，上述SPFA算法必定能求出最小值。证明:每次将点放入队尾，都是经过松弛操作达到的。换言之，每次的优化将会有某个点v的最短路径估计值d[v]变小。所以算法的执行会使d越来越小。由于我们假定图中不存在负权回路，所以每个结点都有最短路径值。因此，算法不会无限执行下去，随着d值的逐渐变小，直到到达最短路径值时，算法结束，这时的最短路径估计值就是对应结点的最短路径值。实际上，如果一个点进入队列达到n次，则表明图中存在负环，没有最短路径。段凡丁论文中的复杂度证明 (O(kE)，k 是小常数)是错误的，在此略去。该算法的最坏复杂度为 O(VE)。对SPFA的一个很直观的理解就是由无权图的BFS转化而来。在无权图中，BFS首先到达的顶点所经历的路径一定是最短路(也就是经过的最少顶点数)，所以此时利用数组记录节点访问可以使每个顶点只进队一次，但在带权图中，最先到达的顶点所计算出来的路径不一定是最短路。一个解决方法是放弃数组，此时所需时间自然就是指数级的，所以我们不能放弃数组，而是在处理一个已经在队列中且当前所得的路径比原来更好的顶点时，直接更新最优解。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 和 LLF 在随机数据上表现优秀，但是在正权图上最坏情况为 O(VE)，在负权图上最坏情况为达到指数级复杂度。
三、实现方法：建立一个队列，初始时队列里只有起始点，在建立一个表格记录起始点到所有点的最短路径（该表格的初始值要赋为极大值，该点到他本身的路径赋为0）。然后执行松弛操作，用队列里有的点去刷新起始点到所有点的最短路，如果刷新成功且被刷新点不在队列中则把该点加入到队列最后。重复执行直到队列为空。此外，SPFA算法还可以判断图中是否有负权环，即一个点入队次数超过N。
四.算法图解
给定一个有向图，求A~E的最短路。

源点A首先入队，并且AB松弛

A	0
B	INF
C	INF
D	INF
E	INF

queue

A
空
空
空
空

扩展与A相连的边，B，C 入队并松弛。

A	0
B	3
C	2
D	INF
E	INF

queue

空
B
C
空
空

B，C分别开始扩展，D入队并松弛

A	0
B	3
C	2
D	4
E	INF

queue

空
空
空
D
空

D出队，E入队并松弛。

A	0
B	3
C	2
D	4
E	8

queue

空
空
空
空
空
E

E出队，此时队列为空，源点到所有点的最短路已被找到，A->E的最短路即为8

A	0
B	3
C	2
D	4
E	8

queue

空
空
空
空
空
空

五、代码

#include<bits/stdc++.h>
using namespace std;
inline int read(){
   
    int X=0,w=1; char ch=0;
    while(ch<'0' || ch>'9') {
   if(ch=='-') w=-1;ch=getchar();}
    while(ch>='0' && ch<='9') X=(X<<3)+(X<<1)+ch-'0',ch=getchar();
    return X*w;
}
const int inf=1e9;
int S,T,n,m;
int head[5005],cnt;
struct E{
   
    int next,to,val,cost;
}tu[150000];
void add(int next,int to,int val,int cost)