結果

問題 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;
}
0