Davenport 圆法 · 专题演示 · 高中生详解版

圆法完整演示:用"五个平方和"从零讲一遍The Hardy–Littlewood circle method, worked end-to-end on sums of five squares (k=2, s=5)

这篇不是翻译某一章,而是把圆法的每一步,针对最简单、最透明的情形 \(k=2,\ s=5\)(把整数写成五个平方之和)从头演示一遍。书里圆法的论证之所以让人觉得"不等式放缩很不自然",是因为一般 \(k\) 的工具(Weyl、Hua)必须照顾最坏情形,把所有具体的东西都藏进了符号里。而 \(k=2\) 时每一个量都能写成你看得见、算得出的具体数:高斯和恰好等于 \(\sqrt q\),"次弧上 \(T\) 很小"是一行因式分解 + 一个等比数列就能看出来的。看完这一篇,再回头看正文那些 \(\le\),就知道每一步在干什么、为什么"刚好够用"。

读这篇需要的全部预备知识加减乘除、平方与平方根、等比数列求和、三角函数 \(\sin\cos\)、复数 \(i^2=-1\)。其余从零讲起。涉及的"放缩"我都会先说为什么要放放掉的是什么会不会放过头

0. 我们到底要证明什么

问题:给一个很大的整数 \(N\),问它能不能写成五个平方之和 \[ N=x_1^2+x_2^2+x_3^2+x_4^2+x_5^2, \qquad x_i\in\mathbb Z. \] 更进一步,我们想知道有多少种写法。把这个写法个数记作 \(r(N)\)(\(x_i\) 允许取负数和 \(0\),不同顺序算不同写法——这样计数最干净)。

我们的目标是证明一条渐近公式

\[ r(N)\;\sim\;\underbrace{\frac{4\pi^2}{3}\,N^{3/2}}_{\text{几何主项}}\;\times\;\underbrace{\mathfrak S(N)}_{\text{算术修正}},\qquad \mathfrak S(N)>0. \]

这里两个因子各有极其朴素的含义,先认个脸熟(后面会一一兑现):

几何主项 \(\frac{4\pi^2}{3}N^{3/2}\)
就是半径 \(\sqrt N\) 的五维球面的"表面积"。如果完全不管整除性、把 \(x_i\) 当成连续变量,落在球面 \(\sum x_i^2=N\) 上的点"有多少",答案就是这个表面积。这是"连续世界"的解的个数。
算术修正 \(\mathfrak S(N)\)(奇异级数,singular series)
真实的 \(x_i\) 是整数,要受"模 \(p\) 同余"的约束。\(\mathfrak S(N)\) 就是把每个素数 \(p\) 处"\(N\) 能否凑成五个平方 \(\bmod\ p^a\)"的难易程度乘起来的修正因子。对五个平方,它永远 \(>0\)——这正是"每个大数都写得出来"的算术原因。
先说清一个常见误会"五个平方"这件事本身,拉格朗日(Lagrange, 1770)早就证明了:每个正整数都是四个平方之和,五个当然更没问题。所以圆法在这里不是为了证明"存在"——存在是已知的。圆法真正给出的是上面那条精确的计数渐近公式,而且它是一套通用机器:同一套流程换成 \(k=3\)(立方)、\(k=4\) 也能转,那些情形就没有拉格朗日式的初等定理了。我们用"五平方和"当教具,是因为它每个零件都透明,能让你看清机器怎么转。

热身:先把机器空转一遍 —— 最简单的 \(k=1,\ s=1\)

在给机器加上 \(k=2,\ s=5\) 的"负载"之前,先让它空转一圈:取 \(k=1,\ s=1\)(把 \(N\) 写成一个一次方)。这是整套圆法最极端的退化情形——所有零件都塌缩成最朴素的样子,正好让你先认全"骨架",再看满载时它们如何一个个变难。(下面要用的记号 \(e(\theta)=\cos2\pi\theta+i\sin2\pi\theta\)、以及"积分挑出傅里叶系数"这件事,第 1 节会从头细讲;这里先借最简单的情形认个脸。)

问题本身:答案显然是 1

\(k=1,\ s=1\) 就是把 \(N\) 写成 \(N=x^1=x\)。解只有 \(x=N\) 这一个,所以表法数 \(r(N)=1\),毫无悬念。有趣的是看圆法怎样把这个 \(1\) 算出来。

圆法在这里 = 一行正交关系

仿照后面的做法,令 \(T(\alpha)=\sum_{x=-P}^{P}e(\alpha x)\)(注意指数里是 \(x\),不是 \(x^2\)),取 \(P\ge N\)。\(s=1\) 时不需要任何五次方、不需要展开,直接:

\[ r(N)=\int_0^1 T(\alpha)e(-N\alpha)\,d\alpha=\sum_{x=-P}^{P}\int_0^1 e\big((x-N)\alpha\big)d\alpha=\sum_{x=-P}^{P}[\,x=N\,]=1. \]

最后一步用的,正是第 1 节将记为 \((1)\) 的正交关系(\(\int_0^1 e(m\alpha)d\alpha\) 在 \(m=0\) 时取 \(1\)、否则取 \(0\))。换句话说:

骨架就这一句\(s=1\) 的圆法 = 正交关系本身。 公式 \(r(N)=\int_0^1 T^s e(-N\alpha)d\alpha\) 在 \(s=1\) 时退回成"提取 \(T(\alpha)\) 的第 \(N\) 个傅里叶系数"。对一般 \(s\),\(r(N)\) 是 \(T(\alpha)^s\) 的第 \(N\) 个系数——难就难在 \(T^s\) 没有封闭式、必须估计;\(s=1\) 时 \(T\) 自己就是答案。
0½1T(α)=Σ e(αx) (k=1,P=10)= Dirichlet 核,峰值 2P+1=21½处仅=1⅓处=0
图 0a \(k=1\) 的指数和 \(T(\alpha)=\sum_{x=-P}^{P}e(\alpha x)\)(这里 \(P=10\))就是 Dirichlet 核 \(\dfrac{\sin((2P+1)\pi\alpha)}{\sin(\pi\alpha)}\)——一个封闭公式。它只有 \(\alpha=0,1\) 一根主峰(高 \(2P+1=21\)),在 \(\tfrac13,\tfrac23\) 处恰为 \(0\)、\(\tfrac12\) 处仅为 \(1\):没有图 1a 那片"尖峰森林"。原因见下表——\(k=1\) 的高斯和在 \(q>1\) 时全为 \(0\)。

每个零件怎么退化

① \(T\) 是等比数列(这正是 \(k\ge2\) 缺的东西)
\(k=1\) 时 \(T(\alpha)\) 有上面那个封闭式。\(k\ge2\) 起 \(\sum e(\alpha x^2)\) 再没有这种封闭式,只能靠高斯和给"大小"、靠 Weyl 给"上界"——文档里所有"难看的放缩"都从这里来。\(k=1\) 一概不需要。
② 高斯和塌成只剩 \(q=1\)
文档里 \(S(q,a)=\sum_{r\bmod q}e(ar^2/q)=\sqrt q\)(第 3 节)。换成 \(k=1\):\(S(q,a)=\sum_{r=0}^{q-1}e(ar/q)=q\) 当 \(q=1\)、\(=0\) 当 \(q>1\)。所以除 \(\alpha=\) 整数外没有任何次级尖峰(正是图 0a 的样子)。没有尖峰森林 ⟹ 没有次弧要对付(第 2、5、6 节那一大套全省了)。
③ 主项也全部 \(=1\)
奇异级数只剩 \(q=1\) 一项:\(\mathfrak S(N)=1\);奇异积分(连续世界里只有 \(x=N\) 一个点)\(=1\)。几何主项的增长指数 \(\dfrac sk-1=\dfrac11-1=0\),故 \(r(N)\sim\text{常数}\cdot N^0=1\)——解数不随 \(N\) 增长,恒为 \(1\),与精确值吻合。

插话:\(k=1\) 的"高斯和"为什么是 \(0\)(而不是 \(\sqrt q\))

上面 ② 那句"\(S(q,a)=0\)(\(q\gt1\))"是 \(k=1\) 全部魔力的来源——它让尖峰森林塌成一根,值得单独讲清。\(k=1\) 的高斯和是

\[ S(q,a)=\sum_{r=0}^{q-1}e\Big(\frac{ar}{q}\Big),\qquad (a,q)=1. \]

当 \((a,q)=1\) 时 \(e(a/q)\) 是一个本原 \(q\) 次单位根,于是 \(r=0,1,\dots,q-1\) 让 \(e(ar/q)\) 恰好取遍全部 \(q\) 个 \(q\) 次单位根各一次。所以 \(S(q,a)\) 就是"全部 \(q\) 次单位根之和":

起点=终点 → 和=0k=1:Σe(ar/q),5 个方向均匀铺开首尾相接成正五边形,回到原点 ⇒ S(5,a)=0合矢量 = √5 ≈ 2.236k=2:Σe(ar²/q),方向聚成 3 类(0°,±72°)部分叠加、不闭合 ⇒ |S(5,a)|=√5
图 0b 把 \(q=5\) 的高斯和当成单位向量首尾相接。左(\(k=1\)):\(e(ar/q)\) 的 \(5\) 个方向均匀分布(\(0^\circ,72^\circ,\dots\)),接成一个正五边形、回到原点,所以 \(S(5,a)=0\)。右(\(k=2\)):指数里是 \(r^2\),方向变成 \(72^\circ\times(r^2\bmod5)=0^\circ,72^\circ,288^\circ,288^\circ,72^\circ\)——只剩 \(3\) 类、还成对重复,于是部分叠加、不闭合,合矢量长 \(\sqrt5\)。这就是 §3.2 那个 \(|S|=\sqrt q\) 的来历。
\(k=1\) 的抵消比 \(k=2\) 的"平方根抵消"还彻底关键差别在指数:\(k=1\) 的 \(ar\) 是 \(r\) 的线性函数,\(r\) 取遍 \(\Rightarrow ar\bmod q\) 也均匀取遍 \(\Rightarrow\) 正多边形 \(\Rightarrow\) 完全抵消为 \(0\);\(k=2\) 的 \(ar^2\) 里 \(r^2\bmod q\) 均匀(有的值重复、有的缺席),向量聚拢、只抵消到剩 \(\sqrt q\)。所以 \(k=1\) 的每个 \(q\gt1\) 尖峰振幅 \(\dfrac{S(q,a)}{q}=\dfrac0q=0\)(连峰都不剩),而 \(k=2\) 是 \(\dfrac{\sqrt q}q=\dfrac1{\sqrt q}\ne0\)(矮峰仍在,构成图 1a 的森林)。"尖峰森林塌成一根"的根,就在这里。

④ 怎么估这条曲线的大小:封闭式直接给界

\(k=1\) 估 \(|T|\) 不需要任何不等式技巧——既然有了封闭式,看分子分母就够了。这正是 \(k=1\) 全部的"轻松"所在。把大小写成分式:

\[ |T(\alpha)|=\frac{\big|\sin\big((2P+1)\pi\alpha\big)\big|}{|\sin(\pi\alpha)|}. \]
两个一眼可得的界
  • 分子永远 \(\le 1\):\(|\sin(\cdot)|\le1\),免费。
  • 分母有下界(Jordan 不等式):在 \([0,\tfrac12]\) 上 \(\sin(\pi x)\) 是凹的、落在弦 \(y=2x\) 之上,所以 \[ |\sin(\pi\alpha)|=\sin(\pi\|\alpha\|)\ \ge\ 2\|\alpha\|, \] 其中 \(\|\alpha\|\) 是 \(\alpha\) 到最近整数的距离(端点对得上:\(\|\alpha\|=0\to0\)、\(\|\alpha\|=\tfrac12\to1\))。
两个界一拼,再配上"总项数 \(2P+1\)"这个平凡上界(三角不等式),得到完整估计
\[ \boxed{\;|T(\alpha)|\ \le\ \min\Big(2P+1,\ \frac{1}{2\|\alpha\|}\Big)\;} \]
0¼½¾1红:上界 min(2P+1, 1/(2‖α‖))  蓝:真实 |T(α)| (P=10)平台 = 2P+1 = 21双曲线 1/(2‖α‖) 衰减½处上界=|T|=1(相切)
图 0c 把上界画出来:红线 是 \(\min\big(2P+1,\tfrac1{2\|\alpha\|}\big)\),蓝线 是真实的 \(|T(\alpha)|\)(\(P=10\))。红线稳稳罩住蓝线:整数附近是高 \(2P+1=21\) 的平台,离开整数后变成双曲线 \(\tfrac1{2\|\alpha\|}\) 衰减,到 \(\alpha=\tfrac12\) 降到 \(1\)、与 \(|T|\) 的局部峰相切。这条与 \(P\) 无关的衰减正是"完全相消"的定量刻画——也就是后面公式 \((7)\) 反复要用的那条界。

怎么读它:

靠近整数(\(\|\alpha\|\lesssim\tfrac1P\))
\(\min\) 取 \(2P+1\),就是那根高 \(2P+1\) 的主峰(图 0a 的尖峰),峰宽约 \(1/P\)。
远离整数
\(\min\) 取 \(\dfrac1{2\|\alpha\|}\),与 \(P\) 无关!哪怕 \(P\) 再大,只要 \(\alpha\) 离整数有固定距离,\(|T|\) 就被一个 \(O(1)\) 常数压住——这就是等比数列的完全相消(单位向量转着圈走、几乎抵消干净),正是图 0a 中间那片矮振荡。

顺带,均方是精确的(不是估计):\(\displaystyle\int_0^1|T(\alpha)|^2 d\alpha=2P+1\)(正交关系一步到位),这是华罗庚不等式在 \(k=1\) 的平凡版。

这条界就是后面的公式 (7)——\(k=1\) 是 \(k=2\) 的内核上面的 \(\min\big(2P+1,\tfrac1{2\|\alpha\|}\big)\) 你在后面会再遇到:它正是第 5 节的公式 \((7)\)。\(k=2\) 做完差分降次 \(x^2-y^2=(x-y)(x+y)\) 后,内层和 \(\sum_y e(2\alpha h\,y)\) 是一个 \(k=1\) 的等比数列,用的就是这同一条界(把 \(\alpha\) 换成 \(2\alpha h\))。换句话说,\(k=1\) 的估计是 \(k=2\) 估计的"内核":\(k=2\) 没有封闭式,于是先用一次差分把二次和降回等比数列,再套用这条 \(\min\) 界。\(k=1\) 之所以轻松,是因为它一上来就是等比数列、无需降次;Weyl 的全部工作,就是"把高次和差分降到能用这条界为止"。

⑤ 完整空跑一遍:按 §1→§7 把机器开一遍(看它在哪卡住)

有了上界,就能把正文 §1→§7 的整套流程在 \(k=1,\ s=1\) 上"空载试转"一遍。结论早知道是 \(r(N)=1\),但重点是看每个齿轮怎么转、又在哪一个齿轮上诚实地卡住——那个卡点正是 §5、§6 存在的理由。设 \(T(\alpha)=\sum_{x=-P}^{P}e(\alpha x)\),\(P\ge N\)。

(对应 §1)写成积分。 \(r(N)=\int_0^1 T(\alpha)e(-N\alpha)\,d\alpha\)(\(s=1\),不用乘方)。

(对应 §2–§3)切弧,尖峰森林塌成一根。 \(k=1\) 的高斯和 \(S(q,a)=0\)(\(q>1\)),所以只剩 \(q=1\) 一根尖峰(图 0a/0b)。于是主弧 \(\mathfrak M=\{\|\alpha\|\le\delta\}\)(取 \(\delta\sim1/P\),只此一段),次弧 \(\mathfrak m\) 是其余全部。

(对应 §3–§4)主弧 → 主项。 在 \(\mathfrak M\) 上把求和换成它的连续模型(这是 \(k=1\) 版的主弧近似 \((4)\),其中高斯和因子是 \(S(1,0)/1=1\))。即用积分逼近求和:

\[ T(\alpha)=\sum_{x=-P}^{P}e(\alpha x)\ \approx\ I(\alpha):=\int_{-P}^{P}e(\alpha t)\,dt. \]
新函数 ①:先把 \(I(\alpha)\) 算到底这是一个初等积分,逐步写全: \[ I(\alpha)=\int_{-P}^{P}e^{2\pi i\alpha t}\,dt=\Big[\frac{e^{2\pi i\alpha t}}{2\pi i\alpha}\Big]_{t=-P}^{t=P}=\frac{e^{2\pi i\alpha P}-e^{-2\pi i\alpha P}}{2\pi i\alpha}=\frac{2i\sin(2\pi P\alpha)}{2\pi i\alpha}=\frac{\sin(2\pi P\alpha)}{\pi\alpha}, \] 末步用 \(e^{i\theta}-e^{-i\theta}=2i\sin\theta\)。它在 \(\alpha\to0\) 的极限是 \(I(0)=2P\)(用 \(\sin u\approx u\)),正是峰高。和 \(T\) 的关系注意 \(T\) 自己也有封闭式(§④ 的 Dirichlet 核 \(T(\alpha)=\sin((2P+1)\pi\alpha)/\sin(\pi\alpha)\))。主弧上 \(\|\alpha\|\to0\),\(\sin(\pi\alpha)=\pi\alpha+O(\alpha^3)\)、且 \((2P+1)\pi\alpha\) 与 \(2\pi P\alpha\) 仅差 \(\pi\alpha\),故 \(T(\alpha)=I(\alpha)\big(1+O(P\alpha)\big)\)——两者主阶一致(这个相对误差 \(O(P\alpha)\) 在 \(\alpha\sim1/P\) 处是 \(O(1)\),是 \(k=1\) 第一处"不干净",记下,后面算总账时再面对)。
-10-50510蓝:精确 T(α)=Σe(αx) (Dirichlet 核,周期 1)红:新函数 I(α)=∫₋ₚᵖ e(αt)dt = sin(2πPα)/(πα)中央峰:I≈T(主弧依据)↑ T 在每个整数都重现高 21 的峰(周期)而 I 已按 1/(π|α|) 衰减(长尾不归零)
图 0d 新函数 \(I(\alpha)\)(红,连续模型)与精确的 \(T(\alpha)\)(蓝,Dirichlet 核),\(P=10\),\(\alpha\in[-10,10]\)。中央 \(\alpha=0\) 附近两条线几乎重合——这正是主弧上"用 \(I\) 代 \(T\)"的依据。但往外立刻分道扬镳:\(T\) 是周期函数(周期 \(1\)),在每个整数处都重现同样高 \(2P+1=21\) 的峰;\(I\) 不周期,只在 \(\alpha=0\) 有一座峰,此后按 \(1/(\pi|\alpha|)\) 缓慢衰减、长尾迟迟不归零。\(I\) 这条 \(1/|\alpha|\) 长尾正是上文"奇异积分只条件收敛"的根源;而在真正的周期 \(T\) 里,这些长尾被卷回(\(\bmod 1\))叠成了相邻整数处的那些峰。两图相差的部分,全在主弧之外。
差这么大,\(I\) 还算"近似"吗?——算,但它是局部(主弧)近似,不是全局近似 看到图 0d 差异这么夸张,先别慌:\(I\) 从来只声称在主弧上近似 \(T\),圆法也只在主弧上用它。
  • 为什么只在主弧成立。 \(I\) 来自"求和换积分" \(\sum_x e(\alpha x)\approx\int e(\alpha t)dt\),而它成立的唯一条件是相位每步 \(2\pi\alpha\) 变化要慢——只有 \(\alpha\) 靠近 \(0\)(\(\bmod 1\))才行,即 \(|\alpha|\lesssim 1/P\)(正是主弧)。图 0d 中央两线重合就是这一小段。
  • 为什么差得夸张。 \(T\) 周期:\(T(m+\beta)=\sum_x e(mx)e(\beta x)=T(\beta)\),于是它在每个整数都重放同一座峰;\(I\) 不周期,是把那一座峰"摊开不折叠"的连续版,只有 \(\alpha=0\) 一座、其余 \(1/|\alpha|\) 衰减。但 \(\alpha=1\equiv0\,(\bmod1)\) 在圆上就是同一条主弧,那里 \(I\approx T\);图把 \(T\) 的一个周期复制了 21 份给你看,而积分只在 \([0,1]\) 一个周期上做。
  • 差异落在哪。 \(I\) 与 \(T\) 差得最凶的地方(整数处的复现峰、中间的振荡背景)全部在次弧——而次弧上圆法根本不用 \(I\),改用 §5 的大小上界 \(1/(2\|\alpha\|)\)。"近似只在主弧、上界只在次弧"正是主弧/次弧二分的全部意义。
所以把 \(I\) 当 \(T\) 的全局近似是范畴错误:它是主弧上的局部主阶模型,圆法的设计就是确保"只在它准的地方用它"。(也正因 \(I\) 准的区域 \(\sim1/P\) 太窄、尾巴又肥又慢,把它延拓到 \(\pm\infty\) 定义奇异积分才只得到条件收敛——见上一框。)
新积分 ②:奇异积分 \(J\) 完整求值把主弧的积分限放宽到 \(\pm\infty\)(这一步本身的误差留到总账),定义并计算奇异积分 \[ J:=\int_{-\infty}^{\infty}I(\alpha)\,e(-N\alpha)\,d\alpha=\int_{-\infty}^{\infty}\frac{\sin(2\pi P\alpha)}{\pi\alpha}\,e^{-2\pi iN\alpha}\,d\alpha. \] 把 \(e^{-2\pi iN\alpha}=\cos(2\pi N\alpha)-i\sin(2\pi N\alpha)\) 代入,分成实、虚两部分: \[ J=\underbrace{\int_{-\infty}^{\infty}\frac{\sin(2\pi P\alpha)\cos(2\pi N\alpha)}{\pi\alpha}\,d\alpha}_{\text{实部}}\ -\ i\underbrace{\int_{-\infty}^{\infty}\frac{\sin(2\pi P\alpha)\sin(2\pi N\alpha)}{\pi\alpha}\,d\alpha}_{\text{虚部}}. \] 虚部 = 0:被积函数 \(=\dfrac{\text{偶}\times\text{奇}}{\text{奇}}\)(\(\sin(2\pi P\alpha)/(\pi\alpha)\) 是偶函数,\(\sin(2\pi N\alpha)\) 是奇函数)\(=\) 奇函数,在对称区间上积分为 \(0\)。
实部:用积化和差 \(\sin A\cos B=\tfrac12[\sin(A+B)+\sin(A-B)]\): \[ \text{实部}=\frac12\int_{-\infty}^{\infty}\frac{\sin\big(2\pi(P+N)\alpha\big)+\sin\big(2\pi(P-N)\alpha\big)}{\pi\alpha}\,d\alpha. \] 再用狄利克雷积分(这是这里唯一要借的"已知积分"):\(\displaystyle\int_{-\infty}^{\infty}\frac{\sin(ux)}{x}\,dx=\pi\,\mathrm{sgn}(u)\),取 \(x=\alpha,\ u=2\pi c\) 并除以 \(\pi\) 得 \(\displaystyle\int_{-\infty}^{\infty}\frac{\sin(2\pi c\alpha)}{\pi\alpha}d\alpha=\mathrm{sgn}(c)\)。于是 \[ J=\frac12\big[\,\mathrm{sgn}(P+N)+\mathrm{sgn}(P-N)\,\big]=\frac12[\,1+1\,]=1\qquad(0\lt N\lt P). \] (若 \(N\gt P\) 则 \(\mathrm{sgn}(P-N)=-1\),\(J=0\)——正确地反映"\(N\) 超出范围则无解"。)所以 奇异积分 \(J=\mathbf 1_{[-P,P]}(N)=1\),含义就是"连续世界里满足 \(t=N\) 的解数 = 一个点"。
先说清:这个"级数"是哪儿冒出来的——因为主弧不止一条 到这里我们只盯着 \(\alpha=0\) 一条主弧,于是"奇异级数"显得突兀。其实每个分母小的有理点 \(a/q\) 都是一条主弧(回想图 1a 的尖峰森林——每座尖峰就是一条弧)。主项是所有这些主弧贡献之和,奇异级数就是这个"对所有有理点求和"的产物。
单条 \(a/q\) 主弧(\(k=1\) 现推):令 \(\alpha=\tfrac aq+\beta\),按 \(x=r+qm\) 分组,因 \(e(\tfrac aq\,qm)=e(am)=1\) 故 \(e(\tfrac aq x)=e(\tfrac aq r)\),于是 \[ T\Big(\tfrac aq+\beta\Big)\approx\Big(\underbrace{\tfrac1q\sum_{r=0}^{q-1}e(\tfrac{ar}{q})}_{=\,S(q,a)/q}\Big)\underbrace{\int e(\beta t)\,dt}_{=\,I(\beta)},\qquad\text{峰高}\ \propto\ \frac{S(q,a)}{q}. \] 再乘 \(e(-N\alpha)=e(-N\tfrac aq)\,e(-N\beta)\) 积分:算术相位 \(e(-Na/q)\) 提出来,解析因子 \(J=\int I(\beta)e(-N\beta)d\beta\)(每条弧都一样)也提出来,得 \[ \int_{a/q\,\text{附近}}T\,e(-N\alpha)\,d\alpha\ \approx\ \frac{S(q,a)}{q}\,e\Big(-\frac{Na}{q}\Big)\cdot J. \] 把所有有理点 \(a/q\) 的算术因子加起来(\(J\) 提到外面),就是下面的奇异级数,主项 \(=\mathfrak S\cdot J\)。所以 \(\mathfrak S\) = "每个有理点那座尖峰(峰高 \(=S(q,a)/q\))贡献的总和":\(J\) 管峰长什么形状(解析),\(\mathfrak S\) 管哪些有理点有峰、各占多重(算术)。
细推:单条 \(a/q\) 主弧上 \(T\) 的峰高 \(=\dfrac{S(q,a)}{q}I(\beta)\) 怎么来的 这就是 §4.1(\(q=1\))那套"求和换积分"的带 \(q\) 版本,多出来的那一项正是高斯和。设 \(\alpha=\tfrac aq+\beta\)(\(\beta\) 小),逐步走:

① 拆相位。 \(e(\alpha x)=e\big(\tfrac aq x\big)\,e(\beta x)\):一个算术相位 \(\times\) 一个解析相位

② 按余数分组。 每个 \(x\) 唯一写成 \(x=r+qm\)(\(r=0,1,\dots,q-1\),\(m\) 为整数,使 \(x\in[-P,P]\))。

③ 算术相位只认 \(r\): \[ e\Big(\tfrac aq x\Big)=e\Big(\tfrac aq(r+qm)\Big)=e\Big(\tfrac aq r\Big)\,\underbrace{e(am)}_{=1}=e\Big(\tfrac{ar}q\Big), \] 因为 \(am\) 是整数 \(\Rightarrow e(am)=1\)。同一余数类里它是常数

④ 解析相位缓变。 \(e(\beta x)=e\big(\beta(r+qm)\big)\),\(\beta\) 小、随 \(m\) 慢慢变。于是整体重排成

\[ T\Big(\tfrac aq+\beta\Big)=\sum_{r=0}^{q-1}e\Big(\tfrac{ar}q\Big)\sum_{m}e\big(\beta(r+qm)\big). \]

⑤ 内层和换成积分(关键,带一个 \(1/q\))。 固定 \(r\),\(m\) 走一步则 \(x=r+qm\) 跳 \(q\) 格,所以内层是把光滑函数 \(e(\beta t)\) 每隔 \(q\) 取一个样本再相加。\(\beta\) 小时它在长度 \(q\) 内几乎不变,黎曼和 \(\approx\dfrac1q\times\)积分(间距为 \(q\) \(\Rightarrow\sum f\approx\frac1{\text{间距}}\int f\)): \[ \sum_{m}e\big(\beta(r+qm)\big)\approx\frac1q\int_{-P}^{P}e(\beta t)\,dt=\frac1q\,I(\beta)\qquad(\text{与 }r\ \text{无关}). \]

⑥ 提出来。 内层与 \(r\) 无关,提到求和号外,外层只剩高斯和: \[ T\Big(\tfrac aq+\beta\Big)\approx\Big(\sum_{r=0}^{q-1}e\big(\tfrac{ar}q\big)\Big)\frac1q\,I(\beta)=\frac{S(q,a)}{q}\,I(\beta). \]

取 \(\beta=0\):峰高 \(=\dfrac{S(q,a)}{q}\,I(0)=\dfrac{S(q,a)}{q}\cdot 2P\),正比于 \(S(q,a)/q\)。检验 \(q=1\):\(S(1,0)=1\),回到 \(T(\beta)\approx I(\beta)\)、峰高 \(2P\),正是 §4.1(新函数①)那一步——可见单条 \(a/q\) 主弧的推导和 \(q=1\) 是同一招,只是按余数分组后多结出一个高斯和 \(S(q,a)\)。
新对象 ③:奇异级数 \(\mathfrak S\) 完整求值定义 \[ \mathfrak S=\sum_{q\ge1}\ \sum_{\substack{a\bmod q\\(a,q)=1}}\frac{S(q,a)}{q}\,e\Big(-\frac{Na}{q}\Big),\qquad S(q,a)=\sum_{r=0}^{q-1}e\Big(\frac{ar}{q}\Big). \] 逐项算 \(S(q,a)\):
  • \(q=1\):唯一项 \(a=0\),\(S(1,0)=\sum_{r=0}^{0}e(0)=1\),贡献 \(\dfrac{1}{1}e(0)=1\)。
  • \(q\ge2\):\(S(q,a)=\sum_{r=0}^{q-1}\big(e(a/q)\big)^r\) 是公比 \(e(a/q)\) 的等比数列,求和 \(=\dfrac{e(a/q)^q-1}{e(a/q)-1}=\dfrac{e(a)-1}{e(a/q)-1}=\dfrac{0}{\ne0}=0\)(分子 \(e(a)=1\) 因 \(a\) 是整数;分母 \(\ne0\) 因 \((a,q)=1,\ q\ge2\))。所以全部 \(q\ge2\) 的项逐个为零
于是 奇异级数 \(\mathfrak S=1\)。注意:这不是"没有奇异级数",而是除 \(q=1\) 外每条主弧的峰高 \(S(q,a)/q=0\)、贡献恰好全为零(图 0b:\(q\ge2\) 的高斯和 = 单位根首尾相接成正多边形 = 0;图 0a 也确实只剩一根峰)。机制完整走了一遍,只是结果退化成一项。
它的意义(同余障碍):\(\mathfrak S=\prod_p\beta_p\),\(\beta_p=\) "模 \(p^a\) 的解密度"。\(k=1\) 的方程是 \(x\equiv N\ (\mathrm{mod}\ p^a)\),永远恰好一个解 \(\Rightarrow\beta_p=1\Rightarrow\mathfrak S=1\):没有任何同余障碍、也没有额外加成,故算术修正为 \(1\)。(对照 §4.4 的 \(k=2\):多条主弧"活着",\(\mathfrak S\) 是它们带符号之和 = 素数局部密度之积 \(\ne1\);图 1c 里 \(\tfrac14,\tfrac34\) 那些朝下的尖峰,就是这些活弧的带符号贡献。)

两半相乘得 主项 \(=\mathfrak S\cdot J=1\times1=1\)。验证增长指数 \(\tfrac sk-1=0\),故主项 \(\sim N^0=\) 常数,与 \(r(N)\equiv1\) 吻合。表面上主弧这半给出了正确答案——但慢着,下面要泼一盆冷水。

慢着:这个"主项"其实也站不稳——奇异积分只条件收敛 \(J=1\) 的数值没错(Fourier 反演:连续世界里 \(t=N\) 只是一个点),但"主弧罩住 \(J\)、其余都是小误差"这套逻辑对 \(s=1\) 站不住。症结在模型 \(I(\alpha)=\dfrac{\sin(2\pi P\alpha)}{\pi\alpha}\) 只按 \(1/|\alpha|\) 衰减
  • 奇异积分 \(J=\int_{-\infty}^{\infty}I\,e(-N\alpha)d\alpha\) 只是条件收敛(因 \(\int|I|\,d\alpha=\infty\));
  • 真正的主弧 \(|\alpha|\le\delta\sim1/P\) 只截到 \(J\) 的一部分。把它老实算出来(虚部仍为奇函数 \(=0\),实部用积化和差 + 正弦积分 \(\mathrm{Si}(x)=\int_0^x\frac{\sin u}{u}du\)): \[ \int_{-\delta}^{\delta}I\,e(-N\alpha)d\alpha=\frac1\pi\Big[\mathrm{Si}\big(2\pi(P{+}N)\delta\big)+\mathrm{Si}\big(2\pi(P{-}N)\delta\big)\Big]=1+O(1), \] 因为 \(\delta\sim1/P\) 时宗量 \(2\pi(P\pm N)\delta\) 只是 \(O(1)\),\(\mathrm{Si}\) 离极限 \(\tfrac\pi2\) 还差 \(O(1)\);补到 \(\pm\infty\) 的尾巴同样是 \(O(1)\)。都不小。
是 Fourier 反演(≈ 正交关系)把它钉成 \(1\) 的,不是"主项主导"推出来的。精确判据:奇异积分 \(\int v(\beta)^s e(-N\beta)d\beta\)(\(v(\beta)=\int_{|t|\le P}e(\beta t^k)dt\sim|\beta|^{-1/k}\),故 \(v^s\sim|\beta|^{-s/k}\))绝对收敛 \(\iff s/k\gt1\iff s\gt k\)。\(k{=}1,s{=}1\) 正卡在 \(s=k\) 临界,只条件收敛;要让主项站得稳必须 \(s\gt k\)(\(k{=}1\) 即 \(s\ge2\),\(k{=}2\) 即 \(s\ge3\))。所以 \(s=1\) 不只是次弧卡住——主弧(奇异积分)同样退化,两半都没立住,这才是 \(r(N)=1\) 只能靠一行正交关系拿到的真正原因。

(对应 §5–§6)次弧 → 误差:机器在这里卡住。 按正文套路,把次弧积分逐步放缩。第一步用三角不等式 + \(|e(-N\alpha)|=1\)(丢掉相位):

\[ \Big|\int_{\mathfrak m}T(\alpha)\,e(-N\alpha)\,d\alpha\Big|\ \le\ \int_{\mathfrak m}|T(\alpha)|\cdot\underbrace{|e(-N\alpha)|}_{=1}\,d\alpha\ =\ \int_{\mathfrak m}|T(\alpha)|\,d\alpha. \]

第二步代入图 0c 的上界 \(|T(\alpha)|\le\dfrac1{2\|\alpha\|}\)。次弧 \(\mathfrak m=\{\delta\le\|\alpha\|\le\tfrac12\}\) 关于 \(\alpha\) 与 \(1-\alpha\) 对称(整数 \(0\) 与 \(1\) 两侧各一份),故

\[ \int_{\mathfrak m}|T|\,d\alpha\ \le\ 2\int_{\delta}^{1/2}\frac{1}{2\alpha}\,d\alpha\ =\ \int_{\delta}^{1/2}\frac{d\alpha}{\alpha}\ =\ \big[\ln\alpha\big]_{\delta}^{1/2}\ =\ \ln\frac{1}{2\delta}. \]

取主弧宽度 \(\delta\sim 1/P\)(峰宽),得

\[ \Big|\int_{\mathfrak m}T\,e(-N\alpha)\,d\alpha\Big|\ \le\ \ln\frac{1}{2\delta}\ \sim\ \ln\frac{P}{2}\ \sim\ \log P. \]
误差上界 \(O(\log P)\) 压不过主项 \(1\)对照 §6:那里 \(\int_{\mathfrak m}|T|^5\le(\sup_{\mathfrak m}|T|)\int_0^1|T|^4\ll P^{3/4}\cdot P^2=P^{11/4}\ll P^3\)(主项)。整台机器要运转(主项站得稳 + 次弧压得下),靠三件 \(k=1,\ s=1\) 都不具备的东西:
取胜要件\(k{=}2,\ s{=}5\)(正文)\(k{=}1,\ s{=}1\)(这里)
主弧:奇异积分绝对收敛(需 \(s\gt k\))\(5\gt2\) ✓\(s=k=1\) 临界,只条件收敛 ✗
次弧:比平凡界 \(P\) 小一个Weyl:\(\sup|T|\ll P^{3/4}\) ✓只有 \(\tfrac1{\|\alpha\|}\),无幂节省(边界 \(\|\alpha\|\sim\tfrac1P\) 处仍 \(\sim P\))✗
次弧:变量数 \(s\) 够多\(s=5\)(§6 门槛)✓\(s=1\) 太少,连 \(\int|T|^s\) 都发散 ✗

(对应 §7)总账。

\[ r(N)=\underbrace{\int_{\mathfrak M}T e(-N\alpha)}_{\text{主项}=1}+\underbrace{\int_{\mathfrak m}T e(-N\alpha)}_{\text{真值小,但上界只给 }O(\log P)}. \]

真相是两块加起来精确 \(=1\),所以次弧的真实贡献确实趋于 \(0\)——但那是 \(e(-N\alpha)\) 在次弧上的相位相消,而 \(|T|\) 的大小上界把相位丢了,于是只看到 \(O(\log P)\)。机器的"主项"半边成立,"误差"半边靠这条上界关不上

空跑的收获\(k=1,\ s=1\) 的主弧给出了正确的主项 \(1\),但那是 Fourier 反演(≈ 正交关系)钉出来的——奇异积分只条件收敛,主弧并没真正"罩住"它;次弧的上界估计则原地卡死。要让机器真正运转,缺的齿轮有三个:① 主项绝对收敛需 \(s\gt k\)(\(s=1\) 不够,主弧这半就先垮了);② 幂节省的次弧界 \(P^{1-\delta}\) 需 \(k\ge2\)(Weyl,§5);③ 足够多的变量 \(s\)(把次弧压到主项之下,§6 算出 \(s=5\))。所以 \(k=1,\ s=1\) 不靠渐近机器、而靠一行正交关系拿到精确的 \(1\)——正因为这几个齿轮全缺位。它像一台空载试转的发动机:你能听见主轴想转,但"主项绝对收敛"和"误差"两档都还没挂上——而那几档,正是正文 §3–§6 在装。
真正的难关是 \(k\),不是 \(s\)固定 \(k=1\)、让 \(s\) 长大,\(r(N)\) 仍然精确(\(T^s\) 是等比数列的幂,照样是有限傅里叶和):例如 \(s=2\) 时 \(N=x_1+x_2\) 给 \(r(N)=N-1\)(增长指数 \(\tfrac sk-1=1\)),依旧一步到位、无需估计。圆法的全部困难,是 \(k\) 从 \(1\) 跨到 \(2\) 那一步——指数和从"等比数列(有封闭式)"变成"二次指数和(只能估计)"。本文档其余部分,就是在 \(k=2\) 这道坎上,把每个零件重新拼起来。

1. 第一步:把"数写法"变成"算积分"

圆法的出发点是一个会让你觉得"凭什么"的操作:把一个计数问题(离散、难)翻译成一个积分问题(连续、可分析)。靠的是一条小学生都能验证的恒等式。

1.1 记号 \(e(\theta)\) 与正交关系

定义. 对实数 \(\theta\),记 \(e(\theta):=\cos(2\pi\theta)+i\sin(2\pi\theta)\)。它是单位圆上的点,满足 \(|e(\theta)|=1\)、\(e(\theta_1)e(\theta_2)=e(\theta_1+\theta_2)\)、对整数 \(m\) 有 \(e(m)=1\)。(它就是 \(e^{2\pi i\theta}\)。)

核心引理只有一条,是整个圆法的"开关":

正交关系. 对整数 \(m\), \[ \int_0^1 e(m\alpha)\,d\alpha=\begin{cases}1,&m=0,\\[2pt]0,&m\neq 0.\end{cases}\tag{1} \]
证. \(m=0\) 时被积函数恒为 \(1\),积分为 \(1\)。\(m\neq0\) 时, \[ \int_0^1 e^{2\pi i m\alpha}d\alpha=\Big[\frac{e^{2\pi i m\alpha}}{2\pi i m}\Big]_0^1=\frac{e^{2\pi i m}-1}{2\pi i m}=\frac{1-1}{2\pi i m}=0, \] 因为 \(m\) 是整数时 \(e^{2\pi i m}=1\)。

这条恒等式是一台"是不是零"的探测器:你给它一个整数 \(m\),它积分出来告诉你 \(m\) 是不是 \(0\)。我们马上要用它来探测"\(x_1^2+\dots+x_5^2-N\) 是不是 \(0\)"。

1.2 指数和 \(T(\alpha)\),以及把 \(r(N)\) 写成积分

设 \(P=\lfloor\sqrt N\rfloor\)(这样每个 \(|x_i|\le\sqrt N\) 自动成立)。定义这一篇的主角——指数和

\[ T(\alpha):=\sum_{x=-P}^{P}e(\alpha x^2). \tag{2} \]

它把"所有 \(\le N\) 的平方数"打包成一个 \(\alpha\) 的函数。现在做关键一步:把它五次方再乘 \(e(-N\alpha)\) 展开:

\[ T(\alpha)^5 e(-N\alpha)=\sum_{x_1,\dots,x_5}e\big(\alpha(x_1^2+\cdots+x_5^2-N)\big). \]

对 \(\alpha\) 在 \([0,1]\) 上积分,用正交关系 \((1)\):每一项里 \(m=x_1^2+\cdots+x_5^2-N\) 是整数,只有当这五个平方真的加起来等于 \(N\) 时那一项才积出 \(1\),否则积出 \(0\)。于是

\[ \boxed{\;r(N)=\int_0^1 T(\alpha)^5\,e(-N\alpha)\,d\alpha\;}\tag{3} \]

左边是我们想要的离散计数,右边是一个积分。圆法的全部工作,就是估计右边这个积分。

"圆"在哪里因为 \(e(\alpha)\) 是单位圆上的点,\(\alpha\) 从 \(0\) 跑到 \(1\) 就是绕圆一圈,所以 \((3)\) 是一个沿单位圆的积分——"圆法"得名于此。

1.3 先看一眼:取定 \(N=101\),把这三条曲线画出来

公式 \((3)\) 把 \(r(N)\) 写成了"一条曲线在 \([0,1]\) 上围出的面积"。在抽象地谈"主弧 / 次弧"之前,不如先对一个具体的 \(N\) 把曲线画出来,眼见为实。取 \(N=101\)(于是 \(P=\lfloor\sqrt{101}\rfloor=10\),\(T\) 是 \(2P+1=21\) 项之和,\(T(0)=21\))。剧透结论:下面图 1c 那条上下振荡的曲线,它在 \([0,1]\) 上的净面积(正负相抵后)精确等于 \(r(101)=\mathbf{10640}\)——我们要的答案,就整个藏在这块面积里。

|T(α)| 峰值 |T(0)|=210¼½¾1↑½处塌成0
图 1a \(|T(\alpha)|\):在分母小的分数 \(a/q\) 处冒尖峰,峰高 \(\approx(2P+1)/\sqrt q\)(第 3 节会精确算)。真实图像比"\(1/\sqrt q\) 单调下降"更有趣——\(q=2\)(\(\alpha=\tfrac12\))处反而塌成 0(因为高斯和 \(S(2,a)=0\)),而 \(q=4\)(\(\tfrac14,\tfrac34\))却高过 \(q=3\)。这些"反常"第 3 节都能算准。
|T(α)|⁵ 峰值 21⁵ ≈ 4.08×10⁶0¼½¾1
图 1b 把它五次方之后(注意纵轴峰值已是 \(21^5\approx4.08\times10^6\)):只有 \(q=1\) 的峰还顶天,\(q=4\) 的峰约剩 \(18\%\)、\(q=3\) 约 \(6\%\),其余几乎贴地。乘方把"信号"猛烈放大、把"噪声"压平——这正是为什么后面只需精算少数几根大尖峰(主弧),其余一律证明它小(次弧)。
0¼½¾1Re[ T(α)⁵·e(−Nα) ] 净面积 = r(101) = 10640
图 1c 真正被积的函数 \(\mathrm{Re}\,[\,T(\alpha)^5 e(-N\alpha)\,]\)(积分 \((3)\) 取实部,虚部自动积成 0)。曲线与横轴围成的净面积 = \(r(101)=\mathbf{10640}\)。最高的正尖峰在 \(\alpha=0\);\(\tfrac14,\tfrac34\) 等小分母处也有可见尖峰,且这里恰好朝下——说明各尖峰贡献可正可负,这正是奇异级数 \(\mathfrak S(N)\) 在"算账"。这几根尖峰合起来就是"主弧"(小分母分数全算);它们之间随 \(e(-N\alpha)\) 剧烈振荡、正负相消的平坦背景,才是"次弧"。
|T(α)| 放大 (α∈[0,0.06]):孤峰宽度 ≈ 1/N ≈ 0.0100.010.020.030.040.050.06≈1/N
图 1d 把 \(\alpha=0\) 的尖峰放大:它的宽度只有约 \(1/N\approx0.01\),是一座极窄的孤峰。主弧就是套在这种尖峰外的、总长极小的小区间——这解释了一件初看反直觉的事:主弧总长很小,却贡献了几乎全部面积
这四张图已经把圆法的策略说尽了把 \([0,1]\) 分成两类地方:(主弧) 少数几根又高又窄的尖峰,位置在小分母分数处、形状由高斯和决定——这里要精算(第 3–4 节);(次弧) 振荡的背景,又矮又乱——这里只需证明它小(第 5–6 节)。下面第 2 节就把"尖峰在哪、为什么高、背景为什么矮"这件事正式说清楚。

2. 全局图景:\(T(\alpha)\) 在哪里大、在哪里小

积分 \((3)\) 看似无从下手,但 \(T(\alpha)\) 的形状其实非常有规律。先把它看明白,后面所有"放缩"就都顺理成章了。

2.1 在整点和简单分数处,\(T\) 会"共振"变大

先看 \(\alpha=0\):每一项 \(e(0)=1\),所以 \(T(0)=2P+1\approx 2\sqrt N\),非常大。这是因为 \(2P+1\) 个单位向量方向完全一致,全加在一起。

再看 \(\alpha\) 接近一个分母小的分数 \(a/q\),比如 \(\alpha=1/4\):此时 \(e(\alpha x^2)=e(x^2/4)\) 只依赖 \(x^2\bmod 4\)。而 \(x^2\bmod4\) 只会是 \(0\) 或 \(1\)(偶数的平方 \(\equiv0\),奇数的平方 \(\equiv1\)),所以一半的项指向 \(e(0)=1\)、一半指向 \(e(1/4)=i\),方向只有两种,又叠加出一个可观的和。分母 \(q\) 越小,"方向种类"越少、叠加越猛,\(T\) 越大。

反过来,\(\alpha\) 若是个"分母很大或干脆无理"的数,\(e(\alpha x^2)\) 的方向杂乱无章地铺满整个圆,正负相消,\(T\) 就很小。这就是平方根抵消(square-root cancellation):\(M\) 个方向随机的单位向量相加,长度典型地只有 \(\sqrt M\) 而不是 \(M\)。这里 \(M\approx 2\sqrt N\),所以"典型"的 \(|T|\) 只有 \(\big(2\sqrt N\big)^{1/2}\approx N^{1/4}\) 量级——比峰值 \(2\sqrt N\) 小得多。

0 1 1/2 1/3 2/3 1/4 3/4 峰 ≈ 2√N (α=0) 其余地方杂乱、很小 ≈ N¼
图 2 \(|T(\alpha)|\) 在 \([0,1]\) 上的形状(理想示意,与图 1a 实算结果对照):在分母小的分数 \(a/q\)(\(0,\tfrac12,\tfrac13,\tfrac14,\dots\))处冒出尖峰,峰高随 \(q\) 增大而降低;其余地方因平方根抵消而很小。

2.2 圆法的中心思想:把积分劈成"信号"和"噪声"

主弧 / 次弧(major / minor arcs) 把积分区间 \([0,1]\) 切成两部分:
  • 主弧 \(\mathfrak M\):所有那些尖峰附近的小区间(\(\alpha\) 离某个分母 \(\le Q\) 的分数 \(a/q\) 很近)。这里 \(T\) 大、有漂亮的近似公式,贡献主项——这是"信号",要算准(用等式)。
  • 次弧 \(\mathfrak m\):剩下的所有 \(\alpha\)。这里 \(T\) 小、毫无封闭公式——这是"噪声",算不出也不想算,只需证明它加起来比主项小(用不等式)。
这就一举回答了你最初的困惑:"为什么主弧用等式、次弧只用不等式?"——因为对次弧,等式既不可能(没有公式)也没必要(我们只要它小)。所有"难看的放缩"都只发生在次弧上,它们唯一的任务就是把噪声压到信号以下。

3. 主弧:在尖峰附近,\(T\) 长什么样(高斯和登场)

这一节是 \(k=2\) 最漂亮的地方:尖峰的高度可以精确算出,而且答案出奇地干净。

3.1 把 \(T(a/q+\beta)\) 拆成"算术部分 × 解析部分"

设 \(\alpha=\dfrac aq+\beta\),其中 \(a/q\) 是既约分数、\(\beta\) 很小。按 \(x\bmod q\) 把求和分组:写 \(x=qy+r\)(\(r=0,1,\dots,q-1\))。则 \(x^2=q^2y^2+2qyr+r^2\),于是

\[ e\Big(\tfrac aq x^2\Big)=e\Big(\tfrac aq r^2\Big)\quad(\text{因为含 }q\text{ 的项让 }e\text{ 取整数幂}=1). \]

而 \(e(\beta x^2)\) 随 \(x\) 缓慢变化(\(\beta\) 小)。把同余类内的和近似成积分,得到主弧上的近似公式

主弧近似 \[ T\Big(\tfrac aq+\beta\Big)\;\approx\;\underbrace{\frac1q\,S(q,a)}_{\text{算术:高斯和}}\;\cdot\;\underbrace{\int_{-P}^{P}e(\beta t^2)\,dt}_{\text{解析:积分}},\qquad S(q,a):=\sum_{r=0}^{q-1}e\Big(\tfrac aq r^2\Big).\tag{4} \]

这就把尖峰的高度分解成两个独立的零件:一个纯算术的高斯和(Gauss sum) \(S(q,a)\),和一个纯解析的积分。下面分别算。

3.2 高斯和恰好是 \(\sqrt q\) —— 你能亲手验证

高斯(Gauss). 若 \(q\) 为奇数且 \(\gcd(a,q)=1\),则 \(\;|S(q,a)|=\sqrt q\)。

这是整篇里最该记住的一个数:\(k=2\) 的高斯和,长度精确等于 \(\sqrt q\)——不多不少,没有任何 \(\le\)、没有任何 \(P^\varepsilon\) 的污染。一般 \(k\) 的"Weyl 和" \(\sum e(ax^k/q)\) 只能放缩到 \(\le C\,q^{1-1/k+\varepsilon}\),正是那些放缩让正文显得"不自然"。\(k=2\) 时它退化成一个干净的等式,所以这里你能看得一清二楚。

例(\(q=5,a=1\),亲手算). \(r=0,1,2,3,4\) 时 \(r^2\bmod5=0,1,4,4,1\)。所以 \[ S(5,1)=e(0)+2e(\tfrac15)+2e(\tfrac45)=1+4\cos72^\circ. \] 代入 \(\cos72^\circ=0.3090\dots\):\(S(5,1)=1+4(0.3090)=2.236\dots=\sqrt5\)。恰好 \(\sqrt5\),分毫不差。
\(r\)01234
\(r^2\bmod5\)01441
\(e(r^2/5)\)方向\(1\)\(72^\circ\)\(288^\circ\)\(288^\circ\)\(72^\circ\)
为什么是 \(\sqrt q\) 而不是 \(q\)注意上面 \(q=5\) 个向量并没有全指一个方向(那样才会得到 \(q\)),而是分散成几个方向、部分抵消,最后长度落在 \(\sqrt q\)——这正是平方根抵消的精确版。它定量地说明了 2.1 节的直觉。验证 \(|S|=\sqrt q\) 的标准技巧是算 \(|S|^2=S\overline S=\sum_{r,r'}e(\tfrac aq(r^2-r'^2))\),用 \(r^2-r'^2=(r-r')(r+r')\) 换元后得到一个等比数列,正好加出 \(q\)。这跟下面 5.1 节压次弧用的是同一招

把 \((4)\) 取绝对值,尖峰高度就是 \[ \Big|T\big(\tfrac aq+\beta\big)\Big|\approx\frac{|S(q,a)|}{q}\Big|\int e(\beta t^2)dt\Big|=\frac{\sqrt q}{q}\,|\cdots|=\frac{1}{\sqrt q}\,|\cdots|. \] 当 \(\beta=0\) 时积分 \(\approx 2P\approx2\sqrt N\),所以 \[ \Big|T\big(\tfrac aq\big)\Big|\approx \frac{2\sqrt N}{\sqrt q}. \tag{5} \] \(q\) 越大峰越矮——和图 1a 完全吻合。\(q=1\)(即 \(\alpha=0\))时回到峰值 \(2\sqrt N\)。

4. 主项:把所有尖峰的贡献加起来

第 3 节给了主弧近似 \((4)\),但只写了"\(\approx\)"、没说怎么来。这一节把每一步都补上:先把 \((4)\) 真正推出来,再代进积分 \((3)\),看主项怎样一段段长出来,最后把两个常数(\(\tfrac{4\pi^2}{3}\) 和 \(\mathfrak S(N)\))亲手算出来

4.1 先补上:主弧近似 \((4)\) 是怎么推出来的

设 \(\alpha=\dfrac aq+\beta\)。把 \((2)\) 的求和按 \(x\) 模 \(q\) 的余数分组:每个 \(x\) 唯一写成 \(x=r+qm\),其中 \(r\in\{0,1,\dots,q-1\}\)、\(m\) 是整数。代入 \(x^2\):

\[ \tfrac aq x^2=\tfrac aq(r+qm)^2=\tfrac aq r^2+\underbrace{2arm+aqm^2}_{\text{整数}},\qquad\text{所以}\quad e\big(\tfrac aq x^2\big)=e\big(\tfrac aq r^2\big). \]

也就是说,算术相位 \(e(\tfrac aq x^2)\) 只认 \(x\) 的余数 \(r\)——同一余数类里它是常数。于是

\[ T\big(\tfrac aq+\beta\big)=\sum_{r=0}^{q-1}\underbrace{e\big(\tfrac aq r^2\big)}_{\text{随 }r\text{ 跳变}}\;\sum_{m:\,r+qm\in[-P,P]}\underbrace{e\big(\beta(r+qm)^2\big)}_{\text{随 }m\text{ 缓变}}. \]
把内层和换成积分固定 \(r\),让 \(m\) 走一步,自变量 \(x=r+qm\) 就跳 \(q\) 格。所以内层和是把光滑函数 \(g(t)=e(\beta t^2)\) 每隔 \(q\) 取一个样本再相加。主弧上 \(\beta\) 极小(我们把主弧宽度取成 \(|\beta|\le\tfrac1{qP}\)),\(g\) 在长度 \(q\) 的一段里几乎不变,于是这串样本和 \(\approx\dfrac1q\times(\text{积分})\)(这正是黎曼和,只不过反过来用): \[ \sum_{m}e\big(\beta(r+qm)^2\big)\;\approx\;\frac1q\int_{-P}^{P}e(\beta t^2)\,dt\qquad(\text{与 }r\text{ 无关}). \]

内层和既然与 \(r\) 无关,就能提到求和号外,外层只剩高斯和:

\[ T\big(\tfrac aq+\beta\big)\approx\Big(\sum_{r=0}^{q-1}e\big(\tfrac aq r^2\big)\Big)\cdot\frac1q\int_{-P}^{P}e(\beta t^2)dt=\frac1q\,S(q,a)\int_{-P}^{P}e(\beta t^2)dt. \]

这正是 \((4)\)。注意它天生就是"算术零件 \(\times\) 解析零件"的形状——这是后面一切"分家"的根源。

4.2 代入积分:算术与解析"自动分家"

主弧 \(\mathfrak M\) 是一堆互不相交的小区间 \(\mathfrak M(q,a)=\{\alpha:|\alpha-\tfrac aq|\le\tfrac1{qP}\}\)(\(q\le Q\),\(\gcd(a,q)=1\))。在每个小区间上写 \(\alpha=\tfrac aq+\beta\),代 \((4)\)。关键观察是那个负指数也会分家

\[ e(-N\alpha)=e\big(-N\tfrac aq\big)\,e(-N\beta). \]

于是单个小区间的积分裂成"一个纯算术常数 \(\times\) 一个纯解析积分":

\[ \int_{\mathfrak M(q,a)}T^5e(-N\alpha)\,d\alpha\approx\underbrace{\Big(\tfrac{S(q,a)}{q}\Big)^5 e\big(-N\tfrac aq\big)}_{\text{只含 }q,a}\cdot\underbrace{\int_{-1/(qP)}^{1/(qP)}\!\Big(\int_{-P}^{P}e(\beta t^2)dt\Big)^5 e(-N\beta)\,d\beta}_{\text{只含 }\beta}. \]

现在对 \(a\)(\(\bmod\,q\))求和、再对 \(q\le Q\) 求和。算术那半攒成一个级数,解析那半(几乎与 \(q,a\) 无关,把积分限放宽到 \(\pm\infty\))提到外面,得

\[ \int_{\mathfrak M}T^5e(-N\alpha)d\alpha\;\approx\;\underbrace{\mathfrak S(N)}_{\text{奇异级数(算术)}}\;\times\;\underbrace{J(N)}_{\text{奇异积分(解析)}},\tag{14} \] \[ \mathfrak S(N)=\sum_{q\ge1}\ \sum_{\substack{a=1\\(a,q)=1}}^{q}\Big(\tfrac{S(q,a)}{q}\Big)^5 e\big(-N\tfrac aq\big),\qquad J(N)=\int_{-\infty}^{\infty}\Big(\int_{-P}^{P}e(\beta t^2)dt\Big)^5 e(-N\beta)\,d\beta. \]

"分家"完成。下面两节分别把 \(J(N)\)(第 4.3 节)和 \(\mathfrak S(N)\)(第 4.4 节)算清楚。

4.3 奇异积分 \(J(N)\):把 \(N\) 抽干,剩一个几何常数

记 \(v(\beta)=\displaystyle\int_{-P}^{P}e(\beta t^2)dt\)。做两步换元,目的就是把所有 \(N\) 都提到积分外面

换元 ①:\(t=\sqrt N\,u\) 则 \(dt=\sqrt N\,du\)、\(\beta t^2=\beta N u^2\)、积分限 \(\pm P\to\pm1\)(因 \(P\approx\sqrt N\)): \[ v(\beta)=\sqrt N\int_{-1}^{1}e(\beta N u^2)\,du. \] 换元 ②:\(\gamma=\beta N\)(即 \(\beta=\gamma/N,\ d\beta=d\gamma/N\)) 代进 \(J(N)\),并用 \(v(\beta)^5=N^{5/2}\big(\int_{-1}^1 e(\gamma u^2)du\big)^5\)、\(e(-N\beta)=e(-\gamma)\): \[ J(N)=\int_{-\infty}^{\infty}N^{5/2}\Big(\int_{-1}^{1}e(\gamma u^2)du\Big)^5 e(-\gamma)\,\frac{d\gamma}{N}=N^{3/2}\underbrace{\int_{-\infty}^{\infty}\Big(\int_{-1}^{1}e(\gamma u^2)du\Big)^5 e(-\gamma)\,d\gamma}_{=:\ \mathcal J\ (\text{纯数字})}. \]

所以 \(J(N)=\mathcal J\cdot N^{3/2}\),其中 \(\mathcal J\) 是一个不含 \(N\) 的常数。要算 \(\mathcal J\) 不必硬碰这个复积分——它有一个一眼能懂的几何意思:

\(\mathcal J\) = "连续世界里 \(\sum x_i^2=N\) 的解密度"把 \(x_i\) 当成连续变量。五维实心球 \(\sum x_i^2\le R^2\) 的体积是 \(V(R)=\dfrac{\pi^{5/2}}{\Gamma(7/2)}R^5\)。用 \(N=R^2\) 换成 \(V=\dfrac{\pi^{5/2}}{\Gamma(7/2)}N^{5/2}\),再对 \(N\) 求导——增量 \(\dfrac{dV}{dN}\) 正是"\(N\) 每涨 \(1\),球面 \(\sum x_i^2=N\) 上新扫过的那层壳里有多少解",也就是解密度: \[ \frac{dV}{dN}=\frac52\cdot\frac{\pi^{5/2}}{\Gamma(7/2)}N^{3/2}=\frac{\pi^{5/2}}{\Gamma(5/2)}N^{3/2}\quad(\text{用 }\Gamma(7/2)=\tfrac52\Gamma(5/2)). \] 圆法的奇异积分 \(J(N)\) 自动就是这个量,所以 \[ \mathcal J=\frac{\pi^{5/2}}{\Gamma(5/2)}=\frac{\pi^{5/2}}{\tfrac34\sqrt\pi}=\frac43\pi^2\approx13.16,\qquad\text{用 }\Gamma(\tfrac52)=\tfrac32\cdot\tfrac12\sqrt\pi=\tfrac34\sqrt\pi. \]
一次免费的对账(验证 \(\tfrac{4\pi^2}{3}\) 没错). 把渐近式 \(r(N)\approx\tfrac{4\pi^2}{3}N^{3/2}\mathfrak S(N)\) 对 \(N\le X\) 求和(\(\mathfrak S\) 平均约为 \(1\)),应当重现"球里的整点总数"。两条路算同一个数:
  • 逐项求和:\(\displaystyle\sum_{N\le X}\tfrac{4\pi^2}{3}N^{3/2}\approx\tfrac{4\pi^2}{3}\int_0^X t^{3/2}dt=\tfrac{4\pi^2}{3}\cdot\tfrac25 X^{5/2}=\tfrac{8\pi^2}{15}X^{5/2}.\)
  • 直接数球:\(\displaystyle\sum_{N\le X}r(N)=\#\{\textstyle\sum x_i^2\le X\}\approx V(\sqrt X)=\tfrac{\pi^{5/2}}{\Gamma(7/2)}X^{5/2}=\tfrac{8\pi^2}{15}X^{5/2}\)(\(\Gamma(\tfrac72)=\tfrac{15}{8}\sqrt\pi\)).
两边分毫不差都是 \(\tfrac{8\pi^2}{15}X^{5/2}\)——常数 \(\tfrac{4\pi^2}{3}\) 经得起检验。✓

4.4 奇异级数 \(\mathfrak S(N)\):拆成每个素数的"局部解密度"

解析那半已经给出几何主项。现在轮到算术那半 \(\mathfrak S(N)\)——它修正"整数不像实数那么自由"。三步走。

第一步:拆成素数乘积(Euler 乘积)。 记 \(A(q)=\sum_{(a,q)=1}\big(\tfrac{S(q,a)}{q}\big)^5 e(-N\tfrac aq)\)。用中国剩余定理可证 \(A\) 是积性的:\(A(q_1q_2)=A(q_1)A(q_2)\)(\(q_1,q_2\) 互素)。积性函数求和就能拆成每个素数的"几何级数":

\[ \mathfrak S(N)=\sum_{q\ge1}A(q)=\prod_{p}\Big(1+A(p)+A(p^2)+A(p^3)+\cdots\Big)=\prod_p\beta_p(N). \]

第二步:把 \(\beta_p\) 翻译成"数同余解"。 用 \((1)\) 的离散版(在 \(\mathbb Z/p^J\mathbb Z\) 上 \(\tfrac1{p^J}\sum_{a\bmod p^J}e(\tfrac aq m)=[\,p^J\mid m\,]\))可证:

\[ 1+A(p)+\cdots+A(p^J)=\frac{M(p^J)}{p^{4J}},\qquad M(p^J):=\#\{\mathbf x\bmod p^J:\ x_1^2+\cdots+x_5^2\equiv N\}. \]

所以 \(\beta_p(N)=\lim_{J\to\infty}\dfrac{M(p^J)}{p^{4J}}\) 就是"模 \(p^J\) 的世界里,五个平方凑出 \(N\) 的解所占的比例"(分母 \(p^{4J}\) 是五个自由变量去掉一个方程约束后的总数 \(p^{5J}/p^{J}\))。这就把抽象的高斯和级数变回了能用手数的同余方程

第三步:奇素数处把 \(M(p)\) 数出来。 对奇素数 \(p\),五个平方的同余方程有一条经典计数公式(用勒让德符号 \(\big(\tfrac{\cdot}{p}\big)\)):

\[ M(p)=\#\{x_1^2+\cdots+x_5^2\equiv N \ (\mathrm{mod}\ p)\}=p^4+\Big(\tfrac Np\Big)\,p^2,\qquad\text{于是}\quad \beta_p(N)=1+\Big(\tfrac Np\Big)\frac1{p^2}. \]

(奇素数处密度从 \(J=1\) 起就稳定,故 \(\beta_p=M(p)/p^4\)。\(p\mid N\) 时 \(\big(\tfrac Np\big)=0\),\(\beta_p=1\)。)这个 \(\pm p^{-2}\) 就是算术对几何主项的毫厘修正

例(\(p=3\),纯手算验证 \(M(3)=81\pm9\)). 模 \(3\) 时平方只有两种值:\(0^2=0\)、\(1^2=2^2=1\)。设 \(k\) 个变量取非零(每个非零变量有 \(2\) 种取值 \(x=1,2\),零变量 \(1\) 种),则平方和 \(\equiv k\pmod 3\),方案数 \(=\binom5k 2^k\)。
\(N\bmod3\)需要 \(k\equiv\)\(k\) 取值\(M(3)=\sum\binom5k2^k\)对照 \(81+(\tfrac N3)9\)
\(1\)\(1\)\(1,4\)\(\binom51 2+\binom54 2^4=10+80=90\)\((\tfrac13)=+1\Rightarrow90\) ✓
\(2\)\(2\)\(2,5\)\(\binom52 2^2+\binom55 2^5=40+32=72\)\((\tfrac23)=-1\Rightarrow72\) ✓
\(0\)\(0\)\(0,3\)\(\binom50+\binom53 2^3=1+80=81\)\(3\mid N\Rightarrow81\) ✓
公式 \(p^4+(\tfrac Np)p^2\) 在 \(p=3\) 上逐格命中。对应 \(\beta_3=\tfrac{90}{81}=\tfrac{10}9\)、\(\tfrac{72}{81}=\tfrac89\) 或 \(1\)。

正性与下界(这才是我们真正要的)。 把上面拼起来:

① 每个因子都 \(>0\)
\(\beta_p=1\pm p^{-2}\),而 \(p\ge3\) 时 \(p^{-2}\le\tfrac19<1\),所以 \(\beta_p\ge1-\tfrac19=\tfrac89>0\)。素数 \(2\) 处要单独算(模 \(8\) 的平方只能是 \(0,1,4\)),结果同样是一个确定的正数。没有任何素数把解数掐成 \(0\)——这正是"五个变量足够多,模每个 \(p^a\) 都解得开"。
② 乘积收敛、有正下界
把奇素数因子连乘,因 \(\big(\tfrac Np\big)=\pm1\),每个因子夹在 \(1-p^{-2}\) 与 \(1+p^{-2}\) 之间,于是 \[ \prod_{p\ \text{odd}}\big(1-\tfrac1{p^2}\big)\ \le\ \prod_{p\ \text{odd}}\beta_p\ \le\ \prod_{p\ \text{odd}}\big(1+\tfrac1{p^2}\big). \] 左边 \(\ge\prod_{p}(1-p^{-2})=\dfrac1{\zeta(2)}=\dfrac6{\pi^2}\approx0.61\),右边 \(\le\prod_p(1+p^{-2})=\dfrac{\zeta(2)}{\zeta(4)}=\dfrac{15}{\pi^2}\approx1.52\)。两个端点都是正常数
结论:存在正常数 \(c_1,c_2\) 使 \(\;c_1\le\mathfrak S(N)\le c_2\;\),特别地 \(\;\mathfrak S(N)>0\;\) 对每个 \(N\) 成立。

注意我们没有用到平方根抵消的精细估计,只用了"每个因子 \(=1\pm p^{-2}\)"——但若只想要"收敛",连这点都不必,高斯和的粗界 \(|S(q,a)|=\sqrt q\) 已给出 \(|A(p)|\le p\cdot p^{-5/2}=p^{-3/2}\),\(\sum_p p^{-3/2}<\infty\) 即可。这里 \(5\) 个变量带来的 \(q^{-5/2}\) 是收敛的关键:变量太少(如 \(s=3\))这个级数就微妙得多(三平方和真的有例外 \(N=4^a(8b+7)\)),这是 \(s\) 不能太小的第一个信号;真正卡死下界的要求来自次弧,见第 6 节。

5. 次弧:为什么 \(T\) 在这里一定很小(全篇核心)

到这里"信号"已经算清楚了。剩下唯一的活,是证明"噪声"——次弧上的积分——比信号小。这就是你觉得"放缩不自然"的那些步骤的家。我们就 \(k=2\) 把它从因式分解一路推到底,让每个 \(\le\) 都有名有姓。

5.1 一行因式分解 + 一个等比数列:看穿 \(|T|^2\)

不要直接碰 \(T\)(它的相位没法控制),而是看 \(|T|^2=T\overline T\):

\[ |T(\alpha)|^2=\sum_{x}\sum_{y}e\big(\alpha(x^2-y^2)\big). \]

关键一步——\(k=2\) 专属的因式分解:令 \(h=x-y\),则 \(x^2-y^2=(x-y)(x+y)=h(2y+h)\)。于是

\[ |T|^2=\sum_{h}e(\alpha h^2)\sum_{y}e(2\alpha h\,y). \tag{6} \]
为什么这一步是"自然"的我们想要的"抵消"藏在相位里,而相位 \(x^2-y^2\) 本身是二次的、难处理。\((x-y)(x+y)\) 把它变成关于 \(y\) 的一次式——一次式的指数和就是等比数列,能精确求和!这不是凑出来的技巧,而是"降次"这个唯一可行方向的必然结果。(一般 \(k\) 的 Weyl 不等式就是把这一步重复 \(k-1\) 次差分来逐步降次,所以才"层层放缩"显得繁琐;\(k=2\) 只需一次,干净见底。)

内层 \(\sum_{y}e(2\alpha h\,y)\) 是公比 \(e(2\alpha h)\) 的等比数列,项数约 \(2P\)。等比数列求和给出

\[ \Big|\sum_{y=-P}^{P}e(2\alpha h\,y)\Big|=\Big|\frac{\sin(2\pi(2P+1)\alpha h)}{\sin(2\pi\alpha h)}\Big|\le\min\Big(2P+1,\ \frac{1}{2\|2\alpha h\|}\Big), \tag{7} \]

其中 \(\|\theta\|\) 表示 \(\theta\) 到最近整数的距离。这个 \(\min\) 的含义非常具体:

对 \((6)\) 取绝对值、用 \(|e(\alpha h^2)|=1\) 扔掉外层相位(这一步只丢相位、不丢大小),得

\[ |T(\alpha)|^2\;\le\;\sum_{|h|\le 2P}\min\Big(2P,\ \frac{1}{2\|2\alpha h\|}\Big). \tag{8} \]

5.2 在尖峰附近"数一数"这个和(把分块过程写全)

现在要估 \((8)\) 右边的和。记 \(\theta=2\alpha\approx\dfrac{2a}{q}\)。被加项是 \(\min\big(2P,\ \tfrac1{2\|\theta h\|}\big)\),关键全在 \(\|\theta h\|\)(\(\theta h\) 到最近整数的距离)何时。这一步在正文里总被一句话带过,下面把它掰成分块来数

把 \(h\) 切成长度 \(q\) 的小段 把 \(1\le h\le 2P\) 切成约 \(\dfrac{2P}{q}+1\) 段,每段是 \(q\) 个连续整数 \(h_0+1,\dots,h_0+q\)。在同一段内,\(\theta h\approx\dfrac{2a}{q}h\) 的小数部分随 \(h\) 走一遍 \(\dfrac{2a}{q},\dfrac{2\cdot2a}{q},\dots\);因为 \(\gcd(2a,q)\) 很小,这 \(q\) 个值基本跑遍 \(\big\{\tfrac0q,\tfrac1q,\dots,\tfrac{q-1}q\big\}\) 各一次。于是同一段里:
  • 最靠近整数的那一两个 \(h\)(即 \(\|\theta h\|\) 几乎为 \(0\))触顶,贡献被 \(\min\) 截在 \(2P\);
  • 其余 \(h\) 的 \(\|\theta h\|\approx\tfrac jq\)(\(j=1,2,\dots\)),贡献 \(\dfrac1{2\|\theta h\|}\approx\dfrac q{2j}\),把它们加起来: \[ \sum_{j=1}^{q/2}\frac{q}{2j}\approx\frac q2\sum_{j=1}^{q/2}\frac1j\approx\frac q2\ln q\ \ll\ q\log q. \]
所以每一段的总贡献 \(\ll\ 2P+q\log q\)。

把约 \(\dfrac{2P}{q}+1\) 段的贡献乘起来("每段 \(\ll 2P+q\log q\)" 乘以"段数 \(\tfrac{2P}q+1\)"):

\[ \sum_{|h|\le 2P}\min\Big(2P,\frac1{2\|\theta h\|}\Big)\ \ll\ \Big(\frac{2P}{q}+1\Big)\big(2P+q\log q\big)\ \ll\ \frac{P^2}{q}+P+q\log q. \]

这三项里中间的 \(P\) 总会被另外两项盖住:\(q\le P\) 时 \(\tfrac{P^2}q\ge P\),\(q>P\) 时 \(q>P\)。把 \(\log\) 吸进 \(P^\varepsilon\),就并成干净的两项,代回 \((8)\):

\[ |T(\alpha)|^2\;\ll\;\underbrace{\frac{P^2}{q}}_{\text{近共振段:}\frac{2P}q\text{ 段}\times 2P}+\underbrace{q}_{\text{每段余项 }q\log q\text{ 累计}}\ \ll\ \Big(\frac{P^2}{q}+q\Big)P^{\varepsilon}. \tag{9} \]
\[ \boxed{\;|T(\alpha)|\;\ll\;\Big(\frac{P}{\sqrt q}+\sqrt q\Big)P^{\varepsilon}\;}\qquad(\text{当 }|\alpha-\tfrac aq|\le\tfrac1{q^2}) \tag{10} \]
读懂这个界——它就是 (5) 的"误差版"主项里 \(\dfrac{P}{\sqrt q}\) 正是第 3 节算出的尖峰高度 \((5)\)(精确等式那里),而 \(\sqrt q\) 是把和"数出来"时不可避免的余项。两件事在这里对上了:主弧近似说"峰高 \(=P/\sqrt q\)",次弧不等式说"除了峰高,再没有更大的东西"。所以 \((10)\) 一点都不神秘——它是 \((5)\) 的诚实加强版,告诉你"\(q\) 一旦不小,\(T\) 就上不去"。

5.3 次弧上 \(q\) 必然不小,于是 \(T\) 必然小

主弧只收了 \(q\le Q\) 的尖峰。取 \(Q=\sqrt P\)。由狄利克雷逼近定理,次弧上的 \(\alpha\) 总能写成 \(|\alpha-a/q|\le1/q^2\) 且 \(q\) 落在 \(Q< q\le P^2/Q\),即 \(\sqrt P< q\le P^{3/2}\)。代进 \((10)\):

\[ \frac{P}{\sqrt q}<\frac{P}{P^{1/4}}=P^{3/4},\qquad \sqrt q\le P^{3/4}, \]
\[ \boxed{\;\sup_{\alpha\in\mathfrak m}|T(\alpha)|\;\ll\;P^{3/4+\varepsilon}\;}\tag{11} \]

记住 \(T\) 的峰值是 \(2P\)。在次弧上它被压到了 \(P^{3/4}\)——从 \(P^1\) 掉到了 \(P^{3/4}\),这 \(P^{1/4}\) 的差距,就是我们要拿去赢比赛的本钱。

6. 决胜局:为什么"刚好"需要 \(s=5\)

这是回答你最初问题的高潮:"那些指数 5、\(3/2\) 不是随便定的,而是'刚好够赢'所需的最小代价"。我们来数一遍指数。

6.1 一个免费的均值估计

除了"次弧上 \(T\) 小",我们还有一个对所有 \(\alpha\) 都成立的平均估计,而且它有纯组合意义:

\[ \int_0^1|T(\alpha)|^4\,d\alpha=\#\{\,x_1^2+x_2^2=x_3^2+x_4^2:\ |x_i|\le P\,\}\ll P^{2+\varepsilon}. \tag{12} \]

(左边按正交关系 \((1)\) 正是这个方程的解数;解数 \(\approx P^2\),主要来自"对角线" \(\{x_1,x_2\}=\{x_3,x_4\}\)。)这步是等式,不丢任何东西,免费。

6.2 把次弧积分"一拆为二"

这是圆法最经典的一招:在次弧上把 \(|T|^5\) 拆成"一个 \(\sup\) × 四个放进平均"——一半用 \((11)\) 的逐点小,一半用 \((12)\) 的平均小:

\[ \Big|\int_{\mathfrak m}T^5 e(-N\alpha)d\alpha\Big|\le\int_{\mathfrak m}|T|^5 d\alpha\le\Big(\sup_{\mathfrak m}|T|\Big)\cdot\int_0^1|T|^4 d\alpha\ll P^{3/4+\varepsilon}\cdot P^{2+\varepsilon}=P^{11/4+\varepsilon}. \tag{13} \]

而主项的大小是 \[ \frac{4\pi^2}{3}N^{3/2}\,\mathfrak S(N)\asymp N^{3/2}=P^{3}. \]

比较两个指数:\(\quad\) 次弧 \(=P^{\,11/4+\varepsilon}\quad\) vs \(\quad\) 主项 \(=P^{\,3}\).
\(\dfrac{11}{4}=2.75<3\),次弧严格小于主项,差距 \(P^{1/4}\)。信号赢了!🎉

6.3 为什么 \(s=4\) 会输,\(s=5\) 才赢

把同样的拆法用到 \(s=4\)(四个平方)上,主项变成 \(N^{4/2-1}=N^{1}=P^{2}\),而次弧

\[ \int_{\mathfrak m}|T|^4\le\Big(\sup_{\mathfrak m}|T|\Big)^2\int_0^1|T|^2 d\alpha\ll \big(P^{3/4}\big)^2\cdot P=P^{5/2}. \]

(用了 \(\int_0^1|T|^2=\#\{x_1^2=x_2^2\}\approx 2P\ll P\)。)于是 \(s=4\) 时

四平方:信号输给噪声次弧 \(=P^{5/2}=P^{2.5}\) 反而大于主项 \(P^{2}\)。噪声淹没信号,这套朴素圆法对 \(s=4\) 失效!要救 \(s=4\) 需要 Kloosterman 的精细改良(把次弧再拆得更细)。而 \(s\ge5\) 时每多一个变量,次弧多乘一个 \(\sup|T|\ll P^{3/4}\)、主项多乘一个 \(P^{1}\),主项每次多涨 \(P^{1/4}\),于是从 \(s=5\) 起信号稳定领先
\(s\)主项 \(=P^{s-2}\)次弧上界谁赢
3\(P^{1}\)\(\gg P^{1}\)噪声(远输)
4\(P^{2}\)\(P^{2.5}\)噪声(输)
5\(P^{3}\)\(P^{2.75}\)信号(赢)
6\(P^{4}\)\(P^{3.75}\)信号(更稳)
这就是 "\(s=5\)" 的全部秘密:它是朴素圆法能让"次弧 \(<\) 主项"成立的第一个 \(s\)。指数 \(\tfrac{11}{4}<3\) 的那 \(\tfrac14\) 富余,就是 \(s\) 从 4 涨到 5 挣来的。

7. 收官:拼出渐近公式

把三件事拼起来:

总账 \[ r(N)=\underbrace{\int_{\mathfrak M}T^5 e(-N\alpha)d\alpha}_{\text{主弧(第 3–4 节,等式算准)}}+\underbrace{\int_{\mathfrak m}T^5 e(-N\alpha)d\alpha}_{\text{次弧(第 5–6 节,不等式压小)}}. \]
  • 主弧 \(=\dfrac{4\pi^2}{3}\,\mathfrak S(N)\,N^{3/2}\,(1+o(1))\),且 \(\mathfrak S(N)\ge c_1>0\);
  • 次弧 \(\ll P^{11/4+\varepsilon}=o\big(N^{3/2}\big)\)。
\[ \boxed{\;r(N)=\frac{4\pi^2}{3}\,\mathfrak S(N)\,N^{3/2}+O\big(N^{11/8+\varepsilon}\big)\;}\;>\;0\quad(N\ \text{充分大}). \]

因为主项是正的、误差严格更小,所以对充分大的 \(N\),\(r(N)>0\):每个充分大的整数都是五个平方之和,而且写法数精确到 \(N^{3/2}\) 的主阶。圆法的机器,转完了一整圈。

8. 回看:那些"不自然的放缩",其实步步必然

把你最初的疑问逐条结清
"为什么从等式退化成不等式?"
只在次弧上退化,因为那里没有公式、也只需要"小"这一个结论。主弧自始至终是等式(第 3–4 节)。
"为什么放掉那么多东西(相位、Cauchy–Schwarz)?"
每次丢的都是相位(角度)——那是对一般 \(\alpha\) 我们唯一控制不了的东西;而且每次只丢得起一个常数或 \(P^\varepsilon\),绝不丢一个 \(P\) 的幂。\((6)\to(8)\) 丢外层相位、\((8)\to(11)\) 数 \(\min\),都是如此。
"为什么指数是 \(\tfrac34,\tfrac{11}4,3\) 这些怪数?"
\(\tfrac34\) 来自高斯和的 \(q^{-1/2}\) 配上次弧的 \(q\) 范围;\(3\) 是几何主项 \(N^{3/2}=P^3\);\(\tfrac{11}4\) 是"一个 \(\sup\) + 四次均值"的和。它们是算出来的,不是凑的——而 \(\tfrac{11}4<3\) 这道窄缝,就是定理成立的全部理由。
"为什么 \(k=2\) 看着比正文清爽?"
因为 \(k=2\) 时高斯和是精确等式 \(|S|=\sqrt q\)、降次一次到位。正文为一般 \(k\) 写的 Weyl/Hua 不等式要把这两件事都换成"放缩版"(\(\le q^{1-1/k+\varepsilon}\)、\(k-1\) 次差分),所以才层层 \(\le\)。本质是同一台机器,只是 \(k=2\) 的零件透明。
一句话总结:圆法把"\(N\) 是不是五个平方和"翻译成"信号(主弧,尖峰,可精确计算)能不能压过噪声(次弧,杂乱,只需估小)"的一场赛跑。\(k=2\) 让每个零件都现出原形——高斯和恰是 \(\sqrt q\)、降次一步到位——于是你能亲眼看到:每一个 \(\le\) 都丢得恰到好处,而 \(s=5\) 正是让信号刚好领先 \(P^{1/4}\) 的最小变量数。
带着这张全景图,再回去读正文第 3–6 章那些一般 \(k\) 的不等式,它们就不再"难看"——只是同一招在看不见具体数字时的样子。