結果
問題 | No.534 フィボナッチフィボナッチ数 |
ユーザー |
![]() |
提出日時 | 2019-04-19 17:31:54 |
言語 | C (gcc 13.3.0) |
結果 |
AC
|
実行時間 | 2 ms / 2,000 ms |
コード長 | 1,121 bytes |
コンパイル時間 | 207 ms |
コンパイル使用メモリ | 30,592 KB |
実行使用メモリ | 6,948 KB |
最終ジャッジ日時 | 2024-09-22 11:37:13 |
合計ジャッジ時間 | 1,310 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 42 |
ソースコード
// yukicoder: No.534 フィボナッチフィボナッチ数// 2019.4.19 bal4u#include <stdio.h>#define MOD0 1000000007int MOD;//// 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] % MOD + f01 * M[1][0] % MOD);if (x >= MOD) x -= MOD;y = (f00 * M[0][1] % MOD + f01 * M[1][1] % MOD);if (y >= MOD) y -= MOD;z = (f10 * M[0][0] % MOD + f11 * M[1][0] % MOD);if (z >= MOD) z -= MOD;w = (f10 * M[0][1] % MOD + f11 * M[1][1] % MOD);if (w >= MOD) w -= 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(){int ans;long long n;scanf("%lld", &n);MOD = (MOD0+1) << 1;ans = fib(n);MOD = MOD0;printf("%d\n", fib(ans));return 0;}