MathLabs

算术与数论

黎曼ζ函数与黎曼猜想

一个复函数把自然数和质数联系在一起;如果它的每一个非平凡零点都恰好落在同一条临界线上,我们就能精确知道质数分布得有多均匀——这是数学中最深刻的未解难题之一。

大学本页假定读者已知的两个经典事实

每个大于 11 的整数都能以唯一的方式(不计因数顺序)分解为质因数。

为什么成立?

正是这种唯一性,使得下面的欧拉乘积得以成立:展开并合并同类项后,每个 n−sn^{-s} 都恰好只出现一次,既不重复也不遗漏。

证明

存在性通过对 nn 的强归纳法证明。基础情形 n=2n=2 已经是质数,故命题成立。假设每个整数 2≤m<n2\le m<n 都能分解为质因数之积,考察 nn。若 nn 本身是质数则完成;否则 n=abn=ab 且 1<a,b<n1<a,b<n,由归纳假设可得 aa 与 bb 的质因数分解。把这两个分解拼接起来即得 nn 的质因数分解,归纳完成。

唯一性依赖于欧几里得引理:若质数 pp 满足 p∣abp\mid ab,则 (p∣a)∨(p∣b)(p\mid a)\lor(p\mid b)。为证明这一点,设 p∤ap\nmid a;此时 gcd⁡(p,a)=1\gcd(p,a)=1,由贝祖等式存在整数 x,yx,y 使得 px+ay=1px+ay=1。两边乘以 bb 得 b=pbx+abyb=pbx+aby。由于 p∣abp\mid ab,质数 pp 整除右边两项,因此 p∣bp\mid b。

现设 nn 有两个质因数分解 n=p1p2⋯pr=q1q2⋯qsn=p_1p_2\cdots p_r=q_1q_2\cdots q_s,每个列表都按不减顺序排列。由于 p1∣q1q2⋯qsp_1\mid q_1q_2\cdots q_s,反复使用欧几里得引理可知 p1p_1 整除某个 qjq_j,又因 qjq_j 是质数,故 p1=qjp_1=q_j。从两个乘积中消去这个公共因子,并对较短的列表重复此论证,经过有限步后可得 r=sr=s,且两个分解逐项相同——因此 nn 的质因数分解在不计因数顺序的意义下是唯一的。

不存在最大的质数;质数永远不会用完。

为什么成立?

欧几里得的经典证明(约公元前300年)从一份有限清单出发,用反证法论证。本页稍后,欧拉利用调和级数在 s=1s=1 处发散,给出了同一事实的第二个解析证明。

证明

反证:假设只有有限多个质数 p1,p2,…,pkp_1,p_2,\ldots,p_k。构造数 N=p1p2⋯pk+1N=p_1p_2\cdots p_k+1。

对每个 ii,用 pip_i 去除 NN 余数为 11,即 N≡1(modpi)N\equiv 1\pmod{p_i},故没有任何 pip_i 整除 NN。

但 N>1N>1,由上面证明的算术基本定理,NN 至少有一个质因数 qq。由于 p1,…,pkp_1,\ldots,p_k 都不整除 NN,这个 qq 满足 q∉{p1,…,pk}q\notin\{p_1,\ldots,p_k\}——出现了一个不在所谓完整清单中的质数,矛盾。因此任何有限列表都不能包含所有质数:质数有无穷多个。

研究看待同一信息的两种方式

1737年,欧拉发现,对实数 s>1s>1,取遍所有自然数的级数恰好等于只取遍质数的一个乘积:

ζ(s)=∑n=1∞1ns=∏p prime(1−1ps)−1,Re⁡(s)>1\zeta(s) = \sum_{n=1}^{\infty} \frac{1}{n^s} = \prod_{p \text{ prime}} \left(1 - \frac{1}{p^s}\right)^{-1}, \qquad \operatorname{Re}(s) > 1

左边只涉及整数,右边只涉及质数。因此,关于质数分布方式的一切信息,都藏在 ζ\zeta 之中的某个角落。

当 Re⁡(s)>1\operatorname{Re}(s) > 1 时,∑n=1∞n−s=∏p(1−p−s)−1\displaystyle\sum_{n=1}^{\infty} n^{-s} = \prod_{p} \left(1-p^{-s}\right)^{-1},乘积取遍所有质数 pp。

为什么成立?

将每个因子 (1−p−s)−1=1+p−s+p−2s+⋯(1-p^{-s})^{-1} = 1+p^{-s}+p^{-2s}+\cdots 展开成等比级数,再全部相乘:质因数分解的唯一性意味着每个 n−sn^{-s} 在展开后的乘积中恰好只出现一次。这个乘积之所以需要无穷多个因子,正是因为质数有无穷多个——这正是欧几里得用初等方法证明的同一事实,如今被包进了一个解析恒等式之中。

证明

固定 Re⁡(s)>1\operatorname{Re}(s)>1,将截断乘积 PN(s)=∏p≤N(1−p−s)−1P_N(s)=\prod_{p\le N}\left(1-p^{-s}\right)^{-1} 的每个因子按等比级数展开:由于 ∣p−s∣<1|p^{-s}|<1,有 (1−p−s)−1=1+p−s+p−2s+⋯(1-p^{-s})^{-1}=1+p^{-s}+p^{-2s}+\cdots。

把这些等比级数相乘并展开后,出现的每一项都是形如 n−sn^{-s} 的项,其中 nn 是所有质因数都不超过 NN 的整数——记这样的整数集合为 SNS_N。由算术基本定理,每个这样的 nn 恰好来自一种质数幂的组合,因此没有任何一项被重复计入或遗漏:PN(s)=∑n∈SNn−sP_N(s)=\sum_{n\in S_N} n^{-s}。

从 11 到 NN 的每个整数都属于 SNS_N,因此完整级数与截断乘积之差只涉及不属于 SNS_N 的整数,而这些整数都大于 NN:ζ(s)−PN(s)=∑n∉SNn−s\zeta(s)-P_N(s)=\sum_{n\notin S_N} n^{-s},故 ∣ζ(s)−PN(s)∣≤∑n>Nn−σ\left|\zeta(s)-P_N(s)\right|\le\sum_{n>N} n^{-\sigma},其中 σ=Re⁡(s)>1\sigma=\operatorname{Re}(s)>1。这个收敛级数的尾项在 N→∞N\to\infty 时趋于 00,因此 PN(s)→ζ(s)P_N(s)\to\zeta(s),无穷乘积等于无穷级数。

例题: 用解析方法重新证明质数无穷多

假设质数只有有限个。在 s=1s=1 处会出现什么矛盾?

解答

如果质数只有有限个 p1,…,pkp_1,\ldots,p_k,那么即使让 s→1s\to1,欧拉乘积 ∏i(1−pi−1)−1\prod_i (1-p_i^{-1})^{-1} 也会是一个有限数。但左边 ∑n≥1n−1\sum_{n\ge1} n^{-1} 是调和级数,是发散的。有限的乘积不可能等于无穷的和,所以质数的清单不可能是有限的——这正是欧拉本人利用调和级数发散、而非直接组合论证给出的欧几里得定理证明。

例题: 密码学应用:需要试多少个候选才能碰到质数?

RSA密钥生成需要找到接近 220482^{2048} 的大质数。利用素数定理 π(x)∼xlog⁡x\pi(x)\sim\dfrac{x}{\log x},平均需要测试多少个接近 220482^{2048} 的随机奇数候选,才能找到一个质数?

解答

素数定理表明,靠近大数 xx 处质数的密度大约为 1/log⁡x1/\log x:靠近 xx 随机选取的整数是质数的概率大约为 1/log⁡x1/\log x。这里 x=22048x=2^{2048},故需要 log⁡(22048)=2048log⁡2≈1420\log(2^{2048})=2048\log2\approx1420。

因此,靠近 220482^{2048} 的一个均匀随机整数是质数的概率约为 1/14201/1420,即平均需要测试 14201420 个随机候选才能找到一个质数。

实际中,候选数被限制为奇数(可免费排除掉那些平凡的偶数合数),使搜索范围减半;再用几千以内的质数做试除,可低成本地排除大部分剩余的合数,然后才运行开销较大的米勒–拉宾素性检验。这样可把预期需要的昂贵素性检验次数降到约 710710 次——这是素数定理在选择RSA密钥长度和估计密钥生成时间方面的一个直接、实用的应用。

大学把 ζ 延拓到级数收敛区域之外

黎曼在1859年的论文中证明,原本只在 Re⁡(s)>1\operatorname{Re}(s)>1 上有定义的 ζ\zeta,可以延拓成整个复平面上的单一全纯函数,只在 s=1s=1 处有一个单极点(在那里它的行为如同发散的调和级数)。一旦引入完备化的 zeta 函数 ξ\xi,这个延拓便满足一个惊人的对称性,即函数方程:

ξ(s)=12s(s−1)π−s/2Γ ⁣(s2)ζ(s),ξ(s)=ξ(1−s)\xi(s) = \tfrac{1}{2}s(s-1)\pi^{-s/2}\Gamma\!\left(\tfrac{s}{2}\right)\zeta(s), \qquad \xi(s) = \xi(1-s)

定义: 平凡零点与非平凡零点

函数方程迫使 ζ\zeta 在负偶数 s=−2,−4,−6,…s=-2,-4,-6,\ldots 处取零值;这些称为平凡零点,来自 Γ(s/2)\Gamma(s/2) 极点被抵消。其余的零点——称为非平凡零点——都必须满足 0<Re⁡(s)<10 < \operatorname{Re}(s) < 1,即临界带,并且总是成群出现:若 ρ\rho 是一个非平凡零点,则 ρˉ\bar\rho、1−ρ1-\rho、1−ρˉ1-\bar\rho 也都是零点,这源于函数方程以及 ζ\zeta 的系数为实数这一事实。偏离直线 Re⁡(s)=12\operatorname{Re}(s)=\tfrac12 的零点从不孤单出现。

研究黎曼猜想

就在那篇1859年的论文——他关于数论的唇一著作——中,黎曼提出了可能最强的对称性猜想:所有非平凡零点都恰好位于临界线 Re⁡(s)=12\operatorname{Re}(s)=\tfrac12 上。

黎曼猜想。 ζ(s)\zeta(s) 的每一个非平凡零点,实部都恰好等于 12\tfrac12。

图书馆把这个问题收录为大问题黎曼猜想:它是克雷数学研究所(2000年)七个千禧年大奖难题之一,也是希尔伯特Ā1900年问题清单中的第8题。尽管近以来找到的每一个零点——数量多达数万亿个——都恰好落在这条直线上,它至仍未解决。

映射 $z\mapsto e^z$ 在复数域上的网格图,用作解析延拓的视觉替代:黎曼 ζ 函数本身并不在该组件的预设函数之中,因此此图并未展示 ζ 的任何真实零点。
仅作示意:全纯映射 z↦ezz\mapsto e^z 下的网格变形,用来代替 ζ 的解析延拓。

研究零点如何支配质数

定义 ψ(x)=∑pk≤xlog⁡p\psi(x) = \sum_{p^k \le x} \log p,即对不超过 xx 的质数幂做加权计数。冯·曼戈尔特的显式公式(1895年)用 ζ\zeta 的零点 ρ\rho 精确地表出了 ψ\psi:

ψ0(x)=x−∑ρxρρ−log⁡2π−12log⁡ ⁣(1−x−2)\psi_0(x) = x - \sum_{\rho} \frac{x^{\rho}}{\rho} - \log 2\pi - \tfrac12\log\!\left(1-x^{-2}\right)

每一对零点 ρ,ρˉ\rho,\bar\rho 都贡献一个大小约为 xβ/∣ρ∣x^{\beta}/|\rho| 的振荡项,其中 β=Re⁡(ρ)\beta=\operatorname{Re}(\rho)。β\beta 越大,这个波动就越大——因此零点的实部直接决定了 ψ(x)\psi(x),以及因此质数计数函数 π(x)\pi(x),能够偏离其光滑预测值多远:

RH  ⟺  ψ(x)=x+O ⁣(x log⁡2x)  ⟺  π(x)=Li⁡(x)+O ⁣(x log⁡x)\text{RH} \iff \psi(x) = x + O\!\left(\sqrt{x}\,\log^2 x\right) \iff \pi(x) = \operatorname{Li}(x) + O\!\left(\sqrt{x}\,\log x\right)
定理: 素数定理

当 x→∞x\to\infty 时,π(x)∼xlog⁡x\pi(x) \sim \dfrac{x}{\log x}:不超过 xx 的质数个数渐近于 x/log⁡xx/\log x。

为什么成立?

阿达马和德拉瓦莱·普桑于1896年各自独立地证明了这一点,两个证明都要经过 ζ\zeta:该定理等价于 ζ(s)\zeta(s) 在直线 Re⁡(s)=1\operatorname{Re}(s)=1 上没有零点。这个已被证明的事实,远弱于黎曼猜想(后者要求直线 Re⁡(s)=12\operatorname{Re}(s)=\tfrac12 之外任何地方都没有零点),这正是该定理已解决而RH尚未解决的原因。RH只会把近似的误差项,从大约 x/log⁡xx/\log x 的精度收紧到上面接近最优的 xlog⁡x\sqrt{x}\log x。

证明

素数定理等价于关于切比雪夫函数 ψ(x)=∑pk≤xlog⁡p\psi(x)=\sum_{p^k\le x}\log p 的一个命题,该函数把不超过 xx 的每个质数幂按其质数底的对数加权:证明 ψ(x)∼x\psi(x)\sim x(通过一个常规的分部求和论证)等价于证明 π(x)∼x/log⁡x\pi(x)\sim x/\log x。

黎曼的显式公式几乎精确地把 ψ(x)\psi(x) 用 ζ\zeta 的非平凡零点 ρ=β+iγ\rho=\beta+i\gamma(都满足 0≤β≤10\le\beta\le1)表示出来:ψ(x)=x−∑ρxρρ−log⁡(2π)−12log⁡(1−x−2)\psi(x)=x-\sum_{\rho}\frac{x^{\rho}}{\rho}-\log(2\pi)-\tfrac12\log(1-x^{-2})。主项 xx 来自 ζ\zeta 在 s=1s=1 处的单极点;其余每一项都以大约 xβx^{\beta} 的量级振荡,因此定理归结为证明这个振荡和为 o(x)o(x)。

这个振荡和为 o(x)o(x) 当且仅当没有零点满足 β=1\beta=1,即 ζ(1+it)≠0\zeta(1+it)\ne0 对所有 t∈Rt\in\mathbb{R} 成立。阿达马和德拉瓦莱·普桑于1896年各自独立地证明了这一非消失性,方法是将初等三角不等式 3+4cos⁡θ+cos⁡2θ≥03+4\cos\theta+\cos2\theta\ge0 应用于 σ→1+\sigma\to1^+ 时 3log⁡ζ(σ)+4log⁡ζ(σ+it)+log⁡ζ(σ+2it)3\log\zeta(\sigma)+4\log\zeta(\sigma+it)+\log\zeta(\sigma+2it) 的实部:若 1+it1+it 处有零点,这个组合就会被迫发散到 −∞-\infty,而不等式不允许这种情况。

把这两步结合起来即得 ψ(x)=x+o(x)\psi(x)=x+o(x),再由联系 ψ\psi 与 π\pi 的分部求和恒等式,把它提升为 π(x)∼x/log⁡x\pi(x)\sim x/\log x,这正是素数定理。

三维涟漪曲面 $z=\sin r/r$ 被用作视觉类比,表示每一对 ζ 零点为 ψ(x) 显式公式所贡献的振荡波;这并非该求和的真实图像,因为 ζ 的零点不在该组件的预设函数之中。
仅作示意:涟漪曲面 z=sin⁡r/rz=\sin r/r,用来代替每一对零点为显式公式添加的振荡波。

研究工具箱:无条件方法能走多远

没有任何一种工具能接近证明RH,但每一种都能把临界带的一部分圏出来。无零点区域是指靠近 Re⁡(s)=1\operatorname{Re}(s)=1、已知不含任何零点的一条带状区域:德拉瓦莱·普桑(1899年)找到了区域 Re⁡(s)≥1−c/log⁡t\operatorname{Re}(s) \ge 1 - c/\log t,而维诺格拉多夫与科罗博夫(1958年)利用指数和估计,将其推进到 Re⁡(s)≥1−c/(log⁡t)2/3(log⁡log⁡t)1/3\operatorname{Re}(s) \ge 1 - c/(\log t)^{2/3}(\log\log t)^{1/3}——这个指数六十多年来基本未被改进。它的宽度在 t→∞t\to\infty 时仍会缩到零;即便是一条固定宽度、不含零点的带 Re⁡(s)>1−ε\operatorname{Re}(s) > 1-\varepsilon(“弱黎曼猜想”)也仍未解决。

如果无法排除零点,那就去数它们。设 N(σ,T)N(\sigma,T) 为实部至少为 σ\sigma、高度不超过 TT 的零点个数。英厄姆(1940年)证明了 N(σ,T)≪T3(1−σ)/(2−σ)N(\sigma,T) \ll T^{3(1-\sigma)/(2-\sigma)};赫克斯利(1972年)把 σ=3/4\sigma=3/4 处的指数改进到 12/512/5,并在此停留了近五十年——直到古斯与梅纳德 2024 年关于狄利德雷多项式大值的估计(2026年发表于《数学年刊》)把它降到 30/1330/13。一个直接推论是:素数定理如今在每个区间 [x,x+x17/30+ε][x, x+x^{17/30+\varepsilon}] 上都无条件地成立。

每种无条件工具的停顿之处
工具最佳结果停顿之处
无零点区域维诺格拉多夫–科罗博夫 1958年,log⁡t\log t 的指数为2/3宽度仍趨于0;尚无固定带 Re⁡(s)>1−ε\operatorname{Re}(s)>1-\varepsilon
零点密度古斯–梅纳德 2024年,A(3/4)=30/13A(3/4)=30/13只计数不排除;即便密度猜想 A=2A=2 也达不到RH
计算验证满足 γ≤3×1012\gamma \le 3\times10^{12} 的所有零点均已验证(普拉特–特鲁贾恩,2021年)只是有限范围的检验;对 t→∞t\to\infty 无法说明任何事
等价判据德布鲁因–纽曼:0≤Λ≤0.220 \le \Lambda \le 0.22(罗杰斯–陶 2018年;Polymath 15)只是把RH(Λ=0\Lambda=0)换了种说法,难度并未降低

德布鲁因–纽曼常数 Λ\Lambda 给出了RH的一种精确表述:随着 tt 增大,ξ\xi 的热流形变 HtH_t 一旦 t≥Λt\ge\Lambda 就只有实零点,而RH等价于 Λ≤0\Lambda\le0。罗杰斯和陶(2018年)证明了 Λ≥0\Lambda\ge0,因此如果RH成立,它也只是恰好成立——Λ=0\Lambda=0 精确成立。此后,Polymath 15 项目又把上界压低到了 Λ≤0.22\Lambda\le0.22。

RH的某个版本被完全证明的唇一场合是在函数域上:对于有限域 Fq\mathbb{F}_q 上的一条曲线,其 zeta 函数是 q−sq^{-s} 的有理函数,而RH的类比断言其所有零点的绝对值都等于 q\sqrt{q}。哈塞(1933年)对椅圆曲线证明了这一点,韦伊(1948年)对所有曲线证明了这一点,德利淅(1974年)则利用作用在上同调群上的弗罗贝尼乌斯算子,对有限域上的所有代数簇证明了这一点。把这套方法照搬到 Z\mathbb{Z} 上的障碍在于:整数上不存在弗罗贝尼乌斯算子,也没有已知的几何对象能扮演 Spec⁡Z×Spec⁡Z\operatorname{Spec}\mathbb{Z}\times\operatorname{Spec}\mathbb{Z} 的角色。

广义黎曼猜想(GRH)把同一个猜想推广到每一个狄利德雷 LL 函数,并在图书馆再往前一步的地方作为实用工具出现:对于大问题哥德巴赫猜想,GRH曾帮助证明了足够大的数满足奇数(三素数)情形,而上面建立的无条件工具——无零点区域、零点密度估计——则在完全不靠任何猜想的情况下,限定了偶数情形能走多远。接下来要讲的筛法,恰好补上了这个故事的另一半。

为什么欧拉乘积 ζ(s)=∏p(1−p−s)−1\zeta(s)=\prod_p(1-p^{-s})^{-1} 需要无穷多个因子?

ζ(s)\zeta(s) 的每一个非平凡零点都满足哪个条件?

下列哪一项已经得到证明(而非猜想)?

数万亿个零点已经被计算验证位于临界线上。为什么这不能证明黎曼猜想?

参考文献

  1. H. M. Edwards (1974). Riemann's Zeta Function · DOI:10.1090/s0273-0979-01-00912-0
  2. E. C. Titchmarsh; revised by D. R. Heath-Brown (1986). The Theory of the Riemann Zeta-Function · DOI:10.1070/im1975v009n03abeh001485
  3. J. B. Conrey (2003). The Riemann Hypothesis
  4. Larry Guth, James Maynard (2026). New large value estimates for Dirichlet polynomials · arXiv:2405.20552
  5. Dave Platt, Tim Trudgian (2021). The Riemann hypothesis is true up to 3×10^12 · arXiv:2004.09765
  6. Brad Rodgers, Terence Tao (2020). The de Bruijn–Newman constant is non-negative · arXiv:1801.05914
  7. D. H. J. Polymath (2019). Effective approximation of heat flow evolution of the Riemann ξ function, and a new upper bound for the de Bruijn–Newman constant · arXiv:1904.12438