hdu2544 最短路 spfa

最短路

Time Limit: 5000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 77752    Accepted Submission(s): 33784

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条商店到赛场的路线。

 

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

 

单源最短路模板题

Bellman-Ford 算法

这个算法主要是构造一个最短路径长度数组的序列：dist[1][u], dist[2][u],···,dist[n−1][u]。

 其中 dist[k][u] 表示从源 s 到 u 至多经过 k 条边的最短路径的长度。 

显然我们可以得到这样的关系：

dist[1][u] = w[s][u]

dist[k][u] = min(dist[k−1][u], min (v,u)∈E (dist[k−1][v] +w[v][u]))

 由于每次计算 dist 数组复杂度都是O(|E|) 的，总的复杂度就是 O(|V||E|)。

事实上，你会发现我们每次进行的计算可以看成是一次“松弛”。假设我 们现在已经得到了 Bellman-Ford 算法某个阶段的 dist 数组，然后我们发现 了一条 s 到 u 的距离比 dist[u] 更加短的路径。我们更新了 dist[u]。 

spfa算法

那么接下来直接受到影响的就是与 u 直接关联的顶点 v，也就是如果 dist[u] +w[u][v] < dist[v] 的话，s 到 v 的最短路就可以利用 s 到 u 的最 短路加上 u 到 v 的边来更新。这样的话与 v 直接关联的顶点又会受到影响 ……不断这样持续下去直到最后没有顶点能被影响。 

那么一个优化就是我们利用队列存储这些需要更新的结点，每次从队列 中取出一个结点，计算是否有结点需要更新，如果有，并且这个结点不在队 列中，那么就将它加入队列。 这样的算法被称为 SPFA——一种优化的 Bellman-Ford 算法。

在竞赛中大多数人会选择用 SPFA 作为单源最段路的算法，主要原因在 于它比较好写，而且通常情况下跑得比较快。但是 SPFA 的复杂度实际上是 没有保证的，最坏情况基本和 Bellman-Ford 相同，但是最好的时候可以到 达O(|V|+|E|)。

非常重要的一点就是 SPFA 的队列需要使用环队列，虽然最多队列里 只会有 n 各点，但是每个点可能会入队多次。

#include <bits/stdc++.h>
using namespace std;
const int N = 1005;
const int INF = INT_MAX;
int n,m;
int dist[N];
struct Node{
    int to,w;
};
vector <Node> G[105];
queue <int> q;
void spfa(){
    while(!q.empty()){
        int u = q.front();
        q.pop();
        for (int i = 0;i < G[u].size();i++){
            int v = G[u][i].to,w = G[u][i].w;
            if (dist[u] + w < dist[v]){
                dist[v] = dist[u] + w;
                q.push(v);
            }
        }
    }
}

int main(){
    while(~scanf("%d%d",&n,&m),n+m){
    	for (int i = 0;i <= n;i++) G[i].clear();
        for (int i = 0;i <= n;i++) dist[i] = INF;
        int point = INF;
        for (int i = 0;i < m;i++){
            int a,b,c;
            scanf("%d%d%d",&a,&b,&c);
            G[a].push_back({b,c});
            G[b].push_back({a,c});
        }
        dist[1] = 0;
        q.push(1);
        spfa();
        printf("%d\n",dist[n]);
    }
    return 0;
}