結果
問題 |
No.827 総神童数
|
ユーザー |
|
提出日時 | 2020-09-16 03:03:57 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 191 ms / 2,000 ms |
コード長 | 1,170 bytes |
コンパイル時間 | 2,444 ms |
コンパイル使用メモリ | 196,060 KB |
最終ジャッジ日時 | 2025-01-14 15:08:43 |
ジャッジサーバーID (参考情報) |
judge1 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 36 |
ソースコード
#include <bits/stdc++.h> using namespace std; #define REP(i,n) for(int i=0; i<(int)(n); i++) const long long MOD = 1e9 + 7; int n; vector<int> edges[200000]; int depth[200000+1]; long long f[200000+1]; long long g[200000+1]; long long modpow(long long x, long long p, long long mod) { long long ret = 1; while (p) { if (p & 1) ret = ret * x % mod; x = x * x % mod; p >>= 1; } return ret; } long long comb(int n, int k) { return f[n] * g[n-k] % MOD * g[k] % MOD; } void dfs(int v, int par = -1, int d = 1) { ++depth[d]; for (int w: edges[v]) { if (w != par) dfs(w, v, d+1); } } int main() { ios_base::sync_with_stdio(0); cin.tie(0); f[0] = g[0] = 1; for (int i = 1; i <= 200000; i++) { f[i] = i * f[i-1] % MOD; g[i] = modpow(f[i], MOD-2, MOD); } cin >> n; REP (i, n-1) { int u, v; cin >> u >> v; --u, --v; edges[u].push_back(v); edges[v].push_back(u); } dfs(0); long long ret = 0; for (int d = 1; d <= n; d++) { long long tmp = depth[d] * comb(n, d) % MOD * f[d-1] % MOD * f[n-d] % MOD; ret = (ret + tmp) % MOD; } cout << ret << endl; return 0; }