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不为完全平方数) x2−dy2=1(d>1且d不为完全平方数)的不定方程称为佩尔方程
设 两 个 解 为 ( 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),那么有{
x12−dy12=1x22−dy22=1得到