莫比乌斯反演学习笔记

为啥学习莫比乌斯反演？是因为最近 CF 连着两次莫反让我感觉到自己没有力量。

大概是复读一下题解的话再加上一点自己的理解。

前置知识和引理：

数论分块：

不想写了，反正都2021年了大家肯定都会了。

二项式定理：

$(a+b{)}^n=\sum _{i=0}^n{C}_n^i{a}^{n-i}{b}^i$。

引理 1：

$$\forall a,b,c\in z\lfloor \frac{a}{bc}\rfloor =\lfloor \frac{\lfloor \frac{a}{b}\rfloor }{c}\rfloor$$

证明：$\lfloor \frac{a}{bc}\rfloor =\lfloor \frac{a}{b}\times \frac{1}{c}\rfloor$。

令 $\lfloor \frac{a}{b}\rfloor =\frac{a}{b}-d$。

$\lfloor \frac{a}{b}\times \frac{1}{c}\rfloor =\lfloor (\lfloor \frac{a}{b}\rfloor +d)\times \frac{1}{c}\rfloor =\lfloor \frac{\lfloor \frac{a}{b}\rfloor }{c}+\frac{d}{c}\rfloor =\lfloor \frac{\lfloor \frac{a}{b}\rfloor }{c}\rfloor$。

$Q.E.D$

引理 2：

对于给定的 $N\in {\mathbb{N}}_{+}$，对于任意的 $d\in {\mathbb{N}}_{+}$，$\lfloor \frac{N}{d}\rfloor$ 最多有 $\lfloor 2\sqrt{N}\rfloor$ 种取值。

证明：当 $d\le \sqrt{N}$ 时，$\lfloor \frac{N}{d}\rfloor$ 显然最多有 $\sqrt{N}$ 种取值。

当 $d>\sqrt{N}$，易得最多有 $\sqrt{N}$ 种取值。

$Q.E.D$

引理 2 其实就是数论分块复杂度的证明。

积性函数

定义

若函数 $f(n)$ 满足 $f(1)=1$，$\forall x,y\in {\mathbb{N}}_{+}(gcd(x,y)=1),f(xy)=f(x)f(y)$，则称其为积性函数。

e.g：

欧拉函数 $\varphi (x)=\sum _{i=1}^n[gcd(i,x)=1]$。

除数函数 ${\sigma }_k(x)=\sum _{d|n}{d}^k$。特别的我们记 ${\sigma }_0(x)=\tau (x)$ 或者 $d(x)$。

若函数 $f(n)$ 满足 $f(1)=1$，$\forall x,y\in {\mathbb{N}}_{+},f(xy)=f(x)f(y)$，则称其为完全积性函数。

备注：其实 $f(1)=1$ 是一个多余的东西，因为 $f(1)=f(1)f(1)$

e.g：

常数函数 $I(n)=1$。

狄利克雷卷积

对于 $f(x),g(x)$ 两个函数（注意不一定是积性函数）它们的狄利克雷卷积写作 $f\ast g$，这个 $f\ast g$ 就是卷出来的新函数的名称。

$(f\ast g)(n)=\sum _{d|n}f(d)g(\frac{x}{d})=\sum _{xy=n}f(x)g(y)$

那么有

交换律：$f\ast g=g\ast f$

结合律：$f\ast (g\ast h)=(f\ast g)\ast h$

然后再介绍一个函数 $ϵ(x)=[x=1]$，这个函数有性质：$ϵ\ast f=f$，可以看做卷积中的单位元（所以又叫做元函数）。

积性函数的性质

若 $f1(x)$ 和 $f2(x)$ 都是积性函数，那么以下函数也都是积性函数：

$h(x)=f1(x^p)$

$h(x)=f{1}^p(x)=f{1}^{p-1}\ast f1$

$h(x)=f1(x)\ast f2(x)$

$h(x)=\sum _{d|n}f1(d)f2(\frac{x}{d})$

$f(x)=f({p_1}^{k_1})\ast f({p_2}^{k_2})\ast ...\ast f({p_m}^{k_m})(x={p_1}^{k_1}{p_2}^{k_2}...{p_m}^{k_m},isprime(p_i)=1)$。

积性函数的逆也是积性函数（什么是函数的逆，对于函数 $f$，若存在函数 $g$ 使得 $f\ast g=ϵ$，就说 $g$ 是 $f$ 的逆（是不是有点像逆元的定义））

证明：不存在的

关于函数的逆，还有一点东西：

• 如何去构造函数的逆：先给出结论，只要 $f(1)\ne 0$，那么 $f$ 存在逆 $g$。

首先，$(f\ast g)(n)=\sum _{d|n}f(d)g(\frac{n}{d})$，把 $d=1$ 拆出来，得 $\sum _{d|n,d\ne 1}f(d)g(\frac{n}{d})+f(1)g(n)$。

然后，我们递归构造 $g(n)=\frac{\left(e-\sum _{d|n,d\ne 1}f(d)g(\frac{n}{d})\right)}{f(1)}$。

带入上式得到 $\sum _{d|n,d\ne 1}f(d)g(\frac{n}{d})+f(1)\frac{\left(e-\sum _{d|n,d\ne 1}f(d)g(\frac{n}{d})\right)}{f(1)}=e$。

即构造出了 $f$ 的逆 $g$。但是注意到 $f(1)\ne 0$。

从某种意义上前置知识已经结束了

---

莫比乌斯反演

莫比乌斯函数

先随便写点东西。

现在有这么两个函数 $F$ 和 $f$。

假定有：

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

如果你理解了卷积的概念：不难得到：

$F(n)=\sum _{d|n}(1\ast f(d))$
$=\sum _{d|n}(I(\frac{n}{d})\ast f(d))$，其中 $I(d)=1$
$=I\ast f$。

于是我们把 $F$ 函数写成了卷积的形式。

有什么好处呢？

如果知道 $f$ 函数，我们就可以得到 $F$ 函数。

那么…如果我们知道 $F$ 函数呢？

那我们就要写成 $f=...$ 的形式，两边同乘 $I^{-1}$，可得
$$f=F\ast I^{-1}$$

所以 $I^{-1}$ 是个什么东西？？？这个时候大佬莫比乌斯站出来说，$I^{-1}=\mu$ 。
大家知道 $I$ 是个积性函数，$I\ast I^{-1}=ϵ$，所以 $I^{-1}$ 也是个积性函数。

研究一个积性函数,先研究其在质数的幂时的表现。——command_block

$\mu (p^0)=1$（显然）

$\mu (p^1)$ 怎么计算呢？注意到 $ϵ(p)=(\mu \ast I)(p)$

$ϵ(p)=\sum _{d|p}\mu (d)I(\frac{p}{d})$
$=\sum _{d|p}\mu (d)$

好像原地踏步了…?注意到 $p$ 是个质数。

$ϵ(p)=\mu (1)+\mu (p)$

由于 $p>1$

$0=1+\mu (p)$

得到 $\mu (p)=-1$

接着我们来看 $\mu (p^k)$ 怎么算。

还是套用刚才的卷积形式：

$ϵ(p^k)=\sum _{d|p}\mu (d)I(\frac{p^k}{d})$

再次去考虑 $d$ 的取值，会发现 $d=p^0,p^1...p^k$ 这些东西。

$ϵ(p^k)=\mu (p^0)I(p^k)+\mu (p^1)I(p^{k-1})...+\mu (p^k)I(p^0)$

$ϵ(p^k)=\sum _{i=0}^k\mu (p^i)=0$

那么可以得到：$\mu (p^k)=0(k>1)$。

当我们研究完了这个东西，注意到积性函数的一个性质：

$f(x)=f({p_1}^{k_1})\ast f({p_2}^{k_2})\ast ...\ast f({p_m}^{k_m})$。

那就不难想到如果有一个 $k_i>1$，$\mu (x)$ 就等于 $0$，而其他情况 $\mu (x)=\pm 1$。具体的：

$$\mu (n)=\{\begin{array}{ll}0& d\in {\mathbb{N}}_{+},d\ne 1,d^2|n\\ (-1{)}^{\sum _{d|n}(isprime(d)=1)}& \text{other}\end{array}$$

说人话就是若 $n$ 存在一个大于 $1$ 的平方因子，那么 $\mu (n)=0$，否则就是 $-1$ 的 $k$ 次方，其中 $k$ 表示 $n$ 的质因数种类数。

既然知道了这些，作为一个臭名昭著著名的积性函数，它和欧拉函数一样可以筛出来：

void Mobius()
{
  e[1]=1;mu[1]=1;
  for (int i=2;i<=n;i++){
    if (!e[i]){p[++tn]=i;mu[i]=-1;}
    for (int j=1;j<=tn;j++){
      if (p[j]*i>n)break;
      mu[p[j]*i]=i%p[j]==0 ? 0 : -mu[i];
      e[p[j]*i]=1;
      if (i%p[j]==0)break;
    }
  }
}

（转自 command_block，有时间再补自己的）

莫比乌斯反演

终于…

嵌入式莫比乌斯反演

由 $\mu \ast I=ϵ$

得 $\sum _{d|n}\mu (d)=[n=1]$

然后：

$[n|m][n/m=1]=[n=m]$

$[n|m]\sum _{d|(n/m)}\mu (d)=[n=m]$

这个东西有什么用呢？来看这样一道题：

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

我们把 $gcd(i,j)$ 看成一个东西。

$\sum _{d|gcd(i,j)}\mu (d)=[gcd(i,j)=1]$

考虑什么是 $d|gcd(i,j)$，也就是说 $d|i$ 且 $d|j$。

所以题目所求变为了$\sum _{i=1}^n\sum _{j=1}^m\sum _{d|i,d|j}\mu (d)$。

考虑变换求和顺序，注意到 $\sum _{d|gcd(i,j)}\mu (d)=ϵ\le 1$，所以可以改成 $\sum _{d=1}^{d=n}\mu (d)\times \sum _{i=1,d|i}^m\sum _{j=1,d|j}^{n}1$。

是不是逐渐可做了？注意到每个 $\lfloor \frac{n}{i}\rfloor$ 都只有一个满足条件（$j$ 同理），这个式子等价于$$\sum _{d=1}^{d=n}\mu (d)\times [n/d][m/d]$$

这个向下取整的形式很类似于整除分块，我们看看怎么用整除分块维护，若确定了一段 $d$ 后面两项的取值是相同的，可以 $(\sqrt{n})$ 计算，但是 $\mu (i)$ 怎么办？记录个前缀和就好了嘛！

也许这道题并不能看出数论分块的优秀复杂度，但是后面就能看出来了。

例题：

洛谷 P1390 公约数的和

题意：求 $\sum _{i=1}^n\sum _{j=i+1}^{m}gcd(i,j)$

这里介绍一个莫凡常用技巧：枚举 gcd，我们设 $gcd(i,j)=g$

当前的问题转化成求：$\sum _{i=1}^n\sum _{j=i+1}^m(gcd(i,j)=d)$

我们先考虑怎么计算 $\sum _{i=1}^n\sum _{j=1}^m(gcd(i,j)=g)$

这个东西是个套路，把这东西换成是 $i,j$ 互质，这样 $gcd(gi,gj)$ 就可以一定是 $g$。原式改成 $\sum _{i=1}^{\lfloor \frac{n}{g}\rfloor }\sum _{j=1}^{\lfloor \frac{m}{g}\rfloor }(gcd(i,j)=1)$这东西大家都很熟了，直接改成重复上面的过程：

$\sum _{i=1}^{\lfloor \frac{n}{g}\rfloor }\sum _{j=1}^{\lfloor \frac{m}{g}\rfloor }\sum _{d|i,d|j}\mu (d)$

$\sum _{d=1}^{min(\lfloor \frac{n}{d}\rfloor ,\lfloor \frac{m}{d}\rfloor )}\mu (d)[\lfloor \frac{\lfloor \frac{n}{g}\rfloor }{d}\rfloor ][\lfloor \frac{\lfloor \frac{m}{g}\rfloor }{d}\rfloor ]$

然后还要删去多算的东西，也就是：

$$\sum _{i=1}^n\sum _{j=1}^{max(i,\lfloor \frac{m}{g}\rfloor )}gcd(i,j)=1$$

同样反演。

然后数论分块一下，预处理一下前缀和就好了（