初探莫比乌斯反演

前言

作为一只数论只会$gcd$的蒟蒻，强行学习了欧拉函数和莫比乌斯反演后，已经暴毙而亡。马上就要$AFO$了，“死”前还是留点遗产吧。

---

简单的定义

• $\sum :$求和
• $\prod :$ 乘积
• $[x]:$ $x$为一个逻辑表达式，若$x$的结果为真，则值为$1$，反之为$0$
• $\lfloor x\rfloor :$ 表示不超过 $x$ 的最大整数。

---

莫比乌斯函数

• $\mu (n)=\{\begin{array}{cc}1& n=1\\ (-1{)}^k& n=\prod _{i=1}^m{p}_i\\ 0& otherwise\end{array}$

---

莫比乌斯反演

（注：以下例题基本出自 $WR$ 的数论博客中，略有删改。）

一个很有用的结论

$\sum _{d\mid n}\mu (d)=[n==1]$
$Proof:$
当 $n=1$ 时，显然成立。
当 $n>1$ 时，设 $n=\prod _{i=1}^m{p}_i^{c_i}$，
若 $d$ 存在一个质因子的次数 $>1$，则 $\mu (d)=0$，
否则，设 $d=\prod _{i=1}^k{p}_i^{\prime }$，则 $\mu (d)=(-1{)}^k$
$\therefore \sum _{d\mid n}=\mu (1)+\sum _{i=1}^m(-1{)}^i\ast (\genfrac{}{}{0ex}{}{m}{i})=\sum _{i=0}^m(-1{)}^i\ast 1^{m-i}\ast \left(\genfrac{}{}{0ex}{}{m}{i}\right)$
根据二项式定理，上式 $=0$。
证毕。
该结论是莫比乌斯反演的核心。

---

以下内容默认 $n\le m$ ，不满足则交换。

$e.g.1$

求 $\sum _{i=1}^n\sum _{j=1}^m[gcd(i,j)==1]$
$Solution:$
根据上面的结论可得：
$\begin{array}{l}\sum _{i=1}^n\sum _{j=1}^m[gcd(i,j)==1]\\ =\sum _{i=1}^n\sum _{j=1}^m\sum _{d\mid gcd(i,j)}\mu (d)\\ =\sum _{d=1}^n\mu (d)\sum _{i=1}^{\lfloor \frac{n}{d}\rfloor }\sum _{j=1}^{\lfloor \frac{m}{d}\rfloor }1\\ =\sum _{d=1}^n\mu (d)\ast \lfloor \frac{n}{d}\rfloor \ast \lfloor \frac{m}{d}\rfloor \end{array}$
对第二行到第三行做出解释：
$d\mid gcd(i,j)\iff d\mid i\wedge d\mid j$
$\therefore$ 满足 $d\mid gcd(i,j)$ 的有序数对 $(i,j)$ 可表示为 $(xd,yd)(x\in [1,\lfloor \frac{n}{d}\rfloor ],y\in [1,\lfloor \frac{m}{d}\rfloor ])$
直接暴力做的话单次询问是 $O(n)$ 的。
于是有了一个叫整除分块的东西。
打表找规律可得 $\lfloor \frac{n}{d}\rfloor$ 只有 $O(\sqrt{n})$ 种取值，
$\therefore \lfloor \frac{n}{d}\rfloor \ast \lfloor \frac{m}{d}\rfloor$ 只有 $O(\sqrt{n}+\sqrt{m})$ 种取值，可以分块。
当 $n=16,m=20$ 时的表如下：

$d$	$1$	$2$	$3$	$4$	$5$	$6$	$7$	$8$	$9$	$10$	$11$	$12$	$13$	$14$	$15$	$16$
$\lfloor \frac{n}{d}\rfloor$	$16$	$8$	$5$	$4$	$3$	$2$	$2$	$2$	$1$	$1$	$1$	$1$	$1$	$1$	$1$	$1$
$\lfloor \frac{m}{d}\rfloor$	$20$	$10$	$6$	$5$	$4$	$3$	$2$	$2$	$2$	$2$	$1$	$1$	$1$	$1$	$1$	$1$

以下给出证明：
$Proof:$
当 $d\in [1,\sqrt{n}]$ 时，$d$ 的取值只有 $O(\sqrt{n})$ 个，
$\therefore \lfloor \frac{n}{d}\rfloor$ 的取值只有 $O(\sqrt{n})$ 个。
当 $d\in (\sqrt{n},n]$ 时， $\lfloor \frac{n}{d}\rfloor \in [1,\sqrt{n}]$，
$\therefore \lfloor \frac{n}{d}\rfloor$ 的取值只有 $O(\sqrt{n})$ 个。
综上所述，$\lfloor \frac{n}{d}\rfloor$ 的取值只有 $O(\sqrt{n})$ 个。
$\therefore \lfloor \frac{n}{d}\rfloor$ 会被划分成 $O(\sqrt{n})$ 段区间，每段区间内 $\lfloor \frac{n}{d}\rfloor$ 的值相等
$\therefore$ 整张表会被划分成 $O(\sqrt{n}+\sqrt{m})$ 段区间，每段区间内 $\lfloor \frac{n}{d}\rfloor$ 和 $\lfloor \frac{n}{d}\rfloor$ 都相等。
\quad 当且仅当上一行的每段区间都与下一行的两段区间有交集时取到最大值。
证毕。
此外，还需要确定当前区间的末尾。
设当前区间的左端点为 $l$ ，右端点为 $r$ 。
则满足 $\lfloor \frac{n}{l}\rfloor =\lfloor \frac{n}{r}\rfloor$，
又 $\because \lfloor \frac{n}{l}\rfloor \le \frac{n}{l}$
$\therefore \lfloor \frac{n}{l}\rfloor \le \frac{n}{r}$
$\therefore r\le \frac{n}{\lfloor \frac{n}{l}\rfloor }$
$\therefore r\le \lfloor \frac{n}{\lfloor \frac{n}{l}\rfloor }\rfloor$
对 $m$ 同理，将两个右端点取 $min$ 即可。

---

$e.g.2$

求 $\sum _{i=1}^n\sum _{j=1}^m[gcd(i,j)==k]$
$Solution:$
显然 $gcd(i,j)==k\iff i\mid k\wedge j\mid k\wedge gcd(\frac{i}{k},\frac{j}{k})==1$
$\therefore \sum _{i=1}^n\sum _{j=1}^m[gcd(i,j)==k]=\sum _{i=1}^{\lfloor \frac{n}{k}\rfloor }\sum _{j=1}^{\lfloor \frac{m}{k}\rfloor }[gcd(i,j)==1]$
然后就和 $e.g.1$一样了。

---

$e.g.3$

求 $\sum _{i=1}^n\sum _{j=1}^{m}i\ast j\ast [gcd(i,j)==k]$
$Solution:$
$\begin{array}{l}\sum _{i=1}^n\sum _{j=1}^{m}i\ast j\ast [gcd(i,j)==k]\\ =\sum _{i=1}^{\lfloor \frac{n}{k}\rfloor }\sum _{j=1}^{\lfloor \frac{m}{k}\rfloor }i\ast k\ast j\ast k\ast [gcd(i,j)==1]\\ =k^2\ast \sum _{i=1}^{\lfloor \frac{n}{k}\rfloor }i\ast \sum _{j=1}^{\lfloor \frac{m}{k}\rfloor }j\ast [gcd(i,j)==1]\\ =k^2\ast \sum _{d=1}^{\lfloor \frac{n}{k}\rfloor }\mu (d)\ast \sum _{i=1}^{\lfloor \frac{n}{k\ast d}\rfloor }i\ast \sum _{j=1}^{\lfloor \frac{m}{k\ast d}\rfloor }j\end{array}$
后面的两个 $\sum$ 都可以预处理，直接整除分块即可 $O(\sqrt{n})$ 求解。

---

$e.g.4$

求 $\sum _{i=1}^n\sum _{j=1}^{m}gcd(i,j)$
$Solution:$
$\begin{array}{l}\sum _{i=1}^n\sum _{j=1}^{m}gcd(i,j)\\ =\sum _{i=1}^n\sum _{j=1}^m\sum _{k=1}^{n}k\ast [gcd(i,j)==k]\\ =\sum _{k=1}^{n}k\ast \sum _{i=1}^n\sum _{j=1}^m[gcd(i,j)==k]\\ =\sum _{k=1}^{n}k\ast \sum _{i=1}^{\lfloor \frac{n}{k}\rfloor }\sum _{j=1}^{\lfloor \frac{m}{k}\rfloor }\sum _{d\mid gcd(i,j)}\mu (d)\\ =\sum _{k=1}^{n}k\ast \sum _{d=1}^{\lfloor \frac{n}{k}\rfloor }\mu (d)\sum _{i=1}^{\lfloor \frac{n}{k}\rfloor }\sum _{j=1}^{\lfloor \frac{m}{k}\rfloor }1\\ =\sum _{k=1}^{n}k\ast \sum _{d=1}^{\lfloor \frac{n}{k}\rfloor }\mu (d)\ast \lfloor \frac{n}{k\ast d}\rfloor \ast \lfloor \frac{m}{k\ast d}\rfloor \end{array}$
做到这一步之后，有一个套路性的东西：
枚举$k\ast d$
设 $T=k\ast d$
$\therefore \sum _{k=1}^{n}k\ast \sum _{d=1}^{\lfloor \frac{n}{k}\rfloor }\mu (d)\ast \lfloor \frac{n}{k\ast d}\rfloor \ast \lfloor \frac{m}{k\ast d}\rfloor =\sum _{T=1}^n\lfloor \frac{n}{T}\rfloor \ast \lfloor \frac{m}{T}\rfloor \ast \sum _{d\mid T}\mu (d)\ast \frac{T}{d}$
根据之前 $\phi$ 函数的性质可得，
$\sum _{d\mid T}\mu (d)\ast \frac{T}{d}=\phi (T)$
$\therefore \sum _{T=1}^n\lfloor \frac{n}{T}\rfloor \ast \lfloor \frac{m}{T}\rfloor \ast \sum _{d\mid T}\mu (d)\ast \frac{T}{d}=\sum _{T=1}^n\lfloor \frac{n}{T}\rfloor \ast \lfloor \frac{m}{T}\rfloor \ast \phi (T)$
预处理 $\phi$ 函数的前缀和即可整除分块， $O(\sqrt{n})$求解。

---

$e.g.5$

求 $\sum _{i=1}^n\sum _{j=1}^{m}lcm(i,j)$
$Solution:$
$\begin{array}{l}\sum _{i=1}^n\sum _{j=1}^{m}lcm(i,j)\\ =\sum _{i=1}^n\sum _{j=1}^m\frac{i\ast j}{gcd(i,j)}\\ =\sum _{i=1}^n\sum _{j=1}^m\sum _{k=1}^n[gcd(i,j)==k]\ast \frac{i\ast j}{k}\\ =\sum _{k=1}^n\sum _{i=1}^{\lfloor \frac{n}{k}\rfloor }\sum _{j=1}^{\lfloor \frac{m}{k}\rfloor }[gcd(i,j)==1]\ast \frac{i\ast k\ast j\ast k}{k}\\ =\sum _{k=1}^{n}k\ast \sum _{i=1}^{\lfloor \frac{n}{k}\rfloor }\sum _{j=1}^{\lfloor \frac{m}{k}\rfloor }[gcd(i,j)==1]\ast i\ast j\\ =\sum _{k=1}^{n}k\ast \sum _{d=1}^{\lfloor \frac{n}{k}\rfloor }\mu (d)\ast d^2\ast \sum _{i=1}^{\lfloor \frac{n}{k\ast d}\rfloor }i\ast \sum _{j=1}^{\lfloor \frac{m}{k\ast d}\rfloor }j\\ =\sum _{T=1}^n\sum _{i=1}^{\lfloor \frac{n}{T}\rfloor }i\ast \sum _{j=1}^{\lfloor \frac{m}{T}\rfloor }j\ast \sum _{d\mid T}\mu (d)\ast d^2\ast \frac{T}{d}\\ =\sum _{T=1}^{n}T\ast \sum _{i=1}^{\lfloor \frac{n}{T}\rfloor }i\ast \sum _{j=1}^{\lfloor \frac{m}{T}\rfloor }j\ast \sum _{d\mid T}\mu (d)\ast d\end{array}$
设 $f(n)=\sum _{i=1}^{n}i$ ，则
$\sum _{T=1}^{n}T\ast \sum _{i=1}^{\lfloor \frac{n}{T}\rfloor }i\ast \sum _{j=1}^{\lfloor \frac{m}{T}\rfloor }j\ast \sum _{d\mid T}\mu (d)\ast d=\sum _{T=1}^{n}f(\lfloor \frac{n}{T}\rfloor )\ast f(\lfloor \frac{m}{T}\rfloor )\ast \sum _{d\mid T}\mu (d)\ast d$
显然 $f$ 是可以预处理的，那么主要问题就是求最后一个 $\sum$ 。
设 $g(n)=\sum _{d\mid n}\mu (d)\ast d$。
然后打个表，找找规律，会发现一个一个性质：$g(n)$是个积性函数！
以下给出证明：
$Proof:$
设 $gcd(a,b)=1$，则
$\begin{array}{l}g(a)\ast g(b)\\ =\sum _{d_1\mid a}\mu (d_1)\ast d_1\ast \sum _{d_2\mid b}\mu (d_2)\ast d_2\\ \because gcd(a,b)=1\\ d_1\mid a\\ d_2\mid b\\ \therefore gcd(d1,d2)=1\\ \therefore \sum _{d_1\mid a}\mu (d_1)\ast d_1\ast \sum _{d_2\mid b}\mu (d_2)\ast d_2\\ =\sum _{(d_1\ast d_2)\mid (a\ast b)}\mu (d_1)\ast \mu (d_2)\ast d_1\ast d_2\\ =\sum _{(d_1\ast d_2)\mid (a\ast b)}\mu (d_1\ast d_2)\ast (d_1\ast d_2)\\ =g(a\ast b)\end{array}$
对两 $\sum$ 合并做出解释：
当 $\mu (d_1)\ast \mu (d_2)=0$ 时，显然对结果无影响。
当 $\mu (d_1)\ast \mu (d_2)̸=0$ 时，
设 S = { x ∣ μ ( x ) ̸ = 0 ∧ x ∣ d 1 } , T = { x ∣ μ ( x ) ̸ = 0 ∧ x ∣ d 2 } S=\{x \mid \mu(x) \not=0 \land x\mid d_1\},T=\{x \mid \mu(x) \not=0 \land x\mid d_2\} S={x∣μ(x)̸​=0∧x∣d1​},T={x∣μ(x)̸​=0∧x∣d2​} ，则 $S\cap T=\varnothing$ 。
显然 $(x_1\ast x_2)\mid (a\ast b)$ 且 $x_1\ast x_2̸=x_1^{\prime }\ast x_2^{\prime },x_1^{\prime }\in S,x_2^{\prime }\in T,x_1̸=x_1^{\prime },x_2̸=x_2^{\prime }$。
即 $a\ast b$ 的每一个 $\mu$ 值不为 $0$ 的约数都唯一对应着一个二元组 $(x_1,x_2)$ 。
证毕。
众所周知，所有的积性函数都可以用线性筛筛出来。
设 $n$ 的最小质因子为 $p$， $x=\frac{n}{p}$
考虑求 $g(n)$ ：
当 $n\in prime$ 时， $g(n)=\mu (1)\ast 1+\mu (n)\ast n=1-n$；
当 $n=p^k$ 时， 显然 $\mu (p^i)=0(i\ge 2)\Rightarrow g(n)=1-p=g(p)$；
当 $p\nmid x$ 时， 则 $gcd(x,p)=1\Rightarrow g(n)=g(x)\ast g(p)$；
当 $p\mid x$ 时， 则 $g(n)=g(\frac{n}{p^k})\ast g(p^k)=g(\frac{n}{；}{p}^k)\ast g(p)=g(\frac{n}{p^k})\ast g(p^{k-1})=g(x)(p^k\mid n\wedge p^{k+1}\nmid n,k\in {\mathbb{N}}^{\ast },k\ge 2)$。
那么就可以 $O(n)$ 预处理 $g(n)$ 了。
$\therefore \sum _{i=1}^n\sum _{j=1}^{m}lcm(i,j)=\sum _{T=1}^{n}f(\lfloor \frac{n}{T}\rfloor )\ast f(\lfloor \frac{m}{T}\rfloor )\ast g(T)$
单次询问即可做到 $O(\sqrt{n})$ 。

---

$e.g.6$

求 $\sum _{i=1}^n\sum _{j=1}^{m}gcd(i,j{)}^x$
$Solution:$
$\begin{array}{l}\sum _{i=1}^n\sum _{j=1}^{m}gcd(i,j{)}^x\\ =\sum _{k=1}^n\sum _{i=1}^n\sum _{j=1}^m[gcd(i,j)==k]\ast k^x\\ =\sum _{k=1}^n{k}^x\ast \sum _{i=1}^{\lfloor \frac{n}{k}\rfloor }\sum _{j=1}^{\lfloor \frac{m}{k}\rfloor }[gcd(i,j)==1]\\ =\sum _{k=1}^n{k}^x\ast \sum _{d=1}^{\lfloor \frac{n}{k}\rfloor }\mu (d)\ast \lfloor \frac{n}{k\ast d}\rfloor \ast \lfloor \frac{m}{k\ast d}\rfloor \\ =\sum _{T=1}^n\lfloor \frac{n}{T}\rfloor \ast \lfloor \frac{m}{T}\rfloor \ast \sum _{d\mid T}{d}^x\ast \mu (\frac{T}{d})\end{array}$
发现 $f(x)=d^x$ 是个完全积性函数，当然也是积性函数。
$\therefore g(n)=\sum _{d\mid n}{d}^x\ast \mu (\frac{n}{d})=(f\ast \mu {)}_{(n)}$ 也是积性函数。
设 $n$ 的最小质因子为 $p$， $x=\frac{n}{p}$
考虑求 $g(n)$ ：
当 $n\in prime$ 时， $g(n)=\mu (1)\ast n^x+\mu (n)\ast 1^x=n^x-1$；
当 $n=p^k$ 时， 显然 $\mu (p^i)=0(i\ge 2)\Rightarrow g(n)=n^x-(\frac{n}{p}{)}^x=(\frac{n}{p}{)}^x\ast g(p)$；
当 $p\nmid x$ 时， 则 $gcd(x,p)=1\Rightarrow g(n)=g(x)\ast g(p)$；
当 $p\mid x$ 时， 则 $g(n)=g(\frac{n}{p^k})\ast g(p^k)=g(\frac{n}{p^k})\ast g(p)\ast p^{(k-1)\ast x}=g(\frac{n}{p^k})\ast g(p)\ast p^{(k-2)\ast x}\ast p^x=g(\frac{n}{p^k})\ast g(p^{k-1})\ast p^k=g(x)\ast p^k$
$(p^k\mid n\wedge p^{k+1}\nmid n,k\in {\mathbb{N}}^{\ast },k\ge 2)$。
和 $e.g.5$ 的筛法大同小异。
之后 $O(\sqrt{n})$ 求解即可。

---

$e.g.7$

求$\sum _{i=1}^n\sum _{j=1}^m\phi (i\ast j)$
$Solution:$
显然 $\phi (i\ast j)$ 不好处理，考虑将其转化。
结论： $\phi (i)\ast \phi (j)\ast gcd(i,j)=\phi (i\ast j)\ast \phi (gcd(i,j))$
以下给出该结论的证明。
$Proof:$
设 S = { x ∣ x ∈ p r i m e ∧ x ∣ i } , T = { x ∣ x ∈ p r i m e ∧ x ∣ j } S=\{x \mid x \in prime \land x \mid i\},T=\{x \mid x \in prime \land x \mid j\} S={x∣x∈prime∧x∣i},T={x∣x∈prime∧x∣j} ，则 S ∪ T = { x ∣ x ∈ p r i m e ∧ x ∣ ( i ∗ j ) } , S ∩ T = { x ∣ x ∈ p r i m e ∧ x ∣ g c d ( i , j ) } S \cup T = \{x \mid x \in prime \land x \mid (i*j)\},S \cap T = \{x \mid x \in prime \land x \mid gcd(i,j)\} S∪T={x∣x∈prime∧x∣(i∗j)},S∩T={x∣x∈prime∧x∣gcd(i,j)}
$\begin{array}{l}\phi (i)\ast \phi (j)\\ =i\ast \prod _{x\in S}\frac{x-1}{x}\ast j\ast \prod _{x\in T}\frac{x-1}{x}\\ =i\ast j\ast \prod _{x\in (S\cup T)}\frac{x-1}{x}\ast \prod _{x\in (S\cap T)}\frac{x-1}{x}\\ =\frac{\phi (i\ast j)\ast \phi (gcd(i,j))}{gcd(i,j)}\end{array}$
$\begin{array}{l}\therefore \sum _{i=1}^n\sum _{j=1}^m\phi (i\ast j)\\ =\sum _{i=1}^n\sum _{j=1}^m\frac{\phi (i)\ast \phi (j)\ast gcd(i,j)}{\phi (gcd(i,j))}\\ =\sum _{k=1}^n\sum _{i=1}^n\sum _{j=1}^m[gcd(i,j)==k]\ast \frac{\phi (i)\ast \phi (j)\ast k}{\phi (k)}\\ =\sum _{k=1}^n\frac{k}{\phi (k)}\sum _{i=1}^{\lfloor \frac{n}{k}\rfloor }\sum _{j=1}^{\lfloor \frac{m}{k}\rfloor }[gcd(i,j)==1]\ast \phi (i\ast k)\ast \phi (j\ast k)\\ =\sum _{k=1}^n\frac{k}{\phi (k)}\sum _{d=1}^{\lfloor \frac{n}{k}\rfloor }\mu (d)\ast \sum _{i=1}^{\lfloor \frac{n}{k\ast d}\rfloor }\sum _{j=1}^{\lfloor \frac{m}{k\ast d}\rfloor }\phi (i\ast k\ast d)\ast \phi (j\ast k\ast d)\\ =\sum _{T=1}^n\sum _{d\mid T}\frac{d}{\phi (d)}\ast \mu (\frac{T}{d})\ast \sum _{i=1}^{\lfloor \frac{n}{T}\rfloor }\phi (i\ast T)\ast \sum _{j=1}^{\lfloor \frac{m}{T}\rfloor }\phi (j\ast T)\end{array}$
式子已经化到最简了，考虑预处理某些 $\sum$ 。
为了便于表达，引入下列函数。

• $f(n)=\sum _{d\mid n}\frac{d}{\phi (d)}\ast \mu (\frac{n}{d})$
• $g(k,n)=\sum _{i=1}^n\phi (i\ast k)$
• $h(x,y,n)=\sum _{T=1}^n\sum _{d\mid T}\frac{d}{\phi (d)}\ast \mu (\frac{T}{d})\ast \sum _{i=1}^x\phi (i\ast T)\ast \sum _{j=1}^y\phi (j\ast T)$
  考虑计算这三个函数。
  对于 $f$ ，可以用类似埃氏筛法的方式，在预处理出 $\phi$ 函数的情况下 $O(n\ln n)$ 求出。
  对于 $g$ ，有递推式 $g(k,n)=g(k,n-1)+\phi (n\ast k)$ 。由于总的有用的状态只有 $O(n\ln n)$ 级别，开 $n$ 个 $vector$ 预处理即可。
  对于 $h$ ，有递推式 $h(x,y,n)=h(x,y,n-1)+f(n)\ast g(n,x)\ast g(n,y)$ 。
  全部预处理空间复杂度过高，考虑设置一个上限 $top$ 。
  预处理所有 $x,y\le top$ 的 $h$ 。
  对于 $\lfloor \frac{n}{T}\rfloor \le top$ 的 $T$ ，可以在 $O(\sqrt{n})$ 的时间内全部计算出来。
  对于 $\lfloor \frac{n}{T}\rfloor >top$ 的 $T$ ，有 $T<\frac{n}{top}$ ，直接用 $f(T)\ast g(T,\lfloor \frac{N}{T}\rfloor )\ast g(T,\lfloor \frac{m}{T}\rfloor )$ 暴力计算，时间复杂度为 $O(\frac{n}{top})$。
  总时间复杂度： $O(n\ln n+to{p}^2\ast n+q\ast (\frac{n}{top}+\sqrt{n}))$。
  当 $top=\sqrt[3]{n}$ 时取到最小值 $O(n^{\frac{5}{3}})$。