結果

問題 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;
}
0