結果
| 問題 |
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 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| other | AC * 1 |
ソースコード
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))