結果

問題 No.1050 Zero (Maximum)
ユーザー 👑 rin204
提出日時 2022-04-17 17:10:17
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 122 ms / 2,000 ms
コード長 754 bytes
コンパイル時間 576 ms
コンパイル使用メモリ 82,304 KB
実行使用メモリ 71,424 KB
最終ジャッジ日時 2024-12-26 10:51:18
合計ジャッジ時間 3,006 ms
ジャッジサーバーID
(参考情報)
judge1 / judge2
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 3
other AC * 15
権限があれば一括ダウンロードができます

ソースコード

diff #

MOD = 10 ** 9 + 7

def matpow(A, B, w):
    l = len(A)
    while w:
        if w & 1:
            C = [0] * l
            for i in range(l):
                for j in range(l):
                    C[i] += A[i][j] * B[j]
                    C[i] %= MOD
            B = C
        C = [[0] * l for _ in range(l)]
        for i in range(l):
            for j in range(l):
                for k in range(l):
                    C[i][j] += A[i][k] * A[k][j]
                    C[i][j] %= MOD
        A = C
        w >>= 1
    return B

m, k = map(int, input().split())
B = [0] * m
B[0] = 1
A = [[0] * m for _ in range(m)]
for i in range(m):
    for j in range(m):
        A[(i + j) % m][j] += 1
        A[(i * j) % m][j] += 1

B = matpow(A, B, k)
print(B[0])

0