結果

問題 No.658 テトラナッチ数列 Hard
コンテスト
ユーザー betit0919
提出日時 2019-08-10 15:43:49
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 55 ms / 2,000 ms
コード長 1,525 bytes
コンパイル時間 952 ms
コンパイル使用メモリ 103,472 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-07-19 17:08:11
合計ジャッジ時間 1,740 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 8
権限があれば一括ダウンロードができます

ソースコード

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{
  int space[4][4];

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

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

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(0);
  ios::sync_with_stdio(false);
  int Q;
  cin>>Q;
  matrix mat={{{1, 1, 1, 1}, {1, 0, 0, 0}, {0, 1, 0, 0}, {0, 0, 1, 0}}};

  REP(i,Q){
    ll n;cin>>n;
    cout<<power(n).space[3][3]<<endl;
  }
}
0