結果

問題 No.534 フィボナッチフィボナッチ数
ユーザー Manuel1024Manuel1024
提出日時 2020-11-29 05:54:20
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 2 ms / 2,000 ms
コード長 1,588 bytes
コンパイル時間 1,084 ms
コンパイル使用メモリ 73,492 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-09-13 01:52:42
合計ジャッジ時間 2,139 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

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

ソースコード

diff # 

#include <iostream>
#include <vector>
using namespace std;
using ll = long long int;
const ll MOD = 1000'000'007;

vector<vector<ll>> update(vector<vector<ll>> &matrix, ll mod){
    vector<vector<ll>> ans(matrix.size(), vector<ll>(matrix.size(), 0));
    for(int i = 0; i < matrix.size(); i++){
        for(int j = 0; j < matrix.size(); j++){
            for(int k = 0; k < matrix.size(); k++){
                ans[i][j] += matrix[i][k]*matrix[k][j];
                ans[i][j] %= mod;
            }
        }
    }
    return ans;
}

vector<ll> times(vector<ll> &memo, vector<vector<ll>> &matrix, ll mod){
    vector<ll> ans(memo.size(), 0);
    for(int i = 0; i < memo.size(); i++){
        for(int j = 0; j < memo.size(); j++){
            ans[i] += memo[j]*matrix[i][j];
            ans[i] %= mod;
        }
    }
    return ans;
}

vector<ll> matpow(vector<ll> &memo, vector<vector<ll>> &matrix, ll n, ll mod){
    vector<ll> ans = memo;
    while(n){
        if(n & 1){
            ans = times(ans, matrix, mod);
        }
        matrix = update(matrix, mod);
        n >>= 1;
    }
    return ans;
}

int main(){
    ll n;
    cin >> n;
    vector<vector<ll>> matrix(2, vector<ll>(2));
    vector<ll> memo(2);
    matrix[0][0] = matrix[0][1] = matrix[1][0] = 1;
    matrix[1][1] = 0;
    memo[0] = 1;
    memo[1] = 0;

    memo = matpow(memo, matrix, n, (MOD+1)*2);
    ll fn = memo[1];
    matrix[0][0] = matrix[0][1] = matrix[1][0] = 1;
    matrix[1][1] = 0;
    memo[0] = 1;
    memo[1] = 0;
    memo = matpow(memo, matrix, fn, MOD);
    cout << memo[1] << endl;
    return 0;
}
0