傲娇少女幽香是一个很萌很萌的妹子，而且她非常非常地有爱心，很喜欢为幻想乡的人们做一些自己力所能及的事情来帮助他们。 这不，幻想乡突然发生了地震，所有的道路都崩塌了。现在的首要任务是尽快让幻想乡的交通体系重新建立起来。幻想乡一共有n个地方，那么最快的方法当然是修复n-1条道路将这n个地方都连接起来。 幻想乡这n个地方本来是连通的，一共有m条边。现在这m条边由于地震的关系，全部都毁坏掉了。每条边都有一个修复它需要花费的时间，第i条边所需要的时间为ei。地震发生以后，由于幽香是一位人生经验丰富，见得多了的长者，她根据以前的经验，知道每次地震以后，每个ei会是一个0到1之间均匀分布的随机实数。并且所有ei都是完全独立的。 现在幽香要出发去帮忙修复道路了，她可以使用一个神奇的大魔法，能够选择需要的那n-1条边，同时开始修复，那么修复完成的时间就是这n-1条边的ei的最大值。当然幽香会先使用一个更加神奇的大魔法来观察出每条边ei的值，然后再选择完成时间最小的方案。 幽香在走之前，她想知道修复完成的时间的期望是多少呢？
#include<cstdio> #include<cctype> #include<queue> #include<cmath> #include<cstring> #include<algorithm> #define rep(i,s,t) for(int i=s;i<=t;i++) #define dwn(i,s,t) for(int i=s;i>=t;i--) #define ren for(int i=first[x];i;i=next[i]) using namespace std; inline int read() { int x=0,f=1;char c=getchar(); for(;!isdigit(c);c=getchar()) if(c==‘-‘) f=-1; for(;isdigit(c);c=getchar()) x=x*10+c-‘0‘; return x*f; } const int maxn=11; const int maxm=50; typedef long long ll; ll f[1<<maxn][maxm],dp[maxm],C[maxm][maxm]; int n,m,e[maxn],g[1<<maxn],cnt[1<<maxn]; int main() { C[0][0]=1; rep(i,1,49) { C[i][0]=C[i][i]=1; rep(j,1,i-1) C[i][j]=C[i-1][j-1]+C[i-1][j]; } n=read();m=read(); rep(i,1,m) { int u=read()-1,v=read()-1; e[u]|=(1<<v);e[v]|=(1<<u); } rep(S,0,(1<<n)-1) rep(i,0,n-1) if(S>>i&1) cnt[S]++; rep(S,0,(1<<n)-1) { rep(i,0,n-1) if(S>>i&1) g[S]+=cnt[S&e[i]]; g[S]>>=1; } rep(S,1,(1<<n)-1) { int c; rep(i,0,n-1) if(S>>i&1) {c=i;break;} if(cnt[S]==1) f[S][0]=1; else for(int S2=(S-1)&S;S2;S2=(S2-1)&S) if(S2>>c&1) { rep(i,0,g[S2]) rep(j,0,g[S^S2]) dp[i+j]+=f[S2][i]*C[g[S^S2]][j]; } rep(i,0,g[S]) f[S][i]=C[g[S]][i]-dp[i],dp[i]=0; } double ans=0.0; rep(i,0,m) ans+=1.0-(double)f[(1<<n)-1][i]/C[m][i]; printf("%.6lf\n",ans/(m+1)); return 0; }
原文：http://www.cnblogs.com/wzj-is-a-juruo/p/5280944.html
