莫比乌斯反演！

在平时，我们经常会需要对带 $\gcd$ 的和式进行处理，然而，$\gcd$ 实在难以处理，所以我们一般退而求其次，改为枚举 $\gcd$，然而，这也不是个很好的办法，因为我们无法很快的求出有哪些数对 $(x,y)$ 满足 $\gcd (x,y)=k$，而处理带 $\gcd$ 的和式的一大杀器便是莫比乌斯反演，所以我们来介绍一下。

莫比乌斯函数

首先定义这样一个函数 $\mu (x)$，它满足：
$\mu (1)=1$
当 $x$ 存在平方数因子的时候，即当存在整数 $d$ 满足 $d^2|x$ 时，$\mu (x)=0$;
否则，$x$ 可以写成若干素数的乘积的形式，假设 $x=p_1{p}_2\dots p_n$，则 $\mu (x)=(-1{)}^n$;

其中 $|$ 符号为整除，表示 $x\%d^2=0$;

不难发现 $\mu (x)$ 其实是一个积性函数，我们可以很容易地通过线性筛来求它。

莫比乌斯函数的重要等式

这样定义的 $\mu (x)$ 满足一个非常良好的性质，即：
$$[x=1]=\sum _{d|x}\mu (d)$$

其中 $[p]$ 表示当 $p$ 为真时结果为 $1$，当 $p$ 为假时结果为 $0$。
后边的 $\sum$ 表示求和符号，即所有满足 $d|x$ 的整数 $d$，对 $\mu (d)$ 进行求和。

这个结果可能并不是那么显然，我们将 $x=1$ 带入验证试一下，化简后得到 $1=\mu (1)$，这似乎是没问题的。

当 $x\ne 1$ 时，我们不妨设 $x=p_1^{r_1}{p}_2^{r_2}\dots p_k^{r_k}$，我们发现，只有当 $d$ 为若干素数的乘积的时候才会对整个和式产生贡献，此时 $d=p_1^{r_1^{\prime }}{p}_2^{r_2^{\prime }}\dots p_k^{r_k^{\prime }}$，其中对任意 $1\le i\le k$，$r_i^{\prime }$ 的值为 $1$ 或 $0$。

那么我们不妨枚举 $r_i^{\prime }$ 中 $1$ 的个数，以另一种形式对这个和式进行求和，得到：

$$[x=1]=\sum _{i=0}^k{C}_k^i(-1{)}^i=(1-1{)}^k=0$$

至此，这个等式就得到了证明。

为什么说这个等式重要呢，因为我们可以通过这个等式推导出许多东西。

莫比乌斯反演公式

对于一个函数 $f(n)$，我们可以利用 $\mu (x)$ 来改写它的形式，具体来说，我们有如下等式：
$$f(n)=\sum _{d|n}[d=n]f(d)$$

这是非常显然的，然后我们进一步对其推导：

$$f(n)=\sum _{d|n}[\frac{n}{d}=1]f(d)$$

直接把 $[\frac{n}{d}=1]$ 用 $\mu$ 函数展开，得到：

$$f(n)=\sum _{d|n}\sum _{k|\frac{n}{d}}\mu (k)f(d)$$

交换 $k$ 和 $d$ 的求和顺序，得到：

$$f(n)=\sum _{k|n}\mu (k)\sum _{d|\frac{n}{k}}f(d)$$

如果说，我们设
$$g(n)=\sum _{d|n}f(d)$$

可以得到：

$$f(n)=\sum _{k|n}\mu (k)g(\frac{n}{k})$$

我们不习惯用 $k$ 来枚举因数，改成用 $d$，就可以得到莫比乌斯反演的结论了，即：
设

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

那么便有：

$$f(n)=\sum _{d|n}\mu (d)g(\frac{n}{d})$$

一定成立。

这里关于 $f(n),g(n)$ 还有一个性质：这两个函数中如果有一个函数是积性函数，那么另一个函数也一定是积性函数。

事实上，设

$$h(n)=\sum _{d|n}f(d)g(\frac{n}{d})$$

如果 $f(n),g(n)$ 都是积性函数，那么 $h(n)$ 一定是积性函数，这一点并不难证明：

假设 $a,b$ 互质，那么他们的任意两个因数一定互质，所以我们进行如下推导：

$$h(a)h(b)=\sum _{d|a}f(d)g(\frac{a}{d})\sum _{k|b}f(k)g(\frac{b}{k})$$

$$=\sum _{d|a}\sum _{k|b}f(d)f(k)g(\frac{a}{d})g(\frac{b}{k})$$

$$=\sum _{dk|ab}f(dk)g(\frac{ab}{dk})=h(ab)$$

我们称这样的函数 $h$ 为函数 $f$ 和函数 $g$ 的狄利克雷卷积，记作 $h=f\ast g$，狄利克雷卷积支持交换律结合律，然而一般情况下，它最大的作用就是两个积性函数的卷积一定是一个积性函数，通常用来证明一个函数是否是积性函数。

莫比乌斯反演的作用

说了这么多，莫比乌斯反演到底怎么处理带 $\gcd$ 的和式呢？
$\gcd$ 有一个非常良好的性质，那就是 $d|i$ 且 $d|j$ 是 $d|\gcd (i,j)$ 的充要条件，也就是说，如果我们把 $\gcd (i,j)$ 弄到条件上，就可以轻松去掉这个 $\gcd$ 了。

我们来看一下如下的和式：

$$\sum _{i=1}^n\sum _{j=1}^m\gcd (i,j)$$

我们可以找到一个函数 $f(n)$ 满足
$$n=\sum _{d|n}f(d)$$

利用莫比乌斯反演，得到

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

由狄利克雷卷积的性质我们知道 $f(n)$ 是一个积性函数，可以用线性筛求。
那么有

$$\sum _{i=1}^n\sum _{j=1}^m\sum _{d|\gcd (i,j)}f(d)$$

这样的话我们就可以把 $d$ 提前，首先枚举 $d$，得到：

$$\sum _{d=1}^{n}f(d)\sum _{i=1}^n[d|i]\sum _{j=1}^m[d|j]$$

$$=\sum _{d=1}^{n}f(d)\lfloor \frac{n}{d}\rfloor \lfloor \frac{m}{d}\rfloor$$

线性筛处理后就可以利用数论分块 $O(\sqrt{n})$ 内求出这个和式了。

事实上，这里 $f(n)=\phi (n)$，这并不难判断，对于积性函数，只要证明这两个函数在 $n=p^k$ 时相等就可以了，其中 $p$ 是质数。

另一个反演法

有时候，我们会考虑枚举最大公约数，而枚举最大公约数是不好处理的，因为你无法很快地找到所有满足 $\gcd (i,j)=k$ 的数对，但是，如果我们枚举的只是一般的约数，那么找这个数对就会变得非常简单，这个时候，我们不妨考虑如下的反演：

对于一个函数 $f(r)$，我们有
$$f(r)=\sum _{i=1}^{\lfloor \frac{n}{r}\rfloor }f(ir)[i=1]$$

我们沿用之前的套路对这个和式进行处理，希望能够找到一个反演公式。

$$f(r)=\sum _{i=1}^{\lfloor \frac{n}{r}\rfloor }f(ir)\sum _{d|i}\mu (d)$$

同样进行交换求和变量得到：

$$f(r)=\sum _{d=1}^{\lfloor \frac{n}{r}\rfloor }\mu (d)\sum _{i=1}^{\lfloor \frac{n}{r}\rfloor }[d|i]f(ir)$$

我们用 $id$ 去替换以前的 $i$，得到：

$$=\sum _{d=1}^{\lfloor \frac{n}{r}\rfloor }\mu (d)\sum _{i=1}^{\lfloor \frac{n}{rd}\rfloor }f(ird)$$

看到这里，我们知道我们的目的达成了，我们设：

$$g(r)=\sum _{i=1}^{\lfloor \frac{n}{r}\rfloor }f(ir)$$

那么一定有：

$$f(r)=\sum _{i=1}^{\lfloor \frac{n}{r}\rfloor }\mu (i)g(ir)$$

我们就得到了另一个反演公式，这个公式在很多时候是非常有用的，我们如果把 $f(r)$ 当成对于所有最大公约数为 $r$ 的数对的函数，那么 $g(r)$ 就代表对于所有公约数为 $r$ 的数对的函数（显而易见，即最大公约数为 $r$ 的倍数），很多时候，$f(r)$ 并不是那么好求，但是 $g(r)$ 很好求，这样我们就可以通过求出 $g(r)$，然后用这个反演公式来得到 $f(r)$。

很多时候，这个反演公式比上面那个反演公式更为直观好用。

另一个反演法的应用

我们可以用这个反演法去解一下上边那个和式，来体会一下这两个反演法的区别。

$$\sum _{i=1}^n\sum _{j=1}^m\gcd (i,j)$$

我们不妨直接枚举最大公约数，假设当前最大公约数为 $r$，那么上式可以转化为：

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

我们设

$$f(r)=\sum _{i=1}^n\sum _{j=1}^m[\gcd (i,j)=r]$$

这个函数并不是很好求，但是我们可以设：

$$g(r)=\sum _{i=1}^{\lfloor \frac{n}{r}\rfloor }f(ir)$$

我们考虑 $f(r)$ 与 $g(r)$ 的实际意义：

$f(r)$ 代表所有满足 $1\le i\le n,1\le j\le m$ 的数对 $(i,j)$ 中最大公约数为 $r$ 的数量。
那么 $g(r)$ 就是代表所有满足 $1\le i\le n,1\le j\le m$ 的数对 $(i,j)$ 中公约数为 $r$ 的数量。

公约数为 $r$ 的数对个数非常好求，只要是满足 $r|i$ 且 $r|j$ 就可以，满足 $r|i$ 的 $i$ 有 $\lfloor \frac{n}{r}\rfloor$ 个，满足 $r|j$ 的 $j$ 有 $\lfloor \frac{m}{r}\rfloor$，所以 $g(r)=\lfloor \frac{n}{r}\rfloor \lfloor \frac{m}{r}\rfloor$，我们借此可以求出：

$$f(r)=\sum _{i=1}^{\lfloor \frac{n}{r}\rfloor }\mu (i)\lfloor \frac{n}{r}\rfloor \lfloor \frac{m}{r}\rfloor$$

直接将 $f(r)$ 代入最开始的式子，得到：

$$\sum _{r=1}^{n}r\sum _{i=1}^{\lfloor \frac{n}{r}\rfloor }\mu (i)\lfloor \frac{n}{r}\rfloor \lfloor \frac{m}{r}\rfloor$$

$$=\sum _{r=1}^{n}r\lfloor \frac{n}{r}\rfloor \lfloor \frac{m}{r}\rfloor \sum _{i=1}^{\lfloor \frac{n}{r}\rfloor }\mu (i)$$

处理一下 $\mu (i)$ 的前缀和后依然可以进行 $O(\sqrt{n})$ 的数论分块。

后记

两种反演法有着不同的效果，但是通常情况下，第二种可能会更好用一点，因为这种反演法允许你无脑枚举 $\gcd$，然后通过这种反演法转化为枚举一般因数，不过对于具体的题目还是要具体分析。