結果
問題 | No.117 組み合わせの数 |
ユーザー | Shirotsume |
提出日時 | 2022-11-23 18:43:55 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 638 ms / 5,000 ms |
コード長 | 1,618 bytes |
コンパイル時間 | 200 ms |
コンパイル使用メモリ | 82,060 KB |
実行使用メモリ | 249,072 KB |
最終ジャッジ日時 | 2024-09-24 23:33:03 |
合計ジャッジ時間 | 1,961 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge1 |
(要ログイン)
ソースコード
import sys from collections import deque, Counter sys.setrecursionlimit(5 * 10 ** 5) from pypyjit import set_param set_param('max_unroll_recursion=-1') input = lambda: sys.stdin.readline().rstrip() ii = lambda: int(input()) mi = lambda: map(int, input().split()) li = lambda: list(mi()) inf = 2 ** 63 - 1 mod = 10 ** 9 + 7 # from Combinatorics import Combinatorics class Combinatorics(): def __init__(self, mod, maxi): self.mod = mod self.maxi = maxi self.facs = [1, 1] self.factinvs = [1, 1] self.invs = [0, 1] for i in range(2, self.maxi + 1): self.facs.append((self.facs[-1] * i) % self.mod) self.invs.append((-self.invs[self.mod % i] * (self.mod // i)) % self.mod) self.factinvs.append((self.factinvs[-1] * self.invs[-1]) % self.mod) def choose(self, n, k) -> int: if k < 0 or k > n: return 0 if k == 0 or k == n: return 1 k = min(k, n - k) return (((self.facs[n] * self.factinvs[k]) % self.mod) * self.factinvs[n-k]) % self.mod def perm(self, n, k) -> int: return (self.choose(n, k) * self.facs[k]) % self.mod def homop(self, n, k) -> int: if n == k == 0: return 1 return self.choose(n + k - 1, k) def factorial(self, n): return self.facs[n] C = Combinatorics(mod, 2 * 10 ** 6 + 10) for _ in range(ii()): s = input() x, y = map(int, s[2:-1].split(',')) if s[0] == 'C': print(C.choose(x, y)) if s[0] == 'P': print(C.perm(x, y)) if s[0] == 'H': print(C.homop(x, y))