結果
問題 |
No.658 テトラナッチ数列 Hard
|
ユーザー |
![]() |
提出日時 | 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; } }