結果

問題 No.1812 Uribo Road
ユーザー MasKoaTS
提出日時 2021-12-30 11:18:40
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 216 ms / 5,000 ms
コード長 2,235 bytes
コンパイル時間 1,020 ms
コンパイル使用メモリ 94,112 KB
最終ジャッジ日時 2025-01-27 07:37:44
ジャッジサーバーID
(参考情報)
judge4 / judge1
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 4
other AC * 30
権限があれば一括ダウンロードができます

ソースコード

diff #

#include <iostream>
#include <vector>
#include <queue>
#include <algorithm>
#define rep(i, l, n) for (int i = (l); i < (n); i++)
using namespace std;
using Pair = pair<int, int>;
template <class T>
using V = vector<T>;
template <class T>
using VV = V<V<T> >;
const int inf = 2000000000;


VV<int> dijkstra(int n, VV<int>& edge, V<int>& R, VV<Pair>& route) {
	VV<int> dist(n, V<int>(0));
	V<int> st = { 0,n - 1 };
	rep(i, 0, R.size()) {
		st.push_back(edge[R[i]][0]);
		st.push_back(edge[R[i]][1]);
	}
	rep(i, 0, st.size()) {
		int k = st[i];
		V<int> di(n, -1);
		priority_queue<Pair, V<Pair>, greater<Pair> > pq;
		pq.push({ 0,k });
		int cnt = 0;
		while (pq.size() && cnt < n) {
			Pair p = pq.top();	pq.pop();
			int d = p.first;	int v = p.second;
			if (di[v] == -1) {
				cnt++;
				di[v] = d;
				rep(j, 0, route[v].size()) {
					Pair p = { d + route[v][j].second,route[v][j].first };
					pq.push(p);
				}
			}
		}
		dist[k] = di;
	}
	return dist;
}


int main(void) {
	int N, M, K;	cin >> N >> M >> K;

	V<int> R(K);
	rep(i, 0, K) {
		cin >> R[i];	R[i] -= 1;
	}

	VV<int> edge(M, V<int>(3));
	VV<Pair> route(N, V<Pair>(0));
	rep(i, 0, M) {
		rep(j, 0, 3) {
			cin >> edge[i][j];
		}
		edge[i][0] -= 1;	edge[i][1] -= 1;
		route[edge[i][0]].push_back({ edge[i][1],edge[i][2] });
		route[edge[i][1]].push_back({ edge[i][0],edge[i][2] });
	}

	int ans = inf;
	int rsum = 0;
	rep(i, 0, R.size()) {
		rsum += edge[R[i]][2];
	}
	VV<int> dist = dijkstra(N, edge, R, route);


	V<VV<int> > dp(1 << K, VV<int>(K, V<int>(2, inf)));
	rep(i, 0, K) {
		rep(j, 0, 2) {
			dp[1 << i][i][j] = rsum + dist[0][edge[R[i]][j]];
		}	
	}
	rep(bit, 1, 1 << K) {
		rep(s, 0, K) {
			if ((bit & (1 << s)) == 0) {
				continue;
			}
			rep(i, 0, 2) {
				rep(t, 0, K) {
					if (bit & (1 << t)) {
						continue;
					}
					int b = bit | (1 << t);
					rep(j, 0, 2) {
						int p = edge[R[t]][j];
						int d = dp[bit][s][i ^ 1] + dist[edge[R[s]][i]][edge[R[t]][j]];
						if (dp[b][t][j] > d) {
							dp[b][t][j] = d;
						}
					}		
				}
			}
		}
	}

	ans = inf;
	rep(i, 0, K) {
		rep(j, 0, 2) {
			int d = dp[(1 << K) - 1][i][j ^ 1] + dist[edge[R[i]][j]][N - 1];
			if (ans > d) {
				ans = d;
			}	
		}	
	}

	cout << ans << endl;
	return 0;
}
0