結果

問題 No.673 カブトムシ
ユーザー 双六
提出日時 2020-08-09 01:37:07
言語 Python3
(3.13.1 + numpy 2.2.1 + scipy 1.14.1)
結果
RE  
(最新)
AC  
(最初)
実行時間 -
コード長 630 bytes
コンパイル時間 97 ms
コンパイル使用メモリ 12,416 KB
実行使用メモリ 44,608 KB
最終ジャッジ日時 2024-10-02 12:34:04
合計ジャッジ時間 12,874 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample RE * 4
other RE * 14
権限があれば一括ダウンロードができます

ソースコード

diff # 

import sys; input = sys.stdin.buffer.readline
sys.setrecursionlimit(10**7)
from collections import defaultdict
con = 10 ** 9 + 7; INF = float("inf")

def getlist():
	return list(map(int, input().split()))

import numpy as np

def matrixPow(A, n, N, mod):
	B = np.identity(N, dtype = np.int)
	while n > 0:
		if n & 1 == 1:
			B = (A @ B) % mod
		A = (A @ A) % mod
		n >>= 1
	return B


#処理内容
def main():
	B, C, D = getlist()
	A = np.array([[C % con, B % con], [0, 1]])
	mat = matrixPow(A, D, 2, con)

	anspre = mat @ np.array([[B % con], [1]])
	ans = (anspre[0, 0] - B) % con
	print(ans)

if __name__ == '__main__':
	main()
0