莫比乌斯反演

莫比乌斯反演

积性定义: φ(ab)=φ(a)φ(b) (a,b互质时)

具有积性的函数叫积性函数。

如果 a,b 不互质时仍然是积性,那就是完全积性函数。

常见积性函数

  • 单位函数:ε(n)=[n=1](完全积性)
  • 恒等函数:idk(n)=nk(完全积性),当 k=1 时,简记为 id(n)=n
  • 常数函数:1(n)=1(完全积性)
  • 除数函数:σk(n)=d|ndkσ0(n) 即因数个数,简记为 d(n)σ1(n) 即因数之和,简记为 σ(n)
  • 欧拉函数:φ(n)=d=1n[gcd(d,n)=1]
  • 莫比乌斯函数:μ(x)={1  (x=1)(1)k  (x没有平方数因子,且x的质因子个数为k)0(x有平方数因子)

(注:莫比乌斯函数的定义是为了满足下文的某些奇妙性质。)

这些积性函数都可以用线性欧拉筛 O(n) 地筛出。

可以证明,积性函数在 1 处的取值为 1

Dirchlet 卷积

定义

(fg)(n)=d|nf(d)g(nd)

奇妙性质

  • φ1=id
  • μ1=ε
  • μid=φ
  • 11=d
  • id1=σ
  • idk1=σk

前两个换一种写法就是:dnφ(d)=n, dnμ(d)=[n=1]。证明比较骚。

第三个相当于第一个式子两边 μ,又因为显然有 φε=φ,化简出来就是第三个式子。

最后三个用定义可得到,比较好理解。

$\small{个人感觉 ε 相当于函数中的单位元,μ 相当于 1 的逆元}$

莫比乌斯反演

g(n)=i|nf(i)f(n)=i|nμ(ni)g(i)

不过这玩意一般不用。

题目与套路

莫反的套路比较死。基本上看到有 gcd,lcm[k=1] 之类的玩意就可以无脑套上莫反。

【例题】 求:

i=1nj=1mlcm(i,j)

(1)i=1nj=1mlcm(i,j)=i=1nj=1mgcd(i,j)igjg(2)=dndin/djm/dij[gcd(i,j)=1]先枚举gcd,并要求i,j互质,表示的是gcd多少倍(3)=dndin/djm/dijk|gcd(i,j)μ(k)奇妙的性质第二条(4)=dnkkdmin(n,m)μ(k)dk2in/kdijm/kdj类比第二步

S(n)=ini=(n+1)n/2,T=kd,那么上式等于:

(5)=dnkmin(n,m)/dμ(k)dk2S(n/kd)S(n/kd)(6)=Tnk|Tμ(k)TkS(n/T)S(n/T)(7)=TnS(n/T)S(n/T)Tk|Tμ(k)k

后面是积性函数,用线性筛 O(n),前面整除分块。

积性证明:

f(pq)=xy|pqμ(xy)xy (p,q互质)

f(p)=x|pμ(x)y,f(q)=y|qμ(y)y

因为 p,q互质,所以 x,y 整除 p,q 中的其中一个。

钦定 x | p, y | q(注意,如果 xyp 则,等价于把 xy 看成新的 x, 把 1 看作 y, 这种情况会在之前统计过)

f(pq)=xy|pqμ(xy)xy=x|pμ(x)xy|pμ(y)y=f(p)f(q)

对于 f(j2p)=f(jp)+μ(j2x)j2x=f(jp),其中 p 是质数。

f(p)=μ(1)1+μ(p)p,其中 p 是质数。

f(1)=μ(1)1=1

5bp94xos.png (300×424)

posted @   花子の水晶植轮daisuki  阅读(9)  评论(0编辑  收藏  举报
相关博文:
阅读排行:
· 阿里最新开源QwQ-32B,效果媲美deepseek-r1满血版,部署成本又又又降低了!
· 开源Multi-agent AI智能体框架aevatar.ai,欢迎大家贡献代码
· Manus重磅发布:全球首款通用AI代理技术深度解析与实战指南
· 被坑几百块钱后,我竟然真的恢复了删除的微信聊天记录!
· AI技术革命,工作效率10个最佳AI工具
https://blog-static.cnblogs.com/files/zouwangblog/mouse-click.js
点击右上角即可分享
微信分享提示