import sys def input(): return sys.stdin.readline()[:-1] MOD = 10**9+7 def mat_dot(one, two, mod): return [[sum([(one[i][k]*two[k][j])%mod for k in range(max(len(two), len(one[0])))])%mod for j in range(len(two[0]))] for i in range(len(one))] def mat_pow(mat, exp, mod): size = len(mat) res = [[0 for _ in range(size)] for _ in range(size)] for i in range(size): res[i][i] = 1 cnt = 0 while (1<>cnt)&1: res = mat_dot(res, mat, mod) mat = mat_dot(mat, mat, mod) cnt += 1 return res k, m, n = map(int, input().split()) pre = [0 for _ in range(k*k)] for i in range(k): pre[i] = 1 mat = [[0 for _ in range(k*k)] for _ in range(k*k)] for _ in range(m): p, q, r = map(int, input().split()) mat[(q-1)*k + r-1][(p-1)*k + q-1] += 1 res = mat_pow(mat, n-2, MOD) ans = 0 for i in range(k): for j in range(k*k): ans += pre[j] * res[i*k][j] ans %= MOD print(ans)