Learning：数论（二） 积性函数

前置知识：狄利克雷卷积、线性筛
积性函数的定义若数论函数$f$满足$\forall a\perp b,f(a)\cdot f(b)=f(ab)$，则$f$是一个积性函数
完全积性函数：若数论函数$f$满足$\forall a,b\in Z_{+},f(a)\cdot f(b)=f(ab)$，则$f$是一个完全积性函数
常见的积性函数有欧拉函数$\varphi$，莫比乌斯函数$\mu$，约数和${\sigma }_0$等
积性函数的一个性质：若$f,g$都是积性函数，则$f$和$g$的狄利克雷卷积$f\ast g$也为积性函数
证明：
记$h=f\ast g$
对于任意的$a\perp b$，有：
$h(ab)=(f\ast g)(ab)$
$={\sum }_{d|ab}f(d)g(\frac{ab}{d})$
$={\sum }_{d|a,k|b}f(dk)g(\frac{ab}{dk})$
$={\sum }_{d|a,k|b}f(d)f(k)g(\frac{a}{d})g(\frac{b}{k})$
$={\sum }_{d|a}f(d)g(\frac{a}{d}){\sum }_{k|b}f(k)g(\frac{b}{k})$
$=h(a)h(b)$
事实上很多积性函数都可以线性筛求解：
如果一个积性函数$f$，我们可以快速求出$f(1),f(p)$和$f(p^k)$，就可以进行线性筛。
方法：
我们需要辅助数组$low$
记$n$的一个素因子为$p$，次数为$\alpha$
$low[n]$即为$p^{\alpha }$
当$n=low[n]$时，$f(n)=f(p^{\alpha })$，直接计算就好了
当$n\ne low[n]$时，$f(n)=f(\frac{n}{low[n]})f(low[n])$
然后再处理一些细节就好了
代码：

void init(){
	f[1]=f在1的值;
	for(int i=2;i<N;i++){
		if(!vis[i]){
			prm[tot++]=i;
			low[i]=i;
			f[i]=f在p的值;
		}
		for(int j=0,p=prm[j],k=p*i;j<tot&&k<N;j++,p=prm[j],k=p*i){
			vis[k]=1;
			if(!(i%prm[j])){
				low[k]=low[i]*p;
				if(i==low[i])
					f[k]=f在p^α的值;
				else f[k]=f[i/low[i]]*f[low[k]];
				break;
			}
			f[k]=f[i]*f[low[k]=p];
		}
	}
	return;
}

有一些特殊数论函数的线性筛不需要$low$数组，在后续会讲到。