結果
問題 | 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;}