【持续更新】莫比乌斯反演简明教程

前言

开始学省选算法了……

感觉莫比乌斯反演好厉害的样子，就先学习一下

一入反演深似海……

相关的东西太多了，以后会不定期更新

前置技能

莫比乌斯函数

莫比乌斯函数$\mu(n)$的定义如下：

设$n=p_1^{k_1}\cdot p_2^{k_2} \dots p_m^{k_m}$

$$\begin{equation} \mu(n)= \begin{cases} 1 & n=1 \\ (-1)^m & \prod_{i=1}^m k_i=1 \\ 0 & \text{otherwise}(k_i\gt 1) \\ \end{cases} \end{equation}$$
显然， $\mu(n)$是积性函数： $\mu(x)\mu(y)=\mu(xy),x\perp y$

那么我们就可以使用线性筛来得到$\mu(n)$：

void prepare(){
    mu[1]=1;
    for (int i=2;i<=N;i++){
        if (!vis[i]) p[++p[0]]=i,mu[i]=-1;
        for (int j=1;j<=p[0]&&i*p[j]<=N;j++){
            vis[i*p[j]]=1;
            if (i%p[j]==0) {mu[i*p[j]]=0;break;}
             else mu[i*p[j]]=-mu[i];
        }
    }
}

狄利克雷卷积

定义：对于数论函数$f(n)$和$g(n)$，定义卷积运算$*$为：

$$(f*g)(n)=\sum_{d|n} f(d)\cdot g\big ( \frac n d \big )$$
狄利克雷卷积的单位元为$e(x)=[x=1]$

任何数论函数$f$满足$f*e=f$

定义函数$I(x)=1$，$id(x)=x$

则有：$I*\mu =e$

这说明$\mu$是$I$的逆元

又有：$\varphi*I= id$

那么可以得到$\varphi =id*\mu$

莫比乌斯反演

对于$f(n),g(n)$满足：

$$f(n)=\sum_{d|n} g(d)$$
则有莫比乌斯反演：
$$g(n)=\sum_{d|n}f(d)\mu \big(\frac n d \big)$$
其实非常显然，把命题换成狄利克雷卷积的形式就是：

已知$f=g*I$，求证$g=f*\mu$

根据$I*\mu=e$直接证明了

应用

求gcd=k的个数

问题：求

$$\sum_{i=1}^n \sum_{j=1}^m [gcd(i,j)=k]$$
这是莫比乌斯反演的入门题，非常经典

推导：

$$\sum _{i=1}^n \sum _{j=1}^m [gcd(i,j)=k] \\ \Rightarrow \sum _{i=1}^{\lfloor \frac n k \rfloor} \sum _{j=1}^{\lfloor \frac m k \rfloor} [gcd(i,j)=1] \\ $$
然后套用莫比乌斯反演：
$$\Rightarrow \sum _{i=1}^{\lfloor \frac n k \rfloor} \sum _{j=1}^{\lfloor \frac m k \rfloor} \sum_{d|(i,j)} \mu(d) \\ \Rightarrow \sum_d \mu(d) \sum _{d|i} \sum _{d|j} 1 \\ \Rightarrow \sum_d \mu(d) \lfloor \frac n {kd} \rfloor \lfloor \frac m {kd} \rfloor$$
然后就可以对 $\lfloor \frac n {kd} \rfloor \lfloor \frac m {kd} \rfloor$分块求和了

然后就可以在$O(\sqrt n)$时间里解决每次询问（$O(n)$预处理）

例题

YY的GCD

BZOJ2820

先考虑枚举质数p，答案就是：

$$\sum_p \sum_d \mu(d) \lfloor \frac n {pd} \rfloor \lfloor \frac m {pd} \rfloor$$
设 $T=pd$，考虑枚举 $T$，则有：
$$\sum_T^{min\{ n,m \}} \lfloor \frac n T \rfloor \lfloor \frac m T \rfloor \sum_{p|T} \mu\Big(\frac T p\Big)$$
如果能够预处理 $\sum_{p|T} \mu\Big(\frac T p\Big)$关于 $T$的前缀和，前面的柿子就可以$O(\sqrt n)$求解

其实可以暴枚$p$

因为均摊每个质数是$ln(n)$的，而质数的个数是$\frac n {ln(n)}$个

所以预处理可以$O(n)$