結果

問題 No.1595 The Final Digit
コンテスト
ユーザー neterukun
提出日時 2021-07-09 21:36:19
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 44 ms / 2,000 ms
コード長 1,122 bytes
コンパイル時間 186 ms
コンパイル使用メモリ 82,324 KB
実行使用メモリ 59,776 KB
最終ジャッジ日時 2024-07-01 15:26:54
合計ジャッジ時間 1,752 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 17
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

def matrix_pow(A, k, MOD):
    def matrix_mul(A, B, MOD):
        C = [[0] * size for _ in range(size)]
        for i in range(size):
            for k in range(size):
                for j in range(size):
                    C[i][j] = (C[i][j] + A[i][k] * B[k][j]) % MOD
        return C

    if len(A) != len(A[0]):
        raise ValueError('square matrix expected')

    size = len(A)
    B = [[0] * size for _ in range(size)]
    for i in range(size):
        B[i][i] = 1
    while k > 0:
        if k & 1:
            B = matrix_mul(A, B, MOD)
        A = matrix_mul(A, A, MOD)
        k = k // 2
    return B


def matvec_mul(A, vec, MOD):
    if len(A[0]) != len(vec):
        raise ValueError('dimension mismatch between matrix and vector')

    h, w = len(A), len(A[0])
    res_vec = [0] * h
    for i in range(h):
        for j in range(w):
            res_vec[i] += A[i][j] * vec[j]
            res_vec[i] %= MOD
    return res_vec


p, q, r, k = map(int, input().split())
MOD = 10


mat = [[1, 1, 1], [1, 0, 0], [0, 1, 0]]


mat = matrix_pow(mat, k - 3, MOD)
res = matvec_mul(mat, [r, q, p], MOD)
print(res[0])
0