結果
問題 |
No.526 フィボナッチ数列の第N項をMで割った余りを求める
|
ユーザー |
![]() |
提出日時 | 2021-04-04 14:14:42 |
言語 | Python3 (3.13.1 + numpy 2.2.1 + scipy 1.14.1) |
結果 |
AC
|
実行時間 | 33 ms / 2,000 ms |
コード長 | 883 bytes |
コンパイル時間 | 262 ms |
コンパイル使用メモリ | 12,672 KB |
実行使用メモリ | 10,624 KB |
最終ジャッジ日時 | 2024-12-27 15:31:48 |
合計ジャッジ時間 | 1,711 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 12 |
ソースコード
N, M = map(int, input().split()) N -= 1 MOD = M def calc_matrix_product(A, B): if not len(A[0]) == len(B): raise Exception("Invalid arguments. A_col: {} B_row: {}".format(len(A), len(B))) C = [[0] * len(B[0]) for _ in range(len(A))] for i in range(len(A)): for j in range(len(B)): for k in range(len(B[0])): C[i][j] = (C[i][j] + A[i][k] * B[k][j]) % MOD return C def calc_matrix_pow(A, n): if not len(A) == len(A[0]): raise Exception("Invalid arguments. A_row: {} A_col: {}".format(len(A), len(A[0]))) B = [[0] * len(A) for _ in range(len(A))] for i in range(len(A)): B[i][i] = 1 while n: if n & 1: B = calc_matrix_product(B, A) A = calc_matrix_product(A, A) n = n >> 1 return B A = [[1, 1], [1, 0]] print(calc_matrix_pow(A, N)[1][0] % MOD)