結果
問題 | No.534 フィボナッチフィボナッチ数 |
ユーザー |
![]() |
提出日時 | 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;}