結果
問題 | No.534 フィボナッチフィボナッチ数 |
ユーザー |
|
提出日時 | 2022-08-26 00:45:32 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 2 ms / 2,000 ms |
コード長 | 2,125 bytes |
コンパイル時間 | 2,439 ms |
コンパイル使用メモリ | 204,696 KB |
最終ジャッジ日時 | 2025-01-31 04:01:04 |
ジャッジサーバーID (参考情報) |
judge4 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 42 |
ソースコード
#include<bits/stdc++.h> using namespace std; #define rep(i,n) for(ll i=0;i<n;i++) #define repl(i,l,r) for(ll i=(l);i<(r);i++) #define per(i,n) for(ll i=(n)-1;i>=0;i--) #define perl(i,r,l) for(ll i=r-1;i>=l;i--) #define fi first #define se second #define ins insert #define pqueue(x) priority_queue<x,vector<x>,greater<x>> #define all(x) (x).begin(),(x).end() #define CST(x) cout<<fixed<<setprecision(x) #define rev(x) reverse(x); using ll=long long; using vl=vector<ll>; using vvl=vector<vector<ll>>; using pl=pair<ll,ll>; using vpl=vector<pl>; using vvpl=vector<vpl>; const ll MOD=1000000007; const ll MOD9=998244353; const int inf=1e9+10; const ll INF=4e18; const ll dy[9]={1,0,-1,0,1,1,-1,-1,0}; const ll dx[9]={0,1,0,-1,1,-1,1,-1,0}; template <typename T> inline bool chmax(T &a, T b) { return ((a < b) ? (a = b, true) : (false)); } template <typename T> inline bool chmin(T &a, T b) { return ((a > b) ? (a = b, true) : (false)); } ll mod; using vec=vector<ll>; using matrix=vector<vector<ll>>; matrix matmul(matrix &left,matrix &right){ int n=left.size(); matrix ret(n,vec(n));//+の単位元で埋める for(int i=0;i<n;i++){ for(int k=0;k<n;k++){//ループ順序の最適化 for(int j=0;j<n;j++){ ret[i][j]=(ret[i][j]+left[i][k]*right[k][j])%mod; } } } return ret; } matrix matpow(matrix &mat,ll k){ int n=mat.size(); matrix ret(n,vec(n));//+の単位元で埋める for(int i=0;i<n;i++)ret[i][i]=1;//*の単位元で埋める while(k){ if(k&1)ret=matmul(mat,ret); k>>=1; if(k)mat=matmul(mat,mat); } return ret; } vec update(matrix &mat,vec &v){ int n=mat.size(); vec ret(n);//+の単位元で埋める for(int i=0;i<n;i++){ for(int j=0;j<n;j++){ ret[i]=(ret[i]+mat[i][j]*v[j])%mod; } } return ret; } int main(){ mod=(MOD+1)*2; matrix mat(2,vl(2,1)); mat[1][1]=0; ll n;cin >> n; mat=matpow(mat,n); ll f=mat[0][1]; //cout << mat[0][1] << endl; mat={{1,1},{1,0}}; mod=MOD; mat=matpow(mat,f); cout << mat[0][1] << endl; }