結果
| 問題 |
No.718 行列のできるフィボナッチ数列道場 (1)
|
| コンテスト | |
| ユーザー |
bal4u
|
| 提出日時 | 2019-04-19 18:18:11 |
| 言語 | C (gcc 13.3.0) |
| 結果 |
AC
|
| 実行時間 | 1 ms / 2,000 ms |
| コード長 | 974 bytes |
| コンパイル時間 | 343 ms |
| コンパイル使用メモリ | 30,464 KB |
| 実行使用メモリ | 6,944 KB |
| 最終ジャッジ日時 | 2024-09-22 11:42:02 |
| 合計ジャッジ時間 | 1,407 ms |
|
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 20 |
ソースコード
// yukicoder: No.718 行列のできるフィボナッチ数列道場 (1)
// 2019.4.19 bal4u
#include <stdio.h>
#define MOD 1000000007
//// fibonacci数の高速計算
void multiply(int F[2][2], int M[2][2])
{
long long f00, f01, f10, f11, 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] = (int)x, F[0][1] = (int)y;
F[1][0] = (int)z, F[1][1] = (int)w;
}
int M[2][2] = { {1,1},{1,0} };
void power(int F[2][2], long long n)
{
if (n == 0 || n == 1) return;
power(F, n >> 1);
multiply(F, F);
if (n & 1) multiply(F, M);
}
int fib(long long n)
{
int F[2][2] = { {1,1},{1,0} };
if (n == 0) return 0;
power(F, n - 1);
return F[0][0];
}
int main()
{
long long N;
scanf("%lld", &N);
printf("%d\n", (int)((long long)fib(N)*fib(N + 1) % MOD));
return 0;
}
bal4u