SRM 573 div1 lev2（状态图上的最短路）

题目链接

从一个神牛博客上无意中发现的题目， 由于TC上的图论题目本来就少， 所以就好奇的想做做， 题意就是给定一个50

个点的图， 每个点有一个高度， 从一个点u到另一个点v的条件是u, v直接相连， 并且h[u] >= h[v]，每个点的高度可以

改变代价为改变前后高度差的绝对值， 求从点0到点n - 1的最小代价。

一开始没什么思路， 后来看了神牛的解法后发现这题其实就是lrj白书上提到的状态图的最短路问题， 定义dp[u][h]表示

顶点u高度为h时最小代价， 由于最优方案中每个点的高度可以是这50个点的起始高度所以可以将高度离散化一下， 由

于转移代价的特殊性， 我们可以看成是一个50 * 50个点的最短路问题， 值得注意的是此题的状态转移不具有拓扑性， 

所以必须用最短路算法求解。。。

#include <vector>
#include <list>
#include <map>
#include <queue>
#include <set>
#include <deque>
#include <stack>
#include <bitset>
#include <algorithm>
#include <functional>
#include <numeric>
#include <utility>
#include <sstream>
#include <iostream>
#include <iomanip>
#include <cstdio>
#include <cmath>
#include <cstdlib>
#include <ctime>

using namespace std;

class SkiResorts {
public:
	long long minCost(vector <string>, vector <int>);
};

#define fi first
#define se second
#define mp make_pair

typedef long long LL;
const int N = 55;
const LL INF = 1LL << 60;


queue<pair<int, int> > Q;
LL dp[N][N];
int tmp[N];
bool inq[N][N];

void init() {
	for (int i = 0; i < N; i++) {
		fill(dp[i], dp[i] + N, INF);
		fill(inq[i], inq[i] + N, 0);
	}
}


long long SkiResorts::minCost(vector<string> adj, vector <int> h) {
	int n = h.size();
	for (int i = 0; i < n; i++)
		tmp[i] = h[i];

	sort(tmp, tmp + n);
	init();
	
	for (int i = 0; i < n; i++) {
		dp[0][i] = abs(tmp[i] - h[0]);
		Q.push(mp(0, i));
		inq[0][i] = 1;
	}

	while (!Q.empty()) {
		pair<int, int> cur = Q.front();
		Q.pop();
		inq[cur.fi][cur.se] = 0;

		int u = cur.fi; 
		for (int v = 0; v < n; v++)
			if (adj[u][v] != 'N') {
				for (int j = cur.se; j >= 0; j--) {
					if (dp[v][j] > dp[u][cur.se] + abs(tmp[j] - h[v])) {
						dp[v][j] = dp[u][cur.se] + abs(tmp[j] - h[v]);
						if (!inq[v][j]) {
							Q.push(mp(v, j));
						}	
					}
				}	
			}
	}

	LL res = INF;

	for (int i = 0; i < n; i++)
		res = min(res, dp[n - 1][i]);

	if (res == INF) 
		return -1;
	else
		return res;
}

<%:testing-code%>
//Powered by [KawigiEdit] 2.0!