結果
問題 | No.140 みんなで旅行 |
ユーザー | _polarbear08 |
提出日時 | 2019-06-18 07:44:03 |
言語 | Python3 (3.12.2 + numpy 1.26.4 + scipy 1.12.0) |
結果 |
AC
|
実行時間 | 491 ms / 5,000 ms |
コード長 | 2,174 bytes |
コンパイル時間 | 87 ms |
コンパイル使用メモリ | 10,916 KB |
実行使用メモリ | 15,688 KB |
最終ジャッジ日時 | 2023-08-20 06:58:14 |
合計ジャッジ時間 | 4,226 ms |
ジャッジサーバーID (参考情報) |
judge11 / judge14 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 17 ms
8,144 KB |
testcase_01 | AC | 17 ms
8,048 KB |
testcase_02 | AC | 37 ms
8,524 KB |
testcase_03 | AC | 17 ms
8,192 KB |
testcase_04 | AC | 17 ms
8,076 KB |
testcase_05 | AC | 16 ms
8,060 KB |
testcase_06 | AC | 16 ms
8,044 KB |
testcase_07 | AC | 17 ms
8,048 KB |
testcase_08 | AC | 16 ms
8,212 KB |
testcase_09 | AC | 17 ms
8,196 KB |
testcase_10 | AC | 16 ms
8,112 KB |
testcase_11 | AC | 488 ms
15,664 KB |
testcase_12 | AC | 29 ms
8,520 KB |
testcase_13 | AC | 20 ms
8,108 KB |
testcase_14 | AC | 486 ms
15,688 KB |
testcase_15 | AC | 491 ms
15,472 KB |
testcase_16 | AC | 187 ms
11,112 KB |
testcase_17 | AC | 93 ms
9,524 KB |
testcase_18 | AC | 368 ms
13,880 KB |
testcase_19 | AC | 419 ms
14,288 KB |
testcase_20 | AC | 76 ms
9,276 KB |
testcase_21 | AC | 19 ms
8,040 KB |
ソースコード
import sys stdin = sys.stdin sys.setrecursionlimit(10 ** 7) def li(): return map(int, stdin.readline().split()) def li_(): return map(lambda x: int(x) - 1, stdin.readline().split()) def lf(): return map(float, stdin.readline().split()) def ls(): return stdin.readline().split() def ns(): return stdin.readline().rstrip() def lc(): return list(ns()) def ni(): return int(stdin.readline()) def nf(): return float(stdin.readline()) class CombTools(object): def __init__(self, cap: int, mod: int): self.cap = cap self.mod = mod self.inv = self._calc_inv() self.fac = self._calc_fac() self.fac_inv = self._calc_fac_inv() def _calc_inv(self): inv = [0, 1] for i in range(2, self.cap+1): inv.append((self.mod - (self.mod // i) * inv[self.mod % i]) % self.mod) return inv def _calc_fac(self): fac = [1] for i in range(1, self.cap+1): fac.append((i * fac[-1]) % self.mod) return fac def _calc_fac_inv(self): fac_inv = [1] for i in range(1, self.cap+1): fac_inv.append((self.inv[i] * fac_inv[-1]) % self.mod) return fac_inv def nCr(self, n: int, r: int): # validation if r > n: raise ValueError("n must be larger than r (n={}, r={})".format(n, r)) # calculation return self.fac[n] * self.fac_inv[n-r] * self.fac_inv[r] % self.mod def nPr(self, n: int, r: int): # validation if r > n: raise ValueError("n must be larger than r (n={}, r={})".format(n, r)) # calculation return self.fac[n] * self.fac_inv[n-r] % self.mod def nHr(self, n: int, r: int): return self.nCr(n+r-1, n) n = ni() MOD = 10**9 + 7 starling = [[0]*(n+1) for _ in range(n+1)] starling[0][0] = 1 for i in range(n): for j in range(n): starling[i+1][j+1] = starling[i][j] + (j+1)*starling[i][j+1] starling[i+1][j+1] %= MOD ct = CombTools(n, MOD) ans = 0 for y in range(1,n+1): for x in range(y, n+1): ans += ct.nCr(n, x) * starling[x][y] * pow(y*(y-1), n-x, MOD) ans %= MOD print(ans)