結果

問題 No.1124 Earthquake Safety
コンテスト
ユーザー QCFium
提出日時 2020-07-03 14:19:03
言語 C++14
(gcc 13.3.0 + boost 1.89.0)
結果
AC  
実行時間 281 ms / 3,000 ms
コード長 1,657 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 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 #
raw source code 

#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;
}
0