从一个神牛博客上无意中发现的题目, 由于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!