数论专题——Dirichlet卷积及积性函数进阶

前置知识

Dirichlet卷积初步

前言

本章节默认$n={\prod }_{i=1}^t{p}_i^{k_i}$

一、欧拉函数

$1.$通项公式
$$\phi (n)=\prod _{i=1}^t(p_i-1)p_i^{k_i-1}=\prod _{i=1}^t(1-\frac{1}{pi})p_i^{k_i}=n\prod _{i=1}^t(1-\frac{1}{pi})$$
证明可以在百度百科中搜到，这里略过。证明

$2.$欧拉定理降幂

欧拉定理
$$a^{p-1}\equiv 1(modp)p为质数$$
证明：
引理：$若ma-na\equiv 0(modp)且gcd(a,p)=1，则(m-n)为p的倍数$

那么 $a,2a,3a,...,(p-1)a\%p$的余数必定各不相同，分别为$1,2,3,...,(p-1)$，($a$不为$p$的倍数)
全部相乘，得
$$(p-1)!a^{p-1}\equiv (p-1)!(modp)$$
$$a^{p-1}\equiv 1(modp)$$
广义欧拉定理
$$a^{\phi (p)}\equiv 1(modp)gcd(a,p)=1$$
证明：
设$1-(p-1)$中与$p$互质的数为$x_1,x_2,...,x_n$，那么$x_{1}a,x_{2}a,...,x_{n}a\%p$的余数也各不相同，且为$x_1,x_2,...,x_n$，证明方式与上面同理。

应用:在快速幂中指数过大，如求$a^{b^c}\%p(p为质数)$时$(a,b,c<=10^{12})$就可以

long long quickpow(long long a,long long b,long long Mod)
int main(){
    int a,b,c,p;
    scanf("%lld%lld%lld%lld",&a,&b,&c,&p);
    printf("%lld\n",quickpow(a,quickpow(b,c,p-1),p);
}

$3.Id=\phi \ast I$

$4.\phi =\mu \ast Id$

$3,4$的证明方法见Dirichlet卷积初步

$5.{\sum }_n^{i-1}[gcd(i,n)=1]\times i=\frac{n\times \phi (n)+[n=1]}{2}$

证明：当$n=1$时，易证等式成立

当$n\ne 1$时，因为$gcd(i,n)=gcd(n-i,n)$，所以所有的$i$必成对出现，一共有$\frac{\phi (n)}{2}对$，所以和为$\frac{n\times \phi (n)}{2}$

$6.\phi (ij)=\frac{\phi (i)\phi (j)gcd(i,j)}{\phi (gcd(i,j))}$

证明：由于$\phi$为积性函数，所以我们只需证明满足质数的整数次幂的情况。

当$i=1或j=1$时显然成立。

设$i=p^s,j=p^t$，则
$$\frac{\phi (i)\phi (j)gcd(i,j)}{\phi (gcd(i,j))}$$
$$=\frac{(p-1{)}^2{p}^{s+t+min(s,t)-2}}{(p-1)p^{min(s,t)-1}}$$
$$=(p-1)p^{s+t-1}$$
$$=\phi (p^{s+t})=\phi (ij)$$

二、莫比乌斯函数

首先，要明白它的定义式：

对于一个数n, { ∑ d ∣ n μ ( d ) } = [ n = 1 ] \{\sum_{d|n}\mu(d)\}=[n=1] {∑d∣n​μ(d)}=[n=1]

当$n$不等于$1$时，$n$所有因子的莫比乌斯函数值的$和$为$0$。

易得，$\mu (1)=1$，取$n=p^k$（$p$为质数，$k>1$），由定义式可得：
$$\mu (n)=\mu (1)+\mu (p)+\mu (p^2)+...+\mu (p^{k-1})=0$$
取$k=2$，得$\mu (p)=-1$

然后可以进一步推出对于$\forall k\ge 2，\mu (p^k)=0$

然后我们就得到了一个结论：当$n$的一个质因子的次数大于等于$2$时，$\mu (n)=0$

由于$\mu$函数为积性函数，所以$\mu (ij)=\mu (i)\mu (j)$

当$n$不存在一个质因子的次数大于$2$时，我们设$n={\prod }_{i=1}^t{p}_i$
$$\therefore \mu (n)=\mu (p_1)\times \mu (p_2)\times ...\times \mu (p_t)=(-1{)}^t$$
综上所述
$$\mu (n)=\{\begin{array}{r}0n有平方因子\\ (-1{)}^{t}otherwise\end{array}$$

void sieve(){  //莫比乌斯函数筛法
	mu[1]=1;
	for(int i=2;i<=n;i++){
		if(!flag[i])  prime[++cnt]=i,mu[i]=-1;
		for(int j=1;j<=cnt&&i*prime[j]<=n;j++){
			flag[i*prime[j]]=1;
			if(i%prime[j]==0)  break;
			mu[i*prime[j]]-=mu[i];
		}
	}
}

三、约数个数函数

一个重要的性质

$$d(ij)=\sum _{k|i}\sum _{l|j}[gcd(k,l)=1]$$
证明：

设$q$为$ij$的一个因子(可以是质因子)，设$q=s\times p^t$（$p$为质数）

设$i=i^{\prime }\times p^a,j=j^{\prime }\times p^b$，由于$q$为$ij$的一个因子，必有$t\le a+b$

如果$t\le a$，就取$k=m\times p^t$以保证$gcd(k,l)=1$

否则取$k$满足$gcd(k,p)=1$，$t=n\times p^{t-a}$，此时$gcd(k,l)=1$。

这样可以保证对于$ij$的每一个约数都存在唯一一种分解方案，该性质成立。