結果
問題 | No.1704 Many Bus Stops (easy) |
ユーザー |
|
提出日時 | 2021-10-08 23:30:40 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 123 ms / 2,000 ms |
コード長 | 2,125 bytes |
コンパイル時間 | 2,294 ms |
コンパイル使用メモリ | 198,152 KB |
最終ジャッジ日時 | 2025-01-24 23:27:47 |
ジャッジサーバーID (参考情報) |
judge5 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 1 |
other | AC * 41 |
ソースコード
#include<bits/stdc++.h>using namespace std;#define rep(i,n) for(ll i=0;i<n;i++)#define repl(i,l,r) for(ll i=(l);i<(r);i++)#define per(i,n) for(ll i=(n)-1;i>=0;i--)#define perl(i,r,l) for(ll i=r-1;i>=l;i--)#define fi first#define se second#define pb push_back#define ins insert#define pqueue(x) priority_queue<x,vector<x>,greater<x>>#define all(x) (x).begin(),(x).end()#define CST(x) cout<<fixed<<setprecision(x)#define vtpl(x,y,z) vector<tuple<x,y,z>>#define rev(x) reverse(x);using ll=long long;using vl=vector<ll>;using vvl=vector<vector<ll>>;using pl=pair<ll,ll>;using vpl=vector<pl>;using vvpl=vector<vpl>;const ll MOD=1000000007;const ll MOD9=998244353;const int inf=1e9+10;const ll INF=4e18;const ll dy[9]={0,1,-1,0,1,1,-1,-1,0};const ll dx[9]={1,0,0,-1,1,-1,1,-1,0};template<class T> inline bool chmin(T& a, T b) {if (a > b) {a = b;return true;}return false;}template<class T> inline bool chmax(T& a, T b) {if (a < b) {a = b;return true;}return false;}ll mod=1000000007;//一次元ベクトルの高速化dpvl matmul(vl &dp, vvl &mt){ll m=dp.size();vl ret(m,0);rep(i,m)rep(j,m){ret[i]+=mt[i][j]*dp[j];ret[i]%=mod;}return ret;}vvl update(vvl &mt){ll m=mt.size();vvl ret(m,vl(m,0));rep(i,m)rep(j,m)rep(k,m){ret[i][j]+=mt[i][k]*mt[k][j];ret[i][j]%=mod;}return ret;}void matpow(vl &dp, vvl &mt, ll k){ll m=dp.size();while(k){if(k&1)dp=matmul(dp,mt);mt=update(mt);k/=2;}}ll modpow(ll a,ll n, ll mod) {ll res = 1;while (n > 0) {if (n & 1) res = res * a % mod;a = a * a % mod;n >>= 1;}return res;}void solve(){ll n;cin >> n;ll inv=modpow(3,mod-2,mod);vvl mat(4,vl(4));rep(i,2)mat[i+2][i]=1;mat[0][0]=inv;mat[0][3]=inv;mat[1][1]=inv;mat[1][2]=inv*2%mod;mat[1][3]=inv;vl dp(4);dp[0]=inv;dp[2]=1;matpow(dp,mat,n);cout << dp[2] << endl;}int main(){ll t;cin >> t;while(t--)solve();}