結果

問題 No.718 行列のできるフィボナッチ数列道場 (1)
ユーザー betit0919betit0919
提出日時 2019-08-18 15:20:40
言語 C++14
(gcc 12.3.0 + boost 1.83.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
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
5,248 KB
testcase_01 AC 2 ms
5,248 KB
testcase_02 AC 1 ms
5,248 KB
testcase_03 AC 2 ms
5,248 KB
testcase_04 AC 1 ms
5,248 KB
testcase_05 AC 2 ms
5,248 KB
testcase_06 AC 1 ms
5,248 KB
testcase_07 AC 1 ms
5,248 KB
testcase_08 AC 2 ms
5,248 KB
testcase_09 AC 2 ms
5,248 KB
testcase_10 AC 1 ms
5,248 KB
testcase_11 AC 2 ms
5,248 KB
testcase_12 AC 1 ms
5,248 KB
testcase_13 AC 2 ms
5,248 KB
testcase_14 AC 1 ms
5,248 KB
testcase_15 AC 2 ms
5,248 KB
testcase_16 AC 1 ms
5,248 KB
testcase_17 AC 2 ms
5,248 KB
testcase_18 AC 2 ms
5,248 KB
testcase_19 AC 1 ms
5,248 KB
testcase_20 AC 1 ms
5,248 KB
testcase_21 AC 2 ms
5,248 KB
testcase_22 AC 2 ms
5,248 KB
権限があれば一括ダウンロードができます

ソースコード

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