[BZOJ3925]地震后的幻想乡 期望 状压DP

不想写题意了……

不错的题解

网上一些大神好像用积分什么的来解……表示蒟蒻看不懂。
上面的题解是用期望的线性性质，要好懂些，不过比较抽象。

以下是我对最后求答案公式的理解，大神可以跳过……

因为有m条边，可以想象一张有m个点的图，节点按编号依次相连（即0连1,1连2......m-1连m），且每个点都跟终点相连。
从0号点开始移动，第i个点有f(i)的概率移动到下一个节点，1-f(i)的概率移动到终点，问到终点的期望步数。

从节点i移动到节点i+1的概率就相当于整张图随机连i条边不连通的概率，直接移动到终点就相当于随机连i条边且连通的概率，
那么到终点的期望步数就是使图连通期望连的边的数量，除以（m+1）就是答案。

感觉这样建张图要好理解一些，原来那个求答案的公式太抽象了。。。也可能是我太弱了。。。

#include <cstdio>
#define N 10

typedef long long ll;

int n,m,A[N][N],sz[1<<N];
ll cnt[1<<N],f[1<<N][N*N],g[1<<N][N*N],C[N*N][N*N];
double Ans;

int main(){
    freopen("mst.in","r",stdin);
    freopen("mst.out","w",stdout);
    scanf("%d%d",&n,&m);
    for(int i=1,u,v;i<=m;i++)scanf("%d%d",&u,&v),A[u][v]=A[v][u]=1;
    for(int i=1;i<=m;i++){
        C[i][0]=C[i][i]=1;
        for(int j=1;j<i;j++) C[i][j]=C[i-1][j]+C[i-1][j-1];
    }
    C[0][0]=1;
    for(int i=1;i<(1<<n);i++){
        for(int j=1;j<=n;j++)
            if(i&(1<<j-1))for(int k=j+1;k<=n;k++)
                if((i&(1<<k-1))&&A[j][k]) cnt[i]++;
    }
    for(int S=1;S<(1<<n);S++){
        sz[S]=sz[S>>1]+(S&1);
        if(sz[S]==1){g[S][0]=1;continue;}
        int lowbit=S&-S;
        for(int T=(S-1)&S;T;T=(T-1)&S)
            if(T&lowbit){
                for(int i=0;i<=cnt[T];i++)
                    for(int j=0;j<=cnt[T^S];j++)
                        f[S][i+j]+=g[T][i]*C[cnt[T^S]][j];
            }
        for(int i=0;i<=cnt[S];i++) g[S][i]=C[cnt[S]][i]-f[S][i];
    }
    for(int i=0;i<=m;i++) Ans+=(double)f[(1<<n)-1][i]/C[cnt[(1<<n)-1]][i];
    return printf("%.6lf",Ans/(m+1)),0;
}