結果

問題 No.1750 ラムドスウイルスの感染拡大-hard
ユーザー inksamurai
提出日時 2021-11-19 21:45:00
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 271 ms / 2,000 ms
コード長 2,636 bytes
コンパイル時間 2,109 ms
コンパイル使用メモリ 184,252 KB
実行使用メモリ 6,820 KB
最終ジャッジ日時 2024-12-31 22:37:37
合計ジャッジ時間 5,877 ms
ジャッジサーバーID
(参考情報)
judge3 / judge1
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 4
other AC * 30
権限があれば一括ダウンロードができます

ソースコード

diff #
プレゼンテーションモードにする

#include <bits/stdc++.h>
#define fi first
#define se second
#define pb push_back
#define sz(a) (int)a.size()
#define all(a) a.begin(),a.end()
#define rep(i,n) for(int i=0;i<n;i++)
#define crep(i,x,n) for(int i=x;i<n;i++)
#define drep(i,n) for(int i=n-1;i>=0;i--)
#define vec(...) vector<__VA_ARGS__>
#define _3HaBFkZ ios_base::sync_with_stdio(false),cin.tie(0),cout.tie(0)
using namespace std;
typedef long long ll;
typedef long double ld;
using pii=pair<int,int>;
using vi=vector<int>;
//snuke's modular int
template <ll mod>
struct modularint{
ll x;
modularint(ll x=0):x(x%mod){}
modularint& operator+=(const modularint a){
if ((x += a.x) >= mod) x -= mod;
return *this;
}
modularint& operator-=(const modularint a){
if ((x += mod-a.x) >= mod) x -= mod;
return *this;
}
modularint& operator*=(const modularint a){
(x *= a.x) %= mod;
return *this;
}
modularint operator+(const modularint a)const{
modularint res(*this);
return res+=a;
}
modularint operator-(const modularint a)const{
modularint res(*this);
return res-=a;
}
modularint operator*(const modularint a)const{
modularint res(*this);
return res*=a;
}
modularint pow(ll n)const{
modularint res=1,x(*this);
while(n){
if(n&1)res*=x;
x*=x;
n>>=1;
}
return res;
}
modularint inv()const{
return pow(mod-2);
}
};
using mint=modularint<998244353>;
template <typename T>
struct matrix{
vec(vec(T)) a;
//initilize matrix here
matrix(){}
matrix(int h,int w){
a.clear();
a.resize(h,vec(T)(w));
}
matrix(vec(vec(T)) nea){
a=nea;
}
//outer vector size here
int size()const{ return a.size(); }
const vector<T>& operator[](int i)const{ return a[i]; }
vector<T>& operator[](int i){ return a[i]; }
matrix<T>& operator *=(const matrix<T>& rhs){
int h=a.size(), w=rhs[0].size(), c=rhs.size();
matrix<T> res(h,w);
rep(i,h){
rep(j,w){
rep(k,c){
res[i][j] += a[i][k] * rhs[k][j];
}
}
}
this->a = res.a;
return *this;
}
matrix<T> operator *(const matrix<T>& rhs){
return (matrix<T>(*this) *= rhs);
}
};
template <typename T>
matrix<T> idenmat(int n){
matrix<T> res(n,n);
rep(i,n) res.a[i][i] = 1;
return res;
};
template <typename T>
matrix<T> pow(matrix<T> mat,ll n){
matrix<T> res=idenmat<T>(mat.size());
while(n){
if(n&1) res*=mat;
mat*=mat;
n>>=1;
}
return res;
};
int main(){
_3HaBFkZ;
ll n,m,t;
cin>>n>>m>>t;
vec(vec(mint)) wys(n,vec(mint)(n,0));
rep(i,m){
int u,v;
cin>>u>>v;
wys[u][v]=wys[u][v]+1;
wys[v][u]=wys[v][u]+1;
}
matrix<mint> mat(wys);
matrix<mint> now=idenmat<mint>(n);
now=now*pow(mat,t);
cout<<now[0][0].x<<"\n";
//
return 0;
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0