結果

問題 No.534 フィボナッチフィボナッチ数
ユーザー chocoruskchocorusk
提出日時 2018-09-22 12:14:53
言語 C++11
(gcc 11.4.0)
結果
AC  
実行時間 2 ms / 2,000 ms
コード長 1,389 bytes
コンパイル時間 535 ms
コンパイル使用メモリ 74,188 KB
実行使用メモリ 4,384 KB
最終ジャッジ日時 2023-09-25 11:00:30
合計ジャッジ時間 2,566 ms
ジャッジサーバーID
（参考情報） 		judge15 / judge12
このコードへのチャレンジ（β）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
4,376 KB
testcase_01 AC 2 ms
4,376 KB
testcase_02 AC 2 ms
4,376 KB
testcase_03 AC 2 ms
4,376 KB
testcase_04 AC 1 ms
4,376 KB
testcase_05 AC 2 ms
4,376 KB
testcase_06 AC 2 ms
4,376 KB
testcase_07 AC 2 ms
4,376 KB
testcase_08 AC 2 ms
4,376 KB
testcase_09 AC 1 ms
4,380 KB
testcase_10 AC 2 ms
4,380 KB
testcase_11 AC 1 ms
4,380 KB
testcase_12 AC 1 ms
4,376 KB
testcase_13 AC 2 ms
4,376 KB
testcase_14 AC 2 ms
4,376 KB
testcase_15 AC 1 ms
4,376 KB
testcase_16 AC 2 ms
4,380 KB
testcase_17 AC 2 ms
4,376 KB
testcase_18 AC 2 ms
4,376 KB
testcase_19 AC 2 ms
4,376 KB
testcase_20 AC 2 ms
4,376 KB
testcase_21 AC 2 ms
4,376 KB
testcase_22 AC 2 ms
4,376 KB
testcase_23 AC 2 ms
4,380 KB
testcase_24 AC 2 ms
4,384 KB
testcase_25 AC 2 ms
4,376 KB
testcase_26 AC 1 ms
4,380 KB
testcase_27 AC 2 ms
4,380 KB
testcase_28 AC 1 ms
4,376 KB
testcase_29 AC 2 ms
4,380 KB
testcase_30 AC 1 ms
4,376 KB
testcase_31 AC 2 ms
4,376 KB
testcase_32 AC 2 ms
4,380 KB
testcase_33 AC 1 ms
4,380 KB
testcase_34 AC 2 ms
4,380 KB
testcase_35 AC 2 ms
4,376 KB
testcase_36 AC 2 ms
4,380 KB
testcase_37 AC 2 ms
4,376 KB
testcase_38 AC 2 ms
4,376 KB
testcase_39 AC 2 ms
4,376 KB
testcase_40 AC 1 ms
4,376 KB
testcase_41 AC 2 ms
4,380 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <cstdio>
#include <cstring>
#include <iostream>
#include <string>
#include <cmath>
#include <bitset>
#include <vector>
#include <map>
#include <set>
#include <queue>
#include <deque>
#include <algorithm>
using namespace std;
typedef long long int ll;
typedef pair<int, int> P;
const ll MOD=1e9+7;

ll mulmod(ll a, ll b, ll m){
	ll x=1000000000;
	if(m<x) return a*b%m;
	ll m1=m/x, m2=m%x, a1=a/x, a2=a%x, b1=b/x, b2=b%x;
	ll q1=a1*b1/m1, r1=a1*b1%m1;
	ll c=r1*x+a1*b2+a2*b1-q1*m2, d=a2*b2;
	ll q2, r2;
	if(c<0) q2=-(-c+m1)/m1;
	else q2=c/m1;
	r2=c-q2*m1;
	ll ans=r2*x+d-q2*m2;
	if(ans<0) ans=(m-(-ans%m))%m;
	else ans%=m;
	return ans;
}

ll fib(ll k, ll m){
	ll a[4]={0, 1, 1, 1}, ans[4]={1, 0, 0, 1};
    while(k>0){
        if(k%2==1){
			ll b[4]={mulmod(ans[0], a[0], m)+mulmod(ans[1], a[2], m), mulmod(ans[0], a[1], m)+mulmod(ans[1], a[3], m), mulmod(ans[2], a[0], m)+mulmod(ans[3], a[2], m), mulmod(ans[2], a[1], m)+mulmod(ans[3], a[3], m)};
            for(int i=0; i<4; i++) ans[i]=b[i]%m;
        }
        ll b[4]={mulmod(a[0], a[0], m)+mulmod(a[1], a[2], m), mulmod(a[0], a[1], m)+mulmod(a[1], a[3], m), mulmod(a[2], a[0], m)+mulmod(a[3], a[2], m), mulmod(a[2], a[1], m)+mulmod(a[3], a[3], m)};
		for(int i=0; i<4; i++) a[i]=b[i]%m;
        k/=2;
    }
    return ans[1]%m;
}

int main()
{
	ll n;
	cin>>n;
	ll n1=fib(n, MOD*MOD-1);
	cout<<fib(n1, MOD)<<endl;
	return 0;
}
0