結果

問題 No.658 テトラナッチ数列 Hard
コンテスト
ユーザー 0x19f
提出日時 2018-03-02 22:39:34
言語 C++14
(gcc 15.2.0 + boost 1.89.0)
コンパイル:
g++-15 -O2 -lm -std=c++14 -Wuninitialized -DONLINE_JUDGE -o a.out _filename_
実行:
./a.out
結果
AC  
実行時間 408 ms / 2,000 ms
コード長 962 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 1,653 ms
コンパイル使用メモリ 191,984 KB
実行使用メモリ 7,716 KB
最終ジャッジ日時 2026-03-06 09:51:25
合計ジャッジ時間 4,287 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 8
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include <bits/stdc++.h>
#define REP(i, a, n) for(ll i = ((ll) a); i < ((ll) n); i++)
using namespace std;
typedef long long ll;

typedef vector<ll> vec;
typedef vector<vec> mat;

mat mul(mat &A, mat &B) {
  mat C(A.size(), vector<ll>(B[0].size()));
  REP(i, 0, A.size())
    REP(j, 0, B[0].size())
      REP(k, 0, B.size())
        C[i][j] = (C[i][j] + A[i][k] * B[k][j]) % 17;
  return C;
}

mat pow(mat A, ll n) {
  mat B(A.size(), vec(A.size()));
  REP(i, 0, 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(void) {
  ll Q;
  cin >> Q;
  vector<ll> N(Q);
  REP(i, 0, Q) cin >> N[i];

  mat A(4, vec(4));
  A[0][0] = 1; A[0][1] = 1; A[0][2] = 1; A[0][3] = 1;
  A[1][0] = 1; A[1][1] = 0; A[1][2] = 0; A[1][3] = 0;
  A[2][0] = 0; A[2][1] = 1; A[2][2] = 0; A[2][3] = 0;
  A[3][0] = 0; A[3][1] = 0; A[3][2] = 1; A[3][3] = 0;

  REP(i, 0, Q) cout << pow(A, N[i] - 1)[3][0] << endl;
}
0