結果

問題 No.526 フィボナッチ数列の第N項をMで割った余りを求める
コンテスト
ユーザー Manuel1024
提出日時 2020-11-28 09:40:40
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 2 ms / 2,000 ms
コード長 1,342 bytes
コンパイル時間 876 ms
コンパイル使用メモリ 74,912 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-09-12 22:11:11
合計ジャッジ時間 1,678 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 12
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <vector>
using namespace std;
using ll = long long int;
ll MOD;
vector<ll> times(vector<ll> &memo, vector<vector<ll>> &matrix){
    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<vector<ll>> update(vector<vector<ll>> &matrix){
    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;
}

void matrixpow(vector<ll> &memo, vector<vector<ll>> &matrix, int n){
    while(n){
        if(n & 1){
            memo = times(memo, matrix);        
        }
        matrix = update(matrix);
        n >>= 1;
    }
}

int main(){
    int n;
    cin >> n >> MOD;
    //MOD = 10000000;
    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;
    matrixpow(memo, matrix, n-1);
    cout << memo[1] << endl;
    return 0;
}
0