結果

問題 No.1595 The Final Digit
ユーザー Tomii9273
提出日時 2021-07-10 03:46:56
言語 Python3
(3.13.1 + numpy 2.2.1 + scipy 1.14.1)
結果
AC  
実行時間 525 ms / 2,000 ms
コード長 667 bytes
コンパイル時間 79 ms
コンパイル使用メモリ 12,672 KB
実行使用メモリ 44,608 KB
最終ジャッジ日時 2024-07-01 23:04:32
合計ジャッジ時間 12,951 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 17
権限があれば一括ダウンロードができます

ソースコード

diff # 

import sys
import numpy as np
# sys.setrecursionlimit(10 ** 9)  # Codeforcesでは350000程度に
def input(): return sys.stdin.readline().strip()


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


mod = 10


def mat_mul(A, B, mod=mod):
    A1, A2 = A >> 15, A & (1 << 15) - 1
    B1, B2 = B >> 15, B & (1 << 15) - 1
    X = (A1 @ B1) % mod
    Y = (A2 @ B2) % mod
    Z = ((A1 + A2) @ (B1 + B2) - X - Y) % mod
    return ((X << 30) + (Z << 15) + Y) % mod


A = np.array([[r % mod], [q % mod], [p % mod]])
M = np.array([[1, 1, 1], [1, 0, 0], [0, 1, 0]])
k -= 3
while k > 0:
    if k % 2 == 1:
        A = mat_mul(M, A)
    M = mat_mul(M, M)
    k //= 2

print(A[0][0] % mod)
0