#include #include #include #include #include #include #include #include #include static const int MOD = 1000000007; using ll = long long; using u32 = unsigned; using u64 = unsigned long long; using namespace std; template constexpr T INF = ::numeric_limits::max() / 32 * 15 + 208; template struct SquareMatrix { using T = typename H::T; using ar = array; using mat = array; mat A; SquareMatrix() = default; static SquareMatrix I(){ SquareMatrix X{}; for (int i = 0; i < SIZE; ++i) { for (int j = 0; j < SIZE; ++j) { if(i == j) X[i][j] = H::one(); else X[i][j] = H::zero(); } } return X; } friend ar operator*=(ar &x, const SquareMatrix &Y) { ar ans{}; for (int i = 0; i < SIZE; ++i) { for (int j = 0; j < SIZE; ++j) { H::add(ans[j], H::mul(x[i], Y[i][j])); } } x.swap(ans); return x; } friend ar operator*(ar x, const SquareMatrix &Y) { return x *= Y; } inline const ar &operator[](int k) const{ return (A.at(k)); } inline ar &operator[](int k) { return (A.at(k)); } SquareMatrix &operator+= (const SquareMatrix &B){ for (int i = 0; i < SIZE; ++i) { for (int j = 0; j < SIZE; ++j) { H::add((*this)[i][j], B[i][j]); } } return (*this); } SquareMatrix &operator-= (const SquareMatrix &B){ for (int i = 0; i < SIZE; ++i) { for (int j = 0; j < SIZE; ++j) { H::add((*this)[i][j], -B[i][j]); } } return (*this); } SquareMatrix &operator*=(const SquareMatrix &B) { SquareMatrix C{}; for (int i = 0; i < SIZE; ++i) { for (int k = 0; k < SIZE; ++k) { for (int j = 0; j < SIZE; ++j) { H::add(C[i][j], H::mul((*this)[i][k], B[k][j])); } } } A.swap(C.A); return (*this); } SquareMatrix pow(ll n) const { SquareMatrix a = (*this), res = I(); while(n > 0){ if(n & 1) res *= a; a *= a; n >>= 1; } return res; } SquareMatrix operator+(const SquareMatrix &B) const {return SquareMatrix(*this) += B;} SquareMatrix operator-(const SquareMatrix &B) const {return SquareMatrix(*this) -= B;} SquareMatrix operator*(const SquareMatrix &B) const {return SquareMatrix(*this) *= B;} }; struct SemiRing { using T = bool; static inline T mul(T x, T y){ return x&y; } static inline void add(T &x, T y){ x |= y; } static inline T one(){ return true; } static inline T zero(){ return false; } }; using ar = array; using mat = SquareMatrix; int main() { int n, m, t; cin >> n >> m >> t; mat X; for (int i = 0; i < m; ++i) { int a, b; scanf("%d %d", &a, &b); X[a][b] = true; } ar p{}; p[0] = true; p *= X.pow(t); int ans = 0; for (int i = 0; i < n; ++i) { ans += p[i]; } cout << ans << "\n"; return 0; }