#include #include #define rep(i,n) for(int i=0;i vi; typedef vector vl; typedef vector> vvi; typedef vector> vvl; typedef long double ld; typedef pair P; ostream& operator<<(ostream& os, const modint& a) {os << a.val(); return os;} template ostream& operator<<(ostream& os, const static_modint& a) {os << a.val(); return os;} template ostream& operator<<(ostream& os, const dynamic_modint& a) {os << a.val(); return os;} template istream& operator>>(istream& is, vector& v){int n = v.size(); assert(n > 0); rep(i, n) is >> v[i]; return is;} template ostream& operator<<(ostream& os, const pair& p){os << p.first << ' ' << p.second; return os;} template ostream& operator<<(ostream& os, const vector& v){int n = v.size(); rep(i, n) os << v[i] << (i == n - 1 ? "\n" : " "); return os;} template ostream& operator<<(ostream& os, const vector>& v){int n = v.size(); rep(i, n) os << v[i] << (i == n - 1 ? "\n" : ""); return os;} template ostream& operator<<(ostream& os, const set& se){for(T x : se) os << x << " "; os << "\n"; return os;} template ostream& operator<<(ostream& os, const unordered_set& se){for(T x : se) os << x << " "; os << "\n"; return os;} template ostream& operator<<(ostream& os, const atcoder::segtree& seg){int n = seg.max_right(0, [](S){return true;}); rep(i, n) os << seg.get(i) << (i == n - 1 ? "\n" : " "); return os;} template ostream& operator<<(ostream& os, const atcoder::lazy_segtree& seg){int n = seg.max_right(0, [](S){return true;}); rep(i, n) os << seg.get(i) << (i == n - 1 ? "\n" : " "); return os;} template void chmin(T& a, T b){a = min(a, b);} template void chmax(T& a, T b){a = max(a, b);} using mint = modint998244353; // https://youtu.be/ylWYSurx10A?t=2352 template struct Matrix : vector> { int h, w; Matrix(int h, int w, T val=0): vector>(h, vector(w, val)), h(h), w(w) {} Matrix(initializer_list> a) : vector>(a.begin(), a.end()){ assert(int(this->size()) >= 1); assert(int((*this)[0].size()) >= 1); h = this->size(); w = (*this)[0].size(); rep(i, h) assert(int((*this)[i].size()) == w); } Matrix& unit() { assert(h == w); rep(i,h) (*this)[i][i] = 1; return *this; } Matrix operator*=(const Matrix& M){ assert(w == M.h); Matrix r(h, M.w); rep(i,h) rep(k,w) rep(j, M.w){ r[i][j] += (*this)[i][k] * M[k][j]; } swap(*this, r); return *this; } Matrix operator*(const Matrix& M) const {return (Matrix(*this) *= M);} Matrix operator*=(const T& a) { for(int i = 0; i < h; i++) for(int j = 0; j < w; j++) (*this)[i][j] *= a; return *this; } Matrix operator*(const T& a) const {return (Matrix(*this) *= a);} Matrix operator+=(const Matrix& M){ assert(h == M.h and w == M.w); for(int i = 0; i < h; i++) for(int j = 0; j < w; j++) (*this)[i][j] += M[i][j]; return *this; } Matrix operator+(const Matrix& M) const {return (Matrix(*this) += M);} Matrix pow(long long t) const { assert(h == w); if (!t) return Matrix(h,h).unit(); if (t == 1) return *this; Matrix r = pow(t>>1);r = r*r; if (t&1) r = r*(*this); return r; } }; int main(){ int n, m, t; cin >> n >> m >> t; vector> path(n, vector(n)); vector path_tot(n); auto f = [&](int u, int v){ return u * n + v; }; rep(i, m){ int u, v; cin >> u >> v; u--; v--; if(u != v){ path[u][v]++; path[v][u]++; path_tot[u]++; path_tot[v]++; }else{ path[u][u]++; path_tot[u]++; } } int k = (n + 1) * n; Matrix mat(k, k); rep(prev, n + 1) rep(from, n){ if(path_tot[from] == 0){ mat[f(n, from)][f(prev, from)] = 1; continue; } if(path_tot[from] == 1 and prev != n){ mat[f(n, from)][f(prev, from)] = 1; continue; } rep(to, n){ mint prob = 1; if(path[from][to] == 0) continue; if(prev == n){ prob *= path[from][to]; prob /= path_tot[from]; }else if(prev != n and prev == to){ prob *= path[from][to] - 1; prob /= path_tot[from] - 1; }else if(prev != n and prev != to){ prob *= path[from][to]; prob /= path_tot[from] - 1; } mat[f(from, to)][f(prev, from)] += prob; } } auto mat_pow = mat.pow(t); vector ans(n); rep(from, n + 1) rep(to, n){ ans[to] += mat_pow[f(from, to)][f(n, 0)]; } rep(to, n){ cout << ans[to] << "\n"; } return 0; }