8.12.31 ACM-ICPC数学 数论 Pell 方程

8.12.31 ACM-ICPC数学 数论 Pell 方程

引言

Pell 方程是一类重要的二次不定方程，其形式为 𝑥2−𝐷𝑦2=1x2−Dy2=1 或 𝑥2−𝐷𝑦2=−1x2−Dy2=−1。这种方程在数论中有着广泛的应用，特别是在研究二次整环和单位数时。本节将详细介绍 Pell 方程的基本概念、性质以及求解方法，并讨论其在 ACM-ICPC 竞赛中的应用。

Pell 方程的基本概念

定义

Pell 方程是形如 𝑥2−𝐷𝑦2=1x2−Dy2=1 或 𝑥2−𝐷𝑦2=−1x2−Dy2=−1 的方程，其中 𝐷D 是一个不含平方因子的正整数。

二次整数

对于二次有理数 𝑎+𝑏𝐷a+bD​，如果 𝑎a 和 𝑏b 满足以下条件之一：

• 𝑎a 与 𝑏b 是整数，且 𝐷≡2(mod4)D≡2(mod4) 或 𝐷≡3(mod4)D≡3(mod4)。
• 𝑎a 与 𝑏b 同时是半整数，且 𝐷≡1(mod4)D≡1(mod4)。

那么 𝑎+𝑏𝐷a+bD​ 就是二次整数。二次整数与首一整系数二次方程的解构成对应关系。

单位数与基本单位数

单位数

如果二次整数 𝑎+𝑏𝐷a+bD​ 的范数 𝑎2−𝐷𝑏2a2−Db2 为 1 或 -1，则它的倒数也是二次整数，恰好是它的共轭或者共轭的相反数。此时称它为整环 𝑍(𝐷)Z(D​) 的单位数，简称单位数。

基本单位数

可以证明，存在基本单位数，使得全体单位数都可以表示成为基本单位数的幂（或幂的相反数）。它也就是对应 Pell 方程的基本解，通解可以表示为基本解的幂（或幂的相反数）。

求解 Pell 方程

Dirichlet 逼近定理

利用 Dirichlet 逼近定理，可以逼近二次根式 𝐷D​，即有无穷个有理数（显然为正有理数）满足： ∣𝑥𝑦−𝐷∣≤1𝑦2∣∣​yx​−D​∣∣​≤y21​

范数估值

根据上面的逼近关系，可以对范数进行估值： ∣𝑁(𝑥+𝑦𝐷)∣=∣𝑥−𝑦𝐷∣∣𝑥+𝑦𝐷∣≤1𝑦(1𝑦+2𝑦𝐷)≤2𝐷+1∣∣​N(x+yD​)∣∣​=∣x−yD​∣∣x+yD​∣≤y1​(y1​+2yD​)≤2D​+1 这说明只要有理数与 𝐷D​ 越接近，范数就越小。

基本解的存在

根据上述估值，可以证明范数为 ±1±1 的单位数存在，并且存在无限个。

利用渐进分数求解

对于所有 𝐷D​ 的渐进分数，配上系数之后得到的二次整数的范数都落在非常小的区间。通过计算每个循环节处前一个渐进分数，可以找到使得范数为 ±1±1 的二次整数。

Pell 方程的求解方法

定理：范数为 ∣𝑠∣<𝐷∣s∣<D​ 的渐进分数

若 𝑥2−𝐷𝑦2=𝑠x2−Dy2=s 且 ∣𝑠∣<𝐷∣s∣<D​，则 𝑥𝑦yx​ 一定是 𝐷D​ 的渐进分数。

定理：渐进分数与 Q_k 的关系

𝑝𝑘2−𝐷𝑞𝑘2=(−1)𝑘+1𝑄𝑘+1pk​2−Dqk​2=(−1)k+1Qk+1​

定理：Q_k 的周期性

当且仅当 𝑘=𝑛𝐿k=nL（其中 𝑛n 是正整数，𝐿L 是最短循环周期）时，有 𝑄𝑘=1Qk​=1。

定理：解的组合

若 (𝑥1,𝑦1)(x1​,y1​) 和 (𝑥2,𝑦2)(x2​,y2​) 都是 𝑥2−𝐷𝑦2=1x2−Dy2=1 的整数解，则： 𝑋=𝑥1𝑥2+𝑦1𝑦2𝐷X=x1​x2​+y1​y2​D 𝑌=𝑥1𝑦2+𝑥2𝑦1Y=x1​y2​+x2​y1​ 也是方程的整数解。

具体例子与应用

例1：求解 Pell 方程 𝑥2−2𝑦2=1x2−2y2=1

表示 22​ 的循环连分数为 [1;2‾][1;2]。通过分析其循环部分，可以找到 Pell 方程的基本解 𝑥=3x=3， 𝑦=2y=2，通解为 (3+22)𝑛(3+22​)n 的形式。

例2：求解 Pell 方程 𝑥2−2𝑦2=−1x2−2y2=−1

当 𝐷D​ 的循环节长度为奇数时，方程 𝑥2−𝐷𝑦2=−1x2−Dy2=−1 有解。

结论

Pell 方程在数论中具有重要地位，其解法涉及连分数、渐进分数和单位数等概念。通过理解这些基本原理和定理，可以高效地求解 Pell 方程，并在 ACM-ICPC 竞赛中解决相关问题。