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𝑄𝑘+1pk2−Dqk2=(−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=x1x2+y1y2D 𝑌=𝑥1𝑦2+𝑥2𝑦1Y=x1y2+x2y1 也是方程的整数解。
具体例子与应用
例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 竞赛中解决相关问题。