新千题计划 10#：[杭电 4624] 无穷动（Endless Spin）

无穷动（Endless Spin¹） n 个白球，每次随机选择一段区间染黑，求使所有球都为黑色的次数期望。

本题解所有内容转自 HDU4624 Endless Spin(概率&&dp) - chanme - 博客园。

容斥。 涉及到一个神奇的转化。

1. 首先要做的是一个题目的转化。如果我定义pi为 恰好i次将区间涂黑的概率，那么显然期望 $E=0p_0+1p_1+2p_2+\cdots$ 换一种角度看这个公式，其实这个公式可以这么写
   $$\begin{gathered} E=p1+p2+p3+p4+p5+\cdots \\ +p2+p3+p4+p5+\cdots \\ +p3+p4+p5+\cdots \end{gathered}$$
   定义 $L_i=p_{i+1}+p_{i+2}+\cdots$ 那么我们可以把Li理解成覆盖了i次都没有将整个区间涂黑的概率。
   所以 $E=L_0+L_1+L_2+\cdots$
2. 理解了上面这一点之后我们可以先考虑一种暴力的做法，那就是枚举 $2^n$ 个最后留下来的点，譬如最后留下来的点是 $v_1,v_2,v_3\cdots v_k$，那么实际上可选的区间只有 $[1,v_1-1],[v_1+1,v_2-1],[v_2+1,v_3-1]\cdots$ 那么这里区间的选法就有 $A$ 种（$A$ 是可以求出来的），要使这些点都能被留下来，那么我们每次选的只能是这 $A$ 个区间，因而概率 $p=\frac{A}{n(n+1)/2}$. 那么对于 $L_i$ 来说 $L_i+=p^i\ast (-1{)}^{k－1}$. 之所以要乘 $(-1{)}^k$ 是因为这里需要有一个容斥。考虑到$p^i$ 的累加和即 $p+p^1+p^2+\cdots =\frac{1}{1-p}$ 所以对于每个 $2^n$，它对最后期望的影响实际上是 $\frac{1}{1-p}\cdot (-1{)}^{k-1}$.
3. 暴力的想法给到了，那么下面就是要优化一下这种做法了，因为我总不可能 $2^{50}$ 枚举所有剩下的点。一个自然的想法是考虑 DP，看上面的式子我们不难发现实际上最后起作用的是 $A$ 和 $k$ 的奇偶性，而 $A$ 的范围是有限的（总的区间个数），那么我们可以定义一个状态 $dp[i][j][k]$ 表示，第 $i$ 个点是白点（未被覆盖），且前 $i$ 个白点的奇偶性为 $j$，可选区间的个数为 $k$ 有多少个子集。 那么对于一个期望值 $E_n(n\ge i)$ 来说这个状态对 $E_n$ 的影响实际上是
   可选的区间 $A=$ 前i个点的可选区间数$+[i..n]$这一段为黑的可选区间数 $＝k+\frac{(n-i+1)(n-i)}{2}$;
   $$\begin{gathered} p=\frac{A}{n(n+1)/2} \\ {E}_n+=dp[i][j][k]\cdot (1-p)\ast (-1{)}^{j-1} \end{gathered}$$

核心代码，亦转自原作者。

dp[0][0][0] = 1;
  for (int i = 0; i <= 50; ++i)
    for (int j = 0; j <= (i + 1)*i / 2; ++j) {
      if (dp[i][0][j])
        for (int k = i + 1; k <= 50; ++k)
          dp[k][1][j + (k - i)*(k - i - 1) / 2] += dp[i][0][j];
      if (dp[i][1][j])
        for (int k = i + 1; k <= 50; ++k)
          dp[k][0][j + (k - i)*(k - i-1) / 2] += dp[i][1][j];
      }
 
for (int i = 1; i <= 50; ++i)
  for (int k = 0; k <= i; ++k)
    for (int j = 0; j <= (k + 1)*k / 2; ++j) {
      if (dp[k][0][j]) {
        long double p = (j + (i - k + 1)*(i - k) / 2 + .0) / ((i + 1)*i / 2);
        if (p == 1.0) continue;
        E[i] -= dp[k][0][j] / (1 - p);
      }
      if (dp[k][1][j]){
        long double p = (j + (i - k + 1)*(i - k) / 2 + .0) / ((i + 1)*i / 2);
        if (p == 1.0) { continue; }
        E[i] += dp[k][1][j] / (1 - p);
      }
    }

---

1. 为防止老实人纠正，此处的 Spin 应理解为“纺纱”而非“自转”。 ↩︎