結果

問題 No.114 遠い未来
ユーザー antaanta
提出日時 2014-12-29 00:51:59
言語 C++11
(gcc 11.4.0)
結果
TLE  
実行時間 -
コード長 3,827 bytes
コンパイル時間 2,494 ms
コンパイル使用メモリ 94,532 KB
実行使用メモリ 5,428 KB
最終ジャッジ日時 2023-09-03 19:09:54
合計ジャッジ時間 36,995 ms
ジャッジサーバーID
（参考情報） 		judge15 / judge11
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 25 ms
4,380 KB
testcase_01 TLE -
testcase_02 AC 2,336 ms
4,384 KB
testcase_03 AC 287 ms
4,384 KB
testcase_04 AC 2 ms
4,384 KB
testcase_05 AC 4 ms
4,380 KB
testcase_06 AC 2,775 ms
4,380 KB
testcase_07 AC 2 ms
4,384 KB
testcase_08 AC 2 ms
4,380 KB
testcase_09 AC 7 ms
4,384 KB
testcase_10 AC 117 ms
4,384 KB
testcase_11 AC 346 ms
5,428 KB
testcase_12 TLE -
testcase_13 AC 4,962 ms
4,384 KB
testcase_14 AC 2,817 ms
4,384 KB
testcase_15 AC 2,506 ms
4,384 KB
testcase_16 AC 1,450 ms
4,384 KB
testcase_17 AC 1,313 ms
4,380 KB
testcase_18 AC 1,314 ms
4,384 KB
testcase_19 AC 703 ms
4,384 KB
testcase_20 AC 163 ms
4,380 KB
testcase_21 AC 47 ms
4,384 KB
testcase_22 AC 43 ms
4,380 KB
testcase_23 AC 3 ms
4,380 KB
testcase_24 AC 3 ms
4,384 KB
testcase_25 AC 2 ms
4,380 KB
testcase_26 AC 2 ms
4,384 KB
testcase_27 AC 1 ms
4,384 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#define _CRT_SECURE_NO_WARNINGS
#include <string>
#include <vector>
#include <algorithm>
#include <numeric>
#include <set>
#include <map>
#include <queue>
#include <iostream>
#include <sstream>
#include <cstdio>
#include <cmath>
#include <ctime>
#include <cstring>
#include <cctype>
#include <cassert>
#include <limits>
#include <functional>
#define rep(i,n) for(int (i)=0;(i)<(int)(n);++(i))
#define rer(i,l,u) for(int (i)=(int)(l);(i)<=(int)(u);++(i))
#define reu(i,l,u) for(int (i)=(int)(l);(i)<(int)(u);++(i))
#if defined(_MSC_VER) || __cplusplus > 199711L
#define aut(r,v) auto r = (v)
#else
#define aut(r,v) typeof(v) r = (v)
#endif
#define each(it,o) for(aut(it, (o).begin()); it != (o).end(); ++ it)
#define all(o) (o).begin(), (o).end()
#define pb(x) push_back(x)
#define mp(x,y) make_pair((x),(y))
#define mset(m,v) memset(m,v,sizeof(m))
#define INF 0x3f3f3f3f
#define INFL 0x3f3f3f3f3f3f3f3fLL
using namespace std;
typedef vector<int> vi; typedef pair<int,int> pii; typedef vector<pair<int,int> > vpii;
typedef long long ll; typedef vector<long long> vl; typedef pair<long long,long long> pll; typedef vector<pair<long long,long long> > vpll;
typedef vector<string> vs; typedef long double ld;
template<typename T, typename U> inline void amin(T &x, U y) { if(y < x) x = y; }
template<typename T, typename U> inline void amax(T &x, U y) { if(x < y) x = y; }

struct UnionFind {
	static const int MAXN = 35;
	typedef char Index;	//適切に書き換えること。符号付きでないといけないので注意
	Index data[MAXN];
	void init(int size_) { memset(data, -1, size_ * sizeof(*data)); }
	bool unionSet(int x, int y) {
		x = root(x); y = root(y);
		if (x != y) {
			if (data[y] < data[x]) swap(x, y);
			data[x] += data[y]; data[y] = x;
		}
		return x != y;
	}
	inline bool findSet(int x, int y) { return root(x) == root(y); }
	inline int root(int x) { return data[x] < 0 ? x : data[x] = root(data[x]); }
	inline int size(int x) { return -data[root(x)]; }
};

int main() {
	int N, M, T;
	scanf("%d%d%d", &N, &M, &T);
	vector<vector<int> > g(N, vector<int>(N, INF));
	rep(i, N)
		g[i][i] = 0;
	rep(i, M) {
		int a, b, c;
		scanf("%d%d%d", &a, &b, &c), -- a, -- b;
		amin(g[a][b], c);
		amin(g[b][a], c);
	}
	rep(k, N) rep(i, N) rep(j, N)
		amin(g[i][j], g[i][k] + g[k][j]);
	vi required(T);
	vector<bool> important(N, false);
	rep(i, T) {
		int v;
		scanf("%d", &v), -- v;
		required[i] = v;
		important[v] = true;
	}
	if(T == 1) {
		puts("0");
		return 0;
	}
	vi steiners;
	rep(i, N) if(!important[i])
		steiners.push_back(i);

	int X = steiners.size();
	int ans = INF;

	if(T >= 15) {
		vector<pair<int,pii> > edges;
		rep(i, N) rep(j, i)
			edges.push_back(mp(g[i][j], mp(i, j)));
		sort(all(edges));
		vector<bool> use(N, false);
		UnionFind uf;
		rep(usesteiners, 1 << X) {
			rep(i, N)
				use[i] = important[i];
			int components = T;
			rep(i, X) if(usesteiners >> i & 1) {
				use[steiners[i]] = true;
				++ components;
			}
			int totalcost = 0;
			uf.init(N);
			rep(i, edges.size()) {
				int x = edges[i].second.first, y = edges[i].second.second;
				if(use[x] && use[y] && uf.unionSet(x, y)) {
					totalcost += edges[i].first;
					-- components;
				}
			}
			if(components == 1)
				amin(ans, totalcost);
		}
	}else {
		vector<int> dp((1 << T) * N, INF);

		rep(p, T) rep(q, N)
			dp[(1 << p) * N + q] = g[required[p]][q];

		rep(S, 1 << T) {
			if(!(S & (S-1))) continue;
			rep(p, N) {
				int x = INF;
				for(int E = (S-1) & S; E > 0; (-- E) &= S)
					amin(x, dp[E * N + p] + dp[(S - E) * N + p]);
				dp[S * N + p] = x;
			}
			rep(p, N) {
				int x = INF;
				rep(q, N)
					amin(x, dp[S * N + q] + g[p][q]);
				dp[S * N + p] = x;
			}
		}

		int U = (1 << T) - 1;
		rep(S, 1 << T) rep(q, N)
			amin(ans, dp[S * N + q] + dp[(U - S) * N + q]);
	}
	printf("%d\n", ans);
	return 0;
}
0