「PKUWC2018」随机算法 DP

题解：

一开始想的是对于每个最大独立集计算它的贡献，但是不行。然后就考虑一位一位填数，$f[i][S]$表示填了$i$个数，当前的最大独立集状态为$S$的方案数。然后实际上点可以分为$3$类，一是最大独立集中的点，二是对最大独立集没有贡献的点，三是会使最大独立集扩大的点，清楚三种点的数量，转移分是否会使最大独立集变大来转移即可。复杂度$O(n^2 2^n)$。（不是最优的）

代码：

#include<bits/stdc++.h>
using namespace std;
#define LL long long
#define pa pair<int,int>
const int Maxn=21;
const int inf=2147483647;
const int mod=998244353;
int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9')x=(x<<3)+(x<<1)+(ch^48),ch=getchar();
    return x*f;
}
void upd(int &x,int y){x+=y;if(x>=mod)x-=mod;}
bool mp[Maxn][Maxn];
int n,m,f[Maxn][1<<Maxn],a[Maxn],bin[Maxn],one[1<<Maxn],mx=0,ans=0,inv[Maxn],fin[Maxn];
int tmp[Maxn],lt;
void pre()
{
    bin[0]=1;for(int i=1;i<=n;i++)bin[i]=(bin[i-1]<<1);
    one[0]=0;for(int i=1;i<(1<<n);i++)one[i]=one[i>>1]+(i&1);
    inv[0]=inv[1]=1;for(int i=2;i<=n;i++)inv[i]=(LL)(mod-mod/i)*inv[mod%i]%mod;
    fin[0]=1;for(int i=1;i<=n;i++)fin[i]=(LL)fin[i-1]*inv[i]%mod;
}
int main()
{
    memset(f,0,sizeof(f));
    n=read(),m=read();
    pre();
    for(int i=1;i<=m;i++)
    {
        int x=read(),y=read();
        a[x]|=bin[y-1];
        a[y]|=bin[x-1];
    }
    for(int i=1;i<=n;i++)a[i]|=bin[i-1];
    f[0][0]=1;
    for(int i=0;i<n;i++)
    for(int j=0;j<(1<<n);j++)
    if(f[i][j])
    {
        int t=f[i][j],S=0;
        for(int k=0;k<n;k++)
        if((1<<k)&j)S|=a[k+1];
        if(S)upd(f[i+1][j],(LL)t*(one[S-j]-i+one[j])%mod);
        for(int k=1;k<=n;k++)
        if(!((1<<(k-1))&S))upd(f[i+1][j|(1<<(k-1))],t);
    }
    for(int i=0;i<(1<<n);i++)
    if(f[n][i])mx=max(mx,one[i]);
    for(int i=0;i<(1<<n);i++)
    if(one[i]==mx)upd(ans,f[n][i]);
    printf("%d",(LL)ans*fin[n]%mod);
}