結果
| 問題 |
No.526 フィボナッチ数列の第N項をMで割った余りを求める
|
| コンテスト | |
| ユーザー |
bal4u
|
| 提出日時 | 2019-04-14 08:48:15 |
| 言語 | C (gcc 13.3.0) |
| 結果 |
AC
|
| 実行時間 | 2 ms / 2,000 ms |
| コード長 | 928 bytes |
| コンパイル時間 | 174 ms |
| コンパイル使用メモリ | 29,696 KB |
| 実行使用メモリ | 5,376 KB |
| 最終ジャッジ日時 | 2024-09-16 08:23:32 |
| 合計ジャッジ時間 | 788 ms |
|
ジャッジサーバーID (参考情報) |
judge4 / judge3 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 12 |
ソースコード
// yukicoder: No.526 フィボナッチ数列の第N項をMで割った余りを求める
// 2019.4.14 bal4u
#include <stdio.h>
int MOD;
//// fibonacci数の高速計算
void multiply(int F[2][2], const int M[2][2])
{
long long f00, f01, f10, f11;
int x, y, z, w;
f00 = F[0][0], f01 = F[0][1], f10 = F[1][0], f11 = F[1][1];
x = (f00 * M[0][0] + f01 * M[1][0]) % MOD;
y = (f00 * M[0][1] + f01 * M[1][1]) % MOD;
z = (f10 * M[0][0] + f11 * M[1][0]) % MOD;
w = (f10 * M[0][1] + f11 * M[1][1]) % MOD;
F[0][0] = x, F[0][1] = y, F[1][0] = z, F[1][1] = w;
}
void power(int F[2][2], int n)
{
const int M[2][2] = { {1,1},{1,0} };
if (n == 0 || n == 1) return;
power(F, n >> 1);
multiply(F, F);
if (n & 1) multiply(F, M);
}
int fib(int n)
{
int F[2][2] = { {1,1},{1,0} };
if (n == 0) return 0;
power(F, n - 1);
return F[0][0];
}
int main()
{
int N;
scanf("%d%d", &N, &MOD);
printf("%d\n", fib(N-1));
return 0;
}
bal4u