結果

問題 No.718 行列のできるフィボナッチ数列道場 (1)
ユーザー betit0919betit0919
提出日時 2019-08-18 15:20:40
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 2 ms / 2,000 ms
コード長 1,362 bytes
コンパイル時間 941 ms
コンパイル使用メモリ 102,520 KB
実行使用メモリ 5,248 KB
最終ジャッジ日時 2024-10-01 14:07:15
合計ジャッジ時間 1,797 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 20
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <vector>
#include <set>
#include <algorithm>
#include <math.h>
#include <string.h>
#include <map>
#include <numeric>
#include <random>
#include <queue>
#include <deque>
#include <tuple>
#include <iomanip>
#include <iterator>
#include <functional>

using namespace std;
typedef long long ll;

const int INF = (1 << 30) - 1;
const ll INFLL= (1LL << 61) - 1;
const int MOD = 1000000007;
#define ALL(a) (a).begin(),(a).end()
#define rALL(a) (a).rbegin(),(a).rend()
#define FOR(i,a,b) for(int i=(a);i<(b);++i)
#define REP(i,n) FOR(i,0,n)

struct matrix{
  ll space[2][2];

  matrix operator *(const matrix &another) const{
    matrix ans;
    for(int i = 0; i < 2; i++){
      for(int j = 0; j < 2; j++){
        ans.space[i][j] = 0;
        for(int k = 0; k < 2; k++){
          ans.space[i][j] += space[i][k] * another.space[k][j];
        }
        ans.space[i][j] %= MOD;
      }
    }
    return ans;
  }
};

const matrix E={{{1,0},{0,1}}};
const matrix A={{{1,1},{1,0}}};

matrix power(ll exponent){
  if(exponent % 2){
    return power(exponent - 1) * A;
  }else if(exponent){
    matrix root_ans = power(exponent / 2);
    return root_ans * root_ans;
  }else{
    return E;
  }
}


int main(){
  cin.tie(nullptr);
  ios::sync_with_stdio(false);
  ll n;
  cin>>n;
  cout<<power(n+1).space[1][1]*power(n+2).space[1][1]%MOD<<endl;
}
0