結果

問題 No.827 総神童数
ユーザー _____TAB__________TAB_____
提出日時 2019-05-03 23:07:04
言語 C++17
(gcc 13.2.0 + boost 1.83.0)
結果
AC  
実行時間 193 ms / 2,000 ms
コード長 1,233 bytes
コンパイル時間 865 ms
コンパイル使用メモリ 86,436 KB
実行使用メモリ 25,248 KB
最終ジャッジ日時 2023-08-30 06:37:30
合計ジャッジ時間 5,965 ms
ジャッジサーバーID
(参考情報)
judge14 / judge13
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
4,376 KB
testcase_01 AC 1 ms
4,376 KB
testcase_02 AC 2 ms
4,380 KB
testcase_03 AC 2 ms
4,380 KB
testcase_04 AC 2 ms
4,376 KB
testcase_05 AC 1 ms
4,376 KB
testcase_06 AC 1 ms
4,376 KB
testcase_07 AC 2 ms
4,376 KB
testcase_08 AC 2 ms
4,376 KB
testcase_09 AC 188 ms
25,248 KB
testcase_10 AC 60 ms
10,872 KB
testcase_11 AC 4 ms
4,380 KB
testcase_12 AC 25 ms
6,380 KB
testcase_13 AC 147 ms
20,960 KB
testcase_14 AC 6 ms
4,380 KB
testcase_15 AC 54 ms
10,092 KB
testcase_16 AC 142 ms
21,228 KB
testcase_17 AC 172 ms
24,280 KB
testcase_18 AC 69 ms
12,400 KB
testcase_19 AC 188 ms
25,228 KB
testcase_20 AC 132 ms
19,420 KB
testcase_21 AC 92 ms
14,868 KB
testcase_22 AC 150 ms
20,624 KB
testcase_23 AC 3 ms
4,376 KB
testcase_24 AC 174 ms
22,420 KB
testcase_25 AC 123 ms
18,044 KB
testcase_26 AC 174 ms
22,744 KB
testcase_27 AC 106 ms
16,264 KB
testcase_28 AC 103 ms
15,916 KB
testcase_29 AC 78 ms
13,560 KB
testcase_30 AC 23 ms
6,172 KB
testcase_31 AC 60 ms
11,184 KB
testcase_32 AC 60 ms
10,980 KB
testcase_33 AC 193 ms
24,520 KB
testcase_34 AC 172 ms
22,788 KB
testcase_35 AC 61 ms
11,240 KB
testcase_36 AC 94 ms
14,740 KB
testcase_37 AC 122 ms
17,980 KB
testcase_38 AC 193 ms
24,056 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

#include <queue>
#include <iostream>
#include <vector>
using namespace std;
const long long mod = 1e9+7;

long long mod_inv(long long x){
  long long k = mod - 2, ret = 1;
  while(k){
    if(k&1) ret *= x, ret %= mod;
    x *= x, x %= mod;
    k >>= 1;
  }
  return ret;
}

int main(){
  int N;
  cin >> N;
  vector<vector<int>> G(N);
  for(int i = 0; i < N-1; ++i){
    int u, v;
    cin >> u >> v;
    --u,--v;
    G[u].push_back(v);
    G[v].push_back(u);
  }
  vector<long long> F(N+2), F_(N+2);
  F[0] = 1;
  for(int i = 0; i <= N; ++i) F[i+1] = (i+1)*F[i]%mod;
  F_[N+1] = mod_inv(F[N+1]);
  for(int i = N+1; i > 0; --i) F_[i-1] = i*F_[i]%mod;
  auto f = [=](int a, int b){
             if(a < b) return 0LL;
             return F[a]*F_[b]%mod*F_[a-b]%mod;
           };
  auto g = [=](int k){
             return F[k]*F[N-k-1]%mod*((f(N+1,k+1)+mod-f(N,k))%mod)%mod;
           };
  queue<int> Q;
  Q.push(0);
  vector<int> D(N,-1);
  D[0] = 0;
  while(not Q.empty()){
    int v = Q.front();
    Q.pop();
    for(auto v_ : G[v]){
      if(D[v_] >= 0) continue;
      D[v_] = D[v]+1;
      Q.push(v_);
    }
  }
  long long ans = 0;
  for(int i = 0; i < N; ++i){
    ans += g(D[i]);
    ans %= mod;
  }
  cout << ans << endl;
}
0