#include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #define llint long long #define inf 1e18 #define rep(x, s, t) for(llint (x) = (s); (x) < (t); (x)++) #define Rep(x, s, t) for(llint (x) = (s); (x) <= (t); (x)++) #define chmin(x, y) (x) = min((x), (y)) #define chmax(x, y) (x) = max((x), (y)) #define mod 1000000007 using namespace std; typedef pair P; struct Matrix{ int h, w; vector mat; Matrix(){ h = w = 0; } Matrix(int h, int w){ this->h = h, this->w = w; mat.resize(h*w); } llint& at(int i, int j){ return mat[w*(i-1)+(j-1)]; } static Matrix ident(int size){ Matrix ret(size, size); for(int i = 1; i <= size; i++) ret.at(i, i) = 1; return ret; } Matrix operator*(Matrix& ope){ Matrix ret(h, ope.w); for(int i = 1; i <= h; i++){ for(int j = 1; j <= ope.w; j++){ for(int k = 1; k <= w; k++){ ret.at(i, j) += at(i, k) * ope.at(k, j) % mod; ret.at(i, j) %= mod; } } } return ret; } Matrix pow(llint n){ if(n == 0) return ident(h); if(n % 2){ return pow(n-1) * (*this); } else{ Matrix tmp = pow(n/2); return tmp * tmp; } } void print(){ for(int i = 1; i <= h; i++){ for(int j = 1; j <= w; j++){ cout << at(i, j) << " "; } cout << endl; } } }; llint k, m, n; Matrix mat; llint enc(llint p, llint q) { return (p-1)*k+(q-1)+1; } int main(void) { ios::sync_with_stdio(0); cin.tie(0); cin >> k >> m >> n; mat = Matrix(k*k, k*k); llint p, q, r; for(int i = 1; i <= m; i++){ cin >> p >> q >> r; mat.at(enc(q, r), enc(p, q)) = 1; } Matrix pmat = mat.pow(n-2); Matrix vec = Matrix(k*k, 1); for(int i = 1; i <= k; i++) vec.at(enc(1, i), 1) = 1; vec = pmat * vec; llint ans = 0; for(int i = 1; i <= k; i++) ans += vec.at(enc(i, 1), 1), ans %= mod; cout << ans << endl; return 0; }