2018.06.29 NOIP模拟 Gcd（容斥原理）

传送门

这是我做过的最可（$du$）做（$liu$）的$NOIP$模拟级别的数论题。考场上打算写随机算法，想了想$30$暴力更稳，于是交了暴力，结果$10$分滚粗了。

这题让我们联想到正难则反的思想，题目上要我们求出所有可能解的方案数，那么我们这样想，我们先把总方案数求出来，显然是$n\ast (2^{n-1}-1)$。

然后我们将不合法的情况去掉，我们设加入的数为$number$，那么我们要去掉的就是在加入$number$之前就已经使得$gcd$为$1$的集合对应的$number$的选择方案数之和，以及在加入$number$之后$gcd$仍然不为$1$的方案数。这东西怎么算啊？变形代数式看看吧。

首先解决第二种情况，即在加入$number$之后$gcd$仍然不为$1$的方案数。这种情况直接统计所有集合中使得$gcd$不为一的方案数即可，即用总方案数减去所有集合中使得$gcd$为$1$的方案数。
另外，我们还知道，情况二其实是让我们求这个东西：${\sum }_{|S|,gcd(S)!=1}|S|$

这样的话，我们似乎将第二种情况转化成了第一种情况的子问题。怎么求出所有集合中使得$gcd$为$1$的方案数呢？

考虑容斥原理，我们先将$gcd$为$1$的倍数的集合数求出，减去$gcd$为$2$的倍数的集合数……一直这样迭代下去，然后有个显然的结论，每个$gcd$的容斥系数就对应着它们的莫比乌斯函数值。于是只需要提前用线性筛筛出所有数的莫比乌斯函数值，每次将情况二的答案加上$mu[i]\ast (2^{sum[i]}-1)$就行了，其中$sum[i]$表示所有数列中所有是$i$的倍数的数的个数。

情况一？考虑用情况二反推情况一，再次将式子变形：我们要求的是在加入$number$之前就已经使得$gcd$为$1$的集合对应的$number$的选择方案数之和，即${\sum }_{S,gcd(S)==1}(n-|S|)$,显然这个式子可以展开成${\sum }_{S,(gcd)==1}1-{\sum }_{S,gcd(S)==1}|S|$,哇，后面一坨不是可以跟情况二合并吗？好像还可以抵掉一坨东西。最后发现只需要容斥计算乱搞一下就行了啊。

代码