結果

問題 No.3635 Probability trip
コンテスト
ユーザー tau1235
提出日時 2026-08-21 22:14:59
言語 C++23
(gcc 15.2.0 + boost 1.90.0)
コンパイル:
g++-15 -O2 -lm -std=c++23 -Wuninitialized -DONLINE_JUDGE -o a.out _filename_
実行:
./a.out
結果
AC  
実行時間 67 ms / 2,000 ms
+ 411µs
コード長 1,499 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 2,264 ms
コンパイル使用メモリ 339,632 KB
実行使用メモリ 5,888 KB
最終ジャッジ日時 2026-08-21 22:15:07
合計ジャッジ時間 7,006 ms
ジャッジサーバーID
（参考情報） 		judge2_0 / judge1_0
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 43
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include<bits/stdc++.h>
#include<atcoder/modint>
using namespace std;

template<typename T,int N>
struct Matrix{
  array<array<T,N>,N> val;
  Matrix(){
    for (int i=0;i<N;i++) for (int j=0;j<N;j++) val[i][j]=T(0);
  }

  array<T,N> &operator[](int i){return val[i];}
	const array<T,N> &operator[](int i)const{return val[i];}
  Matrix operator*(const Matrix &b){
    Matrix a=*this,c;
    for (int i=0;i<N;i++){
      for (int j=0;j<N;j++){
        for (int k=0;k<N;k++){
          c[i][j]+=a[i][k]*b[k][j];
        }
      }
    }
    return c;
  }

  Matrix e(){
    Matrix ret;
    for (int i=0;i<N;i++) ret.val[i][i]=1;
    return ret;
  }

  Matrix pow(long long k){
    Matrix ret=e(),m=*this;
    while (k>0){
      if (k%2==1) ret=ret*m;
      m=m*m;
      k/=2;
    }
    return ret;
  }
};

int main(){
  using ll=long long;
  using mint=atcoder::modint998244353;
  using mat=Matrix<mint,52>;
  int n,m;
  cin>>n>>m;
  vector<int> d(n);
  vector<vector<int>> g(n);
  for (int i=0;i<m;i++){
    int u,v;
    cin>>u>>v;
    u--;v--;
    d[u]++;
    d[v]++;
    g[u].push_back(v);
    g[v].push_back(u);
  }
  mat one;
  for (int v=0;v<n;v++){
    for (int u:g[v]){
      one[v][u]=mint(d[v]).inv();
    }
  }
  auto f=[&](int u,int v,ll k){
    mat pw=one,ret=one.e();
    while (k){
      if (k&1) ret=ret*pw;
      pw=pw*pw;
      k/=2;
    }
    return ret[u][v];
  };
  ll s,t,a,b;
  cin>>s>>t>>a>>b;
  a--;b--;
  mint ans=(f(0,b,t-1)*f(b,a,s-t))/f(0,a,s-1);
  cout<<ans.val()<<endl;
}
0