[JLOI2011]飛行路線 題解

題目TP門

題目描述

Alice和Bob現在要乘飛機旅行，他們選擇了一家相對便宜的航空公司。該航空公司一共在n個城市設有業務，設這些城市分別標記為0到n-1，一共有m種航線，每種航線連線兩個城市，並且航線有一定的價格。Alice和Bob現在要從一個城市沿著航線到達另一個城市，途中可以進行轉機。航空公司對他們這次旅行也推出優惠，他們可以免費在最多k種航線上搭乘飛機。那麼Alice和Bob這次出行最少花費多少？

輸入格式

資料的第一行有三個整數，n,m,k，分別表示城市數，航線數和免費乘坐次數。
第二行有兩個整數，s,t，分別表示他們出行的起點城市編號和終點城市編號。(0<=s,t<n)
接下來有m行，每行三個整數，a,b,c，表示存在一種航線，能從城市a到達城市b，或從城市b到達城市a，價格為c。(0<=a,b<n,a與b不相等，0<=c<=1000)

輸出格式

只有一行，包含一個整數，為最少花費。

樣例

樣例輸入

5 6 1
0 4
0 1 5
1 2 5
2 3 5
3 4 5
2 3 3
0 2 100

樣例輸出

8

資料範圍與提示
對於30%的資料,2<=n<=50,1<=m<=300,k=0;

對於50%的資料,2<=n<=600,1<=m<=6000,0<=k<=1;

對於100%的資料,2<=n<=10000,1<=m<=50000,0<=k<=10.

解析

10分做法：
爆搜刪邊，求出其最小解法。

#include <map>
#include <queue>
#include <cstdio>
#include <vector>
#include <cstring>
#include <iostream>
using namespace std;
const int MAXN = 1e5 + 5;
const int MAXM = 5e5 + 5;
struct Node {
	int v, m;
	Node() {}
	Node(int V, int M) {
		v = V;
		m = M;
	}
	friend bool operator < (Node x, Node y) {
		return x.m > y.m;
	}
};
priority_queue<Node> q;
map<pair<int, int>, int> mp;
map<pair<int, int>, bool> vis;
map<pair<int, int>, bool> vis1;
vector<Node> a[MAXN];
int d[MAXN], x[MAXM], y[MAXM];
bool f[MAXN];
int n, m, S, T, s, minn = 0x7f7f7f7f;
void Read();
void dfs(int, int);
int Dijkstra(int, int);
int main() {
	Read();
	dfs(1, 0);
	cout << minn;
	return 0;
}
void dfs(int now, int last) {
	if(last == s)
		minn = min(minn, Dijkstra(S, T));
	if(now == m + 1)
		return;
	a[x[now]].push_back(Node(y[now], 0));
	a[y[now]].push_back(Node(x[now], 0));
	dfs(now + 1, last + 1);
	a[x[now]].pop_back();
	a[y[now]].pop_back();
	dfs(now + 1, last);
}
void Read() {
	scanf("%d %d %d %d %d", &n, &m, &s, &S, &T);
	for(int i = 1; i <= m; i++) {
		int A, B, C;
		scanf("%d %d %d", &A, &B, &C);
		x[i] = A; y[i] = B;
		if(!vis[make_pair(A, B)]) {
			mp[make_pair(A, B)] = C;
			mp[make_pair(B, A)] = C;
			vis[make_pair(A, B)] = 1;
			vis[make_pair(B, A)] = 1;
		}
		else {
			mp[make_pair(A, B)] = min(mp[make_pair(A, B)], C);
			mp[make_pair(B, A)] = mp[make_pair(A, B)];
		}
	}
	for(int i = 1; i <= m; i++) {
		if(!vis1[make_pair(x[i], y[i])]) {
			a[x[i]].push_back(Node(y[i], mp[make_pair(x[i], y[i])]));
			a[y[i]].push_back(Node(x[i], mp[make_pair(x[i], y[i])]));
			vis1[make_pair(x[i], y[i])] = 1;
			vis1[make_pair(y[i], x[i])] = 1;
		}
	}
}
int Dijkstra(int c, int t) {
	q.push(Node(c, 0));
	memset(d, 0x3f, sizeof(d));
	memset(f, 0, sizeof(f));
	d[c] = 0;
	while(!q.empty()) {
		Node now = q.top();
		q.pop();
		int i = now.v;
		if(f[i]) continue;
		f[i] = 1;
		int SIZ = a[i].size();
		for(int j = 0; j < SIZ; j++) {
			if(d[a[i][j].v] > d[i] + a[i][j].m) {
				d[a[i][j].v] = d[i] + a[i][j].m;
				q.push(Node(a[i][j].v, d[a[i][j].v]));
			}
		}
	}
	return d[t];
}
 

當然，這是考試時實在沒想到正解而做出的無奈之舉。
可以發現k<=10，範圍極小。
所以本題就可以用到分層圖：

由上圖可知：複製k+1份圖連線相鄰兩層的對應點，由題意得：每層圖由長度為0的邊相連，就可以等價於有k次機會能夠免費。再跑一邊最短路Dijkstra即可。
所以說這道題就是一道板題。

正解C++程式碼:

#include <queue>
#include <cstdio>
#include <vector>
#include <cstring>
#include <iostream>
using namespace std;
const int MAXN = 1e6 + 5;
const int MAXM = 1e6 + 5;
struct Node {
	int v, m;
	Node() {}
	Node(int V, int M) {
		v = V;
		m = M;
	}
	friend bool operator < (Node x, Node y) {
		return x.m > y.m;
	}
};
priority_queue<Node> q;
vector<Node> a[MAXN];
int d[MAXM], x[MAXM], y[MAXM];
bool f[MAXN];
int n, m, S, T, k, minn = 0x7f7f7f7f;
void Read();
int Dijkstra(int, int);
int main() {
	Read(); 
	cout << Dijkstra(S, T + k * n);
	return 0;
}
void Read() {
	scanf("%d %d %d %d %d", &n, &m, &k, &S, &T);
	for(int i = 1; i <= m; i++) {
		int u, v, w;
		scanf("%d %d %d", &u, &v, &w);
		a[u].push_back(Node(v, w));
		a[v].push_back(Node(u, w));
		for(int j = 1; j <= k; j++) {
			a[u + (j - 1) * n].push_back(Node(v + j * n, 0));
			a[v + (j - 1) * n].push_back(Node(u + j * n, 0));
			a[u + j * n].push_back(Node(v + j * n, w));
			a[v + j * n].push_back(Node(u + j * n, w));
        }
	}
	for(int i = 1; i <= k; i++)
		a[T + (i - 1) * n].push_back(Node(T + i * n, 0));
}
int Dijkstra(int c, int t) {
	q.push(Node(c, 0));
	memset(d, 0x7f, sizeof(d));
	memset(f, 0, sizeof(f));
	d[c] = 0;
	while(!q.empty()) {
		Node now = q.top();
		q.pop();
		int i = now.v;
		if(f[i]) continue;
		f[i] = 1;
		int SIZ = a[i].size();
		for(int j = 0; j < SIZ; j++) {
			if(d[a[i][j].v] > d[i] + a[i][j].m) {
				d[a[i][j].v] = d[i] + a[i][j].m;
				q.push(Node(a[i][j].v, d[a[i][j].v]));
			}
		}
	}
	return d[t];
}