#pragma GCC target("avx2") #pragma GCC optimize("O3") #pragma GCC optimize("unroll-loops") #include #define For(i, a, b) for(int (i)=(int)(a); (i)<(int)(b); ++(i)) #define rFor(i, a, b) for(int (i)=(int)(a)-1; (i)>=(int)(b); --(i)) #define rep(i, n) For((i), 0, (n)) #define rrep(i, n) rFor((i), (n), 0) #define fi first #define se second using namespace std; typedef long long lint; typedef unsigned long long ulint; typedef pair pii; typedef pair pll; template bool chmax(T &a, const T &b){if(a bool chmin(T &a, const T &b){if(a>b){a=b; return true;} return false;} template T div_floor(T a, T b){ if(b < 0) a *= -1, b *= -1; return a>=0 ? a/b : (a+1)/b-1; } template T div_ceil(T a, T b){ if(b < 0) a *= -1, b *= -1; return a>0 ? (a-1)/b+1 : a/b; } constexpr int mod = 1000000007; //constexpr lint INF = mod * mod; constexpr int MAX = 200010; using mat = array, 36>; mat c; void mul(mat a, mat b){ int n = 36; rep(i, n)rep(j, n) c[i][j] = 0; rep(i, n){ rep(k, n){ rep(j, n){ c[i][j] += (lint)a[i][k] * b[k][j] % mod; if(c[i][j] >= mod) c[i][j] -= mod; } } } } mat pow(mat a, lint n){ int m = 36; mat b; rep(i, m)rep(j, m) b[i][j] = (i == j); while(n){ if(n & 1) mul(a, b), swap(b, c), --n; else mul(a, a), swap(a, c), n >>= 1; } return b; } int main(){ int K, m; lint n; scanf("%d%d%lld", &K, &m, &n); int p[m], q[m], r[m]; rep(i, m){ scanf("%d%d%d", &p[i], &q[i], &r[i]); --p[i]; --q[i]; --r[i]; } mat a; rep(i, 36)rep(j, 36) a[i][j] = 0; int idx[6][6], cur = 0; rep(i, 6)rep(j, 6) idx[i][j] = cur++; rep(i, m){ a[idx[p[i]][q[i]]][idx[q[i]][r[i]]] = 1; } a = pow(a, n-2); int ans = 0; rep(j, K)rep(s, K){ ans += a[idx[0][j]][idx[s][0]]; if(ans >= mod) ans -= mod; } printf("%d\n", ans); }