結果
| 問題 |
No.658 テトラナッチ数列 Hard
|
| コンテスト | |
| ユーザー |
xuzijian629
|
| 提出日時 | 2018-08-30 01:20:29 |
| 言語 | C++11(廃止可能性あり) (gcc 13.3.0) |
| 結果 |
AC
|
| 実行時間 | 138 ms / 2,000 ms |
| コード長 | 1,500 bytes |
| コンパイル時間 | 1,045 ms |
| コンパイル使用メモリ | 92,088 KB |
| 実行使用メモリ | 6,944 KB |
| 最終ジャッジ日時 | 2024-09-13 19:44:15 |
| 合計ジャッジ時間 | 2,130 ms |
|
ジャッジサーバーID (参考情報) |
judge4 / judge5 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 8 |
ソースコード
#include <iostream>
#include <vector>
#include <array>
#include <algorithm>
#include <cmath>
#include <stack>
#include <queue>
#include <map>
#include <unordered_map>
#include <cassert>
#pragma GCC optimize("O3")
#pragma comment(linker, "STACK:36777216")
using namespace std;
using i64 = int64_t;
constexpr i64 mod = 17;
using vi = vector<i64>;
using vvi = vector<vi>;
using ii = pair<i64, i64>;
using vii = vector<ii>;
struct mat {
array<array<int, 4>, 4> a;
mat operator*(const mat m) {
mat ret;
ret.a = {{{0,0,0,0},{0,0,0,0},{0,0,0,0},{0,0,0,0}}};
for (int i = 0; i < 4; i++) {
for (int j = 0; j < 4; j++) {
for (int k = 0; k < 4; k++) {
ret.a[i][j] += a[i][k] * m.a[k][j];
ret.a[i][j] %= mod;
}
}
}
return ret;
}
};
mat matpow(mat m, i64 n) {
if (n == 0) {
mat ret;
ret.a = {{{1,0,0,0},{0,1,0,0},{0,0,1,0},{0,0,0,1}}};
return ret;
}
if (n % 2 == 0) {
mat tmp = matpow(m, n / 2);
return tmp * tmp;
}
return m * matpow(m, n - 1);
}
int main() {
int q;
cin >> q;
mat f;
f.a = {{{1,1,1,1},{1,0,0,0},{0,1,0,0},{0,0,1,0}}};
for (int i = 0; i < q; i++) {
i64 n;
cin >> n;
if (n < 4) {
cout << 0 << endl;
continue;
}
mat ret = matpow(f, n - 4);
cout << (mod + ret.a[0][0]) % mod << endl;
}
}
xuzijian629