#include #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=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; using vi=vector; //snuke's modular int template 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 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& operator[](int i)const{ return a[i]; } vector& operator[](int i){ return a[i]; } matrix& operator *=(const matrix& rhs){ int h=a.size(), w=rhs[0].size(), c=rhs.size(); matrix 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 operator *(const matrix& rhs){ return (matrix(*this) *= rhs); } }; template matrix idenmat(int n){ matrix res(n,n); rep(i,n) res.a[i][i] = 1; return res; }; template matrix pow(matrix mat,ll n){ matrix res=idenmat(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 mat(wys); matrix now=idenmat(n); now=now*pow(mat,t); cout<