[PKUWC2018]随机算法 [概率DP+状压DP]

传送门

看到 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;
}