佩尔方程

最新推荐文章于 2020-02-11 19:09:29 发布

theArcticOcean

最新推荐文章于 2020-02-11 19:09:29 发布

阅读量1k

点赞数

分类专栏： algorithm_数论文章标签：数学

本文链接：https://blog.csdn.net/thearcticocean/article/details/50678650

版权

本文探讨了佩尔方程的性质，重点在于当d不是完全平方数且d>1时，如何寻找其最小特解x1, y1。通过对佩尔方程的深入分析，阐述了解决这类方程的方法及其在数论中的重要应用。" 82025759,1009978,深度学习在人脸检测与对齐的应用,"['深度学习', '计算机视觉', '人脸识别', '图像处理']

摘要由CSDN通过智能技术生成

佩尔方程(Pell Equation)为： $x^2-dy^2=1$ 其中d不为完全平方数且d>1.
如果已知它的最小特解：x1,y1

那么存在迭代公式：

$\\ x_n=x_{n-1}x_1+dy_{n-1}y_1\\ y_n=x_{n-1}y_1+y_{n-1}x_1$
通过简单的证明：
$\\ x_n^2-dy_n^2 \\=x_{n-1}^2x_1^2+d^2y_{n-1}^2y_1^2+2x_{n-1}y_{n-1}x_1y_1d-d(x_{n-1}^2y_1^2+y_{n-1}^2x_1^2+2x_{n-1}y_{n-1}x_1y_1)\\=x_{n-1}^2(x_1^2-dy_1^2)+dy_{n-1}^2(dy_1^2-x_1^2)+2x_{n-1}y_{n-1}x_1y_1d(1-1)\\=x_{n-1}^2-dy_{n-1}^2=1$
由此得到矩阵递推式：
$\\ \begin{bmatrix} x_k \\ y_k \end{bmatrix}=\begin{bmatrix} x_1 & dy_1\\ y_1 & x_1 \end{bmatrix}^{k-1}\begin{bmatrix} x_1 \\ y_1 \end{bmatrix}$

暴力法寻找最小特解：

typedef long long LL;
void search(LL &x,LL &y,LL d){
    y=1;
    while(1>0){
        x=(LL)sqrt(1+d*y*y);
        if(x*x-d*y*y==1) break;
        y++;
    }
}

hdu 3292 No more tricks, Mr Nanguo
http://acm.hdu.edu.cn/showproblem.php?pid=3292

大意：转化题意，即x^2-NY^2=1 求出第K个X解

分析：先暴力法求出特解，然后，Pell方程的解递推矩阵求出答案

#include <iostream>
#include <cstdio>
#in

最低0.47元/天解锁文章

theArcticOcean

关注

0
点赞
踩
1

收藏

觉得还不错? 一键收藏
0
评论
复制链接

分享到 QQ

分享到新浪微博

扫一扫

专栏目录