积性函数前缀和求和的方法

　　将近一周的时间内，我专门学习了数论中有关积性函数求前缀和的一些方法，在被虐心的数论折磨得痛不欲生之后（然而我明明已经逃离了数学系为什么还要被数学虐啊摔！！），我终于基本掌握了这一方法，所以在这里记录一下基本的思想并放上一些题解，以便日后回顾。

一、基本知识

1.积性函数

　　积性函数是指这样一类数论函数，对于$\forall a,b\in{N}^{*}$且$gcd(a,b)=1$有$f(ab)=f(a)*f(b)$.常见的积性函数有：

$$除数函数{\sigma}_{k}(n)，表示n的所有约数的k次幂之和$$
$$欧拉函数\varphi(n)，表示所有小于n且与n互质的正整数的个数$$
$$莫比乌斯函数\mu(n)=\left\{ \begin{aligned} 1　　　　　　　　　　　(n=1) \\ {(-1)}^{k}　　(n为k个不同质数之积)\\ 0　　　　　　　　　　(其它情况) \end{aligned} \right.$$

2.完全积性函数

　　若数论函数$f$对于$\forall a,b\in{N}^{*}$有$f(ab)=f(a)*f(b)$，这样的$f$称为完全积性函数。常见的完全积性函数有：

$$幂函数{Id}_{k}(n)={n}^{k}$$

$$单位函数\epsilon(n)=\left\{ \begin{aligned} 1　　(n=1) \\ 0　　(n>1) \end{aligned} \right.$$

3.Dirichlet卷积

　　$Dirichlet$卷积的定义是，对于两个数论函数$f,g$，$f*g(n)=\sum_{d|n}f(d)g(\frac{n}{d})$，$Dirichlet$卷积有如下的一些性质：卷积运算满足交换律、结合律、分配律，单位函数$\epsilon$是$Dirichlet$卷积运算的单位元，两个积性函数的$Dirichlet$卷积仍然是积性函数。

4.两个重要的公式

　　$\sum_{d|n}\varphi(n)=n$
　　$\sum_{d|n}\mu(n)=[n=1]$$([sth.]意思是括号中的条件成立取1，否则取0，以后还会经常碰到这个记号)$
　　这两个公式在之后的推导中有极大的作用。
　

5.Mobius反演

　　定理：若两个数论函数满足$f(n)=\sum_{d|n}g(d)$，那么$g(n)=\sum_{d|n}\mu(d)f(\frac{n}{d})$。在接下来的讨论中将会看到，这一定理对之后的许多讨论有重要的作用。

6.线性筛

　　线性筛是指一种在线性时间复杂度内求出$n$以内所有正整数的最小质因子的算法。该算法描述如下：从2开始，若一个数没有被筛过，则它为质数，最小质因子为其本身，之后对所有的质数$j\leq min[i]$，有$min[i*j]=j$。这样，由于每个合数仅被其最小质因子筛去一次，所以算法的复杂度是线性的。
　　
　　接下来，就是具体的求积性函数前缀和的算法了，但是我不打算把它详细的写出来，而是在题解中体现这样的思想。
　　

二、题解

　　首先，我们在线性筛的基础上改进一下就能在线性时间内求出积性函数的值。我们增加一个数组记录$n$的最小质因子的次数，在做线性筛的同时，利用积性函数的性质，就可以同时计算出它的值。设$f$是一个积性函数，首先根据定义计算$f(1)$，然后从2开始，先对$i$进行线性筛的操作，然后分两种情况：设$i$的最小质因子为$p$，次数为$k$，若$\frac{i}{{p}^{k}}>1$，根据积性函数的性质有$f(i)=f({p}^{k})*f(\frac{i}{{p}^{k}})$，否则