首页 > 代码库 > bzoj3925: [Zjoi2015]地震后的幻想乡

bzoj3925: [Zjoi2015]地震后的幻想乡

Description

 傲娇少女幽香是一个很萌很萌的妹子,而且她非常非常地有爱心,很喜欢为幻想乡的人们做一些自己力所能及的事情来帮助他们。 这不,幻想乡突然发生了地震,所有的道路都崩塌了。现在的首要任务是尽快让幻想乡的交通体系重新建立起来。幻想乡一共有n个地方,那么最快的方法当然是修复n-1条道路将这n个地方都连接起来。 幻想乡这n个地方本来是连通的,一共有m条边。现在这m条边由于地震的关系,全部都毁坏掉了。每条边都有一个修复它需要花费的时间,第i条边所需要的时间为ei。地震发生以后,由于幽香是一位人生经验丰富,见得多了的长者,她根据以前的经验,知道每次地震以后,每个ei会是一个0到1之间均匀分布的随机实数。并且所有ei都是完全独立的。 现在幽香要出发去帮忙修复道路了,她可以使用一个神奇的大魔法,能够选择需要的那n-1条边,同时开始修复,那么修复完成的时间就是这n-1条边的ei的最大值。当然幽香会先使用一个更加神奇的大魔法来观察出每条边ei的值,然后再选择完成时间最小的方案。 幽香在走之前,她想知道修复完成的时间的期望是多少呢? 

Input

第一行两个数n,m,表示地方的数量和边的数量。其中点从1到n标号。 
接下来m行,每行两个数a,b,表示点a和点b之间原来有一条边。 
这个图不会有重边和自环。 

Output

一行输出答案,四舍五入保留6位小数。 
假设边权分别为1/(m+1) 2/(m+1) ... m/(m+1),答案不变
枚举集合划分,容斥处理出有几种可能的边的大小关系,使前k小的边恰好使图联通,然后可以统计答案
#include<cstdio>typedef long double ld;int n,m,c[1024],cs[16],mx;ld fac[107],fs[107],ans=0,C[107][107];int max(int a,int b){return a>b?a:b;}void dfs(int w,int mx){    if(w==n){        ld f[107]={};        f[0]=1;        ld v=(mx-1&1)?-1:1;        v*=fac[mx-1];        for(int t=0;t<mx;++t){            int k=c[cs[t]];            for(int i=m;i;--i){                for(int j=max(i-k,0);j<i;++j)f[i]+=f[j]*C[k][i-j];            }        }        for(int i=0;i<=m;++i)fs[i]+=f[i]*v;        return;    }    for(int i=0;i<=mx;++i){        cs[i]^=1<<w;        dfs(w+1,mx+(i==mx));        cs[i]^=1<<w;    }}int main(){    scanf("%d%d",&n,&m);    mx=1<<n;    for(int i=fac[0]=1;i<=100;++i)fac[i]=fac[i-1]*i;    C[0][0]=1;    for(int i=0;i<=60;++i){        for(int j=0;j<=i;++j)C[i+1][j]+=C[i][j],C[i+1][j+1]+=C[i][j];    }    for(int i=0,a,b;i<m;++i){        scanf("%d%d",&a,&b);        --a,--b;        for(int s=0;s<mx;++s)if(s>>a&s>>b&1)++c[s];    }    dfs(0,0);    for(int i=0;i<=m;++i)fs[i]*=fac[i]*fac[m-i];    for(int i=m;i;--i)fs[i]-=fs[i-1],ans+=fs[i]*i;    ans/=fac[m+1];    printf("%.6f",(double)ans);    return 0;}

 

bzoj3925: [Zjoi2015]地震后的幻想乡