← 返回 资料库 › 算术与数论 › 解析数论 算术与数论
圆法 利用单位圆上的傅里叶积分来计数加性方程解数的一种解析技巧。
直观 直觉:通过绕圆一周来计数 假设你想计数将整数 n n n 写成取自整数集合 A A A (例如素数,或恰好 k k k 次方数)的 k k k 个元素之和的方法数。将 A A A 编码为绕单位圆恰好一周的"波":对 [ 0 , 1 ) [0,1) [ 0 , 1 ) 中每个实数 α \alpha α ,构造指数和 F A ( α ) = ∑ a ∈ A e ( a α ) F_A(\alpha) = \sum_{a \in A} e(a\alpha) F A ( α ) = ∑ a ∈ A e ( a α ) ,其中 e ( x ) : = e 2 π i x e(x) := e^{2\pi i x} e ( x ) := e 2 π i x 。当 α \alpha α 从 0 0 0 扫到 1 1 1 时,F A ( α ) k F_A(\alpha)^k F A ( α ) k 会振荡;在少数"共振"点附近——分母 q q q 较小的有理数 a / q a/q a / q ——和式各项方向一致而相长叠加(称为主弧),而在几乎所有其余地方,相位方向近乎随机而相互抵消(称为次弧)。圆法把这幅几何图景变成一个精确公式,再逐弧估计。
拖动 θ \theta θ 会使点 e 2 π i θ / 360 e^{2\pi i\theta/360} e 2 π i θ /360 沿单位圆移动;当 θ / 360 = α \theta/360=\alpha θ /360 = α 取遍 [ 0 , 1 ) [0,1) [ 0 , 1 ) 时,圆法就是对 F A ( α ) k F_A(\alpha)^k F A ( α ) k 在这些点上做积分。 大学 定义:圆上的生成函数 定义: 指数和与表示数
对有限集合 A ⊂ Z ≥ 0 A\subset\mathbb Z_{\ge0} A ⊂ Z ≥ 0 ,记 F A ( α ) = ∑ a ∈ A e ( a α ) F_A(\alpha) = \sum_{a \in A} e(a\alpha) F A ( α ) = ∑ a ∈ A e ( a α ) ,其中 e ( x ) : = e 2 π i x e(x) := e^{2\pi i x} e ( x ) := e 2 π i x 。对 n ≥ 0 n\ge0 n ≥ 0 、k ≥ 1 k\ge1 k ≥ 1 ,设 r k ( n ) = # { ( a 1 , … , a k ) ∈ A k : a 1 + ⋯ + a k = n } r_k(n)=\#\{(a_1,\dots,a_k)\in A^k : a_1+\cdots+a_k=n\} r k ( n ) = # {( a 1 , … , a k ) ∈ A k : a 1 + ⋯ + a k = n } 为将 n n n 表示为 A A A 中 k k k 个元素之和的有序 表示个数。
F A ( α ) = ∑ a ∈ A e ( a α ) , e ( x ) : = e 2 π i x F_A(\alpha) = \sum_{a \in A} e(a\alpha), \qquad e(x) := e^{2\pi i x} F A ( α ) = a ∈ A ∑ e ( a α ) , e ( x ) := e 2 π i x 将 F A F_A F A 取 k k k 次幂并展开,e ( n α ) e(n\alpha) e ( n α ) 的系数"应当"恰好是 r k ( n ) r_k(n) r k ( n ) ;下面的提取恒等式(定理1)使这一点严格化:r k ( n ) = ∫ 0 1 F A ( α ) k e ( − n α ) d α r_k(n) = \int_0^1 F_A(\alpha)^k\, e(-n\alpha)\, d\alpha r k ( n ) = ∫ 0 1 F A ( α ) k e ( − n α ) d α 。这里 F A ( α ) k F_A(\alpha)^k F A ( α ) k 正是上图中的那个波,积分恰好挑出我们关心的那一个频率 n n n ——这就是该方法的全部内容:把一个计数 问题换成一个积分估计 问题。
r k ( n ) = ∫ 0 1 F A ( α ) k e ( − n α ) d α r_k(n) = \int_0^1 F_A(\alpha)^k\, e(-n\alpha)\, d\alpha r k ( n ) = ∫ 0 1 F A ( α ) k e ( − n α ) d α 主弧与次弧的对比 特征 主弧(靠近 q q q 较小的 a / q a/q a / q ) 次弧(其余部分) 位置 每个满足 q ≤ Q q\le Q q ≤ Q 的有理数 a / q a/q a / q 周围的一个短区间 去掉所有主弧后 [ 0 , 1 ) [0,1) [ 0 , 1 ) 剩下的部分 F A ( α ) F_A(\alpha) F A ( α ) 的大小接近其平凡最大值 ∣ A ∣ |A| ∣ A ∣ :各项相长 预期远小于 ∣ A ∣ |A| ∣ A ∣ :各项相消 在估计中的作用 给出主项(一个"奇异级数") 必须给出上界并归入误差项
大学 两个基础定理 对任意有限集合 A ⊂ Z ≥ 0 A\subset\mathbb Z_{\ge0} A ⊂ Z ≥ 0 、整数 k ≥ 1 k\ge1 k ≥ 1 及 n ≥ 0 n\ge0 n ≥ 0 ,有 r k ( n ) = ∫ 0 1 F A ( α ) k e ( − n α ) d α r_k(n) = \int_0^1 F_A(\alpha)^k\, e(-n\alpha)\, d\alpha r k ( n ) = ∫ 0 1 F A ( α ) k e ( − n α ) d α ,其中 r k ( n ) r_k(n) r k ( n ) 是取自 A A A 且和为 n n n 的有序 k k k 元组个数。
为什么成立? 这使我们能把纯组合计数问题换成关于某个解析积分大小的问题:只要证明该积分为正,就必定存在一种表示。
证明 首先利用正交关系:∫ 0 1 e ( m α ) d α \int_0^1 e(m\alpha)\, d\alpha ∫ 0 1 e ( m α ) d α = 1 =1 = 1 (当 m = 0 m=0 m = 0 时),= 0 =0 = 0 (当 m ≠ 0 m\neq 0 m = 0 时):写 e ( m α ) = cos ( 2 π m α ) + i sin ( 2 π m α ) e(m\alpha)=\cos(2\pi m\alpha)+i\sin(2\pi m\alpha) e ( m α ) = cos ( 2 π m α ) + i sin ( 2 π m α ) ,当 m ≠ 0 m\neq0 m = 0 时这恰好是正弦/余弦波在 [ 0 , 1 ) [0,1) [ 0 , 1 ) 上整数个周期,积分为 0 0 0 ;当 m = 0 m=0 m = 0 时被积函数是常数 1 1 1 。
现在展开生成函数的 k k k 次幂:F A ( α ) k = ( ∑ a ∈ A e ( a α ) ) k = ∑ ( a 1 , … , a k ) ∈ A k e ( ( a 1 + ⋯ + a k ) α ) 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) F A ( α ) k = ( ∑ a ∈ A e ( a α ) ) k = ∑ ( a 1 , … , a k ) ∈ A k e ( ( a 1 + ⋯ + a k ) α ) ,这是对取自 A A A 的所有有序 k k k 元组按指数 m = a 1 + ⋯ + a k m=a_1+\cdots+a_k m = a 1 + ⋯ + a k 分组后的一个有限和。
两边乘以 e ( − n α ) e(-n\alpha) e ( − n α ) ,再在 [ 0 , 1 ) [0,1) [ 0 , 1 ) 上逐项积分(因为是有限和,积分与求和可交换,故合法):∫ 0 1 F A ( α ) k e ( − n α ) d α = ∑ ( a 1 , … , a k ) ∈ A k ∫ 0 1 e ( ( a 1 + ⋯ + a k − 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 ∫ 0 1 F A ( α ) k e ( − n α ) d α = ∑ ( a 1 , … , a k ) ∈ A k ∫ 0 1 e ( ( a 1 + ⋯ + a k − n ) α ) d α 。
由正交关系,右边每一项恰好当 a 1 + ⋯ + a k = n a_1+\cdots+a_k=n a 1 + ⋯ + a k = n 时为 1 1 1 ,否则为 0 0 0 。于是整个和恰好收缩为满足 a 1 + ⋯ + a k = n a_1+\cdots+a_k=n a 1 + ⋯ + a k = n 的元组个数,即 r k ( n ) r_k(n) r k ( n ) 。这就精确地(无任何近似)证明了该恒等式:它是一个恒等式,而不是渐近公式。
对任意实数 α \alpha α 与任意整数 Q ≥ 1 Q\ge1 Q ≥ 1 ,存在有理数 a / q a/q a / q ,满足 1 ≤ q ≤ Q 1\le q\le Q 1 ≤ q ≤ Q 、gcd ( a , q ) = 1 \gcd(a,q)=1 g cd( a , q ) = 1 ,使得 ∣ α − a q ∣ ≤ 1 q ( Q + 1 ) \left|\alpha - \dfrac{a}{q}\right| \le \dfrac{1}{q(Q+1)} α − q a ≤ q ( Q + 1 ) 1 。
为什么成立? 这保证了圆上每一点都靠近某个分母较小的有理数,这正是我们能把 [ 0 , 1 ) [0,1) [ 0 , 1 ) 划分为主弧(围绕分母 q q q 较小的良好逼近 a / q a/q a / q 的短区间,在此 F A F_A F A 较大)与次弧(其余部分)的原因。
证明 考虑 Q + 2 Q+2 Q + 2 个数 0 , { α } , { 2 α } , … , { Q α } , 1 0,\{\alpha\},\{2\alpha\},\dots,\{Q\alpha\},1 0 , { α } , { 2 α } , … , { Q α } , 1 ,其中 { x } \{x\} { x } 表示 x x x 的小数部分;它们都属于 [ 0 , 1 ] [0,1] [ 0 , 1 ] 。将 [ 0 , 1 ] [0,1] [ 0 , 1 ] 划分为 Q + 1 Q+1 Q + 1 个长度为 1 / ( Q + 1 ) 1/(Q+1) 1/ ( Q + 1 ) 的相等子区间:[ 0 , 1 Q + 1 ) , [ 1 Q + 1 , 2 Q + 1 ) , … [0,\tfrac1{Q+1}), [\tfrac1{Q+1},\tfrac2{Q+1}),\dots [ 0 , Q + 1 1 ) , [ Q + 1 1 , Q + 1 2 ) , … 。
我们有 Q + 2 Q+2 Q + 2 个数却只有 Q + 1 Q+1 Q + 1 个子区间,由鸽笼原理,必有两个数 { j α } \{j\alpha\} { j α } 与 { i α } \{i\alpha\} { i α } (0 ≤ i < j ≤ Q 0\le i<j\le Q 0 ≤ i < j ≤ Q ,允许 i = 0 i=0 i = 0 即 { i α } = 0 \{i\alpha\}=0 { i α } = 0 )落入同一子区间,从而相差小于 1 / ( Q + 1 ) 1/(Q+1) 1/ ( Q + 1 ) :∣ { j α } − { i α } ∣ < 1 Q + 1 |\{j\alpha\}-\{i\alpha\}|<\tfrac{1}{Q+1} ∣ { j α } − { i α } ∣ < Q + 1 1 。
设 q = j − i q=j-i q = j − i ,则 1 ≤ q ≤ Q 1\le q\le Q 1 ≤ q ≤ 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 { j α } − { i α } = ( j α − i α ) − (⌊ j α ⌋ − ⌊ i α ⌋) = q α − a (其中整数 a = ⌊ j α ⌋ − ⌊ i α ⌋ a=\lfloor j\alpha\rfloor-\lfloor i\alpha\rfloor a = ⌊ j α ⌋ − ⌊ i α ⌋ ),故 ∣ q α − a ∣ < 1 Q + 1 |q\alpha - a| < \tfrac{1}{Q+1} ∣ q α − a ∣ < Q + 1 1 ,即 ∣ α − a q ∣ < 1 q ( Q + 1 ) \left|\alpha-\dfrac aq\right| < \dfrac{1}{q(Q+1)} α − q a < q ( Q + 1 ) 1 。
最后,若 gcd ( a , q ) = d > 1 \gcd(a,q)=d>1 g cd( a , q ) = d > 1 ,将 a a a 和 q q q 同除以 d d d :所得分数分母更小,与 α \alpha α 的距离相同或更小,因此不失一般性可设 gcd ( a , q ) = 1 \gcd(a,q)=1 g cd( a , q ) = 1 。证明完成;这是一个有限的、构造性的鸽笼论证,不依赖任何未经证明的估计。
大学 实际应用与典型例题 超出数论范围,"把波相加并寻找共振"的同一想法正是信号处理与电子工程中广泛使用的离散傅里叶变换的工作原理;而哈代与拉马努金用圆法早期版本(应用于分拆函数 p ( n ) p(n) p ( n ) )得到的渐近公式 p ( n ) ∼ 1 4 n 3 exp ( π 2 n 3 ) p(n) \sim \dfrac{1}{4n\sqrt3}\exp\!\left(\pi\sqrt{\dfrac{2n}{3}}\right) p ( n ) ∼ 4 n 3 1 exp ( π 3 2 n ) ,被统计力学用来估计总能量固定为 n n n 的不可区分玻色激发系统的微观态数——也就是熵。
例题: 圆法的一个小规模计算
设 A = { 1 , 2 , … , 8 } A=\{1,2,\dots,8\} A = { 1 , 2 , … , 8 } 。利用提取恒等式,计算满足 a + b = 9 a+b=9 a + b = 9 的有序对 ( a , b ) ∈ A 2 (a,b)\in A^2 ( a , b ) ∈ A 2 的个数 r 2 ( 9 ) r_2(9) r 2 ( 9 ) 。
解答 由定理1,r 2 ( 9 ) = ∫ 0 1 F A ( α ) 2 e ( − 9 α ) d α r_2(9)=\int_0^1 F_A(\alpha)^2 e(-9\alpha)\,d\alpha r 2 ( 9 ) = ∫ 0 1 F A ( α ) 2 e ( − 9 α ) d α ,其中 F A ( α ) = ∑ a = 1 8 e ( a α ) F_A(\alpha)=\sum_{a=1}^{8}e(a\alpha) F A ( α ) = ∑ a = 1 8 e ( a α ) ;无需真的解析计算该积分——恒等式已保证它等于直接计数的结果,因此我们可以用组合方法计数,并信任该恒等式。
列出所有满足 a , b ∈ { 1 , … , 8 } a,b\in\{1,\dots,8\} a , b ∈ { 1 , … , 8 } 且 a + b = 9 a+b=9 a + b = 9 的有序对 ( a , b ) (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) ( 1 , 8 ) , ( 2 , 7 ) , ( 3 , 6 ) , ( 4 , 5 ) , ( 5 , 4 ) , ( 6 , 3 ) , ( 7 , 2 ) , ( 8 , 1 ) 。从 1 1 1 到 8 8 8 的每个值唯一确定 b = 9 − a b=9-a b = 9 − a ,且 b b b 总落回 { 1 , … , 8 } \{1,\dots,8\} { 1 , … , 8 } (因为 1 ≤ a ≤ 8 ⇒ 1 ≤ 9 − a ≤ 8 1\le a\le8 \Rightarrow 1\le 9-a\le8 1 ≤ a ≤ 8 ⇒ 1 ≤ 9 − a ≤ 8 ),所以 a a a 的全部 8 8 8 个取值都有效。
因此 r 2 ( 9 ) r_2(9) r 2 ( 9 ) 等于 8 8 8 。这个小例子正是定理1恒等式的具体体现:解析积分与组合计数按构造是同一个数——该方法的真正内容只有当 A A A 变成一个无限或增长的集合(如素数集)、必须估计而非枚举时才显现出来。
例题: 利用狄利克雷定理求最佳有理逼近
取 α = 2 \alpha=\sqrt2 α = 2 、Q = 5 Q=5 Q = 5 ,应用狄利克雷逼近定理,给出满足定理2界的分数 a / q a/q a / q (1 ≤ q ≤ 5 1\le q\le5 1 ≤ q ≤ 5 ),并做数值验证。
解答 我们要找 1 ≤ q ≤ 5 1\le q\le5 1 ≤ q ≤ 5 且 ∣ 2 − a / q ∣ < 1 / ( q ⋅ 6 ) |\sqrt2-a/q|<1/(q\cdot6) ∣ 2 − a / q ∣ < 1/ ( q ⋅ 6 ) (取定理界中的 Q + 1 = 6 Q+1=6 Q + 1 = 6 )的 a / q a/q a / q 。试 q = 5 q=5 q = 5 :最接近 5 2 ≈ 7.0711 5\sqrt2\approx7.0711 5 2 ≈ 7.0711 的整数是 a = 7 a=7 a = 7 ,给出 7 / 5 = 1.4 7/5=1.4 7/5 = 1.4 。
验证界:∣ 2 − 7 / 5 ∣ = ∣ 1.41421 … − 1.4 ∣ ≈ 0.01421 |\sqrt2-7/5|=|1.41421\ldots-1.4|\approx0.01421 ∣ 2 − 7/5∣ = ∣1.41421 … − 1.4∣ ≈ 0.01421 ,而定理保证 1 / ( 5 ⋅ 6 ) = 1 / 30 ≈ 0.0333 1/(5\cdot6)=1/30\approx0.0333 1/ ( 5 ⋅ 6 ) = 1/30 ≈ 0.0333 ;确实 0.01421 < 0.0333 0.01421<0.0333 0.01421 < 0.0333 ,界成立且留有余量,正如所保证的那样。
这个 7 / 5 7/5 7/5 实际上是 2 = [ 1 ; 2 , 2 , 2 , … ] \sqrt2=[1;2,2,2,\dots] 2 = [ 1 ; 2 , 2 , 2 , … ] 连分数两步后的渐近分数,这正是它逼近得如此好的原因——狄利克雷的鸽笼证明本身不需要连分数即可成立,但总能得到同等好的逼近。因此 7 5 \dfrac{7}{5} 5 7 是一个有效的见证。
常见错误. 定理1的提取恒等式是精确的 ——圆法的困难本身并不在这个恒等式上。真正的困难,也是二元哥德巴赫至今未解而三元哥德巴赫已被证明的原因,完全在于估计 次弧的贡献:需要证明 ∣ F A ( α ) ∣ |F_A(\alpha)| ∣ F A ( α ) ∣ 在次弧上确实很小,而对于像素数集这样刚性、且只有两个加数的情形,至今没有人找到足以超越平凡界的估计。 历史注记
哈代与拉马努金在1918年引入圆法,用以求出分拆函数 p ( n ) p(n) p ( n ) 的渐近公式;此后哈代与李特尔伍德在1920年代将其推广为一般方法("Partitio Numerorum"),用以攻克华林问题——把每个大整数写成有界个数的 k k k 次幂之和。1930年代,维诺格拉多夫改进了次弧上所需的指数和估计;从1930年代到1950年代,华罗庚系统地应用并改进了该方法以研究华林—哥德巴赫问题,即把整数表示为素数幂之和。
斯里尼瓦瑟·拉马努金 华罗庚
研究前沿 截至 2026 年
圆法的瓶颈一直是次弧估计,过去十年出现了真正的突破:2015—2016年,让·布尔甘、西普里安·德梅特尔与拉里·古思(以及独立地由特雷弗·伍利通过"高效同余"方法)完整证明了维诺格拉多夫均值定理猜想,给出了控制华林型问题中次弧的外尔和均值的本质最优界——这是对方法本身的结构性升级,而不仅是针对某一应用。另外,格林与陶哲轩的转移原理(2008年)把圆法风格的傅里叶分析从整数迁移到素数,证明了素数中含有任意长的等差数列;而赫尔弗戈特在2013年解决三元哥德巴赫猜想(建立在相关主题哥德巴赫猜想 之上)仍是该方法在素数加性问题上的旗舰成果;打磨其显式常数、并将类似的主弧/次弧技术推向仍未解决的二元情形,是目前活跃的研究方向。
在提取恒等式 r k ( n ) = ∫ 0 1 F A ( α ) k e ( − n α ) d α r_k(n) = \int_0^1 F_A(\alpha)^k\, e(-n\alpha)\, d\alpha r k ( n ) = ∫ 0 1 F A ( α ) k e ( − n α ) d α 中,为什么这个积分恰好等于 r k ( n ) r_k(n) r k ( n ) ?
因为 ∫ 0 1 e ( m α ) d α \int_0^1 e(m\alpha)d\alpha ∫ 0 1 e ( m α ) d α 在 m = 0 m=0 m = 0 时为 1 1 1 ,否则为 0 0 0 ,所以展开 F A ( α ) k F_A(\alpha)^k F A ( α ) k 并积分后只剩下和为 n n n 的元组 因为该积分无论 n n n 为何总是等于 ∣ A ∣ k |A|^k ∣ A ∣ k 因为 F A ( α ) F_A(\alpha) F A ( α ) 是关于 α \alpha α 的 n n n 次多项式 因为圆法对 r k ( n ) r_k(n) r k ( n ) 只能给出近似值,永远不是精确值 设 A = { 1 , 2 , … , 8 } A=\{1,2,\dots,8\} A = { 1 , 2 , … , 8 } ,满足 a + b = 9 a+b=9 a + b = 9 的有序对 ( a , b ) ∈ A 2 (a,b)\in A^2 ( a , b ) ∈ A 2 的个数 r 2 ( 9 ) r_2(9) r 2 ( 9 ) 是多少?
根据狄利克雷逼近定理,对任意实数 α \alpha α 与整数 Q ≥ 1 Q\ge1 Q ≥ 1 ,保证存在什么?
满足 1 ≤ q ≤ Q 1\le q\le Q 1 ≤ q ≤ Q 且 ∣ α − a q ∣ ≤ 1 q ( Q + 1 ) \left|\alpha - \dfrac{a}{q}\right| \le \dfrac{1}{q(Q+1)} α − q a ≤ q ( Q + 1 ) 1 的有理数 a / q a/q a / q 使得 q α q\alpha q α 本身是整数的整数 q q q 整除 α \alpha α 分子的素数 p ≤ Q p\le Q p ≤ Q 以 α \alpha α 为中心、长度恰为 1 / Q 1/Q 1/ Q 的主弧 在统计力学中,分拆函数的哈代–拉马努金渐近公式 p ( n ) ∼ 1 4 n 3 exp ( π 2 n 3 ) p(n) \sim \dfrac{1}{4n\sqrt3}\exp\!\left(\pi\sqrt{\dfrac{2n}{3}}\right) p ( n ) ∼ 4 n 3 1 exp ( π 3 2 n ) 有助于估计什么?
总能量固定为 n n n 的不可区分玻色激发系统的微观态数(从而是熵) 氢原子基态的精确能量 n n n 的素因子个数n n n 是素数的概率