結果
| 問題 |
No.1388 Less than K
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2020-11-23 00:41:48 |
| 言語 | PyPy3 (7.3.15) |
| 結果 |
RE
|
| 実行時間 | - |
| コード長 | 1,284 bytes |
| コンパイル時間 | 293 ms |
| コンパイル使用メモリ | 82,244 KB |
| 実行使用メモリ | 548,032 KB |
| 最終ジャッジ日時 | 2024-11-26 06:54:20 |
| 合計ジャッジ時間 | 17,072 ms |
|
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 22 RE * 51 MLE * 1 |
ソースコード
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
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)