mobius inversion 莫比乌斯反演例题总结

最新推荐文章于 2024-04-25 12:45:28 发布

winsunboy

最新推荐文章于 2024-04-25 12:45:28 发布

阅读量1.3k

点赞数 28

分类专栏：总结文章标签：算法

本文链接：https://blog.csdn.net/winsunboy/article/details/136082340

版权

总结专栏收录该内容

5 篇文章 0 订阅

订阅专栏

莫比乌斯反演总结

莫比乌斯反演，实际上主要利用莫比乌斯函数的性质： $\epsilon = \mu * 1$ 。一般会结合 整除分块、线性筛 等算法进行考察。下面总结了几道 Luogu 上较为典型的题目。

每道题若涉及 $n$ 和 $m$ 两个量，默认 $\le m$ 。

P3455 [POI2007] ZAP-Queries

$T$ 组数据，对于每组数据，给定 $n$ ， $m$ ， $d$ ，求：

$\sum_{i=1}^n \sum_{j=1}^m [\gcd(i, j) = d]$

其中： $\le T \le 5 \times 10^4$ ， $\le d \le n, m \le 5 \times 10^4$ 。

解：

$\begin{aligned} &\sum_{i=1}^n \sum_{j=1}^m [\gcd(i, j) = d]\\ =&\sum_{i=1}^{\lfloor\frac nd\rfloor}\sum_{j=1}^{\lfloor\frac md\rfloor}[\gcd(i, j) = 1]\\ =&\sum_{i=1}^{\lfloor\frac nd\rfloor}\sum_{j=1}^{\lfloor\frac md\rfloor}\sum_{k \mid \gcd(i, j)} \mu(k)\\ =&\sum_{k=1}^{\lfloor\frac nd\rfloor}\mu(k)\lfloor\frac n{dk} \rfloor \lfloor\frac m{dk} \rfloor \end{aligned}$

$\mu(k)$ 用前缀和预处理，查询用整除分块算。

P1390 公约数的和

给定 $n$ ，求：

$\sum_{i = 1}^n \sum_{j = i + 1}^n \gcd(i, j)$

其中： $\leq n \leq 2 \times 10^6$ 。

解：

原式可以看成：

$\frac{\sum_{i=1}^n \sum_{i=1}^n \gcd(i, j) -\frac{n(n+1)}2}2$

注意到：

$\begin{aligned} &\sum_{i=1}^n\sum_{j=1}^n\gcd(i, j)\\ =&\sum_{d=1}^n d \sum_{i=1}^n \sum_{j=1}^n [\gcd(i, j) = 1]\\ =&\sum_{d=1}^n d \sum_{i=1}^{\lfloor \frac nd \rfloor}\sum_{j=1}^{\lfloor \frac nd \rfloor} \sum_{k \mid \gcd(i, j)} \mu(k)\\ =&\sum_{d=1}^n d \sum_{k=1}^{\lfloor \frac nd \rfloor} \mu(k) \lfloor \frac n {dk} \rfloor ^2\\ =&\sum_{x=1}^n\sum_{d \mid x} d\mu(\frac xd) \lfloor\frac nx \rfloor^2\\ =&\sum_{x=1}^n\varphi(x)\lfloor\frac nx \rfloor^2 \end{aligned}$

前缀和+数论分块，绰绰有余。

上面用到了欧拉函数的一个重要结论： $\mu = \varphi * 1 * \mu = \varphi * \epsilon = \varphi$ 。

双倍经验：P1447 [NOI2010] 能量采集

提示：

Years of struggling in OI, without long long, everything turns out to be a null.
Please remember to use 1ll* to avoid Wrong Answer!

P1829 [国家集训队] Crash的数字表格 / JZPTAB

给定 $n$ ， $m$ ，求：

$\sum_{i=1}^n \sum_{j=1}^m \operatorname{lcm}(i,j) \bmod 20101009$

其中： $\le n, m \le 10^7$ 。

解：

$\begin{aligned} &\sum_{i=1}^n \sum_{j=1}^m \operatorname{lcm}(i, j)\\ =&\sum_{d=1}^n \sum_{i=1}^n \sum_{j=1}^m \frac{ij}d [\gcd(i, j) = d]\\ =&\sum_{d=1}^n \sum_{i=1}^{\lfloor \frac nd \rfloor} \sum_{j=1}^{\lfloor \frac md \rfloor} ijd [\gcd(i, j) = 1]\\ =&\sum_{d=1}^n d\sum_{i=1}^{\lfloor \frac nd \rfloor} i\sum_{j=1}^{\lfloor \frac md \rfloor} j \sum_{k\mid \gcd(i, j)} \mu(k)\\ =&\sum_{d=1}^n d \sum_{k=1}^{\lfloor\frac nd\rfloor} \mu(k) \sum_{i=1}^{\lfloor\frac n{kd}\rfloor} ik \sum_{j=1}^{\lfloor\frac m{kd}\rfloor} jk\\ =&\sum_{d=1}^n d \sum_{k=1}^{\lfloor\frac nd\rfloor} k^2\mu(k) \sum_{i=1}^{\lfloor\frac n{kd}\rfloor} i \sum_{j=1}^{\lfloor\frac m{kd}\rfloor} j\\ =&\sum_{x=1}^n\sum_{k\mid x}\frac xk \cdot k^2\mu(k) \sum_{i=1}^{\lfloor \frac nx \rfloor} i \sum_{j=1}^{\lfloor \frac mx \rfloor} j\\ =&\sum_{x=1}^nx\sum_{k \mid x} k\mu(k) \sum_{i=1}^{\lfloor \frac nx \rfloor} i \sum_{j=1}^{\lfloor \frac mx \rfloor} j \end{aligned}$

$x = k d$ 也是典型的套路啦。

$\sum_{i=1}^{\lfloor \frac nx \rfloor} i$ 是等差数列求和，确定 $\lfloor\frac nx\rfloor$ 之后， $O (1)$ 可求。

$\sum _{d \mid n}d\mu(d)$ 是个积性函数，线性筛 $O (n)$ 可求。（只要是积性函数，都能够线性筛。）

可参考线性筛 $\sum_{d \mid n} d\varphi(d)$ 的文章。

先证明 $f (n)$ 是积性函数：

若 $id(n)\mu(n)$ ，注意到 $i d (n)$ 和 $\mu(n)$ 都是积性函数，那么 $g (n)$ 是积性函数。
而 $f = 1 * g$ ，根据狄利克雷卷积， $f (n)$ 也是积性函数。

再考虑线性筛的方法：

$\sum _{d \mid n}d\mu(d)$

设 $p$ 为质数。

显然， $f (1) = 1$ ， $f (p) = 1 - p$ 。

若 $\nmid i$ ，显然 $f (i p) = f (i) f (p)$ 。

若 $\mid i$ ， $f (i p) = f (i)$ 。因为 $i p$ 的因数中，比 $i$ 多的那一部分 $d$ ，一定含有平方因子。而根据莫比乌斯函数的定义，含有平方因子会导致 $\mu(d) = 0$ ，从而这一部分对结果没有影响。

以上算法预处理 $O (n)$ ，单次询问 $O(\sqrt n)$ ，可以处理多测的情况。

P3327 [SDOI2015] 约数个数和

$T$ 组数据，对于每组数据，给定 $n, m$ ，求：

$\sum_{i=1}^n\sum_{j=1}^m\tau(ij)$

其中 $1\le T,n,m \le 5 \times 10^4$ ， $\tau(n)$ 表示 $n$ 的约数个数。

遇到陌生函数，不要惊慌，把它转化成熟悉的函数。

引理： $\tau(ij) = \sum_{x \mid i} \sum_{y \mid j} [\gcd(x, y) = 1]$ 。

证明：

先复习一个小学生都会的知识：

若 $\prod p_i^{c_i}$ ，则 $\tau(n) = \prod (c_i + 1)$ 。

不难发现，上述结论意味着 $\tau(n)$ 是积性函数。所以我们可以分别考虑 $i$ 和 $j$ 的每一个质因子 $p$ 。

假设 $i$ 和 $j$ 含有的质因子 $p$ 的个数分别为 $\alpha$ 和 $\beta$ 。

$(\alpha + \beta + 1)$ 可以看成前 $\alpha$ 个从 $i$ 中取，后 $\beta$ 个从 $j$ 中取，还有一个就是两边都不取。

那么就不会出现两边都取到质因子 $p$ 的情况，从而得到 $\tau(ij) = \sum_{x \mid i} \sum_{y \mid j} [\gcd(x, y) = 1]$ 。

现在可以直接推式子了。

$\begin{aligned} &\sum_{i=1}^n \sum_{j=1}^m \tau(ij)\\ =&\sum_{i=1}^n \sum_{j=1}^m \sum_{x \mid i} \sum_{y \mid j} [\gcd(x, y) = 1]\\ =&\sum_{x=1}^n \lfloor\frac nx\rfloor \sum_{y=1}^m \lfloor \frac my \rfloor[\gcd(x, y) = 1]\\ =&\sum_{x=1}^n \lfloor\frac nx\rfloor \sum_{y=1}^m \lfloor \frac my \rfloor \sum_{d\mid\gcd(x, y)} \mu(d)\\ =&\sum_{d=1}^n\mu(d) \sum_{i=1}^{\lfloor \frac nd \rfloor}\lfloor \frac{n}{id}\rfloor \sum_{j=1}^{\lfloor \frac md \rfloor}\lfloor \frac{m}{jd}\rfloor \end{aligned}$

若 $\sum_{i=1}^n \lfloor \frac ni \rfloor$ ，则原式等于：

$\sum_{d=1}^n \mu(d) f(\lfloor \frac nd \rfloor) f(\lfloor \frac md \rfloor)$

预处理 $\mu(d)$ 的前缀和，对上述式子进行 $O(\sqrt n)$ 的整除分块，计算 $f (n)$ 的时候注意进行记忆化（不记忆化会TLE）。总时间复杂度 $O(n\sqrt n + T\sqrt n)$ 。

P6810 「MCOI-02」Convex Hull 凸包

本题标题有误导性，怀疑比赛计算几何题被毙掉后，换了莫比乌斯反演题，但标题没有换？

给定 $n$ ， $m$ ， $p$ ，求：

$\sum_{i=1}^n\sum_{j=1}^m \tau(i)\tau(j)\tau(\gcd(i,j)) \bmod p$

其中： $\le n, m \le 2 \times 10^6$ ， $\le p \le 10^9$ 。

解：

$\begin{aligned} &\sum_{i=1}^n\sum_{j=1}^m \tau(i)\tau(j)\tau(\gcd(i,j))\\ =&\sum_{d=1}^n \tau(d) \sum_{i=1}^n\tau(i) \sum_{j=1}^m\tau(j)[\gcd(i, j) = d]\\ =&\sum_{d=1}^n \tau(d) \sum_{i=1}^{\lfloor\frac nd\rfloor}\tau(id) \sum_{j=1}^{\lfloor \frac md \rfloor} \tau(jd)[\gcd(i, j) = 1]\\ =&\sum_{d=1}^n \tau(d) \sum_{i=1}^{\lfloor\frac nd\rfloor}\tau(id) \sum_{j=1}^{\lfloor \frac md \rfloor} \tau(jd)\sum_{k \mid \gcd(i, j)} \mu(k)\\ =&\sum_{d=1}^n \tau(d) \sum_{k=1}^{\lfloor\frac nd\rfloor} \mu(k) \sum_{i=1}^{\lfloor\frac n{dk}\rfloor}\tau(idk) \sum_{j=1}^{\lfloor \frac m{dk} \rfloor} \tau(jdk)\\ =&\sum_{x=1}^n\sum_{d \mid x} \tau(d) \mu(\frac xd)\sum_{i=1}^{\lfloor\frac nx\rfloor}\tau(ix) \sum_{j=1}^{\lfloor \frac mx \rfloor} \tau(jx)\\ =&\sum_{x=1}^n\sum_{i=1}^{\lfloor\frac nx\rfloor}\tau(ix) \sum_{j=1}^{\lfloor \frac mx \rfloor} \tau(jx) \end{aligned}$

最后一步解释：因为 $\tau = 1 * 1$ ， $\mu * 1 = \epsilon$ ，所以 $\tau * \mu = 1 * 1 * \mu = 1 * \epsilon = 1$ 。其他都是典型的套路。

这里可以线性筛预处理 $\tau$ 函数的值，外层数论分块后内层暴力处理。

$\tau(1) = 1$ ， $\tau(p) = 2$ 。对于质数 $\mid i$ ， $\tau(pi) = 2\tau(i) - \tau(i/p)$ 。

总结

归纳总结做莫比乌斯反演的题目的方法：

把陌生的式子化成熟悉的式子，直到可以拿起式子计算。
牢记、抓住 $\mu * 1 = \epsilon$ ，将 $[\gcd(i, j) = 1]$ 化为 $\sum_{k \mid \gcd(i, j)} \mu(k)$ 。
P1390 P1829 P6810 三题，都用了 $x = k d$ 的套路，套了之后可以出现卷积。
掌握线性筛的方法，学会推出陌生积性函数的线性筛公式。