結果
| 問題 |
No.1044 正直者大学
|
| コンテスト | |
| ユーザー |
tktk_snsn
|
| 提出日時 | 2020-12-25 15:38:01 |
| 言語 | PyPy3 (7.3.15) |
| 結果 |
AC
|
| 実行時間 | 66 ms / 2,000 ms |
| コード長 | 1,411 bytes |
| コンパイル時間 | 169 ms |
| コンパイル使用メモリ | 82,132 KB |
| 実行使用メモリ | 71,036 KB |
| 最終ジャッジ日時 | 2024-09-22 08:27:37 |
| 合計ジャッジ時間 | 2,731 ms |
|
ジャッジサーバーID (参考情報) |
judge2 / judge4 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 2 |
| other | AC * 28 |
ソースコード
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 __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
mod = 10 ** 9 + 7
N, M, K = map(int, input().split())
comb = Combination(N + M + 100, mod)
ans = 0
for i in range(1, N + 1): # Mをiこのグループに分ける
pair = N + M - 2 * i
if pair < K:
break
group = comb(M - 1, i - 1)
insert = comb(N, i)
ans += group * insert % mod
ans %= mod
ans *= comb.fact[N - 1]
ans %= mod
ans *= comb.fact[M]
ans %= mod
print(ans)
tktk_snsn