loj2540「PKUWC2018」随机算法（状压dp+组合数学）

设f[i][S]表示答案集合S大小为i，不可选的点的状态为S的方案数。
考虑在答案集合中再加入一个j来转移，本来还剩x个点可以选，选了j以后少了y个点可以选，那么这y个点可以随便放在排列的这x个位置之中，方案数为$A_x^y$。
复杂度$O(2^nn)$

#include <bits/stdc++.h>
using namespace std;
#define ll long long
#define inf 0x3f3f3f3f
#define N 25
#define mod 998244353
inline char gc(){
    static char buf[1<<16],*S,*T;
    if(S==T){T=(S=buf)+fread(buf,1,1<<16,stdin);if(T==S) return EOF;}
    return *S++;
}
inline int read(){
    int x=0,f=1;char ch=gc();
    while(ch<'0'||ch>'9'){if(ch=='-') f=-1;ch=gc();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=gc();
    return x*f;
}
int n,m,w[N],f[N][1100000],fac[N],inv[N],bin[N],cnt[1100000];
inline void inc(int &x,int y){x+=y;x%=mod;}
inline int A(int x,int y){return (ll)fac[x]*inv[x-y]%mod;}
int main(){
//  freopen("a.in","r",stdin);
    n=read();m=read();bin[0]=1;inv[0]=inv[1]=fac[0]=1;
    for(int i=2;i<=n;++i) inv[i]=(ll)inv[mod%i]*(mod-mod/i)%mod;
    for(int i=1;i<=n;++i) bin[i]=bin[i-1]<<1,w[i-1]|=bin[i-1];
    for(int i=1;i<=n;++i) fac[i]=(ll)fac[i-1]*i%mod,inv[i]=(ll)inv[i-1]*inv[i]%mod;
    while(m--){
        int x=read()-1,y=read()-1;
        w[x]|=bin[y];w[y]|=bin[x];
    }f[0][0]=1;
    for(int i=0;i<bin[n]-1;++i)
        for(int j=0;j<n;++j) cnt[i]+=(bool)(i&bin[j]);
    for(int i=0;i<n;++i)
        for(int s=0;s<bin[n]-1;++s){
            if(!f[i][s]) continue;
            for(int j=0;j<n;++j){
                if(s&bin[j]) continue;
                inc(f[i+1][s|w[j]],(ll)f[i][s]*A(n-cnt[s]-1,cnt[w[j]-(w[j]&s)]-1)%mod);
            }
        }
    for(int i=n;i>=0;--i){
        if(!f[i][bin[n]-1]) continue;
        printf("%lld\n",(ll)f[i][bin[n]-1]*inv[n]%mod);return 0;
    }
}