佩尔方程求解器

佩尔方程 x² − Dy² = 1（D 为非平方正整数）是数论中最富传奇的丢番图方程：总有无穷多组解，且全部解由最小解的幂生成。D=2 的最小解 (3,2)：3²−2×2² = 1。

最小解可能大得惊人：D=61 时最小解 x 已达 1766319049——费马 1657 年用它「刁难」英国数学家。求解标准方法是 √D 的连分数展开：周期渐近分数中藏着最小解。

佩尔方程 x²−Dy²=1（D 为正的非平方整数）是不定方程中最富戏剧性的一个：它要求的其实是 √D 的「最佳有理逼近」——x/y 必须逼近 √D 到误差小于 1/(2y²)，恰好落在连分数渐近分数的判定区间内。这就是求解原理：展开 √D 的连分数（其周期部分必然循环），周期为 r 时，第 r−1 个（r 偶）或第 2r−1 个（r 奇）渐近分数给出最小解；一旦有了最小解 (x₁,y₁)，全部解由幂 (x₁+y₁√D)ⁿ 生成。

「佩尔」之名是数学史上最著名的张冠李戴：欧拉误把布龙克尔（Brouncker）的解法记在佩尔名下，佩尔本人与这个方程几乎无关。更早的源流在印度——7 世纪婆罗摩笈多给出组合两个解的「婆罗门恒等式」，12 世纪婆什迦罗的「循环法」（chakravala）已能系统求解 D=61 这样的硬骨头，比欧洲早 500 年。而人类历史上最疯狂的佩尔方程出自阿基米德牛群问题：归约后是 x²−410286423278424y²=1，最小解的 x 超过 20 万位十进制数字——1880 年才有人算出前几位，1965 年计算机才给出完整解。

佩尔方程的深层意义在代数数论：x+y√D 恰是实二次域 ℚ(√D) 中范数为 1 的代数整数，最小解对应域的「基本单位」——狄利克雷单位定理最简单的实例。解的爆炸式增长（D=61 的解为何 17 亿？）由调节子（regulator）控制，它与类数、L 函数纠缠在 BSD 猜想的邻近地带。连分数周期长度、最小解大小、类数之间的关系至今活跃在研究前沿。