結果

問題 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]);
}
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0