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

标签：turn max logs type amp main 答案 void nbsp

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]地震后的幻想乡

标签：turn max logs type amp main 答案 void nbsp

原文地址：http://www.cnblogs.com/ccz181078/p/6607345.html