結果
問題 |
No.1044 正直者大学
|
ユーザー |
|
提出日時 | 2020-05-02 18:44:52 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 92 ms / 2,000 ms |
コード長 | 1,437 bytes |
コンパイル時間 | 323 ms |
コンパイル使用メモリ | 82,472 KB |
実行使用メモリ | 76,168 KB |
最終ジャッジ日時 | 2024-12-31 00:22:36 |
合計ジャッジ時間 | 3,479 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 |
other | AC * 28 |
ソースコード
import sys sys.setrecursionlimit(10 ** 7) rl = sys.stdin.readline class Combination: def __init__(self, n: int, mod: int): self.mod = mod self.fact = [0] * (n + 1) self.factinv = [0] * (n + 1) self.inv = [0] * (n + 1) self.fact[0] = self.fact[1] = 1 self.factinv[0] = self.factinv[1] = 1 self.inv[1] = 1 for i in range(2, n + 1): self.fact[i] = (self.fact[i - 1] * i) % mod self.inv[i] = (-self.inv[mod % i] * (mod // i)) % mod self.factinv[i] = (self.factinv[i - 1] * self.inv[i]) % mod def ncr(self, n: int, r: int): if r < 0 or n < r: return 0 r = min(r, n - r) return self.fact[n] * self.factinv[r] % self.mod * self.factinv[n - r] % self.mod def nhr(self, n: int, r: int): return self.ncr(n + r - 1, r) def npr(self, n: int, r: int): if r < 0 or n < r: return 0 return self.fact[n] * self.factinv[n - r] % self.mod def solve(): N, M, K = map(int, rl().split()) MOD = 10 ** 9 + 7 com = Combination(200010, MOD) ans = 0 for i in range(1, N + 1): if N - i + M - i < K: continue ans = (ans + com.ncr(N, i) * com.nhr(i, M - i)) % MOD ans = ans * com.fact[N - 1] % MOD ans = ans * com.fact[M] % MOD print(ans) if __name__ == '__main__': solve()