\[ f[1] = 0 \] \[ f[i] = 1 + \frac{1}{m} \sum_{j = 1} ^ n f[gcd(i, j)] \ \ \ \ \ \ (i != 1) \] 然后发现后面这一块gcd的个数只可能是i的约数, 那么考虑枚举约数 \[ f[i] = 1 + \frac{1}{m}\sum_{d | i} f[d] cnt(d, i) \] \(cnt(d, i)\)表示和[1,m]内与i的gcd为d的数字个数 考虑这个东西能够怎么算, \(cnt(d, i)\)显然…
题目链接 啊啊啊我在干什么啊.怎么这么颓一道题做这么久.. 又记错莫比乌斯反演式子了(╯‵□′)╯︵┻━┻ \(Description\) 给定\(n\).有一个初始为空的集合\(S\).令\(g\)表示S中所有数的\(\gcd\).每次随机选择一个\([1,n]\)中的数加到集合\(S\)中去,直到\(g=1\).求集合\(S\)的期望大小.(原题目描述为数列长度,\(n\)是指\(m\),我自己都看混了=-=) \(n\leq10^5\). \(Solution\) 首先不要想\(f[i][…
卧槽,这么秀吗??? 暂时留坑...…
求单个莫比乌斯函数忘记算n本身的质数,WA了一发. http://www.fjutacm.com/Problem.jsp?pid=2360 首先,显然随着n增大,与m互质的数不会变少.可以二分来求k,关键是怎么快速计算n以内和m互质的数的个数. 反过来,我们求n以内与m的gcd至少为2的数的个数. 枚举m的因子d,d从2开始到m,那么每个因子d的贡献就是mu(d),不知道怎么解释好. 比如枚举了gcd为因子2,那么gcd为4,8的情况就已经被2包含了,不用计算,直接mu(4)=mu(8)=0.…
蒟蒻数学渣呀,根本不会做. 解法是参考 https://blog.csdn.net/xs18952904/article/details/88785210 这位大佬的. 状态的设计和转移如上面博客一样:dp[i]代表当前序列的gcd为i的期望长度. 那么可以写出状态转移方程:dp[i]=(1+(x/m)∑(j|i,j≠i)dp[j]) / (1-(m/i)/m) (写得有点乱,其实和上面大佬的一样的) 这里要说一下的是 x=∑(t=1,t<=m) [ gcd(t,i)==j ]  就是怎么求1<…
题目链接: http://codeforces.com/contest/1139/problem/D 题意: 在$1$到$m$中选择一个数,加入到一个初始为空的序列中,当序列的$gcd$和为$1$时,停止加入,求序列的期望长度 数据范围: $1 \leq m \leq 10^{9}$ 分析: 定义$f[x$]为$gcd$等于$x$时把序列$gcd和$改变成1的期望长度,定义$G(x,y)$为$i$在1到$n$范围,满足$gcd(x,i)=y$,$i$的数量,得到以下公式: $$f[i]=1+\f…
题意 题目链接 从\([1, M]\)中随机选数,问使得所有数gcd=1的期望步数 Sol 一个很显然的思路是设\(f[i]\)表示当前数为\(i\),期望的操作轮数,转移的时候直接枚举gcd \(f[i] = 1 + \frac{ \sum_{j=1}^N f[gcd(i, j)]}{N}\) 然后移一下项就可以算出\(f[i]\)了. 发现gcd相同的有很多,可以预处理一下. 复杂度\(O(跑的过)\) 还有一种反演做法表示推不出来qwq #include<bits/stdc++.h> #…
https://codeforces.com/contest/1139/problem/D 题意 每次从1,m中选一个数加入队列,假如队列的gcd==1停止,问队列长度的期望 题解 概率正着推,期望反着推 发现每加入一个数,gcd会变为原来gcd的因数 \(dp[x]\) - > \(dp[gcd(x,i)]\) 但是方程却是反方向的 图片 代码 #include<bits/stdc++.h> #define MOD 1000000007 #define MAXN 100005 #def…
stm这是div2的D题……我要对不住我这个紫名了…… 题目链接:CF原网  洛谷 题目大意:有个一开始为空的序列.每次操作会往序列最后加一个 $1$ 到 $m$ 的随机整数.当整个序列的 $\gcd$ 为 $1$ 时停止.问这个序列的期望长度对 $10^9+7$ 取模的值. $1\le m\le 10^5$. 首先很容易想到DP:$f_i$ 表示目前的 $\gcd$ 为 $i$,期望还要多少次才能结束. 那么有 $f_1=0$. 转移,直接枚举即可:$f_i=1+\dfrac{1}{m}\su…
有一种组合方向的考虑有没有dalao肯高抬啊? 题目大意 有一个初始为空的数组$a$,按照以下的流程进行操作: 在$1\cdots m$中等概率选出一个数$x$并添加到$a$的末尾 如果$a$中所有元素的$\gcd=1$则完成这个数组$a$的修改 重复这一过程 求数组$a$的期望长度,$m \le 10^5,\mod 10^9+7$ 题目分析 质因数分解的期望dp题 以下介绍的两个做法中,第一个做法本人不会所以  求助会做的dalao麻烦高抬一手   : 第二个做法是对推得的dp式子质因数分解求…