线性筛五种积性函数

打表五个积性函数（欧拉函数、莫比乌斯函数、因数个数函数、因数个数和函数、某自定义函数

​ 学校专题题目
由于线性筛几乎可以筛所有积性函数（并且本题中的五个积性函数均可以被线性筛），所以我们考虑一次同时筛这五个积性函数。以下考虑每种积性函数各自的筛法。

1、欧拉函数： $\varphi (x)$

（1）当 $x$ 为质数时，显然，小于 $x$ 的数均与 $x$ 互质，所以此时欧拉函数的值即为 $x-1$。

（2）当 $(i\%prime[j]==0)$ 时，$\varphi (i\ast prime[j])=\varphi (i)\ast prime[j]$。

​ 由第一种情况，我们知道 $\phi (n) = n · \prod_{i = 1}^{s} \frac {p_i - 1}{p_i} $（把 $n$ 质因数分解，分别带入各个乘式中，结合第一种情况，且欧拉函数是积性函数，即可证明）。

​ 所以有：
$$\begin{gathered} \varphi (n\ast prime[j])=(n\ast prime[j])\ast \prod _{i=1}^s\frac{p_i-1}{p_i} \\ =i\ast \prod _{i=1}^s\frac{p_i-1}{p_i}\ast prime[j]=\varphi (n)\ast prime[j] \end{gathered}$$
（3）当 $(i\%prime[j]!=0)$ 时，由积性函数性质，可以直接得出：
$$\varphi (i\ast prime[j])=\varphi (i)\ast \varphi (prime[j])$$

2、莫比乌斯函数 : $\mu (x)$

（1）当 $x$ 为质数时，根据莫比乌斯函数定义，显然此时莫比乌斯函数的值为 $-1$。

（2）当 $(i\ % \ prime[j] == 0) $ 时，$\mu (i\ast prime[j])=0$。

​ 因为莫比乌斯函数在质因数分解后只要某个质因子的个数超过 $1$ 个，那么根据定义，莫比乌斯函数的值即为 $0$。

（3）当 $(i\%prime[j]!=0)$ 时，由积性函数性质，可以直接得出：
$$\mu (i\ast prime[j])=\mu (i)\ast \mu (prime[j])=-\mu (i)$$

3、因数个数函数： ${\sigma }_0(x)$

​ 筛因数个数函数时，我们需要一个辅助数组 $num[i]$，表示数字 $i$ 的最小质因子的个数。

​ 为什么要用这个数组？先给出三种情况的结论，然后再解释。

（1）当 $x$ 为质数时，显然，$x$ 的因数只有两个，所以此时因数个数函数的值即为 $2$。

（2）当 $(i\%prime[j]==0)$ 时，${\sigma }_0(i\ast prime[j])={\sigma }_0(i)/(num[i]+1)\ast (num[i]+2)$。

​ $num[i\ast prime[j]]=num[i]+1$。

（3）当 $(i\%prime[j]!=0)$ 时，由积性函数性质，可以直接得出​：
$${\sigma }_0(i\ast prime[j])={\sigma }_0(i)\ast {\sigma }_0(prime[j])$$
​ 同时，我们还可以得到 $num[i\ast prime[j]]=1$ (因为 $i$ 和 $prime[j]$ 互质，所以二者相乘后最小质因数的个数为1)

​ 对 $num[i]$ 数组的解释：我们发现，$num[i]$ 仅在 $(i\%prime[j]==0)$ 时对答案产生贡献。原因是，我们考虑 $n={\prod }_{i=1}^s{p}_i^{a_i}$，那么 ${\sigma }_0(n)={\prod }_{i=1}^s(1+a_i)$，因为对于每个质因子，我们可以不选，可以选 $1$ 个，选 $2$ 个… 选 $a_i$ 个，一共是 $(1+a_i)$ 种选择，根据乘法原理即可得出这个算式。

​ 考虑 $(i\%prime[j]!=0)$，也就是说 $prime[j]$ 是 $i$ 的一个(最小)质因子。那么 ${\sigma }_0(i\ast pirme[j])=({\prod }_{i=1}^s(1+a_i))/(1+a_1)\ast (1+a_1+1)={\sigma }_0(i)/(1+a_1)\ast (1+a_1+1)$，(因为最小的质因子个数加了一个，所以除以原来的个数再乘以原来的个数加 $1$)

​ 这里我们用 $num[i]$ 代表上式中的 $(1+a_1)$，就得到了：
$${\sigma }_0(i\ast prime[j])={\sigma }_0(i)/(num[i]+1)\ast (num[i]+2)$$

4、因数和函数： ${\sigma }_1(x)$

​ 筛因数和函数时我们也需要一个辅助数组 $sp[x]$。对于数 $x$，若其最小质因子为 $p$，且个数为 $a$，

​ 那么 $sp[x]=(1+p+p^2+...+p^a)$。

​ 同样先给出三种情况的结论再做解释。

（1）当 x 为质数时，显然，$x$ 的因数只有 $1$ 和 $x$，所以此时因数和函数的值即为 $1+x$。

（2）当 $(i\%prime[j]==0)$ 时，${\sigma }_1(i\ast prime[j])={\sigma }_1/sp[i]\ast (sp[i]\ast prime[j]+1)$。

​ $sp[i\ast prime[j]]=sp[i]\ast prime[j]+1$

（3）当 $(i\%prime[j]!=0)$ 时，由积性函数性质，可以直接得出：
$${\sigma }_1(i\ast prime[j])={\sigma }_1(i)\ast {\sigma }_1(prime[j])$$
​ 同时，我们还可以得到 $sp[i\ast prime[j]]=1+prime[j]$，因为 $i$ 和 $prime[j]$ 互质，所以二者相乘后最小质因数的几何级数序列就为 $(1,prime[j{]}^1)$，和即为 $(1+prime[j])$。

​ 对 $sp[x]$ 数组的解释：计算某个数字 $N$ 的约数和时，有一个定理：

${\sigma }_1(N)=(1+p_1+p_1^2+...+p_1^{a_1})(1+p_2+p_2^2+...+p_2^{a_2})...(1+p_m+p_m^2+...+p_m^{a_m})$，其中 $p_i$ 指的是 $N$ 的某个质因子，$a_i$ 指的是质因子 $p_i$ 的个数。

​ 有了这个定理，不难发现，当 $prime[j]$ 为 $i$ 的最小质因子时，
$${\sigma }_1(i\ast prime[j])={\sigma }_1(i)/(1+p_1+p_1^2+...+p_1^{a_i})\ast (1+p_1+p_1^2+...+p_1^{a_i}+p_1^{a_i+1})$$
​ 把 $sp[i]$ 数组带入上式可得：
$${\sigma }_1(i\ast prime[j])={\sigma }_1/sp[i]\ast (sp[i]\ast prime[j]+1)$$
​ 同时我们可以得到，$sp[i\ast prime[j]]=sp[i]\ast prime[j]+1$，因为 $prime[j]$ 是 $i$ 的最小质因子，所以 $i\ast prime[j]$ 的最小质因子个数比 $i$ 多一个，结合 $sp[i]$ 的定义即可计算出。

5、自定义函数 $f(x)=(\mu \cdot \varphi )\ast {\sigma }_1$

​ 由于 $\mu (x),\varphi (x)$ 均为积性函数，所以 $(\mu \cdot \varphi )(x)$ 也是积性函数；又 ${\sigma }_1$ 也是积性函数，所以 $f(x)$ 也是积性函数。

​ 利用卷积定义，首先对这个函数化简得到：
$$f(x)=\sum _{d|x}\mu (d)\varphi (d){\sigma }_1(\frac{x}{d})$$
​ 用类似的方法考虑三种情况下 $f(x)$ 的大小。

（1）当 $x$ 为质数时。$d$ 的取值仅有 $1$ 和 $x$，所以有：
$$f(x)=\mu (1)\varphi (1){\sigma }_1(x)+\mu (x)\varphi (x){\sigma }_1(1)={\sigma }_1(x)-x+1$$
（2）当 $(i\%prime[j]==0)$时，$prime[j]$ 为 $i$ 的最小质因子。

​ 假设 $i=p_1^{a_1}{\prod }_{k=2}^s{p}_k^{a_k}$，那么 $i\ast prime[j]=p_1^{a_i+1}{\prod }_{k=2}^s{p}_k^{a_k}$，

​ 由积性函数的性质，有：
$$f(i\ast prime[j])=f(p_1^{a_i+1})\ast f(\prod _{k=2}^s{p}_k^{a_k})$$
​ 所以，计算 $f(i\ast prime[j])$ 时，只需要统计出 $i\ast prime[j]$ 中有多少个质因子 $prime[j]$，把 $i\ast prime[j]$ 拆成 $p_1^{a_i+1}$ 和 ${\prod }_{k=2}^s(p_k^{a_k})$，利用上式就可以求出答案了。

（3）当 $(i\%prime[j]!=0)$时，由积性函数性质，可以直接得出：
$$f(i\ast prime[j])=f(i)\ast f(prime[j])$$

---