一个和莫比乌斯反演有关的东西

一个和莫比乌斯反演有关的东西

数论柿子题的发现 - 博客 - masterhuang的博客

以下会用到两个作用于函数的算子 Σ,Δ 的定义分别是 Σf(n)=Σf(n1)+f(n)(前缀和),Δf(n)=f(n)f(n1)(差分)。这个记号和《具体数学》上有限微积分的符号不是一样的。

对于

i=1nj=1n[gcd(i,j)=1]

这个式子,我们知道一种做法是先折半然后就能写成 φ 的前缀和:(1 是因为 gcd(1,1)=1

2(i=1nj=1i[gcd(i,j)=1])1=2i=1nφ(i)1

还有一种做法是直接用莫比乌斯函数进行反演:

i=1nj=1nd|gcd(i,j)μ(d)=d=1nμ(d)n/d2

但是为什么这两个式子是相等的呢?我们考虑我们有一个经常拿来用的引理:若 fg=h,则 Σh(n)=i=1nf(i)Σg(n/i)

如果我们有一个函数 g(n) 满足其前缀和 Σg(n)=n2,那么我们就有一点它们相等的希望。而显然 g(n)=n2(n1)2=2n1,于是就知道原式等于

Σ(gμ)(n)

又众所周知狄利克雷卷积存在对加法的分配律,以及常数可以自由进入 号,我们将原式改成

(2Σ(Idμ)Σ(Iμ))(n)

这些都是经典结论了,Idμ=φ,Iμ=ϵ,也就是说原式等于

(2ΣφΣϵ)(n)=(2ΣφI)(n)

ϵ 的前缀和恰好为 I,所以这就等于上文那个 φ 的前缀和的形式。其实它们是有联系的!

以上过程成立,都是因为以下这个重要的式子成立,真是太有趣了。还可以用这个技巧优化更多此前没有注意过的题目。

f(n)=i=1ng(i)h(n/i)Δf=gΔh

posted @   caijianhong  阅读(18)  评论(0编辑  收藏  举报
相关博文:
阅读排行:
· 10年+ .NET Coder 心语 ── 封装的思维:从隐藏、稳定开始理解其本质意义
· 地球OL攻略 —— 某应届生求职总结
· 提示词工程——AI应用必不可少的技术
· Open-Sora 2.0 重磅开源!
· 周边上新:园子的第一款马克杯温暖上架
点击右上角即可分享
微信分享提示