結果
問題 | No.534 フィボナッチフィボナッチ数 |
ユーザー |
![]() |
提出日時 | 2017-06-23 23:48:24 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 2 ms / 2,000 ms |
コード長 | 968 bytes |
コンパイル時間 | 1,544 ms |
コンパイル使用メモリ | 174,788 KB |
実行使用メモリ | 6,820 KB |
最終ジャッジ日時 | 2024-10-03 03:51:29 |
合計ジャッジ時間 | 2,720 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 42 |
ソースコード
#include <bits/stdc++.h>using namespace std;#define FOR(i,l,r) for(int i = (int) (l);i < (int) (r);i++)template<typename T> bool chmax(T& a,const T& b){ return a < b ? (a = b,true) : false; }template<typename T> bool chmin(T& a,const T& b){ return b < a ? (a = b,true) : false; }typedef long long ll;ll N;const ll MOD = 1e9 + 7,MOD2 = 2e9 + 16;void mat_mul(vector< vector<ll> >& a,vector< vector<ll> >& b,ll mod){vector< vector<ll> > res(a.size(),vector<ll>(b.front().size()));FOR(i,0,res.size()) FOR(j,0,res.front().size()){FOR(k,0,a.front().size()){(res [i] [j] += a [i] [k] * b [k] [j]) %= mod;}}a.swap(res);}ll fib(ll x,ll mod){vector< vector<ll> > p{{1,0}};vector< vector<ll> > q{{1,1},{1,0}};x--;while(x){if(x & 1) mat_mul(p,q,mod);mat_mul(q,q,mod);x >>= 1;}return p [0] [0];}int main(){scanf("%lld",&N);if(N == 0){puts("0");}else{printf("%lld\n",fib(fib(N,MOD2),MOD));}return 0;}