結果

問題 No.1704 Many Bus Stops (easy)
コンテスト
ユーザー eQe
提出日時 2021-10-13 03:24:31
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 99 ms / 2,000 ms
コード長 2,792 bytes
コンパイル時間 3,978 ms
コンパイル使用メモリ 234,432 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-09-17 15:46:51
合計ジャッジ時間 7,380 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge6
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 1
other AC * 41
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <atcoder/all>
using namespace atcoder;
#include <bits/stdc++.h>
#define itr(x,c) for(auto x=c.begin();x!=c.end();x++)
#define ritr(x,c) for(auto x=c.rbegin();x!=c.rend();x++)
#define all(a) (a.begin()),(a.end())
#define rall(a) (a.rbegin()),(a.rend())
using namespace std;
using ll=long long;
using pll=pair<ll,ll>;
template<class T>bool chmax(T &a, const T &b) {if(a<b) {a=b; return 1;} return 0;}
template<class T>bool chmin(T &a, const T &b) {if(b<a) {a=b; return 1;} return 0;}
template<class T>istream& operator >> (istream& is, vector<T>& v) { for(T& x:v) is >> x; return is; }
template<class T>ostream& operator << (ostream& os, const vector<T>& v) {for(int i=0;i<v.size();i++) {os << v[i] << (i+1 == v.size() ? "":" ");} return os;}
template<class T1, class T2>istream& operator >> (istream& is, pair<T1,T2>& p) { is >> p.first >> p.second; return is; }
//istream& operator >> (istream& is, modint1000000007& x) { unsigned int t; is >> t; x=t; return is; }
//istream& operator >> (istream& is, modint998244353& x) { unsigned int t; is >> t; x=t; return is; }
ostream& operator << (ostream& os, const modint1000000007& x) { os << x.val(); return os; }
ostream& operator << (ostream& os, const modint998244353& x) { os << x.val(); return os; }
template<class... A> void pt()      { std::cout << "\n"; }
template<class... A> void pt_rest() { std::cout << "\n"; }
template<class T, class... A> void pt_rest(const T& first, const A&... rest) { std::cout << " " << first; pt_rest(rest...); }
template<class T, class... A> void pt(const T& first, const A&... rest)      { std::cout << first; pt_rest(rest...); }

static const ll INF=1e18+7;
static const ll MAX=101010;
static const ll MOD=1e9+7;
static const double PI=acos(-1);
static const double EPS=1e-10;
//static const ll MOD=998244353;

using v1d=vector<modint1000000007>;
using v2d=vector<v1d>;
using v3d=vector<v2d>;
using v4d=vector<v3d>;

typedef v2d MAT;

MAT mul(MAT &A,MAT &B){
  ll n=A.size();
  ll m=B[0].size();
  ll l=B.size();
  
  MAT C(n,v1d(m));
  ll i,j,k;
  for(i=0;i<n;i++)for(j=0;j<m;j++)for(k=0;k<l;k++)C[i][j]+=A[i][k]*B[k][j];
  return C;
}

MAT powmat(MAT X,ll n){
  ll R=X.size();
  MAT res(R,v1d(R));
  for(ll i=0;i<R;i++)res[i][i]=1;
  while(n) {
    if(n&1) res=mul(res,X);
    X=mul(X,X);
    n>>=1;
  }
  return res;
}

int main(void){
  cin.tie(nullptr);
  ios::sync_with_stdio(false);
  ll i,j,k;
  
  ll T;cin>>T;
  while(T--){
    ll N;cin>>N;
    MAT A(6,v1d(6,0));
    modint1000000007 three=3;
    
    for(i=0;i<3;i++)A[i][i]=three.inv();
    for(i=0;i<3;i++)for(j=3;j<6;j++)if(j-3!=i)A[i][j]=three.inv();
    for(i=3;i<6;i++)for(j=0;j<3;j++)if(j==i-3)A[i][j]=1;
    
    A=powmat(A,N);
    MAT b(6,v1d(1,0));
    b[0][0]=three.inv();
    b[3][0]=1;
    pt(mul(A,b)[3][0]);
    
  }
  
}
0