算术与数论
佩尔方程
方程 x2−dy2=1,其整数解与连分数相关联。
直观直觉:双曲线上逼近 d 的整数点
寻找满足 x2−2y2=1 的整数 x,y。尝试较小的值:(x,y)=(3,2) 成立,因为 9−8=1。整理得 (yx)2=2+y21,于是 x/y=3/2=1.5 已经是对 2≈1.41421 惊人精确的近似。形如 x2−dy2=1 的方程——称为佩尔方程——与连分数及最佳有理逼近理论密不可分。
将 d 的连分数渐近分数画成路径图:每个节点是一个渐近分数 pk/qk,高亮节点为基本解。进阶定义、基本解与结构
定义: 佩尔方程与基本解
对不是完全平方数的正整数 d,x2−dy2=1 称为佩尔方程。在其正整数解 (x,y) 中,x 最小(等价地 y 最小)的那个称为基本解 (x1,y1)。
x2−dy2=1,d∈Z>0 not a perfect square 将左边分解为 (x−yd)(x+yd)=1,佩尔方程就变成关于环 Z[d] 的命题:解对应于范数为 1 的元素,两个这样的元素相乘仍得到一个。特别地,一旦知道 (x1,y1),其余所有解都可通过取幂得到:xn+ynd=(x1+y1d)n。
(x1+y1d)(x1−y1d)=1 ⟹ xn+ynd=(x1+y1d)n 较小 d 的基本解| d | d 的连分数 | 基本解 (x1,y1) |
|---|
| 2 | [1;2] | (3,2) |
| 3 | [1;1,2] | (2,1) |
| 5 | [2;4] | (9,4) |
| 7 | [2;1,1,1,4] | (8,3) |
进阶定理与证明
对每个不是完全平方数的正整数 d,x2−dy2=1 都存在正整数解 x,y。
为什么成立?
双曲线 x2−dy2=1 显然有实数点,但它必经过某个格点这一点绝非显然——该定理通过对有理逼近的巧妙鸽笼论证,保证对每个非完全平方数 d 都必然如此。
证明
由狄利克雷逼近定理,对任意整数 Q>0,存在整数 p,q,1≤q≤Q,且 d−qp<qQ1。令 Q→∞ 得到无穷多对 (p,q) 满足 d−qp<q21。
对每一对这样的数,∣p−qd∣<q1,故 ∣p2−dq2∣=∣p−qd∣⋅∣p+qd∣<q1(2qd+q1)<2d+1。因此 p2−dq2 只能取区间 (−2d−1,2d+1) 内有限多个整数值,而 (p,q) 却有无穷多对。
由鸽笼原理,该有限范围内某个固定的非零整数 k 使得 p2−dq2=k 对无穷多对 (p,q) 成立。在这无穷多对中,再由鸽笼原理,有无穷多对满足相同的剩余类 p≡p0, q≡q0(mod∣k∣)。
取两对这样不同的 (p1,q1)=(p2,q2),满足 p12−dq12=p22−dq22=k 且模 ∣k∣ 剩余相同。设 x+yd=k(p1+q1d)(p2−q2d);展开可知 x=kp1p2−dq1q2 与 y=kp1q2−p2q1 正因模 ∣k∣ 剩余相同而确为整数,而范数 N(a+bd)=a2−db2 的积性给出 x2−dy2=k2k⋅k=1。由于 (p1,q1)=(p2,q2) 但极限下比值相同,可验证 y=0,将 (x,y) 换成 (∣x∣,∣y∣)(仍是解,因为只出现平方)即得正整数解。
若 (x1,y1) 是 x2−dy2=1 的基本解,则每个正整数解 (x,y) 都等于某个 n≥1 对应的 (xn,yn),其中 xn+ynd=(x1+y1d)n。
为什么成立?
它说明无穷多解并非神秘散乱的集合,而是通过反复“乘以”最小解生成的完全可预测的类几何序列——把无限的搜索归结为只需找到一个数。
证明
首先验证由 xn+ynd=(x1+y1d)n 定义的 (xn,yn) 对每个 n 确实是解:取共轭,xn−ynd=(x1−y1d)n,故 xn2−dyn2=(xn+ynd)(xn−ynd)=[(x1+y1d)(x1−y1d)]n=(x12−dy12)n=1n=1。
现设 (x,y) 是任意不具有此形式的正整数解;由于当 n→∞ 时 xn→∞,存在唯一的 n 使 xn+ynd≤x+yd<xn+1+yn+1d=(xn+ynd)(x1+y1d)。
两边除以 (xn+ynd),即乘以其逆 (xn−ynd)(由 xn2−dyn2=1 保证有效):设 x′+y′d=(x+yd)(xn−ynd)。则 1≤x′+y′d<x1+y1d,且 x′2−dy′2=(x2−dy2)(xn2−dyn2)=1⋅1=1,故 (x′,y′) 也是佩尔方程的解。
利用 x′+y′d≥1 与 x′2−dy′2=1 做简短计算可得 x′≥1 且 y′≥0(若某解满足 x′+y′d≥1 但 y′<0,则会迫使 x′>x1,与 x′+y′d<x1+y1d 连同范数方程矛盾)。若 y′>0,则 (x′,y′) 是满足 x′+y′d<x1+y1d 的正解,与基本解的最小性矛盾。故 y′=0,从而 x′=1,即 x+yd=xn+ynd,故 (x,y)=(xn,yn)——与假设矛盾。
因此每个正解恰好是某个 (xn,yn),证明了基本解生成了整个解集。
进阶实际应用与典型例题
佩尔方程远不止是一个谜题:它支配着机械齿轮设计中所用的最佳有理逼近,其数论支撑着经典的整数分解算法,而对大 d 计算其基本解正是 Hallgren 的量子算法能在多项式时间内求解的计算问题——目前尚无已知的高效经典方法。
例题: 求 d=2 的基本解
利用 2 的连分数求 x2−2y2=1 的基本解,再生成下一个解。
解答
2 的连分数为 [1;2]=1+2+2+⋯11,其渐近分数为 1, 23, 57, 1217, 2941,…
检验渐近分数 23:32−2⋅22=9−8=1。这是满足条件的最小渐近分数,故基本解为 (x1,y1)=(3,2)。
用递推式 xn+ynd=(x1+y1d)n,取 n=2:x2+y22=(3+22)2=9+122+8=17+122,故 (x2,y2)=(17,12)。
验证:172−2⋅122=289−288=1,确认了下一个解,恰与上面求得的渐近分数 1217 相符——正如理论所预测。
例题: 工程逼近:设计接近 2 的齿轮比
某机构需要两个啮合齿轮,其齿数比用较小的整数齿数尽可能逼近 2,使得多次旋转后误差仍然很小。利用佩尔方程的渐近分数来选择齿数并估计误差。
解答
由 x2−2y2=1 的基本解 (3,2),比值 3/2=1.5 给出 ∣3/2−2∣≈0.0858——可用,但对精密机械而言过于粗糙。
利用由 3+22 平方得到的下一个渐近分数 (17,12):齿数为 17 和 12 的齿轮对给出比值 17/12≈1.41667,且 ∣17/12−2∣≈0.00245,仅略微增加齿数就比 3/2 方案精确了近 35 倍。
对连分数任意渐近分数 pk/qk 的一般界为 qkpk−2<qk21,因此当工程师转向下一个佩尔解 (x3,y3)(由 (3+22)3=99+702 即 99/70 得到)时,齿数增至 70,但误差缩小到 1/702≈0.0002 以下。
这体现了佩尔方程中明确的工程权衡:每个后续解 (xn,yn) 都是用制造成本的增加(更多齿数,即更大的 yn)换取可量化、快速缩小的逼近误差,让设计者能在此曲线上精确选取符合预算与精度要求的点。
x2−2y2=1 的基本解是什么?
若 (x1,y1) 是基本解,解 (x2,y2) 如何得到?
x2−3y2=−1 有整数解吗?
Hallgren 的算法用哪种计算在多项式时间内求解佩尔方程?