#include using namespace std; template ostream &operator<<(ostream &s, const pair &v) { s << "(" << v.first << ", " << v.second << ")"; return s; } template >> istream &operator>>(istream &s, T &&v) { for (auto &&x : v) s >> x; return s; } template >> ostream &operator<<(ostream &s, T &&v) { for (auto &&x : v) s << x << ' '; return s; } #ifdef LOCAL template void dbg(T... x) { char e{}; ((cerr << e << x, e = ' '), ...); } #define debug(x...) dbg(#x, '=', x, '\n') #else #pragma GCC optimize("O3,unroll-loops") #pragma GCC target("avx2,bmi,bmi2,lzcnt,popcnt") #define debug(...) ((void)0) #endif #define all(v) (v).begin(), (v).end() #define rall(v) (v).rbegin(), (v).rend() #define ff first #define ss second template inline constexpr T inf = numeric_limits::max() / 2; template bool chmin(T &a, T b) { return (b < a and (a = b, true)); } template bool chmax(T &a, T b) { return (a < b and (a = b, true)); } using u32 = unsigned int; using i64 = long long; using u64 = unsigned long long; using i128 = __int128; constexpr i64 mod = 998244353; void solve() { int n, m, k; cin >> n >> m >> k; vector g(n); for (int i = 0; i < m; i++) { int u, v; cin >> u >> v; u--, v--; g[u] |= (1 << v); g[v] |= (1 << u); } vector pw(n + 1, 1); for (int i = 1; i <= n; i++) { pw[i] = pw[i - 1] * (k + 1); } vector dp(n, vector(pw[n])); for (int i = 0; i < n; i++) { dp[i][pw[i]] = 1; } for (int s = 0; s < pw[n]; s++) { for (int i = 0; i < n; i++) if (dp[i][s]) { for (int j = 0; j < n; j++) if ((g[i] >> j & 1) and (s / pw[j] % (k + 1)) < k) { (dp[j][s + pw[j]] += dp[i][s]) %= mod; } } } i64 ans = 0; for (int i = 0; i < n; i++) (ans += dp[i].back()) %= mod; cout << ans << '\n'; } signed main() { cin.tie(0)->sync_with_stdio(false); cin.exceptions(cin.failbit); int t = 1; // cin >> t; while (t--) { solve(); } return 0; }