結果

問題 No.1124 Earthquake Safety
ユーザー QCFium
提出日時 2020-07-03 14:19:03
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 281 ms / 3,000 ms
コード長 1,657 bytes
コンパイル時間 1,934 ms
コンパイル使用メモリ 173,024 KB
実行使用メモリ 67,840 KB
最終ジャッジ日時 2024-09-16 17:08:52
合計ジャッジ時間 13,092 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 58
権限があれば一括ダウンロードができます

ソースコード

diff #
プレゼンテーションモードにする

#include <bits/stdc++.h>
int ri() {
    int n;
    scanf("%d", &n);
    return n;
}
#define MOD 1000000007
#define inv2 500000004
std::vector<std::vector<int> > hen;
std::vector<int> subtree_dsum;
void dfs0(int i, int prev) {
    subtree_dsum[i] = inv2;
    for (auto j : hen[i]) if (j != prev) {
        dfs0(j, i);
        subtree_dsum[i] = (subtree_dsum[i] + (int64_t) subtree_dsum[j] * inv2) % MOD;
    }
}
int res = 0;
void dfs1(int i, int prev, int prev_dsum) {
    std::vector<int> all;
    if (prev != -1) all.push_back(prev_dsum);
    for (auto j : hen[i]) if (j != prev) all.push_back(subtree_dsum[j]);
    
    int sum = std::accumulate(all.begin(), all.end(), 0LL) % MOD;
    int sum2 = std::accumulate(all.begin(), all.end(), 0, [] (int a, int b) { return (a + (int64_t) b * b) % MOD; });
    int sum3 = std::accumulate(all.begin(), all.end(), 0, [] (int a, int b) { return (a + (int64_t) b * b % MOD * b) % MOD; });
    
    res = (res
         + (int64_t) sum * sum % MOD * sum - (int64_t) sum2 * sum * 3 + sum3 * 2 // choose 3 from different arms
         + ((int64_t) sum * sum - sum2) * 3 // 1 self + 2 from different arms
         + (int64_t) sum * 3 + 1) % MOD; // 2 self + 1, 3 self
    if (res < 0) res += MOD;
    
    for (auto j : hen[i]) if (j != prev) {
        int next = (int64_t) (sum - subtree_dsum[j] + 1) * inv2 % MOD;
        if (next < 0) next += MOD;
        dfs1(j, i, next);
    }
}
int main() {
    int n = ri();
    hen.resize(n);
    for (int i = 1; i < n; i++) {
        int a = ri() - 1;
        int b = ri() - 1;
        hen[a].push_back(b);
        hen[b].push_back(a);
    }
    subtree_dsum.resize(n);
    dfs0(0, -1);
    
    dfs1(0, -1, 0);
    for (int i = 1; i < n; i++) {
        res += res;
        if (res >= MOD) res -= MOD;
    }
    printf("%d\n", res);
    return 0;
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0