感谢这两位巨佬,部分思路和图片来源于他们的启发:
是什么?
Pell \text{Pell} Pell方程主要用来解决下面形式的方程:
x 2 − d y 2 = 1 x^2-dy^2=1 x2−dy2=1其中 d d d是正数且不是完全平方数,若 x , y x,y x,y都为整数的解称为这个方程的解,其中 x + y d x+y\sqrt d x+yd最小的那组解就称为 Pell \text{Pell} Pell方程的基本解,我们要解决两个问题:求出基本解和用它扩展出所有解。
基本解?
本文结论大多不给出证明,但是会讲清楚作法。
首先要知道连分数表示法,是将一个实数用以下方法表示:
连分数的求法就是分离最多的正整数,使剩下的为非负,可以把连分数简记为 x = [ a 0 ; a 1 , a 2 , a 3 . . . . ] x=[a_0;a_1,a_2,a_3....] x=[a0;a1,a2,a3....]
结论1:有理数的连分数是有限的,无理数的连分数是无限的,所以可以用这种方法逼近无理数的值。
举个例子,比如 22 \sqrt{22} 22的连分数:
简写为: 22 = [ 4 ; 1 , 2 , 4 , 2 , 1 , 8 ] \sqrt{22}=[4;1,2,4,2,1,8] 22=[4;1,2,4,2,1,8],其中后面 6 6 6项是循环,也可以表示为 22 = [ 4 ; < 1 , 2 , 4 , 2 , 1 , 8 > ] \sqrt{22}=[4;<1,2,4,2,1,8>] 22=[4;<1,2,4,2,1,8>]
结论2:循环的最后一项是 a 0 a_0 a0的两倍,无理数的连分数都是成周期性的。
我们把 d \sqrt d d写成连分数的形式: d = [ a 0 ; a 1 , a 2 , a 3 . . . . . , 2 a 0 ] \sqrt d=[a_0;a_1,a_2,a_3.....,2a_0] d=