結果
| 問題 |
No.718 行列のできるフィボナッチ数列道場 (1)
|
| コンテスト | |
| ユーザー |
どらら
|
| 提出日時 | 2018-07-27 22:55:22 |
| 言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
AC
|
| 実行時間 | 2 ms / 2,000 ms |
| コード長 | 896 bytes |
| コンパイル時間 | 1,805 ms |
| コンパイル使用メモリ | 174,420 KB |
| 実行使用メモリ | 5,376 KB |
| 最終ジャッジ日時 | 2024-07-05 03:20:33 |
| 合計ジャッジ時間 | 2,416 ms |
|
ジャッジサーバーID (参考情報) |
judge1 / judge3 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 20 |
ソースコード
#include <bits/stdc++.h>
using namespace std;
#define REP(i,a,n) for(int i=(a); i<(int)(n); i++)
#define rep(i,n) REP(i,0,n)
#define FOR(it,c) for(__typeof((c).begin()) it=(c).begin(); it!=(c).end(); ++it)
#define ALLOF(c) (c).begin(), (c).end()
typedef long long ll;
typedef unsigned long long ull;
typedef vector<ll> vec;
typedef vector<vec> mat;
const ll MOD = 1000000007;
mat mul(mat &A, mat &B){ // A*B
mat C(A.size(), vec(B[0].size()));
rep(i,A.size()) rep(k,B.size()){
rep(j,B[0].size()){
C[i][j] = (C[i][j]+A[i][k]*B[k][j])%MOD;
}
}
return C;
}
mat pow(mat A, ll n){ // A^n
mat B(A.size(), vec(A.size()));
rep(i,A.size()) B[i][i] = 1;
while(n>0){
if(n&1) B=mul(B,A);
A=mul(A,A);
n>>=1;
}
return B;
}
int main(){
ll N;
cin >> N;
mat a = {{1, 1},{1, 0}};
mat b = pow(a, N);
cout << (b[0][0] * b[0][1])%MOD << endl;
return 0;
}
どらら