ICode9

精准搜索请尝试: 精确搜索
首页 > 其他分享> 文章详细

作为试验一份

2022-02-22 23:06:03  阅读:197  来源: 互联网

标签:right sum 作为 geqslant1 mu 试验 莫比 left 一份


随便贴一下而已

莫比乌斯函数

不知道是不是受那个带子的影响,莫比乌斯这个名字有种扭曲的美感,莫比乌斯函数也是这样。莫比乌斯函数 μ ( d ) \mu(d) μ(d) 是这样定义的:

∑ d ∣ m μ ( d ) = [ m = 1 ] \sum_{d| m}\mu(d)=[m=1] d∣m∑​μ(d)=[m=1]

其中方括号是《具体数学》中使用的艾弗森约定(我不知道这种形式是不是惯例),如果其中的命题为真的话这个函数为1,否则为0。换句话说也就是只有当 m=1 的时候为 1,其余都是 0。这种递归式并不是很好求值,尤其是累加记号中涉及到整除操作时。我们可以写出前几项看看:

m m m12345678910
μ ( m ) \mu(m) μ(m)1-1-10-11-1001

这东西看上去很随机还没啥用… 不过很快我们就会看到它存在的意义,而且这个东西还有一个很厉害的名字:

莫比乌斯反演

这种说法听起来就像某种奇幻背景的故事中的设定,以至于和这种名字放在一起没有什么违和感:

「神的不在场证明」!

「九天雷霆双脚蹬」!

「莫比乌斯反演」!

(好像也没有那么中二?)总之莫比乌斯函数遵循着如下的反演原理:

g ( m ) = ∑ d ∣ m f ( d ) ⇔ f ( m ) = ∑ d ∣ m μ ( d ) g ( m d ) g(m)=\sum_{d|m}f(d)\Leftrightarrow f(m)=\sum_{d|m}\mu(d)g\left(\frac md\right) g(m)=d∣m∑​f(d)⇔f(m)=d∣m∑​μ(d)g(dm​)

这看起来还是很有用的(虽然除了欧拉函数我也没见过多少这种形式的函数)要证明这个式子我们还是要花点技巧的:

∑ d ∣ m μ ( d ) g ( m d ) = ∑ d ∣ m μ ( m d ) g ( d ) ( ∑ d ∣ m f ( d ) = ∑ d ∣ m f ( m / d ) ) = ∑ d ∣ m μ ( m d ) ∑ k ∣ d f ( k ) = ∑ k ∣ m ∑ l ∣ ( m / k ) μ ( m k l ) f ( k ) ( ∑ d ∣ m ∑ k ∣ d a d , k = ∑ k ∣ m ∑ l ∣ ( m / k ) a k l , k ) = ∑ k ∣ m ∑ l ∣ ( m / k ) μ ( l ) f ( k ) = ∑ k ∣ m [ m / k = 1 ] f ( k ) = f ( m ) \begin{aligned} \sum_{d|m}\mu(d)g\left(\frac md\right)=&\sum_{d|m}\mu\left(\frac md\right)g(d)\qquad\qquad\left(\sum_{d|m}f(d)=\sum_{d|m}f(m/d)\right)\\ =&\sum_{d|m}\mu\left(\frac md\right)\sum_{k|d}f(k)\\ =&\sum_{k|m}\sum_{l|(m/k)}\mu\left(\frac m {kl}\right)f(k)\quad\left(\sum_{d|m}\sum_{k|d}a_{d,k}=\sum_{k|m}\sum_{l|(m/k)}a_{kl,k} \right)\\ =&\sum_{k|m}\sum_{l|(m/k)}\mu(l)f(k)\\ =&\sum_{k|m}[m/k=1]f(k)=f(m) \end{aligned} d∣m∑​μ(d)g(dm​)=====​d∣m∑​μ(dm​)g(d)⎝⎛​d∣m∑​f(d)=d∣m∑​f(m/d)⎠⎞​d∣m∑​μ(dm​)k∣d∑​f(k)k∣m∑​l∣(m/k)∑​μ(klm​)f(k)⎝⎛​d∣m∑​k∣d∑​ad,k​=k∣m∑​l∣(m/k)∑​akl,k​⎠⎞​k∣m∑​l∣(m/k)∑​μ(l)f(k)k∣m∑​[m/k=1]f(k)=f(m)​

可能需要补充的一点:右侧第二个式子的本质就是改名与整除的传递性。不过事实上吸引莫比乌斯发现(发明?用词取决于读者的数学观)这个函数的是另一个反演性质:

g ( x ) = ∑ d ⩾ 1 f ( x / d ) ⇔ f ( x ) = ∑ d ⩾ 1 μ ( d ) g ( x / d ) g(x)=\sum_{d\geqslant 1}f(x/d)\Leftrightarrow f(x)=\sum_{d\geqslant 1}\mu(d)g(x/d) g(x)=d⩾1∑​f(x/d)⇔f(x)=d⩾1∑​μ(d)g(x/d)

这个性质的证明也很简短:

∑ d ⩾ 1 μ ( d ) g ( x / d ) = ∑ d ⩾ 1 μ ( d ) ∑ k ⩾ 1 f ( x / d k ) = ∑ l ⩾ 1 ( x / l ) ∑ d , k ⩾ 1 μ ( d ) [ l = d k ] ( 这 种 换 元 . . . ) = ∑ m ⩾ 1 f ( x / m ) ∑ d ∣ m μ ( d ) = f ( x ) \begin{aligned} \sum_{d\geqslant1}\mu(d)g(x/d)=&\sum_{d\geqslant1}\mu(d)\sum_{k\geqslant1}f(x/dk)\\ =&\sum_{l\geqslant1}(x/l)\sum_{d,k\geqslant1}\mu(d)[l=dk]\quad(这种换元...)\\ =&\sum_{m\geqslant1}f(x/m)\sum_{d|m}\mu(d)=f(x) \end{aligned} d⩾1∑​μ(d)g(x/d)===​d⩾1∑​μ(d)k⩾1∑​f(x/dk)l⩾1∑​(x/l)d,k⩾1∑​μ(d)[l=dk](这种换元...)m⩾1∑​f(x/m)d∣m∑​μ(d)=f(x)​

不过我没什么头绪… 或许是饿了,不过更有可能是太缺乏练习了。习题什么的下回再写,我们下次再见(笑)

2022-02-19

标签:right,sum,作为,geqslant1,mu,试验,莫比,left,一份
来源: https://blog.csdn.net/weixin_44092155/article/details/123079443

本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享;
2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关;
3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关;
4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除;
5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

专注分享技术,共同学习,共同进步。侵权联系[81616952@qq.com]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有