莫比乌斯反演

莫比乌斯反演

莫比乌斯反演是数论中的重要内容。对于一些函数 ，如果很难直接求出它的值，而容易求出其倍数和或约数和 ，那么可以通过莫比乌斯反演简化运算，求得函数的值。

前置知识

整除定理

$$\lfloor \frac{a}{bc}\rfloor =\lfloor \frac{a}{bc}\rfloor =\lfloor \frac{\lfloor \frac{a}{b}\rfloor }{c}\rfloor$$

商定理

$$\forall n\in {\mathbb{N}}_{+},\left|\left\{\lfloor \frac{n}{d}\rfloor \mid d\in {\mathbb{N}}_{+},d\le n\right\}\right|\le \lfloor 2\sqrt{n}\rfloor$$

数论分块

考虑求和式${\sum }_{i=1}^m\lfloor \frac{k}{i}\rfloor$，其中$m$和$k$均为已知常数。

我们发现，$i$在一定的范围内$\lfloor \frac{k}{i}\rfloor$的值相等，因此整个$\lfloor \frac{k}{i}\rfloor$将会分成几个相等的块，加快计算速度。

即，对于任意的$i$，我们需要找到最大的$j$，并且$m\ge j\ge i$，$\lfloor \frac{k}{i}\rfloor =\lfloor \frac{k}{j}\rfloor$。

给出结论：$j=\lfloor \frac{k}{\lfloor \frac{k}{i}\rfloor }\rfloor$

利用上述结论，我们可以进行分块$[i,j]$求和。

P2261

$$ans=\sum _{i=1}^n(k\mathrm{m}\mathrm{o}\mathrm{d}i)=\sum _{i=1}^{n}k-i\lfloor \frac{k}{i}\rfloor$$

这样就可以利用整除分块去做了。

#include <bits/stdc++.h>

using namespace std;

typedef long long ll;

#define FR freopen("in.txt", "r", stdin)
#define FW freopen("out.txt", "w", stdout)

int main()
{
    ll n, k;
    scanf("%lld %lld", &n, &k);

    ll ans = n * k;

    for (ll l = 1, r; l <= n; l = r + 1)
    {
        // 寻找 r
        if (k / l == 0)
        {
            r = n;
        }
        else
        {
            r = min(n, k / (k / l));
        }

        ll t = (k / l) * (l + r) * (r - l + 1) / 2;
        ans -= t;
    }
    printf("%lld", ans);
    return 0;
}

积性函数和狄利克雷卷积

参考我的狄利克雷卷积文章。

莫比乌斯函数

$\mu$为莫比乌斯函数，定义为：

$$\mu (n)=\{\begin{array}{ll}1& n=1\\ 0& n\text{ 含有平方及以上因子}\\ (-1{)}^k& k\text{ 为 }n\text{ 的本质不同质因子个数}\end{array}$$

其中令$k=\omega (n)$，$\omega (n)$为加性函数。

并且:

$$\mu \ast 1=\epsilon$$

这表明，$\mu$是$1$的狄利克雷卷积的逆元。

莫比乌斯函数的求法

如何求$\mu$函数呢，利用线性筛和积性函数的性质，可以参考我的筛法的文章。

莫比乌斯函数与欧拉函数的关系

我们在之前已经证明$\phi \ast 1=Id$，对两边做$\mu$的卷积得到：

$$\phi (n)=\sum _{d\mid n}d\ast \mu (\frac{n}{d})$$

莫比乌斯反演

已知函数$f=g\ast 1$，我们现在要求解$g$函数的表达值，将两边对$\mu$做卷积，得到：

$$f\ast \mu =g$$

这就是莫比乌斯反演。

具体写成和积式为：

如果有$f(n)={\sum }_{d\mid n}g(d)$，那么有$g(n)={\sum }_{d\mid n}\mu (d)f(\frac{n}{d})$。

如果有$f(n)={\sum }_{d\mid n}g(d)$，那么有$g(n)={\sum }_{d\mid n}\mu (\frac{n}{d})f(d)$。

例题

类型一：利用$\phi \ast 1=Id$反向枚举$d$

P2257

答案求：

$$\sum _{p\in primes}\sum _i^N\sum _j^M[\gcd (i,j)=p]$$

对gcd进行套路除法：

$$\sum _{p\in primes}\sum _i^{\lfloor \frac{N}{p}\rfloor }\sum _j^{\lfloor \frac{M}{p}\rfloor }[\gcd (i,j)=1]$$

我们知道$\mu \ast 1=\epsilon$，因此${\sum }_{d\mid x}\mu (d)=[x=1]$，带入gcd，得到：

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

这是莫比乌斯反演很经典的套路代换，带入到上式中：

$$\sum _{p\in primes}\sum _i^{\lfloor \frac{N}{p}\rfloor }\sum _j^{\lfloor \frac{M}{p}\rfloor }\sum _{d\mid \gcd (i,j)}\mu (d)$$

求和太多了，我们改变一下求和次序，我们枚举$d$，从而带入$i,j$，也就是说如果$d\mid \gcd (i,j)$等价于$d\mid i\wedge d\mid j$，这也是是数论常见的处理套路，即为：

$$\sum _{p\in primes}\sum _{d=1}^{\min (\lfloor \frac{N}{p}\rfloor ,\lfloor \frac{M}{p}\rfloor )}\lfloor \frac{N}{pd}\rfloor \times \lfloor \frac{M}{pd}\rfloor \times \mu (d)$$

此时我们可以枚举$p,d$进行计算，此时复杂度是$O(n^2)$的，我们这次枚举$pd$，即令$T=pd$：

$$\sum _{T=1}^{\min (N,M)}\lfloor \frac{N}{T}\rfloor \times \lfloor \frac{M}{T}\rfloor \sum _{p\mid T}\mu (\frac{T}{p})$$

考虑后面的莫比乌斯函数的部分，可以通过元素贡献计算此函数，因为是枚举的是倍数，时间复杂度是调和级的，大概为$O(N\log N)$，前面枚举$T$可以用二维整除分块来做，总时间复杂度为$O(T(\sqrt{N}+N\log N))$。

---

P2522

跟上一题的思路一致，只不过多了一步容斥，假设设$f(n,m)={\sum }_{i=1}^n{\sum }_{j=1}^m[\gcd (i,j)=k]$，那么所求为：$f(b,d)-f(a-1,d)-f(b,c-1)+f(a-1,c-1)$。

---

CF1559E

同样的套路，反向枚举$d$，做背包即可。

---

P3172

同样的套路，反向枚举$d$，注意此题莫比乌斯函数前缀和范围较大，需要用杜教筛等高级筛法求莫比乌斯函数的前缀和。

---

P3327

求两个数乘积约数个数，有公式：

$$d(ij)=\sum _{x\mid i}\sum _{y\mid i}[\gcd (x,y)=1]$$

证明：

考虑$ij$的一个因子$K$，$K$中存在质数因子$p^c$，$i$存在质数因子$p^a$，$j$存在质数因子$p^b$。如果$c\le a$，那么在$i$中选择$c$，如果$c>a$，那么在$b$中选择$c-a$。我们发现对于同一质数的选择，不能在$i$和$j$中选择同一质数，因此$x$和$y$互质。

有了这个公式之后，对这个式子进行莫反即可。

---

P1829

看见对$\gcd (i,j)$计数，应枚举$d=\gcd$，从而构造出$[\gcd (i,1)=1]$。