#pragma GCC optimize("Ofast") #include using namespace std; typedef long long int ll; typedef unsigned long long int ull; mt19937_64 rng(chrono::steady_clock::now().time_since_epoch().count()); ll myRand(ll B) { return (ull)rng() % B; } template struct static_modint { using mint = static_modint; int x; static_modint() : x(0) {} static_modint(int64_t y) : x(y >= 0 ? y % mod : (mod - (-y) % mod) % mod) {} mint &operator+=(const mint &rhs) { if ((x += rhs.x) >= mod) x -= mod; return *this; } mint &operator-=(const mint &rhs) { if ((x += mod - rhs.x) >= mod) x -= mod; return *this; } mint &operator*=(const mint &rhs) { x = (int)(1LL * x * rhs.x % mod); return *this; } mint &operator/=(const mint &rhs) { return *this = *this * rhs.inv(); } mint pow(long long n) const { mint _x = *this, r = 1; while (n) { if (n & 1) r *= _x; _x *= _x; n >>= 1; } return r; } mint inv() const { return pow(mod - 2); } mint operator+() const { return *this; } mint operator-() const { return mint() - *this; } friend mint operator+(const mint &lhs, const mint &rhs) { return mint(lhs) += rhs; } friend mint operator-(const mint &lhs, const mint &rhs) { return mint(lhs) -= rhs; } friend mint operator*(const mint &lhs, const mint &rhs) { return mint(lhs) *= rhs; } friend mint operator/(const mint &lhs, const mint &rhs) { return mint(lhs) /= rhs; } friend bool operator==(const mint &lhs, const mint &rhs) { return lhs.x == rhs.x; } friend bool operator!=(const mint &lhs, const mint &rhs) { return lhs.x != rhs.x; } friend ostream &operator<<(ostream &os, const mint &p) { return os << p.x; } friend istream &operator>>(istream &is, mint &a) { int64_t t; is >> t; a = static_modint(t); return (is); } }; const unsigned int mod = 998244353; using modint = static_modint; modint mod_pow(ll n, ll x) { return modint(n).pow(x); } modint mod_pow(modint n, ll x) { return n.pow(x); } int main() { ios::sync_with_stdio(false); cin.tie(nullptr); int n, m, k; cin >> n >> m >> k; vector> g(n); for (int i = 0; i < m; ++i) { int x, y; cin >> x >> y; x -= 1, y -= 1; g[x].push_back(y); g[y].push_back(x); } using vp = pair, int>; map dp; vector v; for (int i = 0; i < n; ++i) { vector u(n); u[i] = 1; dp[{u, i}] = 1; v.push_back({u, i}); } for (int l = 0; l < v.size(); ++l) { auto vec = v[l].first; int i = v[l].second; auto cost = dp[v[l]]; for (int j : g[i]) { if (vec[j] == k) continue; vec[j] += 1; if (dp.find({vec, j}) == dp.end()) { v.push_back({vec, j}); } dp[{vec, j}] += cost; vec[j] -= 1; } } modint res = 0; vector u(n, k); for (int i = 0; i < n; ++i) { res += dp[{u, i}]; } cout << res << endl; }