結果
問題 | No.117 組み合わせの数 |
ユーザー | Mao-beta |
提出日時 | 2024-03-08 16:14:48 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 308 ms / 5,000 ms |
コード長 | 2,026 bytes |
コンパイル時間 | 198 ms |
コンパイル使用メモリ | 82,432 KB |
実行使用メモリ | 124,032 KB |
最終ジャッジ日時 | 2024-09-29 18:48:50 |
合計ジャッジ時間 | 1,284 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge5 |
(要ログイン)
ソースコード
import sys import math import bisect from heapq import heapify, heappop, heappush from collections import deque, defaultdict, Counter from functools import lru_cache from itertools import accumulate, combinations, permutations, product sys.setrecursionlimit(1000000) MOD = 10 ** 9 + 7 MOD99 = 998244353 input = lambda: sys.stdin.readline().strip() NI = lambda: int(input()) NMI = lambda: map(int, input().split()) NLI = lambda: list(NMI()) SI = lambda: input() SMI = lambda: input().split() SLI = lambda: list(SMI()) EI = lambda m: [NLI() for _ in range(m)] class Comb: """nCrのnもrも10**7くらいまで""" def __init__(self, n, mod): self.mod = mod self.fac = [1] * (n + 1) self.inv = [1] * (n + 1) for i in range(1, n + 1): self.fac[i] = self.fac[i - 1] * i % self.mod self.inv[n] = pow(self.fac[n], self.mod - 2, self.mod) for i in range(n - 1, 0, -1): self.inv[i] = self.inv[i + 1] * (i + 1) % self.mod def C(self, n, r): if n < r: return 0 if n < 0 or r < 0: return 0 return self.fac[n] * self.inv[r] % self.mod * self.inv[n - r] % self.mod def P(self, n, r): if n < r: return 0 if n < 0 or r < 0: return 0 return self.fac[n] * self.inv[n - r] % self.mod def H(self, n, r): """ n個のものから重複を許してr個取り出す """ if n == r == 0: return 1 return self.C(n + r - 1, r) def multi(self, L): res = self.fac[sum(L)] for l in L: res = res * self.inv[l] % self.mod return res def main(): T = NI() com = Comb(10**6*3, MOD) for _ in range(T): S = SI() a = int(S.split("(")[1].split(",")[0]) b = int(S.split(",")[1].split(")")[0]) if S[0] == "C": print(com.C(a, b)) elif S[0] == "P": print(com.P(a, b)) else: print(com.H(a, b)) if __name__ == "__main__": main()