結果
| 問題 |
No.526 フィボナッチ数列の第N項をMで割った余りを求める
|
| コンテスト | |
| ユーザー |
dn6049949
|
| 提出日時 | 2020-03-29 22:24:17 |
| 言語 | PyPy3 (7.3.15) |
| 結果 |
AC
|
| 実行時間 | 41 ms / 2,000 ms |
| コード長 | 1,071 bytes |
| コンパイル時間 | 408 ms |
| コンパイル使用メモリ | 82,596 KB |
| 実行使用メモリ | 54,476 KB |
| 最終ジャッジ日時 | 2025-01-02 15:53:21 |
| 合計ジャッジ時間 | 1,739 ms |
|
ジャッジサーバーID (参考情報) |
judge4 / judge1 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 12 |
ソースコード
# 行列累乗 O(H^3logN) (H := matrix Aの次元)
# 内積
def product(a, b, mod=None, add=lambda x,y:x+y, mul=lambda x,y:x*y):
if mod is None:
res = 0
for i in map(lambda x:mul(x[0],x[1]), zip(a,b)):
res = add(res,i)
return res
else:
res = 0
for i in map(lambda x:mul(x[0],x[1])%mod, zip(a,b)):
res = add(res,i)%mod
return res
# 行列積
def mul_of_matrix(a, b, mod=None, add=lambda x,y:x+y, mul=lambda x,y:x*y):
bt = [[b[i][j] for i in range(len(b))] for j in range(len(b[0]))]
return [[product(ai, bj, mod, add, mul) for bj in bt] for ai in a]
# 行列累乗
def pow_of_matrix(a, n, mod=None, add=lambda x,y:x+y, mul=lambda x,y:x*y, default=1):
res = [[default if i == j else 0 for j in range(len(a))] for i in range(len(a))]
while n:
if n&1:
res = mul_of_matrix(res,a,mod,add,mul)
a = mul_of_matrix(a,a,mod,add,mul)
n >>= 1
return res
n,m = map(int, input().split())
mat = pow_of_matrix([[1,1],[1,0]],n-1,m)
print(mat[1][0])
dn6049949