http://blog.sina.com.cn/s/blog_5d06e2390100ll92.html
佩尔方程实际上并不是佩尔提出的,而是费尔马提出,却被欧拉误记为佩尔提出,因此佩尔方程的名称沿用至今。身为不定方程的特殊一类,佩尔方程与连分数,二次型,代数论等等有着重要的联系,因而是数论中最经典的篇章之一。令 d 为非平方数的正整数,那么佩尔方程(Pell Equation)为:
![数论笔记 <wbr>路 <wbr>佩尔方程 数论笔记 <wbr>路 <wbr>佩尔方程](https://i-blog.csdnimg.cn/blog_migrate/85523b23fb9a8b1b54ef11556fbf8beb.jpeg)
* 连分数 *
在对佩尔方程进行深入了解之前,先来看看与该方程不可分离的连分数:
![数论笔记 <wbr>路 <wbr>佩尔方程 数论笔记 <wbr>路 <wbr>佩尔方程](https://i-blog.csdnimg.cn/blog_migrate/9ff527055198ca59ff2e50ea1e623b68.jpeg)
由于篇幅的原因,通常记为
![数论笔记 <wbr>路 <wbr>佩尔方程 数论笔记 <wbr>路 <wbr>佩尔方程](https://i-blog.csdnimg.cn/blog_migrate/a79eb1651b018d59ae61d36fc53f405d.jpeg)
其中 pn 和 qn 称为连分数之多项式,对于任意的a 均为一次式,它们的比值称为第 n 个渐进值渐进分数。对于渐进值来讲,有着递归关系式:
![数论笔记 <wbr>路 <wbr>佩尔方程 数论笔记 <wbr>路 <wbr>佩尔方程](https://i-blog.csdnimg.cn/blog_migrate/d50f1a9fda20bc9086cd9b90e67d086c.jpeg)
由数学归纳法可以得到关系式:
在华罗庚的《数论导引》中介绍了连分数一些性质定理,例如:有理数必可表示为有限连分数,无理数的连分数表示法唯一等等。
* 佩尔方程 *
注意到在上述多项式的递归表示中,具有决定多项式数值作用的关键的 an 并没有计算公式。此时回来看佩尔方程,根据数论中的定理:
定理1:对于正整数p、q ,如果有
![数论笔记 <wbr>路 <wbr>佩尔方程 数论笔记 <wbr>路 <wbr>佩尔方程](https://i-blog.csdnimg.cn/blog_migrate/ee046efcb45440aac8109e8f7cd05bb0.jpeg)
则比值 p/q 必为 a 的一个渐进值。
因此可以得出推论:对于佩尔方程,其全部的根的集合为:
![数论笔记 <wbr>路 <wbr>佩尔方程 数论笔记 <wbr>路 <wbr>佩尔方程](https://i-blog.csdnimg.cn/blog_migrate/9d02c1d79724f037ed29a94f5b9c7139.jpeg)
由递归方程式
![数论笔记 <wbr>路 <wbr>佩尔方程 数论笔记 <wbr>路 <wbr>佩尔方程](https://i-blog.csdnimg.cn/blog_migrate/945082f8590fa2862e5766c22f2912c7.jpeg)
可以得到
![数论笔记 <wbr>路 <wbr>佩尔方程 数论笔记 <wbr>路 <wbr>佩尔方程](https://i-blog.csdnimg.cn/blog_migrate/a9bed506315beedfbd80978cf10d8953.jpeg)
而判断二次不定方程有解的充分必要条件有下述定理表述:
定理2:二次不定方程
![数论笔记 <wbr>路 <wbr>佩尔方程 数论笔记 <wbr>路 <wbr>佩尔方程](https://i-blog.csdnimg.cn/blog_migrate/cafca4506b5174a0ea770a22a084562b.jpeg)
常有解,二次方程
![数论笔记 <wbr>路 <wbr>佩尔方程 数论笔记 <wbr>路 <wbr>佩尔方程](https://i-blog.csdnimg.cn/blog_migrate/32e3c36dd8f1683bfc85e421070f85cf.jpeg)
不可解的充要条件为:
![数论笔记 <wbr>路 <wbr>佩尔方程 数论笔记 <wbr>路 <wbr>佩尔方程](https://i-blog.csdnimg.cn/blog_migrate/f92f4cfe873537f3242ebb8cd94aa00b.jpeg)
根据上述定理,由数学归纳法可以得到关系式:
![数论笔记 <wbr>路 <wbr>佩尔方程 数论笔记 <wbr>路 <wbr>佩尔方程](https://i-blog.csdnimg.cn/blog_migrate/a521291527906f4672c09ab6e8abc3ed.jpeg)
![数论笔记 <wbr>路 <wbr>佩尔方程 数论笔记 <wbr>路 <wbr>佩尔方程](http://s14.sinaimg.cn/middle/5d06e239g8f5b4e6f7d3d&690)
![数论笔记 <wbr>路 <wbr>佩尔方程 数论笔记 <wbr>路 <wbr>佩尔方程](https://i-blog.csdnimg.cn/blog_migrate/9719c67ff8e0fa8232fe2d2c679cb86b.jpeg)
这样就可以根据 p、q、P、Q 的递归关系式计算出根式形式的无理数的连分数展开。
* 循环连分数 *
循环连分数指的是连分数表达式中存在循环节。根据数论中的基本定理:一个连分数为循环连分数的充分必要条件为此数是一有理系数二次不可化方程的根。因此根式形式的无理数均为循环连分数,且具有独特的性质:
通过循环表示,可以求解得到佩尔方程的最小整数解:
![数论笔记 <wbr>路 <wbr>佩尔方程 数论笔记 <wbr>路 <wbr>佩尔方程](https://i-blog.csdnimg.cn/blog_migrate/71410dab4923d4127205175b3f3946e3.jpeg)
* 佩尔型方程 *
更加复杂的佩尔形式的方程:
假设 p、q 为上述方程的解,考虑其对应的佩尔方程极其基础解系:
则有
![数论笔记 <wbr>路 <wbr>佩尔方程 数论笔记 <wbr>路 <wbr>佩尔方程](https://i-blog.csdnimg.cn/blog_migrate/4dc0bc83b13b87466ec36762ed0bf5f7.jpeg)
当 (p,q) 为包含所有基础解,这样所有的解集由下式递归构成
![数论笔记 <wbr>路 <wbr>佩尔方程 数论笔记 <wbr>路 <wbr>佩尔方程](https://i-blog.csdnimg.cn/blog_migrate/72c1aa4ebf0de89f9f7e07d9449a22a5.jpeg)
计算佩尔型方程的互不等价的基本解集主要依赖 PQA 算法和 Lagrange-Matthews-Mollin 算法。
PQA算法:
对佩尔方程做表
![数论笔记 <wbr>路 <wbr>佩尔方程 数论笔记 <wbr>路 <wbr>佩尔方程](https://i-blog.csdnimg.cn/blog_migrate/750cb5bb2020c0bcb138743808629dbf.jpeg)
各项递归按照佩尔方程递归关系计算得到。其中各项初始值为: a 的初始值 a0 相当于佩尔方程中的 a1 ,p0 = a0 * Q0 - P0 ,q0 = 1 ,p1 = a1 * p0 + Q0 ,q1 = a1 。表的递归计算到 Q 的第一个周期或者 |Q|=1 。值得注意的是,此时 a 的递归方程式中,a0 不能先取整,因为若K是负数,则会导致错误的计算。
这样就可以进行佩尔型方程基本解集的计算:按照 K^2 和 D 的关系有:
CASE K^2<D
P0 = 0 ,Q0 = 1 作为初始值带入 PQA 算法,当 Kn = K 时得到所有基本解集。
CASE K^2>D : LMM Algorithm
When f>0 , list m = K / f^2
For each m , find all z :-|m|/2 < z <= |m|/2,and z^2 = D modulo |m|
For each z , P0 = z ,Q0 = |m| as initial value in PQA algorithm:
If Qn+1 * sign(m) = 1 , then f*(Pn,Qn) is solution.
If Qn+1 * sign(m) = -1 , then f*(Pn,Qn) is solution of x^2 - Dy^2 = -K .
If x^2 - Dy^2 = -1 has primitive solution (u,v) , then f*(uPn+DvQn,vPn+uQn) is solution.
根据递归,该基本解集的每一个元素都可以产生无穷的解集,综合起来就是佩尔型方程的解集。