結果

問題 No.827 総神童数
ユーザー WarToksWarToks
提出日時 2019-05-04 10:19:22
言語 C++17(clang)
(17.0.6 + boost 1.83.0)
結果
AC  
実行時間 132 ms / 2,000 ms
コード長 5,546 bytes
コンパイル時間 4,249 ms
コンパイル使用メモリ 104,240 KB
実行使用メモリ 13,832 KB
最終ジャッジ日時 2023-08-20 12:23:35
合計ジャッジ時間 5,398 ms
ジャッジサーバーID
（参考情報） 		judge11 / judge14
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
4,376 KB
testcase_01 AC 2 ms
4,380 KB
testcase_02 AC 2 ms
4,380 KB
testcase_03 AC 1 ms
4,376 KB
testcase_04 AC 2 ms
4,380 KB
testcase_05 AC 2 ms
4,380 KB
testcase_06 AC 2 ms
4,376 KB
testcase_07 AC 2 ms
4,384 KB
testcase_08 AC 2 ms
4,376 KB
testcase_09 AC 132 ms
13,832 KB
testcase_10 AC 38 ms
7,028 KB
testcase_11 AC 3 ms
4,376 KB
testcase_12 AC 17 ms
4,684 KB
testcase_13 AC 100 ms
11,668 KB
testcase_14 AC 5 ms
4,376 KB
testcase_15 AC 37 ms
6,536 KB
testcase_16 AC 93 ms
11,748 KB
testcase_17 AC 107 ms
13,316 KB
testcase_18 AC 47 ms
7,528 KB
testcase_19 AC 110 ms
13,760 KB
testcase_20 AC 81 ms
11,072 KB
testcase_21 AC 47 ms
8,824 KB
testcase_22 AC 78 ms
11,748 KB
testcase_23 AC 2 ms
4,384 KB
testcase_24 AC 87 ms
12,528 KB
testcase_25 AC 66 ms
10,428 KB
testcase_26 AC 99 ms
12,788 KB
testcase_27 AC 52 ms
9,376 KB
testcase_28 AC 57 ms
9,380 KB
testcase_29 AC 44 ms
8,096 KB
testcase_30 AC 13 ms
4,688 KB
testcase_31 AC 32 ms
6,984 KB
testcase_32 AC 31 ms
6,696 KB
testcase_33 AC 104 ms
13,516 KB
testcase_34 AC 95 ms
12,516 KB
testcase_35 AC 33 ms
7,036 KB
testcase_36 AC 48 ms
8,944 KB
testcase_37 AC 60 ms
10,456 KB
testcase_38 AC 105 ms
13,188 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <list>
#include <set>
#include <deque>
#include <queue>
#include <stack>
#include <vector>
#include <algorithm>
#include <map>
#include <chrono>
#include <math.h>
using namespace std;

using lli = long long int;
using Vint = std::vector<int>;
using Vlli = std::vector<lli>;
using Wint = std::vector<Vint>;
using Wlli = std::vector<Vlli>;
using Vbool = std::vector<bool>;
using pii = std::pair<int, int>;
using pll = std::pair<lli, lli>;
template <class T>
using Vec = std::vector<T>;

constexpr int MOD = 1e9 + 7;
constexpr int INFi = 2e9 + 1;
constexpr lli INFl = (lli)(9e18) + 1;
const vector<pii> DXDY = {std::make_pair(1, 0), std::make_pair(-1, 0), std::make_pair(0, 1), std::make_pair(0, -1)};
constexpr char BR = '\n';

#define FOR(i, a, b) for(int (i) = (a); (i) < (b); (i)++)
#define FOReq(i, a, b) for(int (i) = (a); (i) <= (b); (i)++)
#define rFOR(i, a, b) for(int (i) = (b); (i) >= (a); i--)
#define FORstep(i, a, b, step) for(int (i) = (a); i < (b); i += (step))
#define REP(i, n) FOR(i, 0, n)
#define rREP(i, n) rFOR(i, 0, (n-1))
#define vREP(ele, vec) for(auto &(ele) : (vec))
#define vREPcopy(ele, vec) for(auto (ele) : (vec))
#define SORT(A) std::sort((A).begin(), (A).end())
#define rSORT(A) std::sort((A).rbegin(), (A).rend())
// 座標圧縮 (for vector) : ソートしてから使うのが一般的 ; SORT(A) => COORDINATE_COMPRESSION(A)
#define COORDINATE_COMPRESSION(A) (A).erase(unique((A).begin(),(A).end()),(A).end())



template <class T> inline int argmin(std::vector<T> vec){return min_element(vec.begin(), vec.end()) - vec.begin();}
template <class T> inline int argmax(std::vector<T> vec){return max_element(vec.begin(), vec.end()) - vec.begin();}
template <class T> inline void chmax(T &a, T b){a = max(a, b);}
template <class T> inline void chmin(T &a, T b){a = min(a, b);}
inline int BitI(int k){return 1 << k;}
inline lli BitL(int k){return 1L << k;}
inline int toInt(string &s){int res = 0; for(char a : s) res = 10 * res + (a - '0'); return res;}
inline int toInt(const string s){int res = 0; for(char a : s) res = 10 * res + (a - '0'); return res;}
inline long long int toLong(string &s){lli res = 0; for(char a : s) res = 10 * res + (a - '0'); return res;}
inline long long int toLong(const string s){lli res = 0; for(char a : s) res = 10 * res + (a - '0'); return res;}
template <class T> inline std::string toString(T n){
  if(n == 0) return "0";
  std::string res = "";
  if(n < 0){n = -n;while(n != 0){res += (char)(n % 10 + '0'); n /= 10;}
  std::reverse(res.begin(), res.end()); return '-' + res;}
  while(n != 0){res += (char)(n % 10 + '0'); n /= 10;} std::reverse(res.begin(), res.end()); return res;
}

// ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~


template <long long int mod_number__ = MOD>
struct Modlli{
  using lli = long long int;
  lli number;
  // +, -, *, /, +=, -=, *=, /= が Modについて定義する
  inline Modlli operator +(Modlli that) const {return Modlli((this->number + that.number) % mod_number__);}
  inline void operator +=(Modlli that){this->number = (this->number + that.number) % mod_number__;}
  inline void operator ++(){this->number = (this->number + 1) % mod_number__;}
  inline Modlli operator -(Modlli that) const {return Modlli((this->number - that.number + mod_number__) % mod_number__);}
  inline void operator -=(Modlli that){this->number = (this->number + that.number + mod_number__) % mod_number__;}
  inline void operator --(){this->number = (this->number - 1) % mod_number__;}
  inline Modlli operator *(Modlli that) const {return Modlli((this->number * that.number) % mod_number__); }
  inline void operator *=(const Modlli that){this->number = (this->number * that.number) % mod_number__;}
  inline Modlli power(lli n) const { // 冪乗 : number^nを計算
    lli res = 1LL, waiting = this->number;
    while(n != 0LL){
      if((n & 1LL) != 0LL) res = (res * waiting) % mod_number__;
      waiting = waiting * waiting % mod_number__; n >>= 1;
    }
    return Modlli(res);
  }
  inline Modlli inv(void) const {return this->power(mod_number__ - 2);}
  inline Modlli operator /(Modlli that) const{ return (*this) * that.inv();}
  inline void operator /=(Modlli that) {
    lli n = mod_number__ - 2;
    lli res = 1LL, waiting = that.number;
    while(n != 0LL){
      if((n & 1LL) != 0LL) res = (res * waiting) % mod_number__;
      waiting = waiting * waiting % mod_number__; n >>= 1;
    }
    this->snumber = (this->number * res) % MOD;
  }
  inline Modlli operator %(Modlli that){return Modlli(this->number % that.number);}
  inline Modlli operator %=(Modlli that){this->number %= that.number;}
  Modlli(lli mn){ number = mn % mod_number__;}
  ~Modlli(void){}
  //operator long long int() const {return number;}
  lli toLong(void) const {return number;}
};

using modint = Modlli<MOD>;

int main(void){
  int n; scanf("%d", &n);Wint Tree(n, Vint(0));
  REP(_, n-1){
    int a, b; scanf("%d%d", &a, &b); a--; b--;
    Tree[a].emplace_back(b); Tree[b].emplace_back(a);
  }
  Vbool hasDone(n, false);
  queue<int> A; A.push(0); hasDone[0] = true;
  modint res(0);
  modint kaijou(1);
  int depth = 1;
  for(int i = 2; i <= n; i++) kaijou *= modint(i);
  while(not A.empty()){
    int sz = A.size();
    REP(_, sz){
      int idx = A.front(); A.pop();
      vREP(ele, Tree[idx]){
        if(not hasDone[ele]){
          hasDone[ele] = true; A.push(ele);
        }
      }
    }
    res += (modint(sz) / modint(depth));
    depth++;
  }
  cout << (res * kaijou).toLong() << BR;
  return 0;
}
0