#include #include #define rep(i,n) for(int i=0;i vi; typedef vector vl; typedef vector> vvi; typedef vector> vvl; typedef long double ld; typedef pair P; ostream& operator<<(ostream& os, const modint& a) {os << a.val(); return os;} template ostream& operator<<(ostream& os, const static_modint& a) {os << a.val(); return os;} template ostream& operator<<(ostream& os, const dynamic_modint& a) {os << a.val(); return os;} template istream& operator>>(istream& is, vector& v){int n = v.size(); assert(n > 0); rep(i, n) is >> v[i]; return is;} template ostream& operator<<(ostream& os, const pair& p){os << p.first << ' ' << p.second; return os;} template ostream& operator<<(ostream& os, const vector& v){int n = v.size(); rep(i, n) os << v[i] << (i == n - 1 ? "\n" : " "); return os;} template ostream& operator<<(ostream& os, const vector>& v){int n = v.size(); rep(i, n) os << v[i] << (i == n - 1 ? "\n" : ""); return os;} template void chmin(T& a, T b){a = min(a, b);} template void chmax(T& a, T b){a = max(a, b);} using mint = modint998244353; int main(){ int n, m, k; cin >> n >> m >> k; vector> G(n, vector(n)); rep(i, m){ int u, v; cin >> u >> v; u--; v--; G[u][v] = 1; G[v][u] = 1; } auto f = [&](vector cnt){ int pos = 0; rep(i, n){ pos *= (k + 1); pos += cnt[i]; } return pos; }; auto f_inv = [&](int pos){ vector cnt(n); rep(i, n){ cnt[(n - 1) - i] = pos % (k + 1); pos /= k + 1; } return cnt; }; int x = 1; rep(i, n) x *= k + 1; vector> dp(x, vector(n)); rep(i, n){ vector cnt(n); cnt[i] = 1; int pos = f(cnt); dp[pos][i] = 1; } rep(pos_from, x){ rep(from, n){ rep(to, n){ if(G[from][to] == 1){ vector cnt = f_inv(pos_from); if(cnt[to] >= k) continue; cnt[to]++; dp[f(cnt)][to] += dp[pos_from][from]; } } } } mint ans = 0; rep(i, n) ans += dp[x - 1][i]; cout << ans << "\n"; return 0; }