#include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #define rep(i, n, s) for(int i = (s); i < int(n); i++) #define per(i, n, s) for(int i = (n - 1); i >= int(s); i--) #define MM << " " << #define all(x) x.begin(), x.end() #define rall(x) rbegin(x), rend(x) template using MinHeap = std::priority_queue, std::greater>; template using MaxHeap = std::priority_queue; using ll = long long; using Pii = std::pair; using Pll = std::pair; using Pdd = std::pair; template bool chmin(T &a, const T b) { if (a > b) { a = b; return true; } else { return false; } } template bool chmax(T &a, const T b) { if (a < b) { a = b; return true; } else { return false; } } template void vdeb(const std::vector &da) { auto n = da.size(); for (size_t i = 0; i < n; i++) { if (i == n - 1) std::cout << da[i]; else std::cout << da[i] << " "; } std::cout << std::endl; } template void vdeb(const std::vector> &da) { auto n = da.size(); for (size_t i = 0; i < n; i++) { std::cout << i << " : "; vdeb(da[i]); } std::cout << std::endl; } template <> void vdeb(const std::vector &da) { auto n = da.size(); for (size_t i = 0; i < n; i++) { std::cout << da[i] << std::endl; } } struct modint { long long num; const static long long p = 998244353; constexpr static long long pow(long long n, long long k) {//n^k(mod p) long long ret = 1; while(k) { if(k&1) ret = ret * n % p; n = n * n % p; k >>= 1; } return ret; } // a*A + b*B = 1 constexpr static void euclid(long long &a, long long &b) { // a>=b A*b+B*(a-a/b*b)=1 if (a == 1) { a = 1; } else { long long A = b, B = a % b; euclid(A, B); b = (A - (p + a / b) % p * B % p + p) % p; a = B; } } constexpr static long long rev(const long long n) {// n*x-p*y=1 //long long q = p; //euclid(p, n, p); //return n % q; return pow(n,p-2); } constexpr modint() : num(0) {} constexpr modint(long long x) : num(x%p < 0 ? x%p+p : x%p) {} constexpr modint inv() const {return rev(num);} modint operator-() const {return modint(p-num);} modint& operator+=(const modint &other){ num = (num + other.num) % p; return *this; } modint& operator-=(const modint &other){ num = (num - other.num + p) % p; return *this; } modint& operator*=(const modint &other){ num = (num * other.num) % p; return *this; } modint& operator/=(const modint &other){ (*this) *= other.inv(); return *this; } modint& operator+=(const long long &other){ num = (num + other) % p; return *this; } modint& operator-=(const long long &other){ num = (num - other + p) % p; return *this; } modint& operator*=(const long long &other){ num = (num * other) % p; return *this; } modint& operator/=(const long long &other){ (*this) *= rev(other); return *this; } modint& operator++(){return *this += 1;} modint& operator--(){return *this -= 1;} modint& operator=(const long long &other){return (*this) = modint(other);} modint operator+(const modint &other) const{return modint(*this) += other;} modint operator-(const modint &other) const{return modint(*this) -= other;} modint operator*(const modint &other) const{return modint(*this) *= other;} modint operator/(const modint &other) const{return modint(*this) /= other;} modint operator+(const long long &other) const{return modint(*this) += other;} modint operator-(const long long &other) const{return modint(*this) -= other;} modint operator*(const long long &other) const{return modint(*this) *= other;} modint operator/(const long long &other) const{return modint(*this) /= other;} bool operator==(const modint &other) const{return num == other.num;} }; std::istream& operator>>(std::istream &is, modint x) { is >> x.num; return is; } std::ostream& operator<<(std::ostream &os, const modint &x){ os << x.num; return os; } template class Matrix { T* table; size_t colomn_, row_; struct Plus { static T apply(const T &l, const T &r) { return l + r; } }; struct Minus { static T apply(const T &l, const T &r) { return l - r; } }; template class Expression { const L& l_; const R& r_; public: Expression(L& l, const R& r) : l_(l), r_(r) {} T at(size_t i, size_t j) const { return Op::apply(l_.at(i, j), r_.at(i, j)); } size_t colomn() const { return l_.colomn(); } size_t row() const { return l_.row(); } }; T* reserve_memory() { T* ret = new T[colomn_*row_]; for(size_t i = 0; i < colomn_*row_; i++) ret[i] = T(); return ret; } void delete_memory(T* p) { delete[] p; } public: Matrix (const size_t colomn, const size_t row) : colomn_(colomn), row_(row) { table = reserve_memory(); } Matrix(const size_t size) : Matrix(size, size) {} Matrix(const std::vector> &_ihs) : Matrix(_ihs.size(), _ihs[0].size()) { for(size_t i = 0; i < colomn_; i++) { for(size_t j = 0; j < row_; j++) { at(i, j) = _ihs[i][j]; } } } Matrix(const Matrix &rhs) : Matrix(rhs.colomn(), rhs.row()) { for(size_t i = 0; i < colomn_; i++) { for(size_t j = 0; j < row_; j++) { at(i, j) = rhs.at(i, j); } } } ~Matrix() { delete_memory(table); } template Expression operator+(const R& rhs) const { assert(colomn() == rhs.colomn() && row() == rhs.row()); return Expression(*this, rhs); } template Expression operator-(const R& rhs) const { assert(colomn() == rhs.colomn() && row() == rhs.row()); return Expression(*this, rhs); } template Matrix operator*(const E& rhs) const { return Matrix(*this)*=rhs; } template Matrix& operator*=(const E& rhs) { assert(row() == rhs.colomn()); T* tmp_p = reserve_memory(); for(size_t i = 0; i < colomn_; i++) { for(size_t j = 0; j < row_; j++) { for(size_t k = 0; k < colomn_; k++) { tmp_p[i*row_+k] += at(i, j) * rhs.at(j, k); } } } delete_memory(table); table = tmp_p; return *this; } template Matrix& operator=(const E& rhs) { for(size_t i = 0; i < colomn_; i++) { for(size_t j = 0; j < row_; j++) { at(i, j) = rhs.at(i, j); } } return *this; } Matrix& operator=(const Matrix& rhs) { if(this == &rhs) return *this; row_ = rhs.row(); colomn_ = rhs.colomn(); T* tmp_p = reserve_memory(); for(size_t i = 0; i < colomn_; i++) { for(size_t j = 0; j < row_; j++) { tmp_p[i*row_+j] = rhs.at(i, j); } } delete_memory(table); table = tmp_p; return *this; } size_t colomn() const {return colomn_;} size_t row() const {return row_;} T& at(size_t i, size_t j) { return table[i*row_+j]; } T at(size_t i, size_t j) const { return table[i*row_+j]; } void output() const { for(size_t i = 0; i < colomn_; i++) { for(size_t j = 0; j < row_; j++) { if(j == row_-1) std::cout << at(i,j) << '\n'; else std::cout << at(i, j) << " "; } } } Matrix e() const { assert(row() == colomn()); Matrix ret(row()); for(size_t i = 0; i < row(); i++) { ret.at(i, i) += 1; } return ret; } }; template Matrix power(Matrix &da, int64_t k) { Matrix ret = da.e(); while(k) { if(k&1) ret *= da; da = da*da; k >>= 1; } return ret; } using namespace std; int main() { ll n, m, t; cin >> n >> m >> t; Matrix da(n); rep(i,m,0) { int a, b; cin >> a >> b; da.at(a, b) += 1; da.at(b, a) += 1; } auto ans = power(da, t); cout << ans.at(0, 0) << endl; }