結果

問題 No.1516 simple 門松列 problem Re:MASTER
ユーザー tktk_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
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 19
権限があれば一括ダウンロードができます

ソースコード

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