結果

問題 No.1513 simple 門松列 problem
ユーザー rlangevinrlangevin
提出日時 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
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 37 ms
52,652 KB
testcase_01 AC 39 ms
53,468 KB
testcase_02 AC 812 ms
88,652 KB
testcase_03 AC 54 ms
66,732 KB
testcase_04 AC 38 ms
53,352 KB
testcase_05 AC 56 ms
67,976 KB
testcase_06 AC 61 ms
71,600 KB
testcase_07 AC 49 ms
64,416 KB
testcase_08 AC 61 ms
69,004 KB
testcase_09 AC 37 ms
52,736 KB
testcase_10 AC 38 ms
54,068 KB
testcase_11 AC 38 ms
52,988 KB
testcase_12 AC 38 ms
53,668 KB
testcase_13 AC 37 ms
53,404 KB
testcase_14 AC 804 ms
88,112 KB
testcase_15 AC 811 ms
88,336 KB
testcase_16 AC 821 ms
88,768 KB
testcase_17 AC 550 ms
83,928 KB
testcase_18 AC 322 ms
80,584 KB
testcase_19 AC 94 ms
76,576 KB
testcase_20 AC 190 ms
78,236 KB
権限があれば一括ダウンロードができます

ソースコード

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