結果

問題 No.344 ある無理数の累乗
ユーザー 👑 rin204
提出日時 2022-11-09 21:30:08
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 39 ms / 2,000 ms
コード長 651 bytes
コンパイル時間 167 ms
コンパイル使用メモリ 82,436 KB
実行使用メモリ 53,892 KB
最終ジャッジ日時 2024-07-23 02:37:44
合計ジャッジ時間 2,506 ms
ジャッジサーバーID
(参考情報)
judge4 / judge1
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
other AC * 30
権限があれば一括ダウンロードができます

ソースコード

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

MOD = 1000
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
n = int(input())
A = [[1, 3], [1, 1]]
B = [1, 0]
ans = matpow(A, B, n)[0] * 2
if n % 2 == 0:
ans -= 1
ans %= MOD
print(ans)
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0