#include #include #include #include #include #include #include #include #include #include #include #include #include #include #include using namespace std; int p[220], q[220], r[220]; int main() { long long mod = 1000000007; int k, m, m_size; cin >> k >> m; m_size = k * k + k; vector> mat(50); long long y; cin >> y; y -= 2; vector> ans(50); for (int i = 0; i < m_size; i++) { for (int j = 0; j < m_size; j++) { mat[i].emplace_back(0); } } for (int i = 0; i < m; i++) { cin >> p[i] >> q[i] >> r[i]; p[i]--, q[i]--, r[i]--; mat[q[i] * k + r[i]][p[i] * k + q[i]]++; } for (int i = 0; i < m_size; i++) { for (int j = 0; j < m_size; j++) { ans[i].emplace_back(mat[i][j]); } } y--; vector> ans_mem(220); for (int i = 0; i < m_size; i++) { for (int j = 0; j < m_size; j++) { ans_mem[i].emplace_back(0); } } vector> mat_mem(220); for (int i = 0; i < m_size; i++) { for (int j = 0; j < m_size; j++) { mat_mem[i].emplace_back(0); } } while (y > 0) { if ((y & 1) == 1) { for (int i = 0; i < m_size; i++) { for (int j = 0; j < m_size; j++) { ans_mem[i][j] = 0; } } for (int i = 0; i < m_size; i++) { for (int j = 0; j < m_size; j++) { for (int k = 0; k < m_size; k++) { ans_mem[i][j] += ans[i][k] * mat[k][j] % mod; ans_mem[i][j] %= mod; } } } ans = ans_mem; } for (int i = 0; i < m_size; i++) { for (int j = 0; j < m_size; j++) { mat_mem[i][j] = 0; } } for (int i = 0; i < m_size; i++) { for (int j = 0; j < m_size; j++) { for (int k = 0; k < m_size; k++) { mat_mem[i][j] += mat[i][k] * mat[k][j] % mod; mat_mem[i][j] %= mod; } } } mat = mat_mem; y >>= 1; } long long ans1 = 0; for (int i = 0; i < k; i++) { for (int j = 0; j < k; j++) { ans1 += ans[i * k][j]; } } cout << ans1 << endl; }