数学/数论专题-学习笔记:狄利克雷卷积

1. 前言

狄利克雷卷积,是学习与继续探究 μ \mu μ 函数和 φ \varphi φ 函数的重要前提,因为这两个函数中有一些更好用的性质可以从狄利克雷卷积中得到,或者使用狄利克雷卷积更好的证明。

前置知识:一些基础的数论知识,可以直接看下去,遇到不明白的再查。

若无特殊说明,本文约定:数论函数如 f ( n ) f(n) f(n) 的定义域是正整数集,值域是数集。

2. 一些基础函数

这里给出下面会用到的几个基础函数:

单位函数 ε \varepsilon ε ε ( n ) = [ n = 1 ] \varepsilon(n )=[n=1] ε(n)=[n=1]

[ A ] [A] [A] 表示 A A A 为真时值为 1,否则为 0。

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

恒等函数 i d id id i d ( n ) = n id(n)=n id(n)=n

欧拉函数 φ \varphi φ φ ( n ) = ∑ i = 1 n [ gcd ⁡ ( i , n ) = 1 ] \varphi(n)=\sum_{i=1}^{n}[\gcd(i,n)=1] φ(n)=i=1n[gcd(i,n)=1]

约数和函数 σ 0 = ∑ d ∣ n 1 \sigma_0=\sum_{d \mid n}1 σ0=dn1,可扩展到 σ k = ∑ d ∣ n d k \sigma_k=\sum_{d \mid n}d^k σk=dndk

莫比乌斯函数 μ \mu μ,定义见后文,该函数的一种定义是由狄利克雷卷积得到的(或者说是性质),会简单介绍,另一种是直接定义,与狄利克雷卷积无关。

本文认为您应当对欧拉函数有一定的了解,至少知道欧拉函数的一些性质,如不知道任何性质,可以看我写的博客数学/数论专题-学习笔记:欧拉函数进行学习。

3. 积性函数

定义:对于一个数论函数 f ( n ) f(n) f(n),如果 gcd ⁡ ( i , j ) = 1 \gcd(i,j)=1 gcd(i,j)=1 时有 f ( i ) × f ( j ) = f ( i × j ) f(i) \times f(j)=f(i \times j) f(i)×f(j)=f(i×j),那么 f ( n ) f(n) f(n) 就是积性函数。

易证 ε , 1 , i d \varepsilon,1,id ε,1,id 时积性函数,而证明 φ ( n ) \varphi(n) φ(n) 是积性函数的过程在学习欧拉函数的过程中应该都学过,此处不赘述,莫比乌斯函数也是积性函数,但本篇博文不准备证明,因为我们目前无法得到其定义。

接下来证明 σ 0 \sigma_0 σ0 是积性函数:

gcd ⁡ ( n , m ) = 1 \gcd(n,m)=1 gcd(n,m)=1,则要证明 σ 0 ( n m ) = σ 0 ( n ) × σ 0 ( m ) \sigma_0(nm)=\sigma_0(n) \times \sigma_0(m) σ0(nm)=σ0(n)×σ0(m)

n , m n,m n,m 分解质因数,有 n = ∏ i = 1 s 1 p 1 i k 1 i , m = ∏ i = 1 s 2 p 2 i k 2 i n=\prod_{i=1}^{s1}p1_i^{k1_i},m=\prod_{i=1}^{s2}p2_i^{k2_i} n=i=1s1p1ik1i,m=i=1s2p2ik2i,因为 gcd ⁡ ( n , m ) = 1 \gcd(n,m)=1 gcd(n,m)=1,所以保证所有 p 1 i p1_i p1i p 2 i p2_i p2i 互不相同。

因此 σ 0 ( n ) = ∏ i = 1 s 1 ( 1 + k 1 i ) , σ 0 ( m ) = ∏ i = 1 s 2 ( 1 + k 2 i ) \sigma_0(n)=\prod_{i=1}^{s1}(1+k1_i),\sigma_0(m)=\prod_{i=1}^{s2}(1+k2_i) σ0(n)=i=1s1(1+k1i),σ0(m)=i=1s2(1+k2i)

如果定义 p = p 1 ∪ p 2 , k = k 1 ∪ k 2 p=p1 \cup p2,k=k1 \cup k2 p=p1p2,k=k1k2,并且 p i p_i pi k i k_i ki 在之前的 p 1 , p 2 , k 1 , k 2 p1,p2,k1,k2 p1,p2,k1,k2 中也是一一对应的,那么 n m = ∏ i = 1 s 1 + s 2 p i k i nm=\prod_{i=1}^{s1+s2}p_i^{k_i} nm=i=1s1+s2piki,则 σ 0 ( n m ) = ∏ i = 1 s 1 + s 2 ( 1 + k i ) = ∏ i = 1 s 1 ( 1 + k 1 i ) × ∏ i = 1 s 2 ( 1 + k 2 i ) = σ 0 ( n ) × σ 0 ( m ) \sigma_0(nm)=\prod_{i=1}^{s1+s2}(1+k_i)=\prod_{i=1}^{s1}(1+k1_i) \times \prod_{i=1}^{s2}(1+k2_i)=\sigma_0(n) \times \sigma_0(m) σ0(nm)=i=1s1+s2(1+ki)=i=1s1(1+k1i)×i=1s2(1+k2i)=σ0(n)×σ0(m),证毕。

4. 狄利克雷卷积

先定义两个数论函数 f ( n ) , g ( n ) f(n),g(n) f(n),g(n) 的加法与单个数论函数的数乘如下:

加法为逐项相加,即 ( f + g ) ( i ) = f ( i ) + g ( i ) (f+g)(i)=f(i)+g(i) (f+g)(i)=f(i)+g(i) ( x f ) ( n ) = x f ( n ) (xf)(n)=xf(n) (xf)(n)=xf(n)

其中 f + g f+g f+g 表示加法之后得到的函数, x f xf xf 表示数乘之后得到的函数。

注意,本文将数和函数符号都用小写字母表示。

两个数论函数 f = g f=g f=g 的充要条件是 ∀ n ∈ N + , f ( n ) = g ( n ) \forall n \in \text{N}_+,f(n)=g(n) nN+,f(n)=g(n)

定义两个数论函数的狄利克雷卷积 ∗ * 如下:

f ( n ) , g ( n ) f(n),g(n) f(n),g(n) 是两个数论函数,则其狄利克雷卷积为 t = f ∗ g t=f*g t=fg,满足:

t ( n ) = ∑ d ∣ n f ( n ) g ( n d ) t(n)=\sum_{d \mid n}f(n)g(\dfrac{n}{d}) t(n)=dnf(n)g(dn)

有一个等价形式是 t ( n ) = ∑ i × j = n f ( i ) g ( j ) t(n)=\sum_{i \times j=n}f(i)g(j) t(n)=i×j=nf(i)g(j),这个形式常用于证明某些性质。


狄利克雷卷积有如下的基础性质:

  1. 交换律 f ∗ g = g ∗ f f*g=g*f fg=gf
  2. 结合律 ( f ∗ g ) ∗ h = f ∗ ( g ∗ h ) (f*g)*h=f*(g*h) (fg)h=f(gh)
  3. 分配律 ( f + g ) ∗ h = f ∗ h + g ∗ h (f+g)*h=f*h+g*h (f+g)h=fh+gh
  4. 对数乘而言的结合律(我自己起的名字) ( x f ) ∗ g = x ( f ∗ g ) (xf)*g=x(f*g) (xf)g=x(fg)
  5. 单位函数 ϵ ∗ f = f \epsilon*f=f ϵf=f

这 5 个性质读者自证不难。


接下来看下狄利克雷卷积的逆元概念:对所有 f ( 1 ) ≠ 0 f(1) \ne 0 f(1)=0 的函数 f f f,都存在 f ∗ g = ϵ f * g=\epsilon fg=ϵ,称 g g g f f f 的逆元,记作 f − 1 f^{-1} f1

f ( 1 ) ≠ 0 f(1) \ne 0 f(1)=0 的原因是 ( f ∗ f − 1 ) ( 1 ) = ∑ i j = 1 f ( i ) f − 1 ( j ) = f ( 1 ) f − 1 ( 1 ) = 1 (f*f^{-1})(1)=\sum_{ij=1}f(i)f^{-1}(j)=f(1)f^{-1}(1)=1 (ff1)(1)=ij=1f(i)f1(j)=f(1)f1(1)=1,因此要满足这个式子,则 f ( 1 ) ≠ 0 f(1) \ne 0 f(1)=0

求逆元的式子如下:

n = 1 ⇒ g ( n ) = 1 f ( n ) n=1 \Rightarrow g(n)=\dfrac{1}{f(n)} n=1g(n)=f(n)1

n > 1 ⇒ g ( n ) = − 1 f ( 1 ) ∑ d ∣ n , d > 1 f ( d ) g ( n d ) n>1 \Rightarrow g(n)=-\dfrac{1}{f(1)}\sum_{d \mid n,d>1}f(d)g(\dfrac{n}{d}) n>1g(n)=f(1)1dn,d>1f(d)g(dn)

还有一种合在一起的简化版本:

g ( n ) = [ n = 1 ] − ∑ d ∣ n , d > 1 f ( d ) g ( n d ) g(n)=[n=1]-\sum_{d \mid n,d>1}f(d)g(\dfrac{n}{d}) g(n)=[n=1]dn,d>1f(d)g(dn)

怎么得到的不知道,但是我们可以验证这个 g ( n ) g(n) g(n) 是对的:

( f ∗ g ) ( n ) = ∑ d ∣ n f ( d ) g ( d n ) = f ( 1 ) g ( n ) + ∑ d ∣ n , d > 1 f ( d ) g ( n d ) = [ n = 1 ] (f*g)(n)=\sum_{d \mid n}f(d)g(\dfrac{d}{n})=f(1)g(n)+\sum_{d \mid n,d>1}f(d)g(\dfrac{n}{d})=[n=1] (fg)(n)=dnf(d)g(nd)=f(1)g(n)+dn,d>1f(d)g(dn)=[n=1]

( f ∗ g ) ( n ) = [ n = 1 ] = ϵ ( n ) (f*g)(n)=[n=1]=\epsilon(n) (fg)(n)=[n=1]=ϵ(n),所以 f ∗ g = ϵ f*g=\epsilon fg=ϵ,满足定义。


f ( n ) , g ( n ) f(n),g(n) f(n),g(n) 是两个积性函数,有:

  1. f ∗ g f*g fg 是积性函数。
  2. f − 1 f^{-1} f1 是积性函数。

其中上述性质 2 需要数学归纳法证明,过于繁琐,这里略去,有兴趣的可以参考参考资料。

性质 1 证明:

t = f ∗ g , gcd ⁡ ( n , m ) = 1 t=f*g,\gcd(n,m)=1 t=fg,gcd(n,m)=1,则要证 t ( n m ) = t ( n ) t ( m ) t(nm)=t(n)t(m) t(nm)=t(n)t(m)

t ( n m ) = ∑ d ∣ n m f ( d ) g ( n m d ) = ∑ i ∣ n , j ∣ m f ( i j ) g ( n i × m j ) t(nm)=\sum_{d \mid nm}f(d)g(\dfrac{nm}{d})=\sum_{i \mid n,j \mid m}f(ij)g(\dfrac{n}{i}\times\dfrac{m}{j}) t(nm)=dnmf(d)g(dnm)=in,jmf(ij)g(in×jm)

其中因为 gcd ⁡ ( n , m ) = 1 \gcd(n,m)=1 gcd(n,m)=1,所以 i , j i,j i,j 无共同质因子,且 gcd ⁡ ( n i , m j ) = 1 \gcd(\dfrac{n}{i},\dfrac{m}{j})=1 gcd(in,jm)=1

= ∑ i ∣ n , j ∣ m f ( i ) g ( n i ) f ( j ) g ( n j ) = [ ∑ i ∣ n f ( i ) g ( n i ) ] [ ∑ j ∣ m f ( j ) g ( m j ) ] = t ( n ) t ( m ) =\sum_{i \mid n,j\mid m}f(i)g(\dfrac{n}{i})f(j)g(\dfrac{n}{j})=[\sum_{i \mid n}f(i)g(\dfrac{n}{i})][\sum_{j \mid m}f(j)g(\dfrac{m}{j})]=t(n)t(m) =in,jmf(i)g(in)f(j)g(jn)=[inf(i)g(in)][jmf(j)g(jm)]=t(n)t(m)

证毕。


然后考虑莫比乌斯函数 μ \mu μ

其实接下来的内容会涉及到莫比乌斯函数与莫比乌斯反演的相关性质,这里简单介绍一下。

我们规定:常数函数 1 1 1 的逆元是 μ \mu μ,即 1 ∗ μ = ϵ 1*\mu=\epsilon 1μ=ϵ

假设 g = f ∗ 1 g=f*1 g=f1,因为 f = f ∗ ϵ f=f*\epsilon f=fϵ,那么 f = f ∗ 1 ∗ μ = g ∗ μ f=f*1*\mu=g*\mu f=f1μ=gμ

带入狄利克雷卷积的式子,也就是:

g ( n ) = ∑ d ∣ n f ( d ) g(n)=\sum_{d \mid n}f(d) g(n)=dnf(d),则 f ( n ) = ∑ d ∣ n g ( d ) μ ( n d ) f(n)=\sum_{d \mid n}g(d)\mu(\dfrac{n}{d}) f(n)=dng(d)μ(dn)

上述式子就是根据狄利克雷卷积得到的莫比乌斯反演的式子,另一个莫比乌斯反演的式子是另外一个类似卷积的东西得到的。

关于莫比乌斯反演这里只讲这么多,因为这里不是莫比乌斯反演的博文。

5. 总结

这里总结了本篇的重要定理与式子。


狄利克雷卷积定义: t = f ∗ g , t ( n ) = ∑ d ∣ n f ( n ) g ( n d ) t=f*g,t(n)=\sum_{d \mid n}f(n)g(\dfrac{n}{d}) t=fg,t(n)=dnf(n)g(dn)

狄利克雷卷积满足交换律,结合律,分配律,对数乘分配律, f ∗ ϵ = f f * \epsilon=f fϵ=f

狄利克雷卷积逆元 g ( n ) = [ n = 1 ] − ∑ d ∣ n , d > 1 f ( d ) g ( n d ) g(n)=[n=1]-\sum_{d \mid n,d>1}f(d)g(\dfrac{n}{d}) g(n)=[n=1]dn,d>1f(d)g(dn)

f , g f,g f,g 是两个积性函数,则 f ∗ g , f − 1 f*g,f^{-1} fg,f1 也是积性函数。

莫比乌斯函数 μ \mu μ 定义为 1 ∗ μ = ϵ 1*\mu=\epsilon 1μ=ϵ,即 1 1 1 的逆元。

g = f ∗ 1 ⇒ f = f ∗ 1 ∗ μ = g ∗ μ g=f*1 \Rightarrow f=f*1*\mu=g*\mu g=f1f=f1μ=gμ,即 g ( n ) = ∑ d ∣ n f ( d ) ⇒ f ( n ) = ∑ d ∣ n g ( d ) μ ( n d ) g(n)=\sum_{d \mid n}f(d) \Rightarrow f(n)=\sum_{d \mid n}g(d)\mu(\dfrac{n}{d}) g(n)=dnf(d)f(n)=dng(d)μ(dn)

6. 参考资料

  1. 铃悬的数学小讲堂——狄利克雷卷积与莫比乌斯反演 - 铃悬 的博客 - 洛谷博客
  2. 算法学习笔记(35): 狄利克雷卷积 - 知乎
  • 3
    点赞
  • 4
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值