結果
問題 | No.1388 Less than K |
ユーザー |
|
提出日時 | 2020-11-23 08:28:26 |
言語 | PyPy3 (7.3.15) |
結果 |
WA
|
実行時間 | - |
コード長 | 1,304 bytes |
コンパイル時間 | 426 ms |
コンパイル使用メモリ | 82,240 KB |
実行使用メモリ | 91,524 KB |
最終ジャッジ日時 | 2024-11-26 06:54:03 |
合計ジャッジ時間 | 12,455 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 1 WA * 2 |
other | AC * 1 WA * 21 RE * 52 |
ソースコード
mod = 998244353 max_fact = 2*10**5 + 10 fact = [1] for i in range(1, max_fact): fact.append(fact[-1] * i % mod) fact_inv = [0] * max_fact fact_inv[-1] = pow(fact[-1], mod - 2, mod) for i in reversed(range(1, max_fact)): fact_inv[i - 1] = fact_inv[i] * i % mod def comb(a, b): if a < 0 or b < 0: return 0 elif a < b: return 0 return fact[a] * fact_inv[b] % mod * fact_inv[a - b] % mod h, w, k = map(int, input().split()) n = min(h, w) ans = 0 if k > 300: k2 = (k + 2) // 2 for i in range(n): s = comb(2*i, i) j = 1 while i - j * k2 >= 0: if j % 2 == 0: s += comb(2*i, i - j * k2) * 2 else: s -= comb(2*i, i - j * k2) * 2 j += 1 ans += s * comb(h+w-2, h-1-i) * comb(w+i-1, w-i-1) ans %= mod else: k2 = k // 2 dp1 = [0] * (k2 * 2 + 1) dp1[k2] = 1 for i in range(n): dp2 = [0] * (k2 * 2 + 1) for j in range(k2 * 2 + 1): if j != 0: dp2[j] += dp2[j-1] if i != 0 and j != k2 * 2: dp2[j] += dp1[j+1] dp2[j] %= mod dp1 = dp2 for i in range(n): ans += dp1[k2] * comb(h+w-2, h-1-i) * comb(w+i-1, w-i-1) ans %= mod print(ans)