結果

問題 No.1513 simple 門松列 problem
ユーザー rlangevinrlangevin
提出日時 2023-09-30 12:56:03
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 913 ms / 3,000 ms
コード長 1,843 bytes
コンパイル時間 277 ms
コンパイル使用メモリ 86,968 KB
実行使用メモリ 90,100 KB
最終ジャッジ日時 2023-09-30 12:56:12
合計ジャッジ時間 7,674 ms
ジャッジサーバーID
（参考情報） 		judge15 / judge12
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 72 ms
71,268 KB
testcase_01 AC 74 ms
71,096 KB
testcase_02 AC 876 ms
89,716 KB
testcase_03 AC 87 ms
76,736 KB
testcase_04 AC 78 ms
71,240 KB
testcase_05 AC 89 ms
76,556 KB
testcase_06 AC 100 ms
77,672 KB
testcase_07 AC 85 ms
76,904 KB
testcase_08 AC 93 ms
77,296 KB
testcase_09 AC 69 ms
71,336 KB
testcase_10 AC 71 ms
71,608 KB
testcase_11 AC 72 ms
71,364 KB
testcase_12 AC 71 ms
71,596 KB
testcase_13 AC 69 ms
71,540 KB
testcase_14 AC 885 ms
89,444 KB
testcase_15 AC 854 ms
89,584 KB
testcase_16 AC 913 ms
90,100 KB
testcase_17 AC 593 ms
84,616 KB
testcase_18 AC 376 ms
82,984 KB
testcase_19 AC 128 ms
77,680 KB
testcase_20 AC 230 ms
79,624 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