結果
問題 |
No.1513 simple 門松列 problem
|
ユーザー |
|
提出日時 | 2021-07-20 18:06:47 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 2,186 ms / 3,000 ms |
コード長 | 1,243 bytes |
コンパイル時間 | 390 ms |
コンパイル使用メモリ | 82,304 KB |
実行使用メモリ | 281,392 KB |
最終ジャッジ日時 | 2024-07-17 13:44:58 |
合計ジャッジ時間 | 12,632 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 18 |
ソースコード
N, K = map(int, input().split()) MOD = 998244353 dp = [[[0, 0] for _ in range(K)] for _ in range(K)] for i in range(K): for j in range(K): if i!=j: dp[i][j] = [1, i+j] for _ in range(N-2): acc0 = [[0] for _ in range(K)] acc1 = [[0] for _ in range(K)] for i in range(K): for j in range(K): acc0[i].append((acc0[i][-1]+dp[j][i][0])%MOD) acc1[i].append((acc1[i][-1]+dp[j][i][1])%MOD) ndp = [[[0, 0] for _ in range(K)] for _ in range(K)] for i in range(K): for j in range(K): if i>j: ndp[i][j][0] += acc0[i][i]-dp[j][i][0] ndp[i][j][0] %= MOD ndp[i][j][1] += acc1[i][i]-dp[j][i][1]+j*ndp[i][j][0]%MOD ndp[i][j][1] %= MOD elif i<j: ndp[i][j][0] += acc0[i][-1]-acc0[i][i+1]-dp[j][i][0] ndp[i][j][0] %= MOD ndp[i][j][1] += acc1[i][-1]-acc1[i][i+1]-dp[j][i][1]+j*ndp[i][j][0]%MOD ndp[i][j][1] %= MOD dp = ndp ans0, ans1 = 0, 0 for i in range(K): for j in range(K): ans0 += dp[i][j][0] ans0 %= MOD ans1 += dp[i][j][1] ans1 %= MOD print(ans0, ans1)