結果
問題 |
No.1044 正直者大学
|
ユーザー |
![]() |
提出日時 | 2021-01-14 00:06:16 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 81 ms / 2,000 ms |
コード長 | 1,166 bytes |
コンパイル時間 | 309 ms |
コンパイル使用メモリ | 82,252 KB |
実行使用メモリ | 80,984 KB |
最終ジャッジ日時 | 2024-11-22 20:59:08 |
合計ジャッジ時間 | 2,983 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 |
other | AC * 28 |
ソースコード
import sys sys.setrecursionlimit(10**7) def I(): return int(sys.stdin.readline().rstrip()) def MI(): return map(int,sys.stdin.readline().rstrip().split()) def LI(): return list(map(int,sys.stdin.readline().rstrip().split())) def LI2(): return list(map(int,sys.stdin.readline().rstrip())) def S(): return sys.stdin.readline().rstrip() def LS(): return list(sys.stdin.readline().rstrip().split()) def LS2(): return list(sys.stdin.readline().rstrip()) def factorial(N,p): # 0!~N!(mod p) のリスト res = [0]*(N+1) res[0] = 1 for i in range(1,N+1): res[i] = (res[i-1]*i) % p return res def factorial_inverse(N,p): # 0!~N!(mod p) の逆元のリスト fac = factorial(N,p) res = [0]*(N+1) a = pow(fac[-1],p-2,p) for i in range(N,-1,-1): res[i] = a a *= i a %= p return fac,res N,M,K = MI() mod = 10**9+7 fac,fac_inv = factorial_inverse(max(N,M),mod) def nCr(n,r): if n < r: return 0 return (fac[n]*fac_inv[r]*fac_inv[n-r]) % mod ans = 0 for i in range(1,(N+M-K)//2+1): ans += nCr(N,i)*nCr(M-1,i-1) ans %= mod ans *= fac[N-1] ans *= fac[M] ans %= mod print(ans)