結果

問題 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
権限があれば一括ダウンロードができます

ソースコード

diff #

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)
0