結果
| 問題 |
No.1513 simple 門松列 problem
|
| コンテスト | |
| ユーザー |
rlangevin
|
| 提出日時 | 2023-09-30 12:56:03 |
| 言語 | PyPy3 (7.3.15) |
| 結果 |
AC
|
| 実行時間 | 821 ms / 3,000 ms |
| コード長 | 1,843 bytes |
| コンパイル時間 | 313 ms |
| コンパイル使用メモリ | 82,500 KB |
| 実行使用メモリ | 88,768 KB |
| 最終ジャッジ日時 | 2024-07-23 06:58:05 |
| 合計ジャッジ時間 | 6,192 ms |
|
ジャッジサーバーID (参考情報) |
judge3 / judge1 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 18 |
ソースコード
N, K = map(int, input().split())
mod = 998244353
pre = [[0] * K for _ in range(K)]
pre_sum = [[0] * K for _ in range(K)]
Ac = [[0] * (K + 1) for _ in range(K + 1)]
Ac_sum = [[0] * (K + 1) for _ in range(K + 1)]
for a in range(K):
for b in range(K):
for c in range(K):
if a == b or b == c or a == c:
continue
if a < b < c or c < b < a:
continue
pre[b][c] += 1
pre_sum[b][c] += a + b + c
for a in range(K):
for b in range(K):
Ac[b][a + 1] = Ac[b][a] + pre[a][b]
Ac_sum[b][a + 1] = Ac_sum[b][a] + pre_sum[a][b]
Ac[b][a + 1] %= mod
Ac_sum[b][a + 1] %= mod
for _ in range(N - 3):
dp = [[0] * K for _ in range(K)]
dp_sum = [[0] * K for _ in range(K)]
for a in range(K):
for b in range(K):
if a == b:
continue
if a < b:
dp[a][b] += Ac[a][-1] - Ac[a][a + 1] - pre[b][a]
dp_sum[a][b] += (Ac_sum[a][-1] - Ac_sum[a][a + 1]) + b * (Ac[a][-1] - Ac[a][a + 1] - pre[b][a]) - pre_sum[b][a]
else:
dp[a][b] += Ac[a][a] - pre[b][a]
dp_sum[a][b] += Ac_sum[a][a] + b * (Ac[a][a] - pre[b][a]) - pre_sum[b][a]
dp[a][b] %= mod
dp_sum[a][b] %= mod
dp, pre = pre, dp
dp_sum, pre_sum = pre_sum, dp_sum
for a in range(K):
for b in range(K):
Ac[b][a + 1] = Ac[b][a] + pre[a][b]
Ac_sum[b][a + 1] = Ac_sum[b][a] + pre_sum[a][b]
Ac[b][a + 1] %= mod
Ac_sum[b][a + 1] %= mod
ans1, ans2 = 0, 0
for a in range(K):
for b in range(K):
ans1 += pre[a][b]
ans2 += pre_sum[a][b]
ans1 %= mod
ans2 %= mod
print(ans1, ans2)
rlangevin