#[allow(unused_imports)] use std::{ convert::{Infallible, TryFrom, TryInto as _}, fmt::{self, Debug, Display, Formatter,}, fs::File, hash::{Hash, Hasher, BuildHasherDefault}, iter::{Product, Sum}, marker::PhantomData, ops::{Add, AddAssign, Sub, SubAssign, Div, DivAssign, Mul, MulAssign, Neg, RangeBounds}, str::FromStr, sync::{atomic::{self, AtomicU32, AtomicU64}, Once}, collections::{*, btree_set::Range, btree_map::Range as BTreeRange}, mem::{swap}, cmp::{self, Reverse, Ordering, Eq, PartialEq, PartialOrd}, thread::LocalKey, f64::consts::PI, time::Instant, cell::RefCell, io::{self, stdin, Read, read_to_string, BufWriter, BufReader, stdout, Write}, }; pub mod fxhash { use std::hash::BuildHasherDefault; const K: u64 = 0x517c_c1b7_2722_0a95; #[derive(Default)] pub struct FxHasher { pub hash: u64, } impl FxHasher { #[inline(always)] fn mix_u64(mut h: u64, x: u64) -> u64 { h = h.rotate_left(5) ^ x; h = h.wrapping_mul(K); let x2 = x ^ (x >> 33) ^ (x << 11); h = h.rotate_left(5) ^ x2; h = h.wrapping_mul(K); h } #[inline(always)] fn write_u64_impl(&mut self, x: u64) { self.hash = Self::mix_u64(self.hash, x); } } impl std::hash::Hasher for FxHasher { #[inline(always)] fn finish(&self) -> u64 { self.hash } #[inline(always)] fn write(&mut self, bytes: &[u8]) { let mut h = self.hash; for &b in bytes { h = h.rotate_left(5) ^ (b as u64); h = h.wrapping_mul(K); } self.hash = h; } #[inline(always)] fn write_u64(&mut self, i: u64) { self.write_u64_impl(i); } #[inline(always)] fn write_u32(&mut self, i: u32) { self.write_u64_impl(i as u64); } #[inline(always)] fn write_u16(&mut self, i: u16) { self.write_u64_impl(i as u64); } #[inline(always)] fn write_u8 (&mut self, i: u8 ) { self.write_u64_impl(i as u64); } #[inline(always)] fn write_usize(&mut self, i: usize) { self.write_u64_impl(i as u64); } #[inline(always)] fn write_i64(&mut self, i: i64) { self.write_u64_impl(i as u64); } #[inline(always)] fn write_i32(&mut self, i: i32) { self.write_u64_impl(i as u64); } #[inline(always)] fn write_i16(&mut self, i: i16) { self.write_u64_impl(i as u64); } #[inline(always)] fn write_i8 (&mut self, i: i8 ) { self.write_u64_impl(i as u64); } #[inline(always)] fn write_isize(&mut self, i: isize) { self.write_u64_impl(i as u64); } } pub type FxBuildHasher = BuildHasherDefault; pub type FxMap = std::collections::HashMap; pub type FxSet = std::collections::HashSet; } pub fn gcd(mut a: i64, mut b: i64)->i64{if a==0{return b;}else if b==0{return a;}let l1 = a.trailing_zeros();let l2 = b.trailing_zeros(); a >>= l1; b >>= l2;while a!=b{let x = (a^b).trailing_zeros();if a>x;}a << l1.min(l2)} pub fn factorial_i64(n: usize)->(Vec, Vec){ let mut res = vec![1; n+1];let mut inv = vec![1; n+1];for i in 0..n{ res[i+1] = (res[i]*(i+1)as i64)%MOD; } inv[n] = mod_inverse(res[n], MOD);for i in (0..n).rev(){ inv[i] = inv[i+1]*(i+1) as i64%MOD; }(res, inv) } pub fn floor(a:i64, b:i64)->i64{let res=(a%b+b)%b;(a-res)/b} pub fn modulo(a: i64, b: i64)->i64{(a%b+b)%b} pub fn extended_gcd(a:i64,b:i64)->(i64,i64,i64) {if b==0{(a,1,0)}else{let(g,x,y)=extended_gcd(b,a%b);(g,y,x-floor(a,b)*y)}} pub fn mod_inverse(a:i64,m:i64)->i64{let(_,x,_) =extended_gcd(a,m);(x%m+m)%m} pub fn comb(a: i64, b: i64, f: &Vec<(i64, i64)>)->i64{ if aVec<(i64, i64)>{ let mut f=vec![(1i64,1i64),(1, 1)];let mut z = 1i64; let mut inv = vec![0; x as usize+10];inv[1] = 1; for i in 2..x+1{z=(z*i)%MOD; let w=(MOD-inv[(MOD%i)as usize]*(MOD/i)%MOD)%MOD; inv[i as usize] = w; f.push((z, (f[i as usize-1].1*w)%MOD));}return f;} pub fn fast_mod_pow(mut x: i64,p: usize, m: i64)->i64{ x %= m; let mut res=1;let mut t=x;let mut z=p;while z > 0{ if z%2==1{res = (res*t)%m;}t = (t*t)%m;z /= 2; }res} pub trait SortD{ fn sort_d(&mut self); } impl SortD for Vec{ fn sort_d(&mut self) {self.sort_by(|u, v| v.cmp(&u));} } pub trait Mx{fn max(&self, rhs: Self)->Self;} impl Mx for f64{ fn max(&self, rhs: Self)->Self{if *self < rhs{ rhs } else { *self } }} pub trait Mi{ fn min(&self, rhs: Self)->Self; } impl Mi for f64{ fn min(&self, rhs: Self)->Self{ if *self > rhs{ rhs } else { *self } } } pub trait Chmax: PartialOrd + Copy {fn chmax(&mut self, rhs: Self) {if *self < rhs { *self = rhs; }}} impl Chmax for T {} pub trait Chmin: PartialOrd + Copy {fn chmin(&mut self, rhs: Self) {if *self > rhs { *self = rhs; }}} impl Chmin for T {} #[allow(unused)] use proconio::{*, marker::*}; #[allow(unused)] use fxhash::FxMap; #[allow(dead_code)] const INF: i64 = 1<<60; #[allow(dead_code)] const I: i32 = 1<<30; #[allow(dead_code)] const MOD: i64 = 998244353; #[allow(dead_code)] const D: [(usize, usize); 4] = [(1, 0), (0, 1), (!0, 0), (0, !0)]; #[allow(dead_code)] pub fn c2d(c: u8)->(usize, usize){match c{b'U'=>(!0,0),b'D'=>(1,0),b'L'=>(0,!0),b'R'=>(0,1),_=>unreachable!()}} #[allow(dead_code)] pub fn c2d_i64(c: u8)->(i64, i64){match c{b'U'=>(-1,0),b'D'=>(1,0),b'L'=>(0,-1),b'R'=>(0,1),_=>unreachable!()}} #[allow(dead_code)] const D2: [(usize, usize); 8] = [(1, 0), (1, 1), (0, 1), (!0, 1), (!0, 0), (!0, !0), (0, !0), (1, !0)]; pub trait MatrixMonoid { type S: Clone; fn one()->Self::S; fn mul(a: &Self::S, b: &Self::S)->Self::S; fn zero()->Self::S; fn sum(a: &Self::S, b: &Self::S)->Self::S; } #[derive(Debug)] pub struct DoublingMatrix where M: MatrixMonoid{ n: usize, g: Vec, } impl Clone for DoublingMatrix where M: MatrixMonoid, M::S: Clone, { #[inline] fn clone(&self) -> Self { Self { n: self.n, g: self.g.clone() } } } impl DoublingMatrix where M: MatrixMonoid{ #[inline] pub fn new(n: usize, a: &Vec)->Self{ DoublingMatrix {n, g: a.clone() } } #[inline] pub fn zeros(n: usize)->Self{ DoublingMatrix { n, g: vec![M::zero(); n*n] } } #[inline] pub fn e(n: usize)->Self{ let mut res = Self::zeros(n); for i in 0..n{ res.set(i, i, M::one()); } res } #[inline] pub fn get(&self, i: usize, j: usize)-> &M::S{ &self.g[i*self.n+j] } #[inline] pub fn set(&mut self, i: usize, j: usize, v: M::S){ self.g[i*self.n+j] = v; } #[inline(always)] pub fn prod(&self, rhs: &Self)->Self{ let n = self.n; let mut res = vec![M::zero(); n*n]; for i in 0..n { let a_row = &self.g[i * n .. (i + 1) * n]; let out_row = &mut res[i * n .. (i + 1) * n]; for k in 0..n { let a = &a_row[k]; let b_row = &rhs.g[k * n .. (k + 1) * n]; for j in 0..n { let addend = M::mul(a, &b_row[j]); out_row[j] = M::sum(&out_row[j], &addend); } } } Self {n, g:res} } #[inline] pub fn pow(&self, mut k: usize)->Self { let n = self.n; let mut res = Self::e(n); let mut r = (*self).clone(); while k > 0 { if (k & 1) == 1 { res = res.prod(&r); } k >>= 1; if k > 0 { r = r.prod(&r); } } res } } pub struct AddMulMonoid; impl MatrixMonoid for AddMulMonoid{ type S = i64; fn zero()->Self::S { 0 } fn one()->Self::S { 1 } fn sum(&a: &Self::S, &b: &Self::S)->Self::S { if a+b < MOD{a+b}else{a+b-MOD} } fn mul(a: &Self::S, b: &Self::S)->Self::S { a*b%MOD } } pub struct MinPlusMonoid; impl MatrixMonoid for MinPlusMonoid{ type S = i64; fn zero()->Self::S { 1<<60 } fn one()->Self::S { 0 } fn sum(&a: &Self::S, &b: &Self::S)->Self::S { a.min(b) } fn mul(&a: &Self::S, &b: &Self::S)->Self::S { a+b } } struct MM; impl MatrixMonoid for MM{ type S = i32; fn sum(&a: &Self::S, &b: &Self::S)->Self::S { a|b } fn one()->Self::S { 1 } fn mul(&a: &Self::S, &b: &Self::S)->Self::S { a&b } fn zero()->Self::S { 0 } } const MULTI: bool = false; //#[fastout] fn solve(){ input!{ n: usize, m: usize, t: usize, e: [(usize, usize); m], } let mut base = vec![0;n*n]; for &(u, v) in &e{ base[u*n+v] = 1; } let mut mat = DoublingMatrix::::new(n, &base); mat = mat.pow(t); let mut ans = 0; for j in 0..n{ if *mat.get(0, j) > 0{ans += 1;} } println!("{}", ans); } fn main() { if MULTI{ input!{ t: usize, } for _ in 0..t{ solve(); } } else { solve(); } }