#include using namespace std; using ll = long long; const ll modc = 998244353; class mint { ll x; public: mint(ll x=0) : x((x%modc+modc)%modc) {} mint operator-() const { return mint(-x);} mint& operator+=(const mint& a) { if ((x += a.x) >= modc) x -= modc; return *this; } mint& operator-=(const mint& a) { if ((x += modc-a.x) >= modc) x -= modc; return *this; } mint& operator*=(const mint& a) { (x *= a.x) %= modc; return *this; } mint operator+(const mint& a) const { mint res(*this); return res+=a; } mint operator-(const mint& a) const { mint res(*this); return res-=a; } mint operator*(const mint& a) const { mint res(*this); return res*=a; } mint pow(ll t) const { if (!t) return 1; mint a = pow(t>>1); a *= a; if (t&1) a *= *this; return a; } mint inv() const {return pow(modc-2);} mint& operator/=(const mint& a){ return (*this) *= a.inv();} mint operator/(const mint& a) const { mint res(*this); return res/=a; } bool operator == (const mint& a) const{ return x == a.x;} bool operator != (const mint& a) const{ return x != a.x;} friend ostream& operator<<(ostream& os, const mint& m){ os << m.x; return os; } friend istream& operator>>(istream& ip, mint &m) { ll t; ip >> t; m = mint(t); return ip; } ll val(){ return x;} }; template using matrix = vector>; template void show(matrix &a){ for (int i=0; i matrix dot(matrix &a, matrix &b){ int N = a.size(); matrix c(N, vector(N)); for (int i=0; i T det(matrix a){ int N = a.size(); T res = 1; for (int i=0; i matrix pow(matrix a, ll N){ int M = a.size(); matrix b(M, vector(M)); for (int i=0; i>= 1; a = dot(a, a); } return b; } int main(){ ll N, M, T, a, b, ans=0; cin >> N >> M >> T; matrix G(N, vector(N)); for (int i=0; i> a >> b; G[b][a] = 1; } G = pow(G, T); for (int i=0; i