MathLabs

算术与数论

圆法

利用单位圆上的傅里叶积分来计数加性方程解数的一种解析技巧。

直观直觉:通过绕圆一周来计数

假设你想计数将整数 nn 写成取自整数集合 AA(例如素数,或恰好 kk 次方数)的 kk 个元素之和的方法数。将 AA 编码为绕单位圆恰好一周的"波":对 [0,1)[0,1) 中每个实数 α\alpha,构造指数和 FA(α)=∑a∈Ae(aα)F_A(\alpha) = \sum_{a \in A} e(a\alpha),其中 e(x):=e2πixe(x) := e^{2\pi i x}。当 α\alpha 从 00 扫到 11 时,FA(α)kF_A(\alpha)^k 会振荡;在少数"共振"点附近——分母 qq 较小的有理数 a/qa/q——和式各项方向一致而相长叠加(称为主弧),而在几乎所有其余地方,相位方向近乎随机而相互抵消(称为次弧)。圆法把这幅几何图景变成一个精确公式,再逐弧估计。

交互式单位圆,显示角度为 theta 的点,用于构造圆法中的指数和。
拖动 θ\theta 会使点 e2πiθ/360e^{2\pi i\theta/360} 沿单位圆移动;当 θ/360=α\theta/360=\alpha 取遍 [0,1)[0,1) 时,圆法就是对 FA(α)kF_A(\alpha)^k 在这些点上做积分。

大学定义:圆上的生成函数

定义: 指数和与表示数

对有限集合 A⊂Z≥0A\subset\mathbb Z_{\ge0},记 FA(α)=∑a∈Ae(aα)F_A(\alpha) = \sum_{a \in A} e(a\alpha),其中 e(x):=e2πixe(x) := e^{2\pi i x}。对 n≥0n\ge0、k≥1k\ge1,设 rk(n)=#{(a1,…,ak)∈Ak:a1+⋯+ak=n}r_k(n)=\#\{(a_1,\dots,a_k)\in A^k : a_1+\cdots+a_k=n\} 为将 nn 表示为 AA 中 kk 个元素之和的有序表示个数。

FA(α)=∑a∈Ae(aα),e(x):=e2πixF_A(\alpha) = \sum_{a \in A} e(a\alpha), \qquad e(x) := e^{2\pi i x}

将 FAF_A 取 kk 次幂并展开,e(nα)e(n\alpha) 的系数"应当"恰好是 rk(n)r_k(n);下面的提取恒等式(定理1)使这一点严格化:rk(n)=∫01FA(α)k e(−nα) dαr_k(n) = \int_0^1 F_A(\alpha)^k\, e(-n\alpha)\, d\alpha。这里 FA(α)kF_A(\alpha)^k 正是上图中的那个波,积分恰好挑出我们关心的那一个频率 nn——这就是该方法的全部内容:把一个计数问题换成一个积分估计问题。

rk(n)=∫01FA(α)k e(−nα) dαr_k(n) = \int_0^1 F_A(\alpha)^k\, e(-n\alpha)\, d\alpha
主弧与次弧的对比
特征主弧(靠近 qq 较小的 a/qa/q)次弧(其余部分)
位置每个满足 q≤Qq\le Q 的有理数 a/qa/q 周围的一个短区间去掉所有主弧后 [0,1)[0,1) 剩下的部分
FA(α)F_A(\alpha) 的大小接近其平凡最大值 ∣A∣|A|:各项相长预期远小于 ∣A∣|A|:各项相消
在估计中的作用给出主项(一个"奇异级数")必须给出上界并归入误差项

大学两个基础定理

对任意有限集合 A⊂Z≥0A\subset\mathbb Z_{\ge0}、整数 k≥1k\ge1 及 n≥0n\ge0,有 rk(n)=∫01FA(α)k e(−nα) dαr_k(n) = \int_0^1 F_A(\alpha)^k\, e(-n\alpha)\, d\alpha,其中 rk(n)r_k(n) 是取自 AA 且和为 nn 的有序 kk 元组个数。

为什么成立?

这使我们能把纯组合计数问题换成关于某个解析积分大小的问题:只要证明该积分为正,就必定存在一种表示。

证明

首先利用正交关系:∫01e(mα) dα\int_0^1 e(m\alpha)\, d\alpha=1=1(当 m=0m=0 时),=0=0(当 m≠0m\neq 0 时):写 e(mα)=cos⁡(2πmα)+isin⁡(2πmα)e(m\alpha)=\cos(2\pi m\alpha)+i\sin(2\pi m\alpha),当 m≠0m\neq0 时这恰好是正弦/余弦波在 [0,1)[0,1) 上整数个周期,积分为 00;当 m=0m=0 时被积函数是常数 11。

现在展开生成函数的 kk 次幂:FA(α)k=(∑a∈Ae(aα))k=∑(a1,…,ak)∈Ake((a1+⋯+ak)α)F_A(\alpha)^k=\Big(\sum_{a\in A}e(a\alpha)\Big)^k=\sum_{(a_1,\dots,a_k)\in A^k} e\big((a_1+\cdots+a_k)\alpha\big),这是对取自 AA 的所有有序 kk 元组按指数 m=a1+⋯+akm=a_1+\cdots+a_k 分组后的一个有限和。

两边乘以 e(−nα)e(-n\alpha),再在 [0,1)[0,1) 上逐项积分(因为是有限和,积分与求和可交换,故合法):∫01FA(α)ke(−nα) dα=∑(a1,…,ak)∈Ak∫01e((a1+⋯+ak−n)α) dα\int_0^1 F_A(\alpha)^k e(-n\alpha)\,d\alpha=\sum_{(a_1,\dots,a_k)\in A^k}\int_0^1 e\big((a_1+\cdots+a_k-n)\alpha\big)\,d\alpha。

由正交关系,右边每一项恰好当 a1+⋯+ak=na_1+\cdots+a_k=n 时为 11,否则为 00。于是整个和恰好收缩为满足 a1+⋯+ak=na_1+\cdots+a_k=n 的元组个数,即 rk(n)r_k(n)。这就精确地(无任何近似)证明了该恒等式:它是一个恒等式,而不是渐近公式。

对任意实数 α\alpha 与任意整数 Q≥1Q\ge1,存在有理数 a/qa/q,满足 1≤q≤Q1\le q\le Q、gcd⁡(a,q)=1\gcd(a,q)=1,使得 ∣α−aq∣≤1q(Q+1)\left|\alpha - \dfrac{a}{q}\right| \le \dfrac{1}{q(Q+1)}。

为什么成立?

这保证了圆上每一点都靠近某个分母较小的有理数,这正是我们能把 [0,1)[0,1) 划分为主弧(围绕分母 qq 较小的良好逼近 a/qa/q 的短区间,在此 FAF_A 较大)与次弧(其余部分)的原因。

证明

考虑 Q+2Q+2 个数 0,{α},{2α},…,{Qα},10,\{\alpha\},\{2\alpha\},\dots,\{Q\alpha\},1,其中 {x}\{x\} 表示 xx 的小数部分;它们都属于 [0,1][0,1]。将 [0,1][0,1] 划分为 Q+1Q+1 个长度为 1/(Q+1)1/(Q+1) 的相等子区间:[0,1Q+1),[1Q+1,2Q+1),…[0,\tfrac1{Q+1}), [\tfrac1{Q+1},\tfrac2{Q+1}),\dots。

我们有 Q+2Q+2 个数却只有 Q+1Q+1 个子区间,由鸽笼原理,必有两个数 {jα}\{j\alpha\} 与 {iα}\{i\alpha\}(0≤i<j≤Q0\le i<j\le Q,允许 i=0i=0 即 {iα}=0\{i\alpha\}=0)落入同一子区间,从而相差小于 1/(Q+1)1/(Q+1):∣{jα}−{iα}∣<1Q+1|\{j\alpha\}-\{i\alpha\}|<\tfrac{1}{Q+1}。

设 q=j−iq=j-i,则 1≤q≤Q1\le q\le Q。因为 {jα}−{iα}=(jα−iα)−(⌊jα⌋−⌊iα⌋)=qα−a\{j\alpha\}-\{i\alpha\} = (j\alpha - i\alpha) - (\lfloor j\alpha\rfloor - \lfloor i\alpha\rfloor) = q\alpha - a(其中整数 a=⌊jα⌋−⌊iα⌋a=\lfloor j\alpha\rfloor-\lfloor i\alpha\rfloor),故 ∣qα−a∣<1Q+1|q\alpha - a| < \tfrac{1}{Q+1},即 ∣α−aq∣<1q(Q+1)\left|\alpha-\dfrac aq\right| < \dfrac{1}{q(Q+1)}。

最后,若 gcd⁡(a,q)=d>1\gcd(a,q)=d>1,将 aa 和 qq 同除以 dd:所得分数分母更小,与 α\alpha 的距离相同或更小,因此不失一般性可设 gcd⁡(a,q)=1\gcd(a,q)=1。证明完成;这是一个有限的、构造性的鸽笼论证,不依赖任何未经证明的估计。

大学实际应用与典型例题

超出数论范围,"把波相加并寻找共振"的同一想法正是信号处理与电子工程中广泛使用的离散傅里叶变换的工作原理;而哈代与拉马努金用圆法早期版本(应用于分拆函数 p(n)p(n))得到的渐近公式 p(n)∼14n3exp⁡ ⁣(π2n3)p(n) \sim \dfrac{1}{4n\sqrt3}\exp\!\left(\pi\sqrt{\dfrac{2n}{3}}\right),被统计力学用来估计总能量固定为 nn 的不可区分玻色激发系统的微观态数——也就是熵。

例题: 圆法的一个小规模计算

设 A={1,2,…,8}A=\{1,2,\dots,8\}。利用提取恒等式,计算满足 a+b=9a+b=9 的有序对 (a,b)∈A2(a,b)\in A^2 的个数 r2(9)r_2(9)。

解答

由定理1,r2(9)=∫01FA(α)2e(−9α) dαr_2(9)=\int_0^1 F_A(\alpha)^2 e(-9\alpha)\,d\alpha,其中 FA(α)=∑a=18e(aα)F_A(\alpha)=\sum_{a=1}^{8}e(a\alpha);无需真的解析计算该积分——恒等式已保证它等于直接计数的结果,因此我们可以用组合方法计数,并信任该恒等式。

列出所有满足 a,b∈{1,…,8}a,b\in\{1,\dots,8\} 且 a+b=9a+b=9 的有序对 (a,b)(a,b):(1,8),(2,7),(3,6),(4,5),(5,4),(6,3),(7,2),(8,1)(1,8),(2,7),(3,6),(4,5),(5,4),(6,3),(7,2),(8,1)。从 11 到 88 的每个值唯一确定 b=9−ab=9-a,且 bb 总落回 {1,…,8}\{1,\dots,8\}(因为 1≤a≤8⇒1≤9−a≤81\le a\le8 \Rightarrow 1\le 9-a\le8),所以 aa 的全部 88 个取值都有效。

因此 r2(9)r_2(9) 等于 88。这个小例子正是定理1恒等式的具体体现:解析积分与组合计数按构造是同一个数——该方法的真正内容只有当 AA 变成一个无限或增长的集合(如素数集)、必须估计而非枚举时才显现出来。

例题: 利用狄利克雷定理求最佳有理逼近

取 α=2\alpha=\sqrt2、Q=5Q=5,应用狄利克雷逼近定理,给出满足定理2界的分数 a/qa/q(1≤q≤51\le q\le5),并做数值验证。

解答

我们要找 1≤q≤51\le q\le5 且 ∣2−a/q∣<1/(q⋅6)|\sqrt2-a/q|<1/(q\cdot6)(取定理界中的 Q+1=6Q+1=6)的 a/qa/q。试 q=5q=5:最接近 52≈7.07115\sqrt2\approx7.0711 的整数是 a=7a=7,给出 7/5=1.47/5=1.4。

验证界:∣2−7/5∣=∣1.41421…−1.4∣≈0.01421|\sqrt2-7/5|=|1.41421\ldots-1.4|\approx0.01421,而定理保证 1/(5⋅6)=1/30≈0.03331/(5\cdot6)=1/30\approx0.0333;确实 0.01421<0.03330.01421<0.0333,界成立且留有余量,正如所保证的那样。

这个 7/57/5 实际上是 2=[1;2,2,2,… ]\sqrt2=[1;2,2,2,\dots] 连分数两步后的渐近分数,这正是它逼近得如此好的原因——狄利克雷的鸽笼证明本身不需要连分数即可成立,但总能得到同等好的逼近。因此 75\dfrac{7}{5} 是一个有效的见证。

在提取恒等式 rk(n)=∫01FA(α)k e(−nα) dαr_k(n) = \int_0^1 F_A(\alpha)^k\, e(-n\alpha)\, d\alpha 中,为什么这个积分恰好等于 rk(n)r_k(n)?

设 A={1,2,…,8}A=\{1,2,\dots,8\},满足 a+b=9a+b=9 的有序对 (a,b)∈A2(a,b)\in A^2 的个数 r2(9)r_2(9) 是多少?

根据狄利克雷逼近定理,对任意实数 α\alpha 与整数 Q≥1Q\ge1,保证存在什么?

在统计力学中,分拆函数的哈代–拉马努金渐近公式 p(n)∼14n3exp⁡ ⁣(π2n3)p(n) \sim \dfrac{1}{4n\sqrt3}\exp\!\left(\pi\sqrt{\dfrac{2n}{3}}\right) 有助于估计什么?

参考文献

  1. G. H. Hardy, S. Ramanujan (1918). Asymptotic formulae in combinatory analysis · DOI:10.1112/plms/s2-17.1.75
  2. J. Bourgain, C. Demeter, L. Guth (2016). Proof of the main conjecture in Vinogradov's Mean Value Theorem for degrees higher than three · DOI:10.4007/annals.2016.184.2.7 · arXiv:1512.01565
  3. B. Green, T. Tao (2008). The primes contain arbitrarily long arithmetic progressions · DOI:10.4007/annals.2008.167.481 · arXiv:math/0404188