莫比乌斯反演

狄利克雷卷积

定义:(fg)(n)=d|nf(d)g(n/d)

很显然满足交换律和结合律。

积性函数

为积性函数的有:

I(n) (或1(n) ),恒等于1,所以叫恒等函数

ϵ(n) (或者e(n) ),当且仅当 n=1 时,其值为 1,否则为 0,其满足(ef=f)(因此为狄利克雷卷积的单位元)

id(n)=n 为单位函数。

以上为完全积性函数。

完全积性函数:对于任意整数 ab 满足 f(ab)=f(a)f(b)

以及:

φ(n) ,欧拉函数,小于 n 的整数中与 n 互质的数的个数。

μ(n) ,莫比乌斯函数,接下来我们重点讲,暂且不介绍。

积性函数:对于两个整数 a,b ,满足 (a,b)=1 ,则 f(ab)=f(a)f(b)

虽然没有完全积性函数优美,但是这很好吧,这可以吧。(

然后研究一下这个积性函数的性质。

  • 积性函数 f,总满足 f(1)=1

这个易证了,f(1)=f(1)f(1)

  • 两积性函数之积为积性函数。

这个稍微难一点。

证明:

定义两个积性函数 f,g ,其卷积为 G=fg.

任取两个互质的数 a,b

G(a)G(b)

=d|af(d)g(b/d)t|bf(t)g(b/t)

=d|at|bf(d)g(a/d)f(t)g(b/t)

=dt|abf(dt)g(ab/dt)

=G(ab)

Q.E.D.

  • 积性函数的逆也是积性函数

归纳证明,就不证明了

莫比乌斯函数

引入

对于两个函数 f,F,满足 F(n)=d|n(1f(d))

等价于 F=If,然后有 f=I1F

我们把 I1 称为 μ 莫比乌斯函数。

也就有 f=μF

定义:

o

然后有个性质:

  • (μ1)=e

从定义出发易证。互逆的两个函数卷起来是单位元。

  • φ1=id ,然后 φ=μid

φ1=id,且 μ1=e

φ1μ=idμφ=μid

然后证一下 φ1=id

想了解可以参考 OI wiki

莫比乌斯反演

进入正题。

  • 嵌入式莫比乌斯反演

μ1=ed|nμ(d)=[n=1]

因为 [n|m][n/m=1]=[n=m]

所以有 [n|m]d|(n/m)μ(d)=[n=m] (因为只有当n=1的时候这个玩意才满足)

可以这么转换。

  • d|(i,j)μ(d)=[(i,j)=1]

因为 d|(i,j)μ(d)=e(gcd(i,j)),易证

然后你肯定是要会算莫比乌斯函数的,开筛!

这个我们之前的博客中有,于是不多说了。

变换形式

  • F(n)=d|nf(d)f(n)=d|nμ(d)F(n/d)

本质还是

F=1ff=FI1f=Fμ

数论分块

用来计算形如 i=1nf(i)g(ni) 的和式。

我们再单独来讲这个 数论分块

我们推个式子:

ans=i=1nj=1m[(i,j)=1]

=i=1nj=1md|(i,j)μ(d)

=d=1min{n,m}μ(d)d|ind|jm1

=d=1min{n,m}μ(d)n/dm/d

这个式子我们可以 O(n) 的算。

接下来我们用数论分块处理,达到 O(n)

总之莫反的题就是分为反演和分块,学懂了还是挺套路的。

posted @   cbdsopa  阅读(43)  评论(0编辑  收藏  举报
相关博文:
阅读排行:
· 物流快递公司核心技术能力-地址解析分单基础技术分享
· 单线程的Redis速度为什么快?
· 展开说说关于C#中ORM框架的用法!
· Pantheons:用 TypeScript 打造主流大模型对话的一站式集成库
· SQL Server 2025 AI相关能力初探
点击右上角即可分享
微信分享提示