网上搜莫比乌斯函数大多都没有写莫比乌斯函数的由来,而是直接写的定义式,这里对莫比乌斯函数进行推导。
学习了狄利克雷卷积后,我们来看以下几个常见的积性函数:

常函数: I ( n ) = 1 I(n)=1 I(n)=1

单位元: ξ ( n ) = [ n = 1 ] ( ξ ∗ f = f ) \xi(n)=[n=1] (\xi*f=f) ξ(n)=[n=1](ξ∗f=f)

莫比乌斯函数

对于一个数我们可以对其求逆,那么函数是否也可以求逆呢?答案是肯定的。我们看下面式子:

f − 1 ∗ f = ξ f^{-1}*f=\xi f−1∗f=ξ

f − 1 f^{-1} f−1即是我们要的 f f f的逆。

f − 1 ∗ f = x i f^{-1}*f=xi f−1∗f=xi

⇓ \Downarrow ⇓

∀ n : ∑ d ∣ n f ( d ) ∗ f − 1 ( n d ) = [ n = 1 ] \forall n:\sum_{d|n}{f(d)*f^{-1}(\frac{n}{d})=[n=1]} ∀n:∑d∣n​f(d)∗f−1(dn​)=[n=1]

若n=1

f − 1 ( 1 ) = 1 f ( 1 ) f^{-1}(1)=\frac{1}{f(1)} f−1(1)=f(1)1​

若n = ̸ 1 =\not1 ≠​1

f − 1 ( n ) ∗ f ( n ) = x i ( n ) f^{-1}(n)*f(n)=xi(n) f−1(n)∗f(n)=xi(n)

f − 1 ( n ) ∗ f ( n ) = 0 f^{-1}(n)*f(n)=0 f−1(n)∗f(n)=0

∑ d ∣ n f ( d ) ∗ f − 1 ( n d ) = 0 \sum_{d|n}{f(d)*f^{-1}(\frac{n}{d})}=0 d∣n∑​f(d)∗f−1(dn​)=0

将d=1从公式中提出来:
f ( 1 ) ∗ f − 1 ( n ) + ∑ d ∣ n f ( d ) ∗ f − 1 ( n d ) = 0 ( d ̸ = 1 ) f(1)*f^{-1}(n)+\sum_{d|n}{f(d)*f^{-1}(\frac{n}{d})}=0  (d\not=1) f(1)∗f−1(n)+d∣n∑​f(d)∗f−1(dn​)=0  (d̸​=1)

f − 1 ( n ) = − ∑ d ∣ n f ( d ) ∗ f − 1 ( n d ) f ( 1 ) ( d ̸ = 1 ) f^{-1}(n)=\frac{-\sum_{d|n}{f(d)*f^{-1}(\frac{n}{d})}}{f(1)}  (d\not=1) f−1(n)=f(1)−∑d∣n​f(d)∗f−1(dn​)​  (d̸​=1)

f f f的逆即为这个式子了

那么接下来我们试试把 f = I f=I f=I,n为质数的逆求出来试试:

f − 1 ( p ) = − f ( p ) ∗ f − 1 ( 1 ) f ( 1 ) = − 1 f^{-1}(p)=\frac{-f(p)*f^{-1}(1)}{f(1)}=-1 f−1(p)=f(1)−f(p)∗f−1(1)​=−1

看看 p 2 p^2 p2:

f − 1 ( p 2 ) = − ( f ( p 2 ) ∗ f − 1 ( 1 ) + f ( p ) ∗ f − 1 ( p ) ) = 0 f^{-1}(p^2)=-(f(p^2)*f^{-1}(1)+f(p)*f^{-1}(p))=0 f−1(p2)=−(f(p2)∗f−1(1)+f(p)∗f−1(p))=0

推广到 p k p^k pk

f − 1 ( p k ) = − ∑ i = 1 k f ( p i ) ∗ f − 1 ( p k − i ) f^{-1}(p^k)=-\sum_{i=1}^{k}{f(p^i)*f^{-1}(p^{k-i})} f−1(pk)=−∑i=1k​f(pi)∗f−1(pk−i)

不难看出:当k>1时, f ( p k ) f(p^k) f(pk)就等于0了

所以对于合数 n = p 1 a 1 p 2 a 2 . . . p k a k n=p_1^{a_1}p_2^{a_2}...p_k^{a^k} n=p1a1​​p2a2​​...pkak​ 有

f − 1 ( n ) = f − 1 ( p 1 a 1 ) f − 1 ( p 2 a 2 ) . . . f − 1 ( p k a k ) f^{-1}(n)=f^{-1}(p_1^{a_1})f^{-1}(p_2^{a_2})...f^{-1}(p_k^{a_k}) f−1(n)=f−1(p1a1​​)f−1(p2a2​​)...f−1(pkak​​) ( I − 1 I^{-1} I−1 为积性函数)

(积性函数的逆也是积性函数)
证明:
设 f f f为积性函数,有
f ( i ∗ j ) = f ( i ) ∗ f ( j ) f(i*j)=f(i)*f(j) f(i∗j)=f(i)∗f(j)
f ( i ∗ j ) ∗ f − 1 ( i ∗ j ) = ξ f(i*j)*f^{-1}(i*j)=\xi f(i∗j)∗f−1(i∗j)=ξ
f ( i ) ∗ f − 1 ( i ) ∗ f ( j ) ∗ f − 1 ( j ) = ξ f(i)*f^{-1}(i)*f(j)*f^{-1}(j)=\xi f(i)∗f−1(i)∗f(j)∗f−1(j)=ξ
联立得:
f ( i ∗ j ) ∗ f − 1 ( i ∗ j ) = f ( i ) ∗ f − 1 ( i ) ∗ f ( j ) ∗ f − 1 ( j ) f(i*j)*f^{-1}(i*j)=f(i)*f^{-1}(i)*f(j)*f^{-1}(j) f(i∗j)∗f−1(i∗j)=f(i)∗f−1(i)∗f(j)∗f−1(j)
f − 1 ( i ∗ j ) = f − 1 ( i ) ∗ f − 1 ( j ) f^{-1}(i*j)=f^{-1}(i)*f^{-1}(j) f−1(i∗j)=f−1(i)∗f−1(j)
证毕

如果 a n > 1 a_n>1 an​>1那么 f ( n ) = 0 f(n)=0 f(n)=0

设 f − 1 = μ ( I − 1 = μ ) f^{-1}=\mu(I^{-1}=\mu) f−1=μ(I−1=μ),则有:

μ = { 1 n = 1 ( − 1 ) k a 1 = a 2 = . . . = a k = 1 0 o t h e r w i s e \mu = \begin{cases} 1    n=1\\ (-1)^k   a^1=a^2=...=a^k=1\\ 0    otherwise\\ \end{cases} μ=⎩⎪⎨⎪⎧​1   n=1(−1)k   a1=a2=...=ak=10   otherwise​

这个便是我们所知的莫比乌斯函数了。

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

由我们的推导可知:

μ ∗ I = ξ ξ ∗ f = f \mu*I=\xi \xi*f=f μ∗I=ξ ξ∗f=f

设 f ∗ I = g f*I=g f∗I=g

f ∗ I ∗ μ = g ∗ μ f*I*\mu=g*\mu f∗I∗μ=g∗μ
公式

  • g ( n ) = ∑ d ∣ n f ( d ) ⇐ ⇒ f ( n ) = ∑ d ∣ n μ ( d ) g ( n d ) g(n)=\sum_{d|n}{f(d)}  \Leftarrow\Rightarrow f(n)=\sum_{d|n}{\mu(d)g(\frac{n}{d})} g(n)=∑d∣n​f(d) ⇐⇒ f(n)=∑d∣n​μ(d)g(dn​)
  • g ( n ) = ∑ n ∣ d f ( d ) ⇐ ⇒ f ( n ) = ∑ n ∣ d μ ( d ) g ( d n ) g(n)=\sum_{n|d}{f(d)}  \Leftarrow\Rightarrow f(n)=\sum_{n|d}{\mu(d)g(\frac{d}{n})} g(n)=∑n∣d​f(d) ⇐⇒ f(n)=∑n∣d​μ(d)g(nd​)

证明
方法一:逆推

  • ∑ d ∣ n μ ( d ) g ( n d ) = ∑ d ∣ n μ ( d ) ∑ x ∣ n d f ( x ) = ∑ x ∣ n f ( x ) ∑ d ∣ n x μ ( d ) = ∑ x ∣ n f ( x ) [ n x = 1 ] = f ( n ) \sum_{d|n}{\mu(d)g(\frac{n}{d})}=\sum_{d|n}{\mu(d)\sum_{x|\frac{n}{d}}{f(x)}=\sum_{x|n}{f(x)}\sum_{d|\frac{n}{x}}{\mu(d)}=\sum_{x|n}{f(x)[\frac{n}{x}=1]}}=f(n) d∣n∑​μ(d)g(dn​)=d∣n∑​μ(d)x∣dn​∑​f(x)=x∣n∑​f(x)d∣xn​∑​μ(d)=x∣n∑​f(x)[xn​=1]=f(n)

  • 第二个式子同理

方法二:卷积法(因为好多人不知道 μ \mu μ的性质所以不知道这点)
g = f ∗ I g=f*I g=f∗I
g ∗ μ = f ∗ I ∗ μ g*\mu=f*I*\mu g∗μ=f∗I∗μ
g ∗ μ = f ∗ ξ g*\mu=f*\xi g∗μ=f∗ξ
g ∗ μ = f g*\mu=f g∗μ=f
得证

莫比乌斯函数的由来及推导相关推荐

  1. 读贾志鹏线性筛有感 (莫比乌斯函数的应用)

    先拜大牛.感谢贾志鹏严谨的思维.以及简单清晰的论文描述. 一定要结合论文看.我只是提出我觉得关键的部分.论文在网上随处可见.贾志鹏线性筛. 开头两种线性筛的比较. 一种是传统的线性筛.时间复杂度为N* ...

  2. 正态分布的由来及推导

    正态分布的由来及推导 一.正态分布 二.二项分布的近似计算 三.De Moivre-Laplace中心极限定理 四.最小二乘法与正态分布 五.基于独立性和旋转对称性的推导 六.Lindeberg-Lé ...

  3. 莫比乌斯函数和莫比乌斯反演

    莫比乌斯函数 概念 莫比乌斯函数的定义如下: 若 d=1d=1d=1 ,则 μ(d)=1\mu(d)=1μ(d)=1 若 d=p1⋅p2⋯pk,pid=p_1\cdot p_2 \cdots p_k, ...

  4. 数论函数 - 莫比乌斯函数与莫比乌斯反演 - 基础杜教筛

    原文链接http://www.cnblogs.com/zhouzhendong/p/8627380.html 省选后发现我数学好差.于是先从数论开始学习. 如果发现本文有任何错误,欢迎留言指正. 本文 ...

  5. 欧拉函数+狄利克雷卷积+莫比乌斯函数+莫比乌斯反演+整除分块+杜教筛

    Powered by:NEFU AB-IN 文章目录 欧拉函数 狄利克雷卷积 莫比乌斯函数 莫比乌斯反演 P3455 [POI2007]ZAP-Queries 整除分块 P2522 [HAOI2011 ...

  6. 狄利克雷卷积与莫比乌斯函数

    狄利克雷卷积与莫比乌斯函数 数论函数与积性函数 狄利克雷卷积 莫比乌斯函数 相关性质 基础问题与公式 若想使用莫比乌斯反演,熟练掌握狄利克雷卷积包括定义.记号以及相关的性质.证明等是非常有好处的. 数 ...

  7. 容斥原理与欧拉函数与莫比乌斯函数,狄利克雷卷积与莫比乌斯变换,反演

    莫比乌斯函数可以看成是一种被内化了的容斥原理,许多数论上的结论定理根据容斥原理和数学归纳法可以推导出来,但是有关容斥原理的表达式的构造往往并不容易,运气不好很难找到,而莫比乌斯函数则是巧妙的把容斥原理 ...

  8. c莫比乌斯函数_莫比乌斯函数总结

    莫比乌斯函数总结 性质:\(\sum_{d|n}\mu(d)=[n==1]\) 这个可以用组合数的性质来证,形象点的话就是杨辉三角. 因为恒等式:\(\sum_{i=0}^{n}(-1)^nC_{n} ...

  9. 莫比乌斯函数(数论)(转载+原创)

    一.莫比乌斯函数(Möbius function)是指以下的函数: 在这里,λ(n)是刘维尔函数 莫比乌斯函数是一个数论函数,它同时也是一个积性函数(i.e.μ(ab) =μ(a)μ(b), a,b互 ...

最新文章

  1. 【Live555】live555源码详解(七):GenericMediaServer、RTSPServer、RTSPClient
  2. 网站设计师必备50教程
  3. windows自带的压缩,解压缩命令
  4. 打开word2016总是出现很抱歉,此功能看似中断需要修复。。问题解决办法
  5. hql可以使用distinct吗_输送食品可以使用白色PVC输送带吗?
  6. 【转】Azure基本核心服务计费
  7. 赞!苏州大学95后硕士一作发《Nature》!
  8. 云服务器,Linux上安装xampp来搭建node,api
  9. DEDE留言板调用导航的方法
  10. kafka jar包_和同事交流不会kafka怎么行,API奉上,不是大神也能编
  11. 浙江大学 PTA习题3.6 一元多项式的乘法与加法运算 (20分)
  12. Android 缓存处理和图片处理
  13. 机器学习中的数学系列
  14. c语言谭浩强第五版复习梳理1
  15. Win7/Win10如何设置局域网文件共享
  16. (附源码)ssm教培管理系统 毕业设计 230932
  17. GB28181协议错误码返回码整理
  18. matlab中画花瓣,matlab花瓣图的编程原理是什么,向天下大侠求解!!!!给力的? 爱问知识人...
  19. Mac 中如何解压.bin文件
  20. Windows开机加速实战之开机设置优…

热门文章

  1. java 新浪短链接_新浪短链接/腾讯短链接的API接口分享(含调用代码)
  2. python计算均方误差_在python中寻找线性回归的均方误差(使用scikit learn)
  3. Java人民币转大写工具类
  4. 今天的爱情多少钱一斤
  5. wed后端和java的区别_web前端和web后端的区别详细分析
  6. SDUT-程序设计基础-实验2-选择结构
  7. 优化问题-半正定规划(Semi-Definite Program, SDP)
  8. IP摄像头与传统模拟摄像头相比
  9. 海康相机 连接相机失败 未安装USB驱动或驱动无效
  10. Vscode 中vetur设置,让代码变得更美观