有n个城镇,这n个城镇之间由m条道路连接,每条道路有一个长度,求出指定两个点间的最短距离。典型的单源最短路径(Single-Source Shortest Paths,SSSP),即从单个源点出发,到所有结点的最短路。算法适用于正权的有向图或无向图图,因为存在负权不一定有最短路,但在有最短路的情况下还是可以求出最短路,这要用到Bellman-Ford算法。
算法从起点开始进行扩展,直到所有的点都求出来了最短路。
描述如下:
设起点为s,终点为t,用d[i]表示从起点到结点i的路径长度,算法完成时,这个长度就是最短路径长度。如果要记录路径的话可以用一个pre[i]记录从起点到结点i的最短路径上的前一个结点最后就可以用回退打印出路径。
清除所有点的标号,将d[s]设为0,其他结点d值设为无穷大(即为INF,这个无穷大的选取要注意,不能溢出,也不能使一个有效的最短路径取到这个值)。
循环n次
在所有没有标号的结点中选出来d值最小的一个结点u,标记u结点。
对于从v出发的所有边(u,v),当d[u]+w[u][v]小于d[v]时更新d[v],并记录pre[v]为u。
循环结束
为了简单,用w[u][v]==INF表示边(u,v)不存在。值得一提的是,当结点u与结点v之间存在多条路径时只需记录最短的那个路径即可。
代码如下;
/*************************************************************************
> File Name: hdu1874.cpp
> Author: gwq
> Mail: 457781132@qq.com
> Created Time: 2014年10月14日 星期二 14时18分33秒
************************************************************************/
#include <cmath>
#include <ctime>
#include <cctype>
#include <climits>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <map>
#include <set>
#include <queue>
#include <stack>
#include <vector>
#include <iostream>
#include <algorithm>
#define INF (INT_MAX / 10)
#define SQR(x) ((x) * (x))
#define rep(i, n) for (int i = 0; i < (n); ++i)
#define repf(i, a, b) for (int i = (a); i <= (b); ++i)
#define repd(i, a, b) for (int i = (a); i >= (b); --i)
#define clr(arr, val) memset(arr, val, sizeof(arr))
#define pb push_back
#define sz(a) ((int)(a).size())
#define mid(x, y) ((x + y) / 2)
using namespace std;
typedef set<int> si;
typedef vector<int> vi;
typedef map<int, int> mii;
typedef long long ll;
#define N 210
#define M 1010
int n, m, w[N][N], v[N], d[N];
int dijkstra(int s, int t)
{
clr(v, 0);
for (int i = 0; i < n; ++i) {
d[i] = ((i == s) ? 0 : INF);
}
for (int i = 0; i < n; ++i) {
int x;
int m = INF;
for (int j = 0; j < n; ++j) {
if (v[j] == 0 && m >= d[j]) {
m = d[j];
x = j;
}
}
v[x] = 1;
for (int j = 0; j < n; ++j) {
d[j] = min(d[j], d[x] + w[x][j]);
}
}
if (d[t] == INF) {
return -1;
} else {
return d[t];
}
}
int main(int argc, char *argv[])
{
while (scanf("%d%d", &n, &m) != EOF) {
//INF表示不存在边
for (int i = 0; i < n; ++i) {
for (int j = 0; j < n; ++j) {
w[i][j] = INF;
}
}
for (int i = 0; i < m; ++i) {
int u, v, a;
scanf("%d%d%d", &u, &v, &a);
//处理多重边
if (a <= w[u][v]) {
w[u][v] = w[v][u] = a;
}
}
int s, t;
scanf("%d%d", &s, &t);
int ans = dijkstra(s, t);
printf("%d\n", ans);
}
return 0;
}