积性函数笔记

莫比乌斯反演的前置知识

定义

如果一个数论函数$f$满足：当$n\perp m$时，有

$$f(nm)=f(n)f(m)$$

则称其为积性函数

当不需要满足互质条件$(n\perp m)$就能满足$f(nm)=f(n)f(m)$时，称其为完全积性函数

若$f$是积性函数，且$n=p_1^{{\alpha }_1}{p}_2^{{\alpha }_2}\dots p_s^{{\alpha }_s}$是$n$的标准分解，则有

$$f(n)=f(p_1^{{\alpha }_1})f(p_2^{{\alpha }_2})\dots f(p_s^{{\alpha }_s})$$

因此研究积性函数$f$可以转化为研究$f(p^{\alpha })$，即$f$在素数和素数的幂上的取值。

积性函数求值

设$f$是积性函数，为求$f(n)$，可以对$n$分解素因子，然后计算所有的$f(p^{\alpha })$乘起来，因此积性函数求值都可以用线性筛来求

如果要对$1$到$n$之间的所有数求出$f$，注意到$\text{Euler}$筛法的过程中可以求出每个数的最小素因子和最小素因子的幂次，利用此就能在线性时间内计算出所需的$f$的值

欧拉筛如何求出最小素因子的幂次

void sieve() {
    ip[0] = ip[1] = 1;
    alpha[1] = 0;
    for (int i = 2; i < XR; i++) {
        if (!ip[i]) p[c++] = i, alpha[i] = 1, p1[i] = i;
        for (int j = 0; j < c && i * p[j] < XR; j++) {
            ip[i * p[j]] = 1;
            p1[i * p[j]] = p[j];
            if (i % p[j]) alpha[i * p[j]] = 1;
            else {
                alpha[i * p[j]] = alpha[i] + 1;
                break;
            }
        }
        //最小素因子：p1[i] ^ alpha[i]
    }
}

常见积性函数

单位函数

单位函数$ϵ(n)$定义为

$$ϵ(n)=[n=1]=\{\begin{array}{rl}1,& n=1\\ 0,& n\ne 1\end{array}$$

其中$[\text{condition}]$表示当$\text{condition}$为真时取值为$1$，否则为$0$的函数。

除数函数

除数函数${\sigma }_k(n)$用来表示$n$的因子的$k$次方之和：

$${\sigma }_k(n)=\sum _{d|n}{d}^k$$

约数个数${\sigma }_0(n)$常记为$d(n)$，约数和${\sigma }_1(n)$常记为$\sigma (n)$。

可以证明除数函数都是积性函数。

证明

每个素因子对函数的贡献都是独立的，它们之间的贡献可以用乘积来衡量，举一个例子：要算$n$的约数个数（即$d(n)$），就要枚举$n$的约数$d$，考虑$n$的标准分解，也就是$n=p_1^{{\alpha }_1}{p}_2^{{\alpha }_2}\dots p_s^{{\alpha }_s}$,可以想象出$d$的素因子也跑不出这些个素因子，只不过上面的指数不一样，假设是$\beta$，那么也就是$d=p_1^{{\beta }_1}{p}_2^{{\beta }_2}\dots p_s^{{\beta }_s}$，${\beta }_1$显然是要在$0$到${\alpha }_1$之间，否则$d$就不是$n$的约数了，同理每个$\beta$都是对应的每个$0$到$\alpha$之间的，那么${\beta }_1$有$0$到${\alpha }_1$共${\alpha }_1+1$种取值，每个$\beta$的选取是独立的，所以他们的贡献的乘积就是$(1+{\alpha }_1)\ast (1+{\alpha }_2)\ast \dots \ast (1+{\alpha }_s)$，这就是${\sigma }_0$的算法，其他的${\sigma }_k$是一样的，只不过$(1+\alpha )$这一步要换成一个求和的形式。

素因子贡献独立是一个常见的思路。

Euler函数（欧拉函数）

$\text{Euler}$函数$\phi (n)$表示不超过$n$且与$n$互质的正整数的个数，即

$$\phi (n)=\sum _{i=1}^n[\gcd (i,n)=1]$$

比如说$\phi (6)=2$，两个数是$1,5$；$\phi (10)=4$，四个数分别为$1,3,5,7$。

由$n$的标准分解并结合容斥原理，我们可以得到欧拉函数的显示表达式：

$$\phi (n)=n\cdot \prod _{i=1}^s(1-\frac{1}{p_i})$$

其中$p_i$是质数，$n=p_1^{{\alpha }_1}{p}_2^{{\alpha }_2}\dots p_s^{{\alpha }_s}$是$n$的标准分解。

性质

对于任$n$，$\text{Euler}$函数具有如下性质：

$$n=\sum _{d|n}\phi (d)$$

要证明这个等式，我们将$1$到$n$中的所有整数按与$n$的最大公约数分类。

若$gcd(n,i)=d$，那么$gcd(\frac{n}{d},\frac{d}{i})=1$。而又$\frac{i}{d}$是不超过$\frac{n}{d}$的整数，故这样的$i$有$\phi (\frac{n}{d})$个。

考虑所有$d|n$，我们就考虑到了所有$1$到$n$之间的$n$个整数，因此有

$$n=\sum _{d|n}\phi (\frac{n}{d})=\sum _{d|n}\phi (d)$$

证明积性

以下是从大佬那里淘来的证明

同样的，$t\perp nm\iff t\perp n,t\perp m\iff (t\mathrm{m}\mathrm{o}\mathrm{d}n)\perp n,(t\mathrm{m}\mathrm{o}\mathrm{d}m)\perp m$，所以每个 $[1,nm]$之间的与$nm$互质的数$t$都可以对应到一个$[1,n]$的与$n$互质的数 $t\mathrm{m}\mathrm{o}\mathrm{d}n$和一个$[1,m]$的与$m$互质的数$t\mathrm{m}\mathrm{o}\mathrm{d}m$。

并且根据中国剩余定理，这种对应是一一对应的（即已知 $a\perp n,b\perp m$后可以唯一确定一个$[1,nm]$之间的$t$使得$t\mathrm{m}\mathrm{o}\mathrm{d}n=a,t\mathrm{m}\mathrm{o}\mathrm{d}m=b$，且$t\perp nm$）。因此 $\phi (nm)=\phi (n)\phi (m)$ 。

---

然而我看不懂……还是从定义上来证明吧！

假设有两个互质的正整数$n,m$，则

$$\phi (n)=n\prod (1-\frac{1}{p_i})$$

$$\phi (m)=m\prod (1-\frac{1}{p_{i^{\prime }}})$$

$$\phi (n)\phi (m)=n\prod (1-\frac{1}{p_i})m\prod (1-\frac{1}{p_{i^{\prime }}})=nm\prod (1-\frac{1}{p_i})\prod (1-\frac{1}{p_{i^{\prime }}})$$

因为$n,m$互质，所以$p_i$和$p_{i^{\prime }}$各各都不相同，且都是$nm$的质因子

因此就可以推出$\phi (nm)=\phi (n)\phi (m)$

至此，积性函数的性质得证。但是由上面的证明可知，$n,m$必须要互质才可以满足欧拉函数是积性函数，由此可见欧拉函数不是完全积性函数

恒等函数

${\mathbf{i}\mathbf{d}}^k=n^k$,$\mathbf{i}\mathbf{d}(n)=n$

（就酱紫一笔带过吧……）

性质

若$\mathbf{f}(x)$和$g(x)$均为积性函数，则以下函数也为积性函数

$$\begin{array}{rl}h(x)& =f(x^p)\\ h(x)& =f^p(x)\\ h(x)& =f(x)g(x)\\ h(x)& =\sum _{d\mid x}f(d)g(\frac{x}{d})\end{array}$$

部分性质的证明

来自铃悬的数学讲堂（抄这个东西也很累……）

两个积性函数的狄利克雷卷积是积性函数

考虑到上面的性质，即“若$n\perp m$则每个$nm$的约数都可以分解成一个$n$的约数和一个$m$的约数的积”，并且有另一个性质：若$n\perp m,a\mid n,b\mid m$则$a\perp b$。

于是若$n\perp m$，我们就有

$$\begin{array}{rl}t(nm)& =\sum _{d\mid nm}f(d)g\left(\frac{nm}{d}\right)\\ & =\sum _{a\mid n,b\mid m}f(ab)g\left(\frac{nm}{ab}\right)\\ & =\sum _{a\mid n,b\mid m}f(a)f(b)g\left(\frac{n}{a}\right)g\left(\frac{m}{b}\right)\\ & =\left(\sum _{a\mid n}f(a)g\left(\frac{n}{a}\right)\right)\left(\sum _{b\mid m}f(b)g\left(\frac{m}{b}\right)\right)\\ & =t(n)t(m)\end{array}$$

积性函数的逆是积性函数

第二个重要的结论：积性函数的逆是积性函数。
对于一个积性函数 $f$ ，如何证明其逆 $g(n)=[n=1]-{\sum }_{i\mid n,i\ne 1}f(i)g\left(\frac{n}{i}\right)$ （注意，积性函数一定满足 $f(1)=1$ ，因为 $f(1)=f(1)f(1)$ ，并且如果 $f(1)=0$ 则 $f(n)\equiv 0$ ，这种情况我们不考虑）也满足积性呢？
对 $nm$ 的大小进行归纳：

$nm=1$ 时， $g(1)=1$ ，结论显然成立；

假设 $nm>1$ ，当 $n^{\prime }{m}^{\prime }<nm$ 的时候有 $g(n^{\prime }{m}^{\prime })=g(n^{\prime })g(m^{\prime })$ ，
$$\begin{array}{rl}g(nm)& =-\sum _{d\mid nm,d\ne 1}f(d)g\left(\frac{nm}{d}\right)\\ & =-\sum _{a\mid n,b\mid m,ab\ne 1}f(ab)g\left(\frac{nm}{ab}\right)\\ & =-\sum _{a\mid n,b\mid m,ab\ne 1}f(a)f(b)g\left(\frac{n}{a}\right)g\left(\frac{m}{b}\right)\\ & =f(1)f(1)g(n)g(m)-\sum _{a\mid n,b\mid m}f(a)f(b)g\left(\frac{n}{a}\right)g\left(\frac{m}{b}\right)\\ & =g(n)g(m)-\left(\sum _{a\mid n}f(a)g\left(\frac{n}{a}\right)\right)\left(\sum _{b\mid m}f(b)g\left(\frac{m}{b}\right)\right)\\ & =g(n)g(m)-ϵ(n)ϵ(m)\\ & =g(n)g(m)\end{array}$$

注意前面几步把 $g(\frac{nm}{ab})$ 拆成 $g(n/a)g(m/b)$ 的时候 $\frac{nm}{ab}<nm$ ，可以运用归纳条件。
最后一步是因为 $nm>1$ 所以 $n,m$ 之间至少一个不为 $1$ ，则 $ϵ(n)ϵ(m)=[n=1][m=1]=0$ 。

就酱紫吧