結果
問題 | No.117 組み合わせの数 |
ユーザー | 草苺奶昔 |
提出日時 | 2023-02-24 19:46:36 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 240 ms / 5,000 ms |
コード長 | 1,473 bytes |
コンパイル時間 | 326 ms |
コンパイル使用メモリ | 82,372 KB |
実行使用メモリ | 107,856 KB |
最終ジャッジ日時 | 2024-09-13 03:53:36 |
合計ジャッジ時間 | 1,024 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge2 |
(要ログイン)
ソースコード
import sys input = sys.stdin.readline sys.setrecursionlimit(10 ** 7) class Combination: """ SIZEが10**6程度以下の二項係数を何回も呼び出したいときに使う 使い方: comb = Combination(SIZE, MOD) comb(10, 3) => 120 """ def __init__(self, N, MOD=10 ** 9 + 7): self.MOD = MOD self.__make_factorial_list(N) def __call__(self, n, k): if k < 0 or k > n: return 0 res = self.fact[n] * self.inv[k] % self.MOD res = res * self.inv[n - k] % self.MOD return res def nPk(self, n, k): if k < 0 or k > n: return 0 return self.fact[n] * self.inv[n - k] % self.MOD def nHk(self, n, k): if k == 0: return 1 return self.__call__(n + k - 1, k) def __make_factorial_list(self, N): self.fact = [1] * (N + 1) self.inv = [1] * (N + 1) MOD = self.MOD for i in range(1, N + 1): self.fact[i] = (self.fact[i - 1] * i) % MOD self.inv[N] = pow(self.fact[N], MOD - 2, MOD) for i in range(N, 0, -1): self.inv[i - 1] = (self.inv[i] * i) % MOD return U = 10 ** 6 comb = Combination(2 * U + 10) T = int(input()) for _ in range(T): s = input().rstrip() n, k = map(int, s[2:-1].split(",")) if s[0] == "C": print(comb(n, k)) elif s[0] == "P": print(comb.nPk(n, k)) else: print(comb.nHk(n, k))