結果
問題 | No.534 フィボナッチフィボナッチ数 |
ユーザー |
![]() |
提出日時 | 2018-06-19 00:53:29 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 2 ms / 2,000 ms |
コード長 | 1,383 bytes |
コンパイル時間 | 1,690 ms |
コンパイル使用メモリ | 174,980 KB |
実行使用メモリ | 5,376 KB |
最終ジャッジ日時 | 2024-06-30 17:09:07 |
合計ジャッジ時間 | 3,035 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 42 |
コンパイルメッセージ
main.cpp:45:1: warning: ISO C++ forbids declaration of 'main' with no type [-Wreturn-type] 45 | main(){ | ^~~~
ソースコード
#include<bits/stdc++.h> #define rep(i,a,b) for(int i=int(a);i<int(b);++i) using namespace std; typedef long long ll; typedef vector<vector<ll> > MAT; int INF = (1LL << 30) - 1; ll MOD = 1e9+7; MAT mulmat(MAT a,MAT b,ll m){ //ll mat[2][2] = {}; auto mat = MAT (2, vector<ll>(2)); mat[0][0] = ((a[0][0] * b[0][0]) % m + (a[0][1] * b[1][0]) % m) % m; mat[0][1] = ((a[0][0] * b[0][1]) % m + (a[0][1] * b[1][1]) % m) % m; mat[1][0] = ((a[1][0] * b[0][0]) % m + (a[1][1] * b[1][0]) % m) % m; mat[1][1] = ((a[1][0] * b[0][1]) % m + (a[1][1] * b[1][1]) % m) % m; return mat; } ll fib(ll N,ll m){ /* ll a[2][2] = { {1 , 1}, {1 , 0} }; ll b[2][2] = { {1 , 0}, 0 , 1} }; */ auto a = MAT(2, vector<ll>(2)); auto b = MAT(2, vector<ll>(2)); a[0][0] = 1; a[0][1] = 1; a[1][0] = 1; a[1][1] = 0; b[0][0] = 1; b[0][1] = 0; b[1][0] = 0; b[1][1] = 1; ll ans = 0; while(N > 0){ if(N & 1){ b = mulmat(a,b,m); } a = mulmat(a,a,m); N /= 2; } //cout << b[0][0] << " " << b[0][1] << endl << b[1][0] << " " << b[1][1] << endl; return b[1][0]; } main(){ ll loop = 2000000016; ll N; cin >> N; //cout << fib(N,MOD) << endl; //cout << fib(N % loop,MOD) << endl; cout << fib(fib(N,loop),MOD) << endl; }