莫比乌斯函数的由来及推导
网上搜莫比乌斯函数大多都没有写莫比乌斯函数的由来,而是直接写的定义式,这里对莫比乌斯函数进行推导。
学习了狄利克雷卷积后,我们来看以下几个常见的积性函数:
常函数: 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∣nf(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∣nf(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=1kf(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=p1a1p2a2...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∣nf(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∣df(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
得证
莫比乌斯函数的由来及推导相关推荐
- 读贾志鹏线性筛有感 (莫比乌斯函数的应用)
先拜大牛.感谢贾志鹏严谨的思维.以及简单清晰的论文描述. 一定要结合论文看.我只是提出我觉得关键的部分.论文在网上随处可见.贾志鹏线性筛. 开头两种线性筛的比较. 一种是传统的线性筛.时间复杂度为N* ...
- 正态分布的由来及推导
正态分布的由来及推导 一.正态分布 二.二项分布的近似计算 三.De Moivre-Laplace中心极限定理 四.最小二乘法与正态分布 五.基于独立性和旋转对称性的推导 六.Lindeberg-Lé ...
- 莫比乌斯函数和莫比乌斯反演
莫比乌斯函数 概念 莫比乌斯函数的定义如下: 若 d=1d=1d=1 ,则 μ(d)=1\mu(d)=1μ(d)=1 若 d=p1⋅p2⋯pk,pid=p_1\cdot p_2 \cdots p_k, ...
- 数论函数 - 莫比乌斯函数与莫比乌斯反演 - 基础杜教筛
原文链接http://www.cnblogs.com/zhouzhendong/p/8627380.html 省选后发现我数学好差.于是先从数论开始学习. 如果发现本文有任何错误,欢迎留言指正. 本文 ...
- 欧拉函数+狄利克雷卷积+莫比乌斯函数+莫比乌斯反演+整除分块+杜教筛
Powered by:NEFU AB-IN 文章目录 欧拉函数 狄利克雷卷积 莫比乌斯函数 莫比乌斯反演 P3455 [POI2007]ZAP-Queries 整除分块 P2522 [HAOI2011 ...
- 狄利克雷卷积与莫比乌斯函数
狄利克雷卷积与莫比乌斯函数 数论函数与积性函数 狄利克雷卷积 莫比乌斯函数 相关性质 基础问题与公式 若想使用莫比乌斯反演,熟练掌握狄利克雷卷积包括定义.记号以及相关的性质.证明等是非常有好处的. 数 ...
- 容斥原理与欧拉函数与莫比乌斯函数,狄利克雷卷积与莫比乌斯变换,反演
莫比乌斯函数可以看成是一种被内化了的容斥原理,许多数论上的结论定理根据容斥原理和数学归纳法可以推导出来,但是有关容斥原理的表达式的构造往往并不容易,运气不好很难找到,而莫比乌斯函数则是巧妙的把容斥原理 ...
- c莫比乌斯函数_莫比乌斯函数总结
莫比乌斯函数总结 性质:\(\sum_{d|n}\mu(d)=[n==1]\) 这个可以用组合数的性质来证,形象点的话就是杨辉三角. 因为恒等式:\(\sum_{i=0}^{n}(-1)^nC_{n} ...
- 莫比乌斯函数(数论)(转载+原创)
一.莫比乌斯函数(Möbius function)是指以下的函数: 在这里,λ(n)是刘维尔函数 莫比乌斯函数是一个数论函数,它同时也是一个积性函数(i.e.μ(ab) =μ(a)μ(b), a,b互 ...
最新文章
- 【Live555】live555源码详解(七):GenericMediaServer、RTSPServer、RTSPClient
- 网站设计师必备50教程
- windows自带的压缩,解压缩命令
- 打开word2016总是出现很抱歉,此功能看似中断需要修复。。问题解决办法
- hql可以使用distinct吗_输送食品可以使用白色PVC输送带吗?
- 【转】Azure基本核心服务计费
- 赞!苏州大学95后硕士一作发《Nature》!
- 云服务器,Linux上安装xampp来搭建node,api
- DEDE留言板调用导航的方法
- kafka jar包_和同事交流不会kafka怎么行,API奉上,不是大神也能编
- 浙江大学 PTA习题3.6 一元多项式的乘法与加法运算 (20分)
- Android 缓存处理和图片处理
- 机器学习中的数学系列
- c语言谭浩强第五版复习梳理1
- Win7/Win10如何设置局域网文件共享
- (附源码)ssm教培管理系统 毕业设计 230932
- GB28181协议错误码返回码整理
- matlab中画花瓣,matlab花瓣图的编程原理是什么,向天下大侠求解!!!!给力的? 爱问知识人...
- Mac 中如何解压.bin文件
- Windows开机加速实战之开机设置优…
热门文章
- java 新浪短链接_新浪短链接/腾讯短链接的API接口分享(含调用代码)
- python计算均方误差_在python中寻找线性回归的均方误差(使用scikit learn)
- Java人民币转大写工具类
- 今天的爱情多少钱一斤
- wed后端和java的区别_web前端和web后端的区别详细分析
- SDUT-程序设计基础-实验2-选择结构
- 优化问题-半正定规划(Semi-Definite Program, SDP)
- IP摄像头与传统模拟摄像头相比
- 海康相机 连接相机失败 未安装USB驱动或驱动无效
- Vscode 中vetur设置,让代码变得更美观