loj 2540「PKUWC 2018」随机算法

http://www.elijahqi.win/archives/3568
题目描述

我们知道，求任意图的最大独立集是一类NP完全问题，目前还没有准确的多项式算法，但是有许多多项式复杂度的近似算法。

例如，小C常用的一种算法是：

对于一个n个点的无向图，先等概率随机一个 1..n1..n1..n 的排列p[1..n]p[1..n]p[1..n]。
维护答案集合 SSS ，一开始 SSS 为空集，之后按照 i=1..ni=1..ni=1..n 的顺序，检查 {p[i]p[i]p[i]}∪S\cup S∪S 是否是一个独立集，如果是的话就令 S=S=S={p[i]p[i]p[i]}∪S\cup S∪S。
最后得到一个独立集 SSS 作为答案。
小C现在想知道，对于给定的一张图，这个算法的正确率，输出答案对 998244353998244353998244353 取模

输入格式

第一行两个非负整数 n,mn,mn,m 表示给定的图的点数和边数。

接下来 mmm 行，每行有两个正整数 (u,v)(u≠v)(u,v) (u\neq v)(u,v)(u≠v) 描述这张图的一条无向边。

输出格式

输出正确率，答案对 998244353998244353998244353 取模。

样例

样例输入

3 2
1 2
2 3
样例输出

665496236
样例解释

这张图的最大独立集显然为2，可以发现只有p[1]=2p[1]=2p[1]=2时会得出S={2}，否则都是S={1,3}，所以答案是23\frac{2}{3}​3​​2​​

数据范围与提示

对于10%10\%10%的数据，有1≤n≤91\leq n\leq 91≤n≤9

对于30%30\%30%的数据，有1≤n≤131\leq n\leq 131≤n≤13

对于50%50\%50%的数据，有1≤n≤171\leq n\leq 171≤n≤17

另有10%10\%10%的数据，满足给定的图是一条链。

另有10%10\%10%的数据，满足给定的图是一棵树。

对于100%100\%100%的数据，有1≤n≤20,0≤m≤n∗(n−1)21\leq n\leq 20,0\leq m\leq \frac{n*(n-1)}{2}1≤n≤20,0≤m≤​2​​n∗(n−1)​​，保证给定的图没有重边和自环。

设dp[i][s]表示s状态的集合一定不能选 i表示这种状态的最大独立集大小

转移yy一下 或者看代码 可以这么想 选了这个不在集合内的点他会贡献和他相邻的点不能选 那么这些贡献出来的点体现在原序列中怎么选就需要应用排列的公式去计算一下即可

#include<cstdio>
#include<cctype>
#include<algorithm>
#define ll long long
using namespace std;
inline char gc(){
    static char now[1<<16],*S,*T;
    if (T==S){T=(S=now)+fread(now,1,1<<16,stdin);if (T==S) return EOF;}
    return *S++;
}
inline int read(){
    int x=0,f=1;char ch=gc();
    while(!isdigit(ch)) {if(ch=='-') f=-1;ch=gc();}
    while(isdigit(ch)) x=x*10+ch-'0',ch=gc();
    return x*f;
}
const int mod=998244353;
int dp[21][1100000],n,m,jc[22],inv[22],cnt[1100000],w[22];
inline int ksm(ll b,int t){static ll tmp;
    for (tmp=1;t;b=b*b%mod,t>>=1) if (t&1) tmp=tmp*b%mod;return tmp;
}
inline int A(int x,int y){
    return (ll)jc[x]*inv[x-y]%mod;
}
inline void inc(int &x,int v){
    x=x+v>=mod?x+v-mod:x+v;
}
int main(){
    freopen("loj2540.in","r",stdin);
    n=read();m=read();jc[0]=1;
    for (int i=1;i<=n;++i) jc[i]=(ll)jc[i-1]*i%mod;
    inv[n]=ksm(jc[n],mod-2);
    for (int i=n-1;~i;--i) inv[i]=(ll)inv[i+1]*(i+1)%mod;
    for (int i=1;i<=m;++i){
        int x=read()-1,y=read()-1;
        w[x]|=1<<x;w[x]|=1<<y;
        w[y]|=1<<x;w[y]|=1<<y;
    }int S=(1<<n)-1;dp[0][0]=1;
    for (int s=0;s<=S;++s)
        for (int i=0;i<n;++i) if ((1<<i)&s) ++cnt[s];
    for (int i=0;i<n;++i){
        for (int s=0;s<=S;++s){
            if (!dp[i][s]) continue;
            for (int j=0;j<n;++j){
                if ((1<<j)&s)continue;
                inc(dp[i+1][s|w[j]],(ll)dp[i][s]*A(n-cnt[s]-1,cnt[w[j]-(w[j]&s)]-1)%mod);
            }
        }
    }
    for (int i=n;i;--i){
        if (dp[i][S]) {
            printf("%lld\n",(ll)dp[i][S]*inv[n]%mod);return 0;
        }
    }
    return 0;
}