結果
| 問題 |
No.534 フィボナッチフィボナッチ数
|
| コンテスト | |
| ユーザー |
Bantako
|
| 提出日時 | 2018-06-19 00:53:29 |
| 言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
AC
|
| 実行時間 | 2 ms / 2,000 ms |
| コード長 | 1,383 bytes |
| コンパイル時間 | 1,690 ms |
| コンパイル使用メモリ | 174,980 KB |
| 実行使用メモリ | 5,376 KB |
| 最終ジャッジ日時 | 2024-06-30 17:09:07 |
| 合計ジャッジ時間 | 3,035 ms |
|
ジャッジサーバーID (参考情報) |
judge1 / judge3 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| other | AC * 42 |
コンパイルメッセージ
main.cpp:45:1: warning: ISO C++ forbids declaration of 'main' with no type [-Wreturn-type]
45 | main(){
| ^~~~
ソースコード
#include<bits/stdc++.h>
#define rep(i,a,b) for(int i=int(a);i<int(b);++i)
using namespace std;
typedef long long ll;
typedef vector<vector<ll> > MAT;
int INF = (1LL << 30) - 1;
ll MOD = 1e9+7;
MAT mulmat(MAT a,MAT b,ll m){
//ll mat[2][2] = {};
auto mat = MAT (2, vector<ll>(2));
mat[0][0] = ((a[0][0] * b[0][0]) % m + (a[0][1] * b[1][0]) % m) % m;
mat[0][1] = ((a[0][0] * b[0][1]) % m + (a[0][1] * b[1][1]) % m) % m;
mat[1][0] = ((a[1][0] * b[0][0]) % m + (a[1][1] * b[1][0]) % m) % m;
mat[1][1] = ((a[1][0] * b[0][1]) % m + (a[1][1] * b[1][1]) % m) % m;
return mat;
}
ll fib(ll N,ll m){
/*
ll a[2][2] = { {1 , 1},
{1 , 0} };
ll b[2][2] = { {1 , 0},
0 , 1} };
*/
auto a = MAT(2, vector<ll>(2));
auto b = MAT(2, vector<ll>(2));
a[0][0] = 1; a[0][1] = 1;
a[1][0] = 1; a[1][1] = 0;
b[0][0] = 1; b[0][1] = 0;
b[1][0] = 0; b[1][1] = 1;
ll ans = 0;
while(N > 0){
if(N & 1){
b = mulmat(a,b,m);
}
a = mulmat(a,a,m);
N /= 2;
}
//cout << b[0][0] << " " << b[0][1] << endl << b[1][0] << " " << b[1][1] << endl;
return b[1][0];
}
main(){
ll loop = 2000000016;
ll N;
cin >> N;
//cout << fib(N,MOD) << endl;
//cout << fib(N % loop,MOD) << endl;
cout << fib(fib(N,loop),MOD) << endl;
}
Bantako