結果
問題 |
No.526 フィボナッチ数列の第N項をMで割った余りを求める
|
ユーザー |
|
提出日時 | 2024-10-05 10:19:29 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 43 ms / 2,000 ms |
コード長 | 801 bytes |
コンパイル時間 | 407 ms |
コンパイル使用メモリ | 82,308 KB |
実行使用メモリ | 52,736 KB |
最終ジャッジ日時 | 2024-10-05 10:19:31 |
合計ジャッジ時間 | 1,666 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 12 |
ソースコード
import sys sys.setrecursionlimit(10**8) def mat_dot_mod(mat1, mat2, mod): if len(set(map(len, mat1))) != 1 or len(set(map(len, mat2))) != 1: raise ValueError("wrong matrix") if len(mat1[0]) != len(mat2): raise ValueError("matrix size dose not match") return [[sum(v1 * v2 for v1, v2 in zip(r1, r2)) % mod for r2 in zip(*mat2)] for r1 in mat1] def mat_pow_mod(mat,p,mod): c = mat len_mat = len(mat) ans = [[1 if i == j else 0 for j in range(len(mat))]for i in range(len(mat))] for i in range(len(bin(p))-2): if 1 << i & p: ans = mat_dot_mod(ans,c,mod) c = mat_dot_mod(c,c,mod) return ans N, M = map(int, input().split()) mat = [[1, 1], [1, 0]] mt = mat_pow_mod(mat,N-2, M) print(mt[0][0]*1 + mt[0][1]*0)