cf

题目大意：有个一开始为空的序列。每次操作会往序列最后加一个 $1$ 到 $m$ 的随机整数。当整个序列的 $\gcd$ 为 $1$ 时停止。问这个序列的期望长度对 $10^9+7$ 取模的值。

1.对于 $50\%$ 的数据 : $1\le m\le 1000$

令 $f[i]$ 表示当前的 $gcd$ 为 $i$ (第一个位置的数为 $i$ 开始)，到序列所有的 $gcd$ 为 $1$ 的状态的期望步数，初值 $f[1]=0$，那么它的转移方程就是：

$$f[i]=\frac{1}{m}\ast \sum _{j=1}^m(f[\gcd (i,j)]+1)=1+\frac{1}{m}\ast \sum _{j=1}^{m}f[\gcd (i,j)]$$，

然后我们把 $i==\gcd (i,j)$ 这种提出来：

得到： $$f[i]=1+\frac{1}{m}\ast \sum _{j=1,i\nmid j}^{m}f[\gcd (i,j)]+\frac{1}{m}\ast f[i]\ast \lfloor \frac{m}{i}\rfloor$$

然后移项得： $$f[i]=\frac{1+\frac{1}{m}\ast \sum _{j=1,i\nmid j}^{m}f[\gcd (i,j)]}{1-\lfloor \frac{m}{i}\rfloor \ast \frac{1}{m}}$$

然后求出 $f[1-m]$ ，第一步每个情况都有 $\frac{1}{m}$ 的概率走到,

所以最终答案为：$$ans=\frac{\sum _{i=1}^{m}f[i]}{m}+1$$

对于求每个 $f[i]$ ，这里暴力求 ，复杂度 $O(m^{2}lo{g}_{2}m)$

2.对于 $100\%$ 的数据 : $1\le m\le 100000$

考虑快速求出 $\gcd$ ，我们套路地枚举 $\gcd$，设 $g(i,j)$ 表示有多少个 $1\le x\le m$ 满足 $\gcd (i,x)=j$。
那么就有：

$$f[i]=\frac{1+\frac{1}{m}\ast \sum _{j|i}f[j]\ast g(i,j)}{1-\lfloor \frac{m}{i}\rfloor \ast \frac{1}{m}}$$

然后考虑求 $g(i,j)(j|i)$。

$$g(i,j)=\sum _{x=1}^m[\gcd (i,x)=j]$$

$$g(i,j)=\sum _{x=1}^{\lfloor \frac{m}{j}\rfloor }[\gcd (\frac{i}{j},x)=1]$$

因为
$$\sum _{d|n}u(d)=[n=1]$$

这里把 $n$ 换成 $\gcd (\frac{i}{j},x)$ 得到

$$g(i,j)=\sum _{x=1}^{\lfloor \frac{m}{j}\rfloor }\sum _{d|\gcd (\frac{i}{j},x)}\mu (d)$$

$$g(i,j)=\sum _{d|\frac{i}{j}}\mu (d)\sum _{d|x}^{\lfloor \frac{m}{j}\rfloor }1$$

$$g(i,j)=\sum _{d|\frac{i}{j}}\mu (d)\lfloor \frac{m}{jd}\rfloor$$

然后最终就等于：

$$f[i]=\frac{1+\frac{1}{m}\ast \sum _{j|i}f[j]\ast \sum _{d|\frac{i}{j}}\mu (d)\lfloor \frac{m}{jd}\rfloor }{1-\lfloor \frac{m}{i}\rfloor \ast \frac{1}{m}}$$

复杂度大概 $O(m\ast \sqrt{m})$