結果

問題 No.658 テトラナッチ数列 Hard
ユーザー hiro71687khiro71687k
提出日時 2023-03-08 12:03:53
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 931 ms / 2,000 ms
コード長 1,451 bytes
コンパイル時間 4,387 ms
コンパイル使用メモリ 267,092 KB
実行使用メモリ 4,348 KB
最終ジャッジ日時 2023-10-18 05:53:55
合計ジャッジ時間 9,446 ms
ジャッジサーバーID
(参考情報)
judge14 / judge15
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
4,348 KB
testcase_01 AC 2 ms
4,348 KB
testcase_02 AC 2 ms
4,348 KB
testcase_03 AC 5 ms
4,348 KB
testcase_04 AC 317 ms
4,348 KB
testcase_05 AC 369 ms
4,348 KB
testcase_06 AC 468 ms
4,348 KB
testcase_07 AC 510 ms
4,348 KB
testcase_08 AC 604 ms
4,348 KB
testcase_09 AC 931 ms
4,348 KB
testcase_10 AC 918 ms
4,348 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

#include <bits/stdc++.h>
#include <atcoder/all>
using namespace std;
using namespace atcoder;
using ll=long long;
using ld=long double;
ld pie=3.141592653589793;
ll inf=10010010010010010;
ll mod=17;
vector<vector<ll>>mul(vector<vector<ll>>&a,vector<vector<ll>>&b){
    vector<vector<ll>>c(a.size(),vector<ll>(b[0].size(),0));
    for (ll i = 0; i < a.size(); i++)
    {
        for (ll j = 0; j < b[0].size(); j++)
        {
            for (ll k = 0; k < a[0].size(); k++)
            {
                c[i][j]+=a[i][k]*b[k][j];
                c[i][j]%=mod;
            }
        }
    }
    return c;
}
vector<vector<ll>>matpow(vector<vector<ll>>a,ll n){
    vector<vector<ll>>b(a.size(),vector<ll>(a.size(),0));
    for (ll i = 0; i < a.size(); i++)
    {
        b[i][i]=1;
    }
    while (n>0)
    {
        if (n%2)
        {
            b=mul(a,b);
            n-=1;
        }else{
            a=mul(a,a);
            n/=2;
        }
    }
    return b;
}
int main(){
    ll q;
    cin >> q;
    for (ll i = 0; i < q; i++)
    {
        ll n;
        cin >> n;
        vector<vector<ll>>a(4,vector<ll>(4,0));
        a[0][0]=1,a[0][1]=1,a[0][2]=1,a[0][3]=1;
        a[1][0]=1,a[2][1]=1,a[3][2]=1;
        if (n<=3)
        {
            cout << 0 << endl;
        }else if (n==4)
        {
            cout << 1 << endl;
        }else{
            vector<vector<ll>>d=matpow(a,n-4);
            cout << d[0][0] << endl;
        }
    }
    
}
0