#include "bits/stdc++.h" using namespace std; #define int long long #define FOR(i, a, b) for(int i=(a);i<(b);i++) #define RFOR(i, a, b) for(int i=(b-1);i>=(a);i--) #define REP(i, n) for(int i=0; i<(n); i++) #define RREP(i, n) for(int i=(n-1); i>=0; i--) #define REP1(i, n) for(int i=1; i<=(n); i++) #define RREP1(i, n) for(int i=(n); i>=1; i--) #define ALL(a) (a).begin(),(a).end() #define UNIQUE_SORT(l) sort(ALL(l)); l.erase(unique(ALL(l)), l.end()); #define CONTAIN(a, b) find(ALL(a), (b)) != (a).end() #define out(...) printf(__VA_ARGS__) #if DEBUG #define debug(...) printf(__VA_ARGS__) #else #define debug(...) /* ... */ #endif templatebool chmax(T &a,const T &b){if(abool chmin(T &a,const T &b){if(b>K>>M>>N; int p,q,r; vector tups; vector s0; REP(i,M) { cin>>p>>q>>r; --p;--q;--r; if (!p) s0.push_back(tups.size()); tups.push_back({p,q,r}); E[p][q][r]=1; } int size = tups.size(); REP(i,size)REP(j,size) { auto [a,b,c] = tups[i]; auto [d,e,f] = tups[j]; if (E[b][c][d] && E[c][d][e]) { DP2[0][i][j] += 1; } } N -= 3; int T = N/3; // i番目,-3,-2,-1 vector ans(size), tmp(size); for (int i:s0) { ans[i] = 1; } REP(i,61) { REP(j,size)REP(k,size)REP(l,size) { (DP2[i+1][j][l] += DP2[i][j][k]*DP2[i][k][l]) %= MOD; } if (T & (1LL<