結果

問題 No.1073 無限すごろく
コンテスト
ユーザー hs0319boy
提出日時 2020-06-05 23:52:03
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 2 ms / 2,000 ms
コード長 1,268 bytes
コンパイル時間 1,741 ms
コンパイル使用メモリ 174,292 KB
実行使用メモリ 6,824 KB
最終ジャッジ日時 2024-12-17 19:14:48
合計ジャッジ時間 2,684 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 30
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
using ll=long long;
using vin=vector<int>;
using vll=vector<long long>;
using vvin=vector<vector<int>>;
using vvll=vector<vector<long long>>;
using vstr=vector<string>;
using vvstr=vector<vector<string>>;
using vch=vector<char>;
using vvch=vector<vector<char>>;
using vbo=vector<bool>;
using vvbo=vector<vector<bool>>;
using vpii=vector<pair<int,int>>;
using pqsin=priority_queue<int,vector<int>,greater<int>>;
#define mp make_pair
#define rep(i,n) for(int i=0;i<(int)(n);i++)
#define rep2(i,s,n) for(int i=(s);i<(int)(n);i++)
#define all(v) v.begin(),v.end()
#define decp(n) cout<<fixed<<setprecision((int)n)
const ll inf=1e9+7;
const ll INF=1e18;

vvll matrixpro(vvll a,vvll b){
    ll s=a.size();
    vvll c(s,vll(s));
    rep(i,s)rep(j,s)rep(k,s){
        c[i][j]+=(a[i][k]*b[k][j])%inf;
        c[i][j]%=inf;
    }
    return c;
}

vvll matrixexp(vvll a,ll n){
    ll s=a.size();
    vvll ans(s,vll(s));
    rep(i,s)ans[i][i]=1;
    while(n){
        if(n&1)ans=matrixpro(ans,a);
        a=matrixpro(a,a);
        n>>=1;
    }
    return ans;
}

int main(){
    ll n;cin>>n;
    const ll tmp=166666668;
    vvll a(6,vll(6));
    rep(i,6)a[0][i]=tmp;
    rep2(i,1,6)a[i][i-1]=1;
    cout<<matrixexp(a,n)[0][0]<<endl;
}
0