結果
問題 | No.1124 Earthquake Safety |
ユーザー |
![]() |
提出日時 | 2020-07-22 23:26:34 |
言語 | C++17(clang) (17.0.6 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 446 ms / 3,000 ms |
コード長 | 3,125 bytes |
コンパイル時間 | 7,493 ms |
コンパイル使用メモリ | 163,456 KB |
実行使用メモリ | 50,176 KB |
最終ジャッジ日時 | 2024-06-23 01:19:19 |
合計ジャッジ時間 | 20,025 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge2 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 55 ms
15,104 KB |
testcase_01 | AC | 53 ms
15,104 KB |
testcase_02 | AC | 57 ms
15,104 KB |
testcase_03 | AC | 56 ms
15,104 KB |
testcase_04 | AC | 54 ms
15,104 KB |
testcase_05 | AC | 55 ms
15,232 KB |
testcase_06 | AC | 56 ms
15,104 KB |
testcase_07 | AC | 60 ms
15,744 KB |
testcase_08 | AC | 112 ms
21,120 KB |
testcase_09 | AC | 293 ms
33,280 KB |
testcase_10 | AC | 194 ms
50,176 KB |
testcase_11 | AC | 52 ms
15,048 KB |
testcase_12 | AC | 52 ms
15,104 KB |
testcase_13 | AC | 53 ms
15,104 KB |
testcase_14 | AC | 324 ms
33,280 KB |
testcase_15 | AC | 354 ms
33,280 KB |
testcase_16 | AC | 320 ms
33,280 KB |
testcase_17 | AC | 345 ms
33,280 KB |
testcase_18 | AC | 329 ms
33,408 KB |
testcase_19 | AC | 375 ms
33,280 KB |
testcase_20 | AC | 334 ms
33,280 KB |
testcase_21 | AC | 333 ms
33,280 KB |
testcase_22 | AC | 375 ms
33,408 KB |
testcase_23 | AC | 344 ms
32,896 KB |
testcase_24 | AC | 319 ms
32,896 KB |
testcase_25 | AC | 308 ms
33,024 KB |
testcase_26 | AC | 323 ms
33,024 KB |
testcase_27 | AC | 332 ms
34,048 KB |
testcase_28 | AC | 384 ms
33,920 KB |
testcase_29 | AC | 446 ms
37,248 KB |
testcase_30 | AC | 341 ms
37,248 KB |
testcase_31 | AC | 404 ms
44,160 KB |
testcase_32 | AC | 340 ms
42,880 KB |
testcase_33 | AC | 358 ms
41,216 KB |
testcase_34 | AC | 400 ms
42,752 KB |
testcase_35 | AC | 293 ms
33,920 KB |
testcase_36 | AC | 285 ms
33,152 KB |
testcase_37 | AC | 288 ms
34,944 KB |
testcase_38 | AC | 279 ms
34,304 KB |
testcase_39 | AC | 293 ms
33,408 KB |
testcase_40 | AC | 267 ms
34,816 KB |
testcase_41 | AC | 245 ms
33,664 KB |
testcase_42 | AC | 248 ms
34,432 KB |
testcase_43 | AC | 280 ms
34,944 KB |
testcase_44 | AC | 240 ms
34,816 KB |
testcase_45 | AC | 246 ms
34,048 KB |
testcase_46 | AC | 251 ms
34,560 KB |
testcase_47 | AC | 244 ms
34,048 KB |
testcase_48 | AC | 236 ms
34,052 KB |
testcase_49 | AC | 263 ms
39,728 KB |
testcase_50 | AC | 55 ms
14,976 KB |
testcase_51 | AC | 55 ms
15,232 KB |
testcase_52 | AC | 55 ms
15,104 KB |
testcase_53 | AC | 55 ms
15,104 KB |
testcase_54 | AC | 54 ms
15,232 KB |
testcase_55 | AC | 54 ms
15,104 KB |
testcase_56 | AC | 54 ms
15,104 KB |
testcase_57 | AC | 54 ms
15,104 KB |
ソースコード
#include<bits/stdc++.h> using namespace std; #define int long long #define ii pair <int, int> #define app push_back #define all(a) a.begin(), a.end() #define bp __builtin_popcountll #define ll long long #define mp make_pair #define f first #define s second #define Time (double)clock()/CLOCKS_PER_SEC #define debug(x) std::cout << #x << ": " << x << '\n'; const int MOD = 1000 * 1000 * 1000 + 7; const int N = 3e5+7; int mod(int n) { n %= MOD; if (n < 0) return n + MOD; else return n; } int fp(int a, int p) { int ans = 1, c = a; for (int i = 0; (1ll << i) <= p; ++i) { if ((p >> i) & 1) ans = mod(ans * c); c = mod(c * c); } return ans; } int dv(int a, int b) { return mod(a * fp(b, MOD - 2)); } int pw[N], inv[N], inv2; vector <int> g[N]; int up[N]; int sum[N]; int par[N]; void add(int &a, int b) { a = mod(a + b); } void dfs1(int u, int p) { par[u] = p; sum[u] = 1; for (int v : g[u]) { if (v != p) { dfs1(v, u); add(sum[u], sum[v] * inv2); } } } void dfs2(int u, int p) { int s = 0; for (int v : g[u]) { if (v != p) { add(s, sum[v]); } } for (int v : g[u]) { if (v != p) { int oth = mod(s - sum[v]); up[v] = mod(mod(up[u] + oth) * inv2 + 1); dfs2(v, u); } } } int n; int solve(int root) { vector <int> d; for (int v : g[root]) { if (v != par[root]) { d.app(sum[v]); } } if (root != 1) { d.app(up[root]); } for (int &e : d) e = mod(e * inv2); /* debug(root); for (auto e : d) cout << e << ' '; cout << endl; */ int t = 0; for (int e : d) { t += e; t %= MOD; } int ans = ((t * t) % MOD) * t % MOD; ans = mod(ans + t * t * 3); for (int e : d) { ans -= 3 * mod(mod(e * e) * t); ans += 2 * mod(mod(e * e) * e); ans = mod(ans); } for (int e : d) { ans -= mod(e * e * 3); ans = mod(ans); } for (int e : d) { ans = mod(ans + e * 3); } ans = mod(ans + 1); return ans; } signed main() { #ifdef HOME freopen("input.txt", "r", stdin); #else #define endl '\n' ios_base::sync_with_stdio(0); cin.tie(0); #endif pw[0] = 1; for (int i = 1; i < N; ++i) pw[i] = pw[i - 1] * 2 % MOD; for (int i = 0; i < N; ++i) inv[i] = fp(pw[i], MOD - 2); inv2 = fp(2, MOD - 2); cin >> n; for (int i = 0; i < n - 1; ++i) { int u, v; cin >> u >> v; g[u].app(v); g[v].app(u); } dfs1(1, 1); dfs2(1, 1); /* for (int i = 1; i <= n; ++i) cout << sum[i] << ' '; cout << endl; for (int i = 1; i <= n; ++i) cout << up[i] << ' '; cout << endl; debug(dv(1, 2)); */ int ans = 0; for (int i = 1; i <= n; ++i) { ans += solve(i); ans %= MOD; } cout << mod(ans * pw[n - 1]) << endl; }