loj 2540. 「PKUWC2018」随机算法

题意：

传送门

题解：

算出所有合法方案数然后算期望。
$f[i][S]$表示已经用了$i$个点，当前最大独立集是$S$，直接转移即可。
code：

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#define LL long long
using namespace std;
const LL mod=998244353;
LL ans1=0,ans2=0;
LL f[21][1<<20],num[1<<20],bin[1<<21];
bool check[21][1<<20],map[21][21];
LL n,m;
LL pow(LL a,LL b)
{
    LL ans=1;
    while(b)
    {
        if(b&1) ans=ans*a%mod;
        a=a*a%mod;b>>=1;
    }
    return ans;
}
int main()
{
    scanf("%lld %lld",&n,&m);
    memset(map,false,sizeof(map));
    for(LL i=1;i<=m;i++)
    {
        LL x,y;scanf("%lld %lld",&x,&y);x--;y--;
        map[x][y]=true;map[y][x]=true;
    }
    memset(check,false,sizeof(check));
    for(LL i=0;i<n;i++) check[i][0]=true;
    for(LL i=1;i<1<<n;i++) num[i]=num[i>>1]+(i&1);
    for(LL i=0;i<n;i++) bin[1<<i]=i;
    for(LL i=0;i<n;i++)
        for(LL s=1;s<1<<n;s++)
        {
            if((1<<i)&s) continue;
            LL t=s&(-s);
            if(check[i][s^t]&&(!map[i][bin[t]])) check[i][s]=true;
        }
    f[0][0]=1;
    for(LL i=0;i<n;i++)
        for(LL s=0;s<1<<n;s++)
            if(f[i][s])
            {
                LL k=0;
                for(LL j=0;j<n;j++)
                {
                    if(!check[j][s]) continue;
                    k++;(f[i+1][s|(1<<j)]+=f[i][s])%=mod;
                }
                (f[i+1][s]+=f[i][s]*(n-i-k)%mod)%=mod;
            }
    for(LL s=0;s<1<<n;s++)
        if(f[n][s])
        {
            if(num[s]>ans1) ans1=num[s],ans2=0;
            if(num[s]==ans1) (ans2+=f[n][s])%=mod;
        }
    LL sum=1;for(LL i=1;i<=n;i++) (sum*=i)%=mod;
    printf("%lld",ans2*pow(sum,mod-2)%mod);
}