佩尔方程

讲解佩尔方程几个不错的博客地址：http://blog.csdn.net/acdreamers/article/details/8529686

http://hi.baidu.com/aekdycoin/item/a45f7c37850e5b9db80c03d1

http://hi.baidu.com/shouzhewei/item/5ff25ee1624c3419585dd832

ax^2 - by^2 = c       Pell 方程一般解法

以下涉及到定理的证明的部分全部略过.

开题自然少不了介绍,以上的公式就是Pell方程的一般形态.
显然如果告诉你a,b,c,一开始想到的只可能是暴力,可是接下来介绍的纯数学的方法可以很快速的求解几乎大部分解.
1.首先构造一个系数矩阵,显然为了构造这个矩阵，我们需要先得到下面方程的一个最小特解(x,y>0)

至于如何得到，可以使用暴力(当某些情况下暴力几乎求不到最小解)或者使用连分数的方法来求
假设我们得到了以上方程的特解: x0 y0 (x0,y0>0,并是最小的满足条件的解)
2.继续求

的一个最小特解.假设是x1,y1(x1,y1>0)
3.
假设你要求第k个解,那么有

例子:
1.求 x^2 - 3y^2 = 1的解
由于这里a=1,b=3,而c=1,所以我们可以知道x0=x1,y0=y1;
不难解得一个最小特解(2,1)，于是有

假设现在要知道第2个解
那么套用上面的公式得到
x2=7
y2=4
即
49-48=1
其他解类似.

2.求 x^2 - 3y^2 = 13的解
显然x^2 - 3y^2 =1的最小特解在上面已经求出来了
即
x0=2;
y0=1;
现在我们需要知道的是x^2 - 3y^2 = 13的最小特解，显然应该是
(4,1)
于是如果继续套用上面的解，可以得到:

那么得到
(x1,y1) = (4,1)
(x2,y2) = (11,6)

在这儿，特殊的佩尔方程x^2-n*y^2=1的求解