結果
| 問題 |
No.1044 正直者大学
|
| コンテスト | |
| ユーザー |
👑 SPD_9X2
|
| 提出日時 | 2020-05-01 22:56:29 |
| 言語 | PyPy3 (7.3.15) |
| 結果 |
AC
|
| 実行時間 | 111 ms / 2,000 ms |
| コード長 | 1,332 bytes |
| コンパイル時間 | 247 ms |
| コンパイル使用メモリ | 82,180 KB |
| 実行使用メモリ | 92,192 KB |
| 最終ジャッジ日時 | 2024-12-25 13:59:44 |
| 合計ジャッジ時間 | 2,835 ms |
|
ジャッジサーバーID (参考情報) |
judge3 / judge1 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 2 |
| other | AC * 28 |
ソースコード
"""
ループ内でのう→し し→うが N+M-K以下ならよい
高橋君基準で考える
ブロック数がN+M-K以下ならよい
Xこのブロックに分けるとする
X//2このブロックはう
X-X//2このブロックはし
N,Mを少なくとも一つ分ける条件で上に分配
あとは (N-1)! * M! を掛ける
"""
def modfac(n, MOD):
f = 1
factorials = [1]
for m in range(1, n + 1):
f *= m
f %= MOD
factorials.append(f)
inv = pow(f, MOD - 2, MOD)
invs = [1] * (n + 1)
invs[n] = inv
for m in range(n, 1, -1):
inv *= m
inv %= MOD
invs[m - 1] = inv
return factorials, invs
def modnCr(n,r,mod,fac,inv): #上で求めたfacとinvsを引数に入れるべし(上の関数で与えたnが計算できる最大のnになる)
return fac[n] * inv[n-r] * inv[r] % mod
N,M,K = map(int,input().split())
ans = 0
mod = 10**9+7
fac,inv = modfac(N+M+10,mod)
for i in range(2,N+M+1):
B = i // 2
W = i - B
#print (B,W,i)
if B > M or W > N:
continue
if N+M-2*B < K:
continue
ans += modnCr(N-W + W-1 , W-1 , mod,fac,inv) * modnCr(M-B + B-1,B-1,mod,fac,inv)
ans %= mod
for i in range(1,N):
ans *= i
ans %= mod
for i in range(1,M+1):
ans *= i
ans %= mod
print (ans)
SPD_9X2