Codeforces 1139D Steps to One 期望,容斥,莫比乌斯反演

题意

一个序列会不断随机产生$[1,m]$之间的一个正整数添加在序列的最后,求整个序列中所有数的最大公约数为$1$时序列长度的期望.

题解

我们先思考转移.很显然$$dp[k]=\frac{{\sum }_{i=1}^{m}dp[gcd(i,k)]}{m}+1$$
但是这样做我们可以成功得到优越的$m^{2}log(m)$超时算法.
考虑使用容斥进行优化.
一般的数论容斥完全看不懂,我就学了一下莫比乌斯反演.
首先运用各种方法筛出$\mu$函数,然后开始思考$[2,m]$的每一个数字的贡献.
令$f_x=\frac{[\frac{m}{x}]}{m-[\frac{m}{x}]}(m以内x的倍数的个数/m以内非x的倍数的数的个数)$
则可得答案为$-{\sum }_{i=2}^m{\mu }_i{f}_i+1$.
写出公式即可.

const int yuzu=2e5,mod=1e9+7;
typedef int fuko[yuzu|10];
fuko mu={0,1};
ll kasumi(ll a,ll b=mod-2) {
  ll zw=1;
  for (;b;b>>=1,a=a*a%mod) if (b&1) zw=zw*a%mod;
  return zw;
}
int main() {
  int i,j,n=read();
  for (i=1;i<=n;++i) 
    for (j=i<<1;j<=n;j+=i) mu[j]-=mu[i];
  ll zw=1;
  for (i=2;i<=n;++i) 
    zw=(zw-mu[i]*(n/i)%mod*kasumi(n-n/i)%mod+mod)%mod; 
  printf("%lld\n",zw);
}

谢谢大家.