結果
問題 | No.95 Alice and Graph |
ユーザー | Demystify |
提出日時 | 2022-05-07 18:10:14 |
言語 | C++17 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 662 ms / 5,000 ms |
コード長 | 6,287 bytes |
コンパイル時間 | 2,415 ms |
コンパイル使用メモリ | 215,272 KB |
実行使用メモリ | 7,616 KB |
最終ジャッジ日時 | 2024-07-06 23:29:05 |
合計ジャッジ時間 | 4,875 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge5 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 253 ms
7,580 KB |
testcase_01 | AC | 35 ms
7,616 KB |
testcase_02 | AC | 29 ms
7,612 KB |
testcase_03 | AC | 2 ms
6,944 KB |
testcase_04 | AC | 366 ms
7,528 KB |
testcase_05 | AC | 289 ms
7,576 KB |
testcase_06 | AC | 356 ms
7,552 KB |
testcase_07 | AC | 35 ms
7,584 KB |
testcase_08 | AC | 2 ms
6,940 KB |
testcase_09 | AC | 2 ms
6,944 KB |
testcase_10 | AC | 1 ms
6,940 KB |
testcase_11 | AC | 2 ms
6,944 KB |
testcase_12 | AC | 2 ms
6,940 KB |
testcase_13 | AC | 662 ms
7,488 KB |
ソースコード
#include <bits/stdc++.h> using namespace std; using ll = long long; using ld = long double; // -------------------------------------------------------- #define FOR(i,l,r) for (int i = (l); i < (r); ++i) #define RFOR(i,l,r) for (int i = (r)-1; (l) <= i; --i) #define REP(i,n) FOR(i,0,n) #define RREP(i,n) RFOR(i,0,n) #define ALL(c) (c).begin(), (c).end() #define RALL(c) (c).rbegin(), (c).rend() #define SORT(c) sort(ALL(c)) #define RSORT(c) sort(RALL(c)) #define MIN(c) *min_element(ALL(c)) #define MAX(c) *max_element(ALL(c)) #define SUM(c) accumulate(ALL(c), 0LL) #define COUNT(c,v) count(ALL(c),(v)) #define SZ(c) ((ll)(c).size()) #define BIT(b,i) (((b)>>(i)) & 1) #define PCNT(b) __builtin_popcountll(b) #define P0(i) (((i) & 1) == 0) #define P1(i) (((i) & 1) == 1) #ifdef _LOCAL #define debug_bar cerr << "--------------------\n"; #define debug(x) cerr << "l." << __LINE__ << " : " << #x << " = " << (x) << '\n' #define debug_pair(x) cerr << "l." << __LINE__ << " : " << #x << " = (" << x.first << "," << x.second << ")\n"; template<class T> void debug_line(const vector<T>& ans, int l, int r, int L = 0) { cerr << "l." << L << " :"; for (int i = l; i < r; i++) { cerr << ' ' << ans[i]; } cerr << '\n'; } #else #define cerr if (false) cerr #define debug_bar #define debug(x) #define debug_pair(x) template<class T> void debug_line([[maybe_unused]] const vector<T>& ans, [[maybe_unused]] int l, [[maybe_unused]] int r, [[maybe_unused]] int L = 0) {} #endif template<class... T> void input(T&... a) { (cin >> ... >> a); } void print() { cout << '\n'; } template<class T> void print(const T& a) { cout << a << '\n'; } template<class T, class... Ts> void print(const T& a, const Ts&... b) { cout << a; (cout << ... << (cout << ' ', b)); cout << '\n'; } template<class T> void cout_line(const vector<T>& ans, int l, int r) { for (int i = l; i < r; i++) { if (i != l) { cout << ' '; } cout << ans[i]; } cout << '\n'; } template<class T> bool chmin(T& a, const T b) { if (b < a) { a = b; return 1; } return 0; } template<class T> bool chmax(T& a, const T b) { if (a < b) { a = b; return 1; } return 0; } template<class T> pair<T,T> divmod(T a, T b) { assert(a >= 0 && b > 0); return make_pair(T(a / b), T(a % b)); } template<class T> T mod(T x, T m) { assert(m != 0); return T((x % m + m) % m); } template<class T> T ceil(T a, T b) { if (b < 0) { return ceil(T(-a), T(-b)); } assert(b > 0); return (a < 0 ? T(a / b) : T((a + b - 1) / b)); } template<class T> T floor(T a, T b) { if (b < 0) { return floor(T(-a), T(-b)); } assert(b > 0); return (a > 0 ? T(a / b) : T((a - b + 1) / b)); } template<class T> T powint(T x, T n) { assert(n >= 0); if (n == 0) { return T(1); }; T res = powint(x, T(n>>1)); res *= res; if (n & 1) { res *= x; } return res; } ll bitlen(ll b) { if (b <= 0) { return 0; } return (64LL - __builtin_clzll(b)); } ll digit_len(ll n) { assert(n >= 0); if (n == 0) { return 1; } ll sum = 0; while (n > 0) { sum++; n /= 10; } return sum; } ll digit_sum(ll n) { assert(n >= 0); ll sum = 0; while (n > 0) { sum += n % 10; n /= 10; } return sum; } ll digit_prod(ll n) { assert(n >= 0); if (n == 0) { return 0; } ll prod = 1; while (n > 0) { prod *= n % 10; n /= 10; } return prod; } ll xor_sum(ll x) { assert(0 <= x); switch (x % 4) { case 0: return x; case 1: return 1; case 2: return x ^ 1; case 3: return 0; } assert(false); } string toupper(const string& S) { string T(S); for (int i = 0; i < (int)T.size(); i++) { T[i] = toupper(T[i]); } return T; } string tolower(const string& S) { string T(S); for (int i = 0; i < (int)T.size(); i++) { T[i] = tolower(T[i]); } return T; } int a2i(const char& c) { assert(islower(c)); return (c - 'a'); } int A2i(const char& c) { assert(isupper(c)); return (c - 'A'); } int d2i(const char& d) { assert(isdigit(d)); return (d - '0'); } char i2a(const int& i) { assert(0 <= i && i < 26); return ('a' + i); } char i2A(const int& i) { assert(0 <= i && i < 26); return ('A' + i); } char i2d(const int& i) { assert(0 <= i && i <= 9); return ('0' + i); } using P = pair<ll,ll>; using VP = vector<P>; using VVP = vector<VP>; using VS = vector<string>; using VVS = vector<VS>; using VI = vector<int>; using VVI = vector<VI>; using VVVI = vector<VVI>; using VLL = vector<ll>; using VVLL = vector<VLL>; using VVVLL = vector<VVLL>; using VB = vector<bool>; using VVB = vector<VB>; using VVVB = vector<VVB>; using VD = vector<double>; using VVD = vector<VD>; using VVVD = vector<VVD>; using VLD = vector<ld>; using VVLD = vector<VLD>; using VVVLD = vector<VVLD>; const ld EPS = 1e-10; const ld PI = acosl(-1.0); constexpr ll MOD = 1000000007; // constexpr ll MOD = 998244353; constexpr int inf = (1 << 30) - 1; // 1073741824 - 1 constexpr ll INF = (1LL << 62) - 1; // 4611686018427387904 - 1 // -------------------------------------------------------- // #include <atcoder/all> // using namespace atcoder; // ワーシャルフロイド : O(N^3) void warshall_floyd(VVLL& G) { int N = G.size(); for (int k = 0; k < N; k++) { for (int i = 0; i < N; i++) if (G[i][k] != INF) { for (int j = 0; j < N; j++) if (G[k][j] != INF) { chmin(G[i][j], G[i][k] + G[k][j]); } } } } int main() { ios::sync_with_stdio(false); cin.tie(nullptr); cout << fixed << setprecision(15); ll N, M, K; cin >> N >> M >> K; VVLL G(N); REP(_,M) { ll u, v; cin >> u >> v; u--; v--; G[u].push_back(v); G[v].push_back(u); } VLL C(N); REP(i,N) C[i] = powint<ll>(2, i) - 1; VVLL dist(N,VLL(N,INF)); REP(u,N) { for (auto v : G[u]) dist[u][v] = 1; } warshall_floyd(dist); VI X; RFOR(s,1,N) { VI Y(X); Y.push_back(s); int n = SZ(Y); VVLL dp(n, VLL(1<<n, INF)); REP(u,n) dp[u][1<<u] = dist[0][Y[u]]; REP(S,1<<n) { REP(u,n) if (BIT(S,u)) { REP(v,n) if (not (BIT(S,v))) { chmin(dp[v][S^(1<<v)], dp[u][S] + dist[Y[u]][Y[v]]); } } } ll mn = INF; REP(u,n) chmin(mn, dp[u][(1<<n)-1]); if (mn <= K) X = Y; if (SZ(X) == K) break; } ll ans = 0; for (const auto& x : X) ans += C[x]; print(ans); return 0; }