結果

問題 No.534 フィボナッチフィボナッチ数
ユーザー BantakoBantako
提出日時 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(){
      | ^~~~

ソースコード

diff #

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