結果

問題 No.1275 綺麗な式
ユーザー ygd.
提出日時 2020-11-25 00:38:38
言語 Python3
(3.13.1 + numpy 2.2.1 + scipy 1.14.1)
結果
AC  
実行時間 512 ms / 2,000 ms
コード長 395 bytes
コンパイル時間 88 ms
コンパイル使用メモリ 12,544 KB
実行使用メモリ 44,476 KB
最終ジャッジ日時 2024-07-23 18:59:21
合計ジャッジ時間 34,368 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 60
権限があれば一括ダウンロードができます

ソースコード

diff #
プレゼンテーションモードにする

import numpy as np
MOD = pow(10,9)+7
a,b = map(int,input().split())
n = int(input())
A = np.array([[a, b], [1, a]])
#print(A)
MAX = 70
L = [A]
for i in range(MAX):
  A = np.dot(A,A)
  A = np.mod(A, MOD)
  #print(A)
  L.append(A)
#print(L)
#print(B)
X = np.array([[1, 0], [0, 1]])
for i in range(MAX):
  if (n>>i)&1 == 1:
    X = np.dot(X,L[i])
    X = np.mod(X,MOD)
ans = X[0,0]*2
print(ans%MOD)
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0