851. spfa求最短路

给定一个n个点m条边的有向图，图中可能存在重边和自环， 边权可能为负数。

请你求出1号点到n号点的最短距离，如果无法从1号点走到n号点，则输出impossible。

数据保证不存在负权回路。

输入格式
第一行包含整数n和m。

接下来m行每行包含三个整数x，y，z，表示存在一条从点x到点y的有向边，边长为z。

输出格式
输出一个整数，表示1号点到n号点的最短距离。

如果路径不存在，则输出”impossible”。

数据范围
1≤n,m≤10⁵,
图中涉及边长绝对值均不超过10000。

输入样例：

3 3
1 2 5
2 3 -3
1 3 4

输出样例：

2

SPFA 算法是 Bellman-Ford算法的队列优化算法的别称，通常用于求含负权边的单源最短路径，以及判负权环。SPFA一般情况复杂度是O(m) 最坏情况下复杂度和朴素 Bellman-Ford 相同，为O(nm)。（n为顶点数，m为边数）

bellman-ford算法操作如下：

for n次
	for 所有边 u,v,w (松弛操作)
		dist[v] = min(dist[v],dist[u] + w)

spfa算法对第二层循环中所有边进行松弛操作进行了优化，原因是在bellman—ford算法中，即使该点的最短距离尚未更新过，也还是需要用尚未更新过的值去更新其他点，由此可知，该操作是不必要的，我们只需要找到更新过的值去更新其他点即可。

代码

#include <iostream>
#include<queue>
#include<cstring>
using namespace std;
int n,m;
/*h数组是邻接表的头它的下标是当前结点的编号，
值是当前结点第一条边的编号（其实是最后加入的那一条边），
e数组是边的集合，它的下标是当前边的编号，数值是当前边的终点。
ne是nextedge，如果ne是-1表示当前结点没有下一条边，
ne的下标是当前边的编号，数值是当前结点的下一条边的编号，
idx记录下一个要分配的边的编号，如果还是不明白，
先学习链式前向星后再看本文效果会更佳
*/
int h[100005],e[100005],ne[100005],w[100005],idx;
int dist[100005];//dist[i]表示1号点走到i号点的最短距离
bool visit[100005];//visit[i]表示第i个结点是否在队列中

//加入一条边
void add(int u,int v,int p){//注意：当前边的编号为idx
   e[idx]=v;//边的终点为v
   w[idx]=p;//边的权值为p
   ne[idx]=h[u];//当前边的下一条边为旧邻接表的头
   h[u]=idx++;//邻接表的头存储当前边的编号
}
int spfa(){
    //将1号点走到所有点的最短距离都初始化为0x3f3f3f3f
    //因为是按字节初始化的，即每个字节都初始化为0x3f
    memset(dist,0x3f,sizeof dist);
    dist[1]=0;//1号点到自己的距离为0
    queue<int> q;
    q.push(1);//将1号点压入队中
    visit[1]=1;//标记1号点已经在队列内
    while(q.size()){//队列不为空则循环
        int u=q.front();//取出队头元素
        q.pop();//弹出队头元素
        visit[u]=0;//标记u号点已出队
        //遍历以u为起点的每一条边
        for(int i=h[u];i!=-1;i=ne[i]){
            int v = e[i];//获得当前边的终点
            //如果借助<u,v>这条边能使到达点v的最短距离变小
            if(dist[v]>dist[u]+w[i]){
                dist[v]=dist[u]+w[i];//更新到达v点的最短距离
                //因为到达v点的最短距离已经变小，其相邻的点的
                //最短距离也可能因此而变小
                //所以将其压入队中去更新其相邻的点的最短距离
                if(!visit[v]){
                    q.push(v);
                    visit[v]=1;//标志v号点已经在队列内
                }
            }
        }
    }
    if(dist[n]==0x3f3f3f3f)//如果1号点走到n号点的最短距离还为初始值
        return -1;//即该点不可达，则返回-1
    return dist[n];//否则返回1号点走到n号点的最短距离
}
int main()
{
    cin>>n>>m;
    //邻接表的头全部初始化为-1，-1会随着插入的边而往后移
    //这也是上面为什么用-1作循环条件的原因for(int i=h[u];i!=-1;i=ne[i])
    memset(h,-1,sizeof h);
    while(m--){//输入m条边
        int a,b,c;
        cin>>a>>b>>c;
        add(a,b,c);
    }
    int r=spfa();
    if(r==-1){//输出结果
        cout<<"impossible"<<endl;
    }else{
        cout<<r<<endl;
    }
    return 0;
}