【算法竞赛学习笔记】佩尔方程-数学提升计划

本文探讨了佩尔方程的性质及其解法,包括如何使用矩阵快速幂迭代法解决佩尔方程,并通过实例HDU3292介绍了第二类佩尔方程的解题思路。通过迭代公式和矩阵运算,可以高效地求解佩尔方程的通解,适合ACM竞赛训练和数学爱好者学习。
摘要由CSDN通过智能技术生成

title : 佩尔方程
date : 2021-10-31
tags : ACM,数学
author : Linno

在这里插入图片描述

佩尔方程

形如 x 2 − d y 2 = 1 ( d > 1 且 d 不 为 完 全 平 方 数 ) x^2-dy^2=1(d>1且d不为完全平方数) x2dy2=1(d>1d)的不定方程称为佩尔方程
设 两 个 解 为 ( x 1 , y 1 ) 和 ( x 2 , y 2 ) , 那 么 有 { x 1 2 − d y 1 2 = 1 x 2 2 − d y 2 2 = 1 得 到 ( x 1 2 − d y 1 2 ) ( x 2 2 − d y 2 2 ) = 1 展 开 得 到 x 1 2 x 2 2 − d x 1 2 y 2 2 − d y 1 2 x 2 2 + d 2 y 1 2 y 2 2 = 1 ⇒ ( x 1 2 x 2 2 + d 2 y 1 2 y 2 2 ) − d ( d x 1 2 y 2 2 + y 1 2 x 2 2 ) = 1 ⇒ ( x 1 2 x 2 2 + d 2 y 1 2 y 2 2 + 2 d x 1 x 2 y 1 y 2 ) − d ( d x 1 2 y 2 2 + y 1 2 x 2 2 − 2 d x 1 x 2 y 1 y 2 ) = 1 ⇒ ( x 1 x 2 + d y 1 y 2 ) − d ( d x 1 2 y 2 2 + y 1 2 x 2 2 − 2 d x 1 x 2 y 1 y 2 ) = 1 所 以 x 3 = x 2 x 1 + d y 2 y 1 , y 3 = x 2 y 1 + x 1 y 2 根 据 迭 代 公 式 , 通 解 为 : x n = x n − 1 x 1 + d y n − 1 y 1 , y n = x n − 1 y 1 + y n − 1 x 1 设两个解为(x_1,y_1)和(x_2,y_2),那么有\\ \begin{cases} x_1^2-dy_1^2=1\\ x_2^2-dy_2^2=1\\ \end{cases} 得到(x_1^2-dy_1^2)(x_2^2-dy_2^2)=1\\ 展开得到x_1^2x_2^2-dx_1^2y_2^2-dy_1^2x_2^2+d^2y_1^2y_2^2=1\\ \Rightarrow (x_1^2x_2^2+d^2y_1^2y_2^2)-d(dx_1^2y_2^2+y_1^2x_2^2)=1\\ \Rightarrow (x_1^2x_2^2+d^2y_1^2y_2^2+2dx_1x_2y_1y_2)-d(dx_1^2y_2^2+y_1^2x_2^2-2dx_1x_2y_1y_2)=1\\ \Rightarrow (x_1x_2+dy_1y_2)-d(dx_1^2y_2^2+y_1^2x_2^2-2dx_1x_2y_1y_2)=1\\ 所以x_3=x_2x_1+dy_2y_1,y_3=x_2y_1+x_1y_2\\ 根据迭代公式,通解为:\\ x_n=x_{n-1}x_1+dy_{n-1}y_1,y_n=x_{n-1}y_1+y_{n-1}x_1\\ (x1,y1)(x2,y2){ x12dy12=1x22dy22=1

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

RWLinno

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值