莫比乌斯反演

最新推荐文章于 2022-04-15 10:14:32 发布

置顶有钱哥哥家的

最新推荐文章于 2022-04-15 10:14:32 发布

阅读量284

点赞数

分类专栏：知识点算法文章标签：莫比乌斯反演

知识点同时被 2 个专栏收录

13 篇文章 0 订阅

订阅专栏

算法

7 篇文章 0 订阅

订阅专栏

莫比乌斯反演

（PS：在评论区中众多dalao的催促下，我认真的写了三天三夜写完了这篇杜教筛，保证是精品！）

前言

（这大概是我第一次写学习笔记吧OvO）
可能每一个刚开始接触莫比乌斯反演的OIer，起初都会厌恶这个神奇的东西。(我也一样233)每一个人厌恶的原因有许多，可能是这个烦人的式子，也可能仅仅只是因为不理解μμ函数而感到不爽。当然，莫比乌斯反演有一个小小的预备知识：整除分块
那么我们先从莫比乌斯反演中最基础的莫比乌斯函数μ开始说起：

莫比乌斯函数

首先，我们可以先明确一点，莫比乌斯函数并不是什么很高大上的东西，它其实只是一个由容斥系数所构成的函数。μ(d)μ(d)的定义是:

当d=1d=1时，μ(d)=1μ(d)=1；
当d=Πki=1pid=Πi=1kpi且pipi为互异素数时，μ(d)=(−1)kμ(d)=(−1)k。(说直白点，就是dd分解质因数后，没有幂次大于平方的质因子，此时函数值根据分解的个数决定)；
只要当dd含有任何质因子的幂次大于等于2，则函数值为0.

当然，莫比乌斯函数也有很多有趣的性质：

对于任意正整数nn，∑d|nμ(d)=[n=1]∑d|nμ(d)=[n=1]。（[n=1][n=1]表示只有当n=1n=1成立时，返回值为11；否则，值为00；(这个就是用μμ是容斥系数的性质可以证明)(PS:这一条性质是莫比乌斯反演中最常用的)
对于任意正整数nn，∑d|nμ(d)d=ϕ(n)n∑d|nμ(d)d=ϕ(n)n。（这个性质很奇妙，它把欧拉函数和莫比乌斯函数结合起来，或许我之后写杜教筛的学习笔记时会去证明吧）

程序实现并不难，我们可以在线性筛素数的程序上略作修改，便可以筛出μμ函数。
那我还是给一段线筛的代码吧

void get_mu(int n)
{
    mu[1]=1;
    for(int i=2;i<=n;i++)
    {
        if(!vis[i]){prim[++cnt]=i;mu[i]=-1;}
        for(int j=1;j<=cnt&&prim[j]*i<=n;j++)
        {
            vis[prim[j]*i]=1;
            if(i%prim[j]==0)break;
            else mu[i*prim[j]]=-mu[i];
        }
    }
 }

那么，莫比乌斯函数就这么告一段落了。

莫比乌斯反演

解决完莫比乌斯函数的问题后，我们便迎来了重头戏莫比乌斯反演
定理：F(n)和f(n)是定义在非负整数集合上的两个函数，并且满足条件：
F(n)=∑d|nf(d)F(n)=∑d|nf(d)

那么存在一个结论：
f(n)=∑d|nμ(d)F(⌊nd⌋)f(n)=∑d|nμ(d)F(⌊nd⌋)

这个定理就称作莫比乌斯反演定理。
莫比乌斯反演的证明主要有两种方式，其中一种就是通过定义来证明；另外一种，我则是会在杜教筛中提到(利用狄利克雷卷积)。那么我先来说一说第一种证明方法：
∑d|nμ(d)F(⌊nd⌋)=∑d|nμ(d)∑i|⌊nd⌋f(i)∑d|nμ(d)F(⌊nd⌋)=∑d|nμ(d)∑i|⌊nd⌋f(i)

=∑i|nf(i)∑d|⌊ni⌋μ(d)=f(n)=∑i|nf(i)∑d|⌊ni⌋μ(d)=f(n)

(PS:如果不知道最后一步怎么来的，可以再去看性质一，至于和式的变换，就自己脑补一下吧)
当然，莫比乌斯反演有另外的一种形式，当F(n)F(n)和f(n)f(n)满足：
F(n)=∑n|df(d)F(n)=∑n|df(d)

可以推出：
f(n)=∑n|dμ(dn)F(d)f(n)=∑n|dμ(dn)F(d)
感觉这个式子，可能在莫比乌斯反演中更加好用。