結果
問題 | No.1750 ラムドスウイルスの感染拡大-hard |
ユーザー |
|
提出日時 | 2023-05-13 13:32:12 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 281 ms / 2,000 ms |
コード長 | 1,021 bytes |
コンパイル時間 | 1,830 ms |
コンパイル使用メモリ | 174,656 KB |
実行使用メモリ | 8,832 KB |
最終ジャッジ日時 | 2024-11-29 06:38:31 |
合計ジャッジ時間 | 6,718 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 30 |
ソースコード
#include<bits/stdc++.h> using namespace std; constexpr int mod=998244353; int N,M; vector<vector<vector<long long>>> dp(102,vector<vector<long long>>(102,vector<long long>(61,0ll))); long long T; int main(){ scanf("%d%d%lld",&N,&M,&T); for(int i=0;i<M;i++){ int u,v; scanf("%d%d",&u,&v); dp[u][v][0]=dp[v][u][0]=1ll; } for(int bit=1;bit<=60;bit++){ //2^bit回の移動 for(int k=0;k<N;k++){ //中継地点k for(int i=0;i<N;i++)for(int j=0;j<N;j++){ dp[i][j][bit]+=(dp[i][k][bit-1]*dp[k][j][bit-1])%mod; dp[i][j][bit]%=mod; } } } vector<long long> ans(N,0ll);//ある時点でiにいる ans[0]=1ll; for(int bit=0;bit<=60;bit++)if(T&1ll<<bit){ vector<long long> ans_(N,0ll); for(int i=0;i<N;i++)for(int j=0;j<N;j++){ ans_[j]+=(ans[i]*dp[i][j][bit])%mod; ans_[j]%=mod; } ans=ans_; } printf("%lld",ans.front()); }