結果

問題 No.526 フィボナッチ数列の第N項をMで割った余りを求める
コンテスト
ユーザー yurukumo
提出日時 2021-05-18 00:36:46
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 2 ms / 2,000 ms
コード長 1,922 bytes
コンパイル時間 4,385 ms
コンパイル使用メモリ 232,836 KB
実行使用メモリ 5,248 KB
最終ジャッジ日時 2024-10-07 06:26:13
合計ジャッジ時間 4,621 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 12
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
#include <atcoder/all>
#define int long long
#define endl "\n"
#define all(n) n.begin(),n.end()
#define rall(n) n.rbegin(),n.rend()
#define rep(i, s, n) for (int i = s; i < (int)(n); i++)
#define floatset(n) fixed<<setprecision(n)
#define rangeout(x,y,h,w) (x<0||y<0||h-1<x||w-1<y)
using namespace std;
using namespace atcoder;
template<class T>using vec=vector<T>;
template<class T>using min_queue=priority_queue<T,vector<T>,greater<T>>;
using ll = long long;
using ld = long double;
using pll = pair<ll,ll>;
using Graph = vector<vector<int>>;
struct Edge1{ll to,cost;};
using WGraph = vec<vec<Edge1>>;
struct BellEdge{ll from,to,cost;};
using BGraph=vec<BellEdge>;
constexpr ll INF = 1e18;
constexpr ll mod = 1000000007;
constexpr ll MOD = 998244353;
constexpr int dx[4]={1,0,-1,0};
constexpr int dy[4]={0,1,0,-1};
template<typename T>
bool chmin(T&a,T b){return a>b?a=b,true:false;}
template<typename T>
bool chmax(T&a,T b){return a<b?a=b,true:false;}
bool Flag(int bit,int i){return (bit&(1LL<<i));}
template<typename A, size_t N, typename T>
void Fill(A (&array)[N], const T &val){
    fill( (T*)array, (T*)(array+N), val );
}
vec<vec<int>>Mul(vec<vec<int>>&A,vec<vec<int>>&B,const int m){
    vec<vec<int>>C(A.size(),vec<int>(B.front().size()));
    for(int i=0;i<A.size();i++){
        for(int k=0;k<B.size();k++){
            for(int j=0;j<B.front().size();j++){
                (C[i][j]+=(A[i][k]*B[k][j]))%=m;
            }
        }
    }
    return C;
}
vec<vec<int>>Pow(vec<vec<int>>A,int x,const int m){
    vec<vec<int>>res(A.size(),vec<int>(A.size()));
    for(int i=0;i<A.size();i++)res[i][i]=1;
    while(x>0){
        if(x&1)res=Mul(res,A,m);
        A=Mul(A,A,m);
        x>>=1;
    }
    return res;
}
signed main(){
    int N,M;
    cin>>N>>M;
    vec<vec<int>>A(2,vec<int>(2));
    A[0][0]=1;A[0][1]=1;
    A[1][0]=1;A[1][1]=0;
    A=Pow(A,N-1,M);
    cout<<A[1][0]<<endl;
}
0