佩尔方程

w20810

于 2015-02-05 14:13:53 发布

阅读量1.3k

点赞数

分类专栏： ACM-数学文章标签：佩尔方程

本文链接：https://blog.csdn.net/w20810/article/details/43527357

版权

73 篇文章 0 订阅

订阅专栏

参见：维基百科

定义：若一个不定方程具有这样的形式： $http://upload.wikimedia.org/math/3/6/9/369132e4b800099cf4c2a41e4a954c55.png$ 则称此二元二次不定方程为佩尔方程

若n是完全平方数，则这个方程式只有平凡解 $(\pm 1, 0)$ 。

解佩尔方程：若佩尔方程 $http://upload.wikimedia.org/math/3/6/9/369132e4b800099cf4c2a41e4a954c55.png$ 的最小特解（最小正整数解）是（x1,y1），那么可有迭代公式

$\displaystyle x_{i+1} = x_1 x_i + n y_1 y_i,$

$\displaystyle y_{i+1} = x_1 y_i + y_1 x_i$ 。

求出所有正整数解（xk,xk）,用矩阵表示如下

已知最小特解可以用矩阵快速幂求出任何一项解。

暴力法求特解代码：

void search(int d,int &x,int &y) //x^2 - d*(y^2) = 1
{
	y=1;
	while(1)
	{
		x=(long long )sqrt(d*y*y+1);
		if(x*x-d*y*y==1)
			break;
		y++;
	}
}

题目1：http://poj.org/problem?id=1320
题目2：acm.hdu.edu.cn/showproblem.php?pid=3292

关注

专栏目录