【算法&题解】莫比乌斯反演及一些经典套路

首先声明，由于我太菜了，本文很多式子可能并不会有详细的证明。若要了解莫比乌斯反演最基本的证明，可以去这里，本文主要讲解一些套路。

莫比乌斯函数

我们定义莫比乌斯函数如下：
$$\begin{gathered} 设n=\prod _{i=1}^m{p}_i^{a_i},k=\prod _{i=1}^m{a}_i \\ \mu (n)=\{\begin{array}{ll}1& (n=1)\\ (-1{)}^m& (k=1)\\ 0& (k>1)\end{array} \end{gathered}$$
然后我们会发现它是一个积性函数，可以用线性筛$O(n)$求出$1-n$的莫比乌斯函数。代码如下：

mu[1] = 1;
for(int i = 2; i <= n; i++){
   
    if(!mark[i]){
   
        p[++cnt] = i;
        mu[i] = -1;		//质数显然k=1
    }
    for(int j = 1; j <= cnt && p[j]*i <= n; j++){
   
        mark[p[j]*i] = 1;
        if(i%p[j] == 0) break;
        mu[i*p[j]] = -mu[i];	//由于欧拉筛中，每个数只会被最小的质因子筛到一次，相当于这个数多了一个因子p[j]，就相当于-mu[i](可以通过分类讨论理解)。
    }
}

同时，莫比乌斯函数还满足：
$$\sum _{i=1}^n\mu (i)=[n=1]$$
这也是一个美妙的性质（可惜我不会证）。

---

莫比乌斯反演

依然直接扔公式（是的，我也不会证）：

若函数$f(n)$和$g(n)$是数论函数，且满足：
$$g(n)=\sum _{d|n}f(d)$$
那么有如下结论：
$$f(n)=\sum _{d|n}g(d)\times \mu (\frac{n}{d})$$
实际上，在做题过程中也会遇到上述式子中的$d|n$变为$n|d$，式子依然成立，但此时$d$往往有明确的上限。

---

经典套路

以下默认$n\le m$。

1.求(i,j)=x的个数

这是莫比乌斯反演最经典的应用。

给定$n,m$，求$f(x)=\sum _{i=1}^n\sum _{j=1}^m[(i,j)=x]$。

我们设$g(x)=\sum _{i=1}^n\sum _{j=1}^m[x|(i,j)]$。那么我们发现显然$g(x)=\sum _{x|d}f(d)$。此处$d$的上限为$n$。

于是我们套公式，直接反演一波：
$$f(x)=\sum _{x|d}g(d)\times \mu (\frac{d}{x})。$$
我们为啥这么麻烦的反演呢？当然是因为$g(x)$比较好求。我们观察$x|(i,j)$这个条件，发现$i,j$均可以表示成$x\times k$，于是在$n$和