結果
問題 | No.1388 Less than K |
ユーザー |
|
提出日時 | 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 = 998244353max_fact = 4*10**5 + 10fact = [1]for i in range(1, max_fact):fact.append(fact[-1] * i % mod)fact_inv = [0] * max_factfact_inv[-1] = pow(fact[-1], mod - 2, mod)for i in reversed(range(1, max_fact)):fact_inv[i - 1] = fact_inv[i] * i % moddef comb(a, b):if a < 0 or b < 0:return 0elif a < b:return 0return fact[a] * fact_inv[b] % mod * fact_inv[a - b] % modh, w, k = map(int, input().split())n = min(h, w)ans = 0if k > 300:k2 = (k + 2) // 2for i in range(n):s = comb(2*i, i)j = 1while i - j * k2>= 0:if j % 2 == 0:s += comb(2*i, i - j * k2) * 2else:s -= comb(2*i, i - j * k2) * 2j += 1ans += s * comb(h+w-2, h-1-i) * comb(w+i-1, w-i-1)ans %= modelse:k2 = k // 2dp = [[0] * n for _ in range(k2 * 2 + 1)]dp[k2][0] = 1for 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] %= modfor i in range(n):ans += dp[k2][i] * comb(h+w-2, h-1-i) * comb(w+i-1, w-i-1)ans %= modprint(ans)