結果

問題 No.534 フィボナッチフィボナッチ数
ユーザー gigime
提出日時 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
権限があれば一括ダウンロードができます

ソースコード

diff #
プレゼンテーションモードにする

#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;
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0