bzoj3925 ZJOI2015 地震后的幻想乡

题意

给定$n$个点，$m$条边的图$G$，每条边$e_i$是0~1的实数，求$G$的最小生成树中边的最大值的期望。

$n\le 10$, $m\le n(n-1)/2$

---

题解

对于$n$个$[0,1]$之间的随机变量$x_1,x_2,\dots ,x_n$，第$k$小的期望值是$\frac{k}{n+1}$。

设$P_{S,i}$表示点集$S$选择$i$条边不连通的概率。

刚好$i$条边连通的概率是$(1-P_{S,i})-(1-P_{S,i-1})=P_{S,i-1}-P_{S,i}$

令总点集为$all$

$\begin{array}{rl}ans& =\frac{1}{m+1}(P_{all,0}-P_{all,1})+\frac{2}{m+1}(P_{all,1}-P_{all,2})+...+\frac{m+1}{m+1}{P}_{all,m}\\ & =\frac{1}{m+1}\sum _{i=0}^m{P}_{all,m}\end{array}$

如何计算$P$：

求$P$即求不连通的方案数除以总方案数。

设$f[S][i]$表示点集$S$选择$i$条边不连通的方案数，$g[S][i]$表示点集$S$选择$i$条边连通的方案数。

$P_{S,i}=\frac{f[S][i]}{\left(\genfrac{}{}{0ex}{}{i}{cnt[S]}\right)}$

设$cnt[S]$表示点集$S$形成的子图的边数。

显然$f[S][i]+g[S][i]=\left(\genfrac{}{}{0ex}{}{i}{cnt[S]}\right)$

转移枚举包含某定点$S$的真子集：$f[S][i+j]={\sum }_{T\subset S}g[T][i]\ast \left(\genfrac{}{}{0ex}{}{j}{cnt[S-T]}\right),i\in [0,cnt[T]],j\in [0,cnt[S-T]]$

---

代码

#include <bits/stdc++.h>
using namespace std;
typedef long long ll ;
int e[15], sz[5010], cnt[5010] ;
ll c[50][50], f[5010][50], g[5010][50] ;
int main() {
    int n, m ;
    cin >> n >> m ;
    for (int i = 1; i <= m; i ++) {
        int u, v ;
        scanf("%d%d", &u, &v) ;
        u --; v -- ;
        e[u] |= (1 << v) ;
        e[v] |= (1 << u) ;
    }
    c[0][0] = 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] ;
    }
    for (int s = 1; s < (1 << n); s ++) {
        sz[s] = sz[s >> 1] + (s & 1) ;
        if (sz[s] == 1) {
            g[s][0] = 1; continue ;
        }
        for (int i = 0; i < n; i ++)
            if (s & (1 << i))
                cnt[s] += sz[e[i] & s] ;
        cnt[s] >>= 1 ;
        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[s ^ t]; j ++)
                        f[s][i + j] += g[t][i] * c[cnt[s ^ t]][j] ;
        for (int i = 0; i <= cnt[s]; i ++) g[s][i] = c[cnt[s]][i] - f[s][i] ;
    }
    double ans = 0 ;
    for (int i = 0; i <= m; i ++)
        ans += (double)f[(1 << n) - 1][i] / c[cnt[(1 << n) - 1]][i] ;
    ans /= (m + 1) ;
    printf("%.6f\n", ans) ;
    return 0 ;
}