結果
問題 | No.1044 正直者大学 |
ユーザー | 👑 SPD_9X2 |
提出日時 | 2020-05-01 22:56:29 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 131 ms / 2,000 ms |
コード長 | 1,332 bytes |
コンパイル時間 | 250 ms |
コンパイル使用メモリ | 87,096 KB |
実行使用メモリ | 94,092 KB |
最終ジャッジ日時 | 2023-08-26 15:37:56 |
合計ジャッジ時間 | 3,799 ms |
ジャッジサーバーID (参考情報) |
judge12 / judge14 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 66 ms
71,344 KB |
testcase_01 | AC | 64 ms
71,388 KB |
testcase_02 | AC | 74 ms
76,476 KB |
testcase_03 | AC | 67 ms
71,404 KB |
testcase_04 | AC | 62 ms
71,256 KB |
testcase_05 | AC | 97 ms
91,588 KB |
testcase_06 | AC | 131 ms
93,588 KB |
testcase_07 | AC | 64 ms
71,532 KB |
testcase_08 | AC | 63 ms
71,332 KB |
testcase_09 | AC | 64 ms
71,344 KB |
testcase_10 | AC | 64 ms
71,504 KB |
testcase_11 | AC | 64 ms
71,552 KB |
testcase_12 | AC | 93 ms
84,940 KB |
testcase_13 | AC | 115 ms
87,196 KB |
testcase_14 | AC | 98 ms
88,816 KB |
testcase_15 | AC | 106 ms
94,092 KB |
testcase_16 | AC | 113 ms
88,748 KB |
testcase_17 | AC | 92 ms
82,268 KB |
testcase_18 | AC | 97 ms
83,380 KB |
testcase_19 | AC | 85 ms
81,864 KB |
testcase_20 | AC | 74 ms
78,100 KB |
testcase_21 | AC | 84 ms
80,700 KB |
testcase_22 | AC | 79 ms
81,112 KB |
testcase_23 | AC | 86 ms
88,428 KB |
testcase_24 | AC | 79 ms
83,592 KB |
testcase_25 | AC | 82 ms
83,440 KB |
testcase_26 | AC | 83 ms
86,900 KB |
testcase_27 | AC | 64 ms
71,352 KB |
testcase_28 | AC | 65 ms
71,344 KB |
testcase_29 | AC | 63 ms
71,460 KB |
ソースコード
""" ループ内でのう→し し→うが 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)