#717777 (C++14) No.1750 ラムドスウイルスの感染拡大-hard

• 問題: No.1750 ラムドスウイルスの感染拡大-hard
• ユーザー: inksamurai
• 提出日時: 2021-11-19 21:45:00
• 言語: C++14 (gcc 12.3.0 + boost 1.83.0)
• 結果: AC
• 実行時間: 230 ms / 2,000 ms
• コード長: 2,636 bytes
• コンパイル時間: 1,664 ms
• コンパイル使用メモリ: 182,112 KB
• 実行使用メモリ: 5,376 KB
• 最終ジャッジ日時: 2024-06-10 08:25:39
• 合計ジャッジ時間: 4,911 ms
• ジャッジサーバーID （参考情報）: judge5 / judge1

入力	結果	実行時間; 実行使用メモリ
testcase_00	AC	2 ms; 5,248 KB
testcase_01	AC	1 ms; 5,376 KB
testcase_02	AC	2 ms; 5,376 KB
testcase_03	AC	1 ms; 5,376 KB
testcase_04	AC	7 ms; 5,376 KB
testcase_05	AC	1 ms; 5,376 KB
testcase_06	AC	2 ms; 5,376 KB
testcase_07	AC	2 ms; 5,376 KB
testcase_08	AC	44 ms; 5,376 KB
testcase_09	AC	45 ms; 5,376 KB
testcase_10	AC	43 ms; 5,376 KB
testcase_11	AC	34 ms; 5,376 KB
testcase_12	AC	46 ms; 5,376 KB
testcase_13	AC	44 ms; 5,376 KB
testcase_14	AC	215 ms; 5,376 KB
testcase_15	AC	220 ms; 5,376 KB
testcase_16	AC	223 ms; 5,376 KB
testcase_17	AC	226 ms; 5,376 KB
testcase_18	AC	230 ms; 5,376 KB
testcase_19	AC	222 ms; 5,376 KB
testcase_20	AC	152 ms; 5,376 KB
testcase_21	AC	174 ms; 5,376 KB
testcase_22	AC	30 ms; 5,376 KB
testcase_23	AC	207 ms; 5,376 KB
testcase_24	AC	24 ms; 5,376 KB
testcase_25	AC	57 ms; 5,376 KB
testcase_26	AC	19 ms; 5,376 KB
testcase_27	AC	2 ms; 5,376 KB
testcase_28	AC	2 ms; 5,376 KB
testcase_29	AC	1 ms; 5,376 KB
testcase_30	AC	30 ms; 5,376 KB
testcase_31	AC	27 ms; 5,376 KB
testcase_32	AC	27 ms; 5,376 KB
testcase_33	AC	27 ms; 5,376 KB

ソースコード

#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;
}