#include #ifdef DEBUG #include #else #define dump(...) ((void)0) #endif template bool chmin(T &a, const U &b){ return (a > b ? a = b, true : false); } template bool chmax(T &a, const U &b){ return (a < b ? a = b, true : false); } template void fill_array(T (&a)[N], const U &v){ std::fill((U*)a, (U*)(a + N), v); } template auto make_vector(const std::array &a, T value = T()){ static_assert(I >= 1); static_assert(N >= 1); if constexpr (I == 1){ return std::vector(a[N - I], value); }else{ return std::vector(a[N - I], make_vector(a, value)); } } template std::ostream& operator<<(std::ostream &s, const std::vector &a){ for(auto it = a.begin(); it != a.end(); ++it){ if(it != a.begin()) s << " "; s << *it; } return s; } template std::istream& operator>>(std::istream &s, std::vector &a){ for(auto &x : a) s >> x; return s; } std::string YesNo(bool value){return value ? "Yes" : "No";} std::string YESNO(bool value){return value ? "YES" : "NO";} std::string yesno(bool value){return value ? "yes" : "no";} template void putl(const T &value){ std::cout << value << "\n"; } template void putl(const Head head, const Tail &... tail){ std::cout << head << " "; putl(tail ...); } namespace haar_lib { template class vector_dyn { public: using value_type = T; private: std::vector data_; public: vector_dyn(): data_(N){} vector_dyn(T value): data_(N, value){} vector_dyn(std::initializer_list list): data_(N){ int i = 0; for(auto it = list.begin(); it != list.end(); ++it) data_[i++] = *it; } vector_dyn(const vector_dyn &that): data_(that.data_){} template vector_dyn(const std::vector &that): data_(that.begin(), that.end()){} bool operator==(const vector_dyn &that){return data_ == that.data_;} bool operator!=(const vector_dyn &that){return !(*this == that);} auto& operator=(const vector_dyn &that){ data_ = that.data_; return *this; } auto& operator+=(const vector_dyn &that){ for(int i = 0; i < N; ++i) data_[i] += that.data_[i]; return *this; } auto& operator-=(const vector_dyn &that){ for(int i = 0; i < N; ++i) data_[i] -= that.data_[i]; return *this; } friend auto dot(const vector_dyn &a, const vector_dyn &b){ T ret = 0; for(int i = 0; i < N; ++i) ret += a.data_[i] * b.data_[i]; return ret; } auto operator+(const vector_dyn &that) const { return vector(*this) += that; } auto operator-(const vector_dyn &that) const { return vector(*this) -= that; } auto& operator[](int i){return data_[i];} const auto& operator[](int i) const {return data_[i];} auto begin() const {return data_.begin();} auto end() const {return data_.end();} int size() const {return N;} friend std::ostream& operator<<(std::ostream &s, const vector_dyn &a){ s << "{"; for(auto it = a.data_.begin(); it != a.data_.end(); ++it){ if(it != a.data_.begin()) s << ","; s << *it; } s << "}"; return s; } }; template class square_matrix_dyn { public: using value_type = T; using vector_type = vector_dyn; private: std::vector data_; public: square_matrix_dyn(): data_(N, vector_type()){} square_matrix_dyn(const T &val): data_(N, vector_type(val)){} square_matrix_dyn(std::initializer_list> list): data_(N){ int i = 0; for(auto it = list.begin(); it != list.end(); ++it){ data_[i++] = vector_type(*it); } } square_matrix_dyn(const square_matrix_dyn &that): data_(that.data_){} square_matrix_dyn(const std::vector> &that): data_(N){ for(int i = 0; i < N; ++i) data_[i] = that[i]; } bool operator==(const square_matrix_dyn &that) const {return data_ == that.data_;} bool operator!=(const square_matrix_dyn &that) const {return !(*this == that);} auto& operator=(const square_matrix_dyn &that){ data_ = that.data_; return *this; } auto& operator+=(const square_matrix_dyn &that){ for(int i = 0; i < N; ++i) data_[i] += that.data_[i]; return *this; } auto& operator-=(const square_matrix_dyn &that){ for(int i = 0; i < N; ++i) data_[i] -= that.data_[i]; return *this; } auto& operator*=(const square_matrix_dyn &that){ square_matrix_dyn ret; for(int i = 0; i < N; ++i) for(int j = 0; j < N; ++j) for(int k = 0; k < N; ++k) ret[i][j] |= data_[i][k] * that.data_[k][j]; //ret[i][j] += data_[i][k] * that.data_[k][j]; return *this = ret; } const auto& operator[](int i) const {return data_[i];} auto& operator[](int i){return data_[i];} int size() const {return N;} static auto unit(){ square_matrix_dyn ret; for(int i = 0; i < N; ++i) ret[i][i] = 1; return ret; } auto operator+(const square_matrix_dyn &that){ return square_matrix_dyn(*this) += that; } auto operator-(const square_matrix_dyn &that){ return square_matrix_dyn(*this) -= that; } auto operator*(const square_matrix_dyn &that){ return square_matrix_dyn(*this) *= that; } auto pow(uint64_t p) const { auto ret = unit(); auto a = *this; while(p > 0){ if(p & 1) ret *= a; a *= a; p >>= 1; } return ret; } auto operator*(const vector_type &that){ vector_type ret; for(int i = 0; i < N; ++i) ret[i] = dot(data_[i], that); return ret; } }; } namespace haar_lib {} namespace solver { using namespace haar_lib; constexpr int m1000000007 = 1000000007; constexpr int m998244353 = 998244353; void init(){ std::cin.tie(0); std::ios::sync_with_stdio(false); std::cout << std::fixed << std::setprecision(12); std::cerr << std::fixed << std::setprecision(12); std::cin.exceptions(std::ios_base::failbit); } static int N; void solve(){ std::cin >> N; int M; std::cin >> M; int64_t T; std::cin >> T; square_matrix_dyn mat; for(int i = 0; i < M; ++i){ int a, b; std::cin >> a >> b; mat[b][a] = 1; } mat = mat.pow(T); // for(int i = 0; i < N; ++i){ // for(int j = 0; j < N; ++j){ // std::cerr << mat[i][j] << " "; // } // std::cerr << "\n"; // } int ans = 0; for(int i = 0; i < N; ++i){ if(mat[0][i]) ++ans; } std::cout << ans << "\n"; } } int main(){ solver::init(); while(true){ try{ solver::solve(); std::cout << std::flush; std::cerr << std::flush; }catch(const std::istream::failure &e){ break; }catch(...){ break; } } return 0; }