hdu1874 畅通工程续 单源最短路径（SSSP）Dijkstra算法

hdu1874 畅通工程续

有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;
}