結果
問題 | No.1388 Less than K |
ユーザー | Mitarushi |
提出日時 | 2020-11-23 08:31:50 |
言語 | PyPy3 (7.3.15) |
結果 |
MLE
(最新)
AC
(最初)
|
実行時間 | - |
コード長 | 1,284 bytes |
コンパイル時間 | 221 ms |
コンパイル使用メモリ | 82,048 KB |
実行使用メモリ | 553,384 KB |
最終ジャッジ日時 | 2024-11-26 06:56:56 |
合計ジャッジ時間 | 33,262 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 73 MLE * 1 |
ソースコード
mod = 998244353 max_fact = 4*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 dp = [[0] * n for _ in range(k2 * 2 + 1)] dp[k2][0] = 1 for i in range(n): for j in range(k2 * 2 + 1): if j != 0: dp[j][i] += dp[j-1][i] if i != 0 and j != k2 * 2: dp[j][i] += dp[j+1][i-1] dp[j][i] %= mod for i in range(n): ans += dp[k2][i] * comb(h+w-2, h-1-i) * comb(w+i-1, w-i-1) ans %= mod print(ans)