結果

問題 No.1516 simple 門松列 problem Re:MASTER
ユーザー tktk_snsntktk_snsn
提出日時 2021-05-24 22:07:53
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 389 ms / 6,000 ms
コード長 1,169 bytes
コンパイル時間 295 ms
コンパイル使用メモリ 82,176 KB
実行使用メモリ 76,416 KB
最終ジャッジ日時 2024-10-13 06:23:09
合計ジャッジ時間 3,778 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 44 ms
58,368 KB
testcase_01 AC 76 ms
63,872 KB
testcase_02 AC 202 ms
76,032 KB
testcase_03 AC 47 ms
60,288 KB
testcase_04 AC 47 ms
60,544 KB
testcase_05 AC 49 ms
61,440 KB
testcase_06 AC 56 ms
62,464 KB
testcase_07 AC 76 ms
64,768 KB
testcase_08 AC 88 ms
65,792 KB
testcase_09 AC 48 ms
58,496 KB
testcase_10 AC 60 ms
62,208 KB
testcase_11 AC 49 ms
59,904 KB
testcase_12 AC 49 ms
60,800 KB
testcase_13 AC 49 ms
59,904 KB
testcase_14 AC 334 ms
76,032 KB
testcase_15 AC 193 ms
75,904 KB
testcase_16 AC 126 ms
72,704 KB
testcase_17 AC 83 ms
67,328 KB
testcase_18 AC 57 ms
62,336 KB
testcase_19 AC 56 ms
64,896 KB
testcase_20 AC 389 ms
75,904 KB
testcase_21 AC 381 ms
76,416 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

N, K = map(int, input().split())
mod = 998244353


def mat_mul(A, B):
    res = [[0] * len(B[0]) for _ in range(len(A))]
    for i in range(len(A)):
        for k in range(len(A[0])):
            for j in range(len(B[0])):
                res[i][j] += A[i][k] * B[k][j]
                res[i][j] %= mod
    return res


def mat_pow(A, n):
    size = len(A)
    res = [[0] * size for _ in range(size)]
    for i in range(size):
        res[i][i] = 1
    while n:
        if n & 1:
            res = mat_mul(res, A)
        A = mat_mul(A, A)
        n >>= 1
    return res


A = [[0]*(K*K) for _ in range(K*K)]
B = [0] * (K*K)
for x in range(K):
    for y in range(K):
        if x == y:
            continue
        if x < y:
            for z in range(y):
                if x != z:
                    A[x*K+y][y*K+z] = 1
        if x > y:
            for z in range(y+1, K):
                if x != z:
                    A[x*K+y][y*K+z] = 1
        B[x*K+y] = 1

P = mat_pow(A, N-2)
ans1 = 0
for i in range(K*K):
    for j in range(K*K):
        ans1 += P[i][j] * B[j] % mod
        ans1 %= mod

ans2 = ans1 * N * (K - 1) * pow(2, mod-2, mod) % mod
print(ans1, ans2)
0