結果
問題 | No.718 行列のできるフィボナッチ数列道場 (1) |
ユーザー |
![]() |
提出日時 | 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; }