看到 n <= 20 应该想到状压 !
每个点有3个状态, 被考虑过在集合中, 被考虑过不在集合中, 没有被考虑过
很明显状态太多, 我们看怎么压成两个状态
首先在集合中必须要, 状态设为1, 那么不在集合中, 不管有没有考虑我们都设成0
换句话说, 不能放进集合的点我们不需要知道它是谁, 只需要知道个数就可以了
我们设 f[i][S] 表示考虑了 i 个点, 集合为 S 的方案数
我们枚举下一个点, 如果能放进集合, 就更新 f[i+1][S'], 并将 cnt ++
剩下有一部分是不能放进集合的, 有多少个呢, n - cnt - i
于是 f[i][S] 对 f[i+1][S] 有 n - cnt - i 的贡献, 加上即可
#include<bits/stdc++.h>
#define N 25
#define M (1<<20)
using namespace std;
typedef long long ll;
int n, m, g[N], bin[M]; ll f[N][M];
const int Mod = 998244353;
ll mul(ll a, ll b){ return (a*b) % Mod;}
void Add(ll &a, ll b){ a = (a+b) % Mod;}
ll power(ll a, ll b){ ll ans = 1; for(;b;b>>=1){if(b&1) ans = mul(ans, a); a = mul(a, a);} return ans;}
int main(){
scanf("%d%d", &n, &m);
for(int i=1; i<=m; i++){
int x, y; scanf("%d%d", &x, &y);
g[x] |= (1<<(y-1)); g[y] |= (1<<(x-1));
}
for(int i=0; i<(1<<n); i++) bin[i] = bin[i>>1] + (i&1);
f[0][0] = 1;
for(int i=0; i<n; i++){
for(int S = 0; S < (1<<n); S++){
if(f[i][S]){
int cnt = 0;
for(int j = 1; j <= n; j++){
if(!(g[j] & S) && !((1 << (j-1)) & S)){
Add(f[i+1][S | (1<<(j-1))], f[i][S]);
cnt++;
}
}
Add(f[i+1][S], mul((n - cnt - i), f[i][S]));
}
}
}
ll ans = 0, inv = 1; for(int i=2; i<=n; i++) inv = mul(inv, i);
inv = power(inv, Mod - 2); int mx = 0;
for(int S = 0; S < (1<<n); S++){
if(!f[n][S]) continue;
if(bin[S] > mx){ mx = bin[S]; ans = f[n][S];}
else if(bin[S] == mx) Add(ans, f[n][S]);
} printf("%lld", mul(ans, inv)); return 0;
}