当前位置: 首页 > news >正文

狄利克雷卷积

lgのhhxxjy

一、定义

狄利克雷卷积是定义在数论函数(即定义域为正整数,值域为复数域的函数)之间的一种二元运算。

f(n)f(n)f(n)g(n)g(n)g(n) 是两个数论函数,它们的狄利克雷卷积是一个新的数论函数,记作 (f∗g)(n)(f * g)(n)(fg)(n),定义为:

(f∗g)(n)=∑d∣nf(d)⋅g(nd)(f * g)(n) = \sum_{d \mid n} f(d) \cdot g\left(\frac{n}{d}\right)(fg)(n)=dnf(d)g(dn)

其中,求和是对 nnn 的所有正因数 ddd 进行的。

这个定义也可以等价地写为(令 k=ndk = \frac{n}{d}k=dn):
(f∗g)(n)=∑ab=nf(a)⋅g(b)(f * g)(n) = \sum_{ab = n} f(a) \cdot g(b)(fg)(n)=ab=nf(a)g(b)
即对所有满足 a×b=na \times b = na×b=n 的正整数对 (a,b)(a, b)(a,b) 求和。


二、为什么叫“卷积”?

这个名字来源于卷积(Convolution) 的通用概念。在数学中,卷积是一种通过两个函数生成第三个函数的数学算子,它表示一个函数在经过另一个函数“加权”和“翻转平移”后的叠加效果。

狄利克雷卷积的具体形式 ∑d∣nf(d)g(n/d)\sum_{d|n} f(d)g(n/d)dnf(d)g(n/d) 完美地符合了这种“加权求和”的思想:

  • 它“翻转”了除数 ddd(将其变为 n/dn/dn/d)。
  • 它对所有可能的“平移”(即所有除数对)进行了“加权”(f(d)f(d)f(d)g(n/d)g(n/d)g(n/d) 的乘积)求和。

因此,它是离散情形下(定义在正整数格点上)的一种典型卷积。


三、基本性质

狄利克雷卷积赋予所有数论函数的集合一个非常优美的代数结构。

1. 交换律 (Commutativity)

f∗g=g∗ff * g = g * ffg=gf
证明
(f∗g)(n)=∑ab=nf(a)g(b)=∑ba=ng(b)f(a)=(g∗f)(n)(f * g)(n) = \sum_{ab=n} f(a)g(b) = \sum_{ba=n} g(b)f(a) = (g * f)(n)(fg)(n)=ab=nf(a)g(b)=ba=ng(b)f(a)=(gf)(n)
本质上是因为乘法 ab=nab=nab=n 具有交换性。

2. 结合律 (Associativity)

(f∗g)∗h=f∗(g∗h)(f * g) * h = f * (g * h)(fg)h=f(gh)
证明
((f∗g)∗h)(n)=∑ab=n(f∗g)(a)⋅h(b)=∑ab=n(∑cd=af(c)g(d))h(b)=∑cdb=nf(c)g(d)h(b)((f * g) * h)(n) = \sum_{ab=n} (f*g)(a) \cdot h(b) = \sum_{ab=n} \left( \sum_{cd=a} f(c)g(d) \right) h(b) = \sum_{cdb=n} f(c)g(d)h(b)((fg)h)(n)=ab=n(fg)(a)h(b)=ab=n(cd=af(c)g(d))h(b)=cdb=nf(c)g(d)h(b)

(f∗(g∗h))(n)=∑ac=nf(a)⋅(g∗h)(c)=∑ac=nf(a)(∑bd=cg(b)h(d))=∑abd=nf(a)g(b)h(d)(f * (g * h))(n) = \sum_{ac=n} f(a) \cdot (g*h)(c) = \sum_{ac=n} f(a) \left( \sum_{bd=c} g(b)h(d) \right) = \sum_{abd=n} f(a)g(b)h(d)(f(gh))(n)=ac=nf(a)(gh)(c)=ac=nf(a)(bd=cg(b)h(d))=abd=nf(a)g(b)h(d)

由于 c,d,bc, d, bc,d,b 都是求和变量,两种形式等价。

3. 分配律 (Distributivity)

f∗(g+h)=(f∗g)+(f∗h)f * (g + h) = (f * g) + (f * h)f(g+h)=(fg)+(fh)
证明:根据求和的线性性质,显然成立。

4. 单位元 (Identity Element)

存在一个特殊的数论函数 ε(n)\varepsilon(n)ε(n),称为卷积单位元,使得对于任何数论函数 fff,都有:
f∗ε=ε∗f=ff * \varepsilon = \varepsilon * f = ffε=εf=f

这个函数就是:
ε(n)={1if n=10if n>1\varepsilon(n) = \begin{cases} 1 & \text{if } n = 1 \\ 0 & \text{if } n > 1 \end{cases}ε(n)={10if n=1if n>1
它也被称为单位函数

证明
(f∗ε)(n)=∑d∣nf(d)ε(n/d)(f * \varepsilon)(n) = \sum_{d|n} f(d) \varepsilon(n/d)(fε)(n)=dnf(d)ε(n/d)
只有当 n/d=1n/d = 1n/d=1,即 d=nd = nd=n 时,ε(n/d)=1\varepsilon(n/d)=1ε(n/d)=1,其他项均为 0。
所以 (f∗ε)(n)=f(n)⋅1=f(n)(f * \varepsilon)(n) = f(n) \cdot 1 = f(n)(fε)(n)=f(n)1=f(n)


四、逆元 (Inverse)

并非所有数论函数都有逆元。如果一个数论函数 f(n)f(n)f(n) 满足 f(1)≠0f(1) \neq 0f(1)=0,那么存在唯一的数论函数 g(n)g(n)g(n),使得:
f∗g=εf * g = \varepsilonfg=ε
这个函数 ggg 称为 fff狄利克雷逆元(Dirichlet Inverse),记作 f−1f^{-1}f1

逆元可以通过递归定义来求出:

  1. g(1)=1f(1)g(1) = \frac{1}{f(1)}g(1)=f(1)1 (这就是为什么要求 f(1)≠0f(1) \neq 0f(1)=0
  2. 对于 n>1n > 1n>1
    g(n)=−1f(1)∑d∣nd>1f(d)g(nd)g(n) = -\frac{1}{f(1)} \sum_{\substack{d|n \\ d>1}} f(d) g\left(\frac{n}{d}\right)g(n)=f(1)1dnd>1f(d)g(dn)

例子:莫比乌斯函数 μ(n)\mu(n)μ(n) 就是常数函数 1(n)=1{\bf 1}(n)=11(n)=1 的逆元,即 μ=1−1\mu = {\bf 1}^{-1}μ=11。这直接导致了著名的莫比乌斯反演公式


五、几个重要的例子

理解狄利克雷卷积的最佳方式就是看例子。我们定义几个常见的数论函数:

  • 常数函数1(n)=1{\bf 1}(n) = 11(n)=1 (对于所有 nnn
  • 恒等函数id(n)=n{\bf id}(n) = nid(n)=n
  • 单位函数ε(n)\varepsilon(n)ε(n) (见上文)
  • 欧拉函数φ(n)\varphi(n)φ(n) (小于等于 nnn 且与 nnn 互质的正整数的个数)
  • 除数函数d(n)d(n)d(n)nnn 的正因数的个数)
  • 莫比乌斯函数μ(n)\mu(n)μ(n)

现在来看一些重要的卷积恒等式:

  1. 1{\bf 1}11{\bf 1}1 的卷积
    (1∗1)(n)=∑d∣n1⋅1=∑d∣n1=d(n)({\bf 1} * {\bf 1})(n) = \sum_{d|n} 1 \cdot 1 = \sum_{d|n} 1 = d(n)(11)(n)=dn11=dn1=d(n)
    结论:除数函数 d=1∗1d = {\bf 1} * {\bf 1}d=11

  2. 1{\bf 1}1id{\bf id}id 的卷积
    (1∗id)(n)=∑d∣n1⋅nd=∑d∣nnd=∑d∣nd=σ(n)({\bf 1} * {\bf id})(n) = \sum_{d|n} 1 \cdot \frac{n}{d} = \sum_{d|n} \frac{n}{d} = \sum_{d|n} d = \sigma(n)(1id)(n)=dn1dn=dndn=dnd=σ(n)
    结论:因数和函数 σ=1∗id\sigma = {\bf 1} * {\bf id}σ=1id

  3. 1{\bf 1}1φ\varphiφ 的卷积(最著名的恒等式之一)
    (1∗φ)(n)=∑d∣nφ(d)({\bf 1} * \varphi)(n) = \sum_{d|n} \varphi(d)(1φ)(n)=dnφ(d)
    这个和有一个美妙的结论:它等于 nnn 本身。
    ∑d∣nφ(d)=n\sum_{d|n} \varphi(d) = ndnφ(d)=n
    证明思路:考虑分母为 nnn 的所有最简分数,将其分组,每一组的分母是 nnn 的某个约数 ddd
    结论1∗φ=id{\bf 1} * \varphi = {\bf id}1φ=id

  4. 1{\bf 1}1μ\muμ 的卷积(逆元的定义)
    (1∗μ)(n)=∑d∣n1⋅μ(d)=∑d∣nμ(d)({\bf 1} * \mu)(n) = \sum_{d|n} 1 \cdot \mu(d) = \sum_{d|n} \mu(d)(1μ)(n)=dn1μ(d)=dnμ(d)
    根据逆元的定义,因为 μ=1−1\mu = {\bf 1}^{-1}μ=11,所以:
    1∗μ=ε{\bf 1} * \mu = \varepsilon1μ=ε
    ∑d∣nμ(d)=ε(n)={1n=10n>1\sum_{d|n} \mu(d) = \varepsilon(n) = \begin{cases} 1 & n=1 \\ 0 & n>1 \end{cases}dnμ(d)=ε(n)={10n=1n>1
    这个性质是莫比乌斯反演公式的核心。


六、莫比乌斯反演公式 (Möbius Inversion Formula)

狄利克雷卷积最著名的应用就是导出莫比乌斯反演公式。

定理:假设 f(n)f(n)f(n)g(n)g(n)g(n) 是两个数论函数,如果有:
f(n)=∑d∣ng(d)f(n) = \sum_{d|n} g(d)f(n)=dng(d)
那么可以有:
g(n)=∑d∣nμ(d)f(nd)g(n) = \sum_{d|n} \mu(d) f\left(\frac{n}{d}\right)g(n)=dnμ(d)f(dn)

用卷积语言证明,极其简洁
原式 f(n)=∑d∣ng(d)f(n) = \sum_{d|n} g(d)f(n)=dng(d) 其实就是 f=1∗gf = {\bf 1} * gf=1g
两边同时卷积上 1{\bf 1}1 的逆元 μ\muμ
μ∗f=μ∗(1∗g)\mu * f = \mu * ({\bf 1} * g)μf=μ(1g)
根据结合律:
μ∗f=(μ∗1)∗g\mu * f = (\mu * {\bf 1}) * gμf=(μ1)g
因为 μ∗1=ε\mu * {\bf 1} = \varepsilonμ1=ε,所以:
μ∗f=ε∗g=g\mu * f = \varepsilon * g = gμf=εg=g
g=μ∗fg = \mu * fg=μf,写开就是:
g(n)=∑d∣nμ(d)f(nd)g(n) = \sum_{d|n} \mu(d) f(\frac{n}{d})g(n)=dnμ(d)f(dn)
证毕。

卷积的代数性质让这个证明变得清晰而强大。


七、总结与重要性

狄利克雷卷积之所以如此重要,是因为它:

  1. 提供了强大的代数工具:它将数论函数的研究转化为一个代数系统(实际上形成一个交换环),可以使用交换律、结合律、逆元等工具进行运算和证明。
  2. 统一了数论函数的关系:许多看似不相关的数论函数和恒等式(如 ∑d∣nφ(d)=n\sum_{d|n}\varphi(d)=ndnφ(d)=n, ∑d∣nμ(d)=0\sum_{d|n}\mu(d)=0dnμ(d)=0)都可以用卷积优雅地表示和推导。
  3. 是更高级理论的基石:它是理解狄利克雷级数(特别是黎曼ζ函数)的基础。卷积关系在级数上对应着乘法关系:∑n=1∞(f∗g)(n)ns=(∑n=1∞f(n)ns)(∑n=1∞g(n)ns)\sum_{n=1}^{\infty}\frac{(f*g)(n)}{n^s} = \left( \sum_{n=1}^{\infty}\frac{f(n)}{n^s} \right) \left( \sum_{n=1}^{\infty}\frac{g(n)}{n^s} \right)n=1ns(fg)(n)=(n=1nsf(n))(n=1nsg(n))
  4. 简化了计算和证明:如莫比乌斯反演的证明所示,卷积语言将复杂的求和问题转化为简洁的代数操作。

希望这份详解能帮助你彻底理解狄利克雷卷积这个优美而有力的数学工具!

http://www.dtcms.com/a/410997.html

相关文章:

  • 沈阳 教育 公司 网站建设种子搜索网站怎么做的
  • 矩阵置零--leetcode
  • 删除iCloud中的照片但保留iPhone上的照片的两种方法
  • 巴中市住房和城乡建设局官方网站广告牌子设计图片
  • 找学校的网站塔城网站seo
  • oracle linux 10 +pg18 源码安装要点
  • 如何快速排查服务器宕机故障
  • 分sheet写入excel
  • 怎么把网站黑了传统文化网站建设
  • 【开题答辩全过程】以 SpringBootVue的旅游租车管理系统为例,包含答辩的问题和答案
  • Windows配置PicGo+Gitee图床——解决你的MarkDown笔记分享无图片的问题
  • 网站排名alexawordpress 商城安全
  • 【开题答辩过程】以《分布式菌菇销售系统》为例,不会开题答辩的可以进来看看
  • 广州优化网站关键词静态网站用什么做
  • 油棕种植密度控制:单位面积产值与光照利用优化
  • 紫金桥软件:深耕工业数据价值,赋能智造未来
  • Unity / C# 闭包详解 —— 按钮回调、协程、事件中的坑与修复
  • 使用Windbg分析dmp文件的方法以及实战分析实例分享
  • 什么是m3u8协议
  • 富文本编辑器Tinymce的使用、语言包配置、工具栏与工具栏组添加
  • 西安网站优化seo网站ico在后台哪里找到
  • 生成式人工智能教育应用的治理策略:构建“法律-伦理-标准”三位一体治理体系
  • 哪些行业对做网站的需求大qux wordpress
  • 智能中控终端:多系统联动管控中枢-LKONE立控信息
  • Mac M2安装VUE3
  • clang编译器 abseil-cpp中的ABSL_MUST_USE_RESULT
  • 网站营销怎么做做公众号用什么网站吗
  • ReactUse 是一个为 React 开发者设计的实用工具库(陀螺仪、地理位置、虚拟列表、标签页通讯函数)
  • ORB-SLAM2之ORB特征提取
  • 【Anaconda】anaconda卸载重安装