CF1139D Steps to One

Description

给一个数列,每次随机从 [1,m][1,m] 中选择一个数添加到数列末尾,直至数列的 gcd=1\gcd =1 时停止。求期望长度。
1m1051\leq m\leq 10^5

Solution

考虑求 E(X)E(X) 的一种公式:

E(X)=i1P(Xi)E(X) = \sum_{i\ge 1}P(X\ge i)

考虑 P(Xi)P(X\ge i) 的意义:前 i1i-1 个数的 gcd1\gcd \ne 1
fif_i 表示前 ii 个数,gcd1\gcd \ne 1 的方案数,则有

fi=d=2mμ(d)(md)if_i=-\sum_{d=2}^{m}\mu(d)(\lfloor\frac{m}{d}\rfloor)^i

可以通过考虑枚举最后的 gcd\gcd 并容斥得到上式。
然后将 ff 带回到 E(X)E(X) 的式子里,可以得到

E(X)=i1P(Xi)=i1fi1mi=i1d=2mμ(d)(md)i1mi\begin{aligned} E(X)&=\sum_{i\ge 1}P(X\ge i)\\ &=\sum_{i\ge 1}\frac{f_{i-1}}{m^i}\\ &=-\sum_{i\ge 1}\frac{\sum\limits_{d=2}^{m}\mu(d)(\lfloor\frac{m}{d}\rfloor)^{i-1}}{m^i} \end{aligned}

分数线上下的指数不统一不方便,我们考虑一下怎么统一。事实上只要让 PP 中的 >= 变成 > 即可。
因为长度至少为 11 ,所以有

E(X)=1+i1P(X>i)=1i1d=2mμ(d)(md)imi=1d=2mμ(d)i1(mdm)i\begin{aligned} E(X)&=1+\sum\limits_{i\ge 1}P(X>i)\\ &=1-\sum_{i\ge 1}\frac{\sum\limits_{d=2}^{m}\mu(d)(\lfloor\frac{m}{d}\rfloor)^i}{m^i}\\ &=1-\sum_{d=2}^{m}\mu(d)\sum_{i\ge 1}(\frac{\lfloor\frac{m}{d}\rfloor}{m})^i \end{aligned}

然后对于等比数列求和,我们有

i0xi=1x+1x=11x     (x(1,1))\sum_{i\ge 0}x^i = \frac{1-x^{+\infin}}{1-x} = \frac{1}{1-x}\ \ \ \ \ (x\in (-1,1))

所以原式可以转化为

E(X)=1d=2mμ(d)mdmmd\begin{aligned} E(X) = 1-\sum_{d=2}^{m}\mu(d)\frac{\lfloor\frac{m}{d}\rfloor}{m-\lfloor\frac{m}{d}\rfloor} \end{aligned}

复杂度 O(m)\mathcal{O}(m)

总结一下:在 OI 中出现了无穷项,应该果断往「等比数列」的形式上去靠。