結果
問題 | No.1513 simple 門松列 problem |
ユーザー | ayaoni |
提出日時 | 2021-05-22 05:57:43 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 615 ms / 3,000 ms |
コード長 | 1,564 bytes |
コンパイル時間 | 278 ms |
コンパイル使用メモリ | 82,228 KB |
実行使用メモリ | 205,024 KB |
最終ジャッジ日時 | 2024-10-10 10:40:16 |
合計ジャッジ時間 | 4,780 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge2 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 37 ms
53,560 KB |
testcase_01 | AC | 38 ms
53,024 KB |
testcase_02 | AC | 615 ms
205,024 KB |
testcase_03 | AC | 48 ms
63,336 KB |
testcase_04 | AC | 38 ms
53,904 KB |
testcase_05 | AC | 50 ms
65,304 KB |
testcase_06 | AC | 53 ms
64,828 KB |
testcase_07 | AC | 45 ms
63,516 KB |
testcase_08 | AC | 56 ms
67,124 KB |
testcase_09 | AC | 37 ms
53,288 KB |
testcase_10 | AC | 37 ms
53,776 KB |
testcase_11 | AC | 37 ms
54,068 KB |
testcase_12 | AC | 37 ms
52,956 KB |
testcase_13 | AC | 37 ms
52,808 KB |
testcase_14 | AC | 585 ms
204,312 KB |
testcase_15 | AC | 593 ms
204,432 KB |
testcase_16 | AC | 599 ms
204,792 KB |
testcase_17 | AC | 401 ms
159,452 KB |
testcase_18 | AC | 182 ms
105,020 KB |
testcase_19 | AC | 65 ms
71,460 KB |
testcase_20 | AC | 138 ms
91,456 KB |
ソースコード
import sys sys.setrecursionlimit(10**7) def I(): return int(sys.stdin.readline().rstrip()) def MI(): return map(int,sys.stdin.readline().rstrip().split()) def LI(): return list(map(int,sys.stdin.readline().rstrip().split())) def LI2(): return list(map(int,sys.stdin.readline().rstrip())) def S(): return sys.stdin.readline().rstrip() def LS(): return list(sys.stdin.readline().rstrip().split()) def LS2(): return list(sys.stdin.readline().rstrip()) N,K = MI() mod = 998244353 dp = [[[0]*K for _ in range(K)] for _ in range(N+1)] # dp[i][j][k] = 長さiで最後がk,jである門松列列の個数 S_dp = [[[0]*K for _ in range(K)] for _ in range(N+1)] # S_dp[i][j][k] = dp[i][j][0]+…+dp[i][j][k] for j in range(K): for k in range(K): if j != k: dp[2][j][k] = 1 if k == 0: S_dp[2][j][k] = dp[2][j][k] else: S_dp[2][j][k] = S_dp[2][j][k-1]+dp[2][j][k] S_dp[2][j][k] %= mod for i in range(3,N+1): for j in range(K): for k in range(K): if j < k: dp[i][j][k] = S_dp[i-1][k][k-1]-dp[i-1][k][j] elif j > k: dp[i][j][k] = S_dp[i-1][k][K-1]-S_dp[i-1][k][k]-dp[i-1][k][j] dp[i][j][k] %= mod if k == 0: S_dp[i][j][k] = dp[i][j][k] else: S_dp[i][j][k] = S_dp[i][j][k-1]+dp[i][j][k] S_dp[i][j][k] %= mod ans1 = 0 for j in range(K): ans1 += S_dp[-1][j][-1] ans1 %= mod ans2 = ans1*N*(K-1)*(mod+1)//2 ans2 %= mod print(ans1,ans2)