【数学】不定方程

b不定方程的相关内容,包括n元一次不定方程、毕达哥拉斯定理、费马大定理、佩尔方程。

一次不定方程

转化为一元线性同余方程

nn 元一次不定方程

转化为多元线性同余方程

毕达哥拉斯定理

x2+y2=z2x^2+y^2=z^2 ,当gcd(x,y,z)=1\gcd(x,y,z)=1 时被称为本原的毕达哥拉斯三元组。

本原的毕达哥拉斯三元组(x,y,z)y为偶数m,n(m>n,m,n互素,不同奇偶),x=m2n2,y=2mn,z=m2+n2本原的毕达哥拉斯三元组(x,y,z)且y为偶数\Leftrightarrow\exist m,n(m>n,m,n互素,不同奇偶),x=m^2-n^2,y=2mn,z=m^2+n^2

费马大定理

xn+yn=zn,n3,nNx^n+y^n=z^n,n\geq 3,n\in N 无非00 整数解

佩尔方程

第一类佩尔方程

形如:x2dy2=1,d>1x^2-dy^2=1,d>1

dd 是完全平方数\Rightarrow 无解

解有迭代公式:

xn=xn1x1+dyn1y1yn=xn1y1+yn1x1x_{n}=x_{n-1} x_{1}+d y_{n-1} y_{1}\\ y_{n}=x_{n-1} y_{1}+y_{n-1} x_{1}

推导:设特解(x1,y1),(x2,y2)(x_1,y_1),(x_2,y_2) ,则有x12dy12=1,x22dy22=1x_{1}^{2}-d y_{1}^{2}=1,x_{2}^{2}-d y_{2}^{2}=1(x12dy12)(x22dy22)=1\left(x_{1}^{2}-d y_{1}^{2}\right)\left(x_{2}^{2}-d y_{2}^{2}\right)=1
展开,有

x12x22dx12y22dy12x22+d2y12y22=(x12x22+d2y12y22)d(x12y22+y12x22)=(x1x2+dy1y2)2d(x1y2+x2y1)2=1\begin{aligned} &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}\\=&(x_{1}^{2} x_{2}^{2}+d^{2} y_{1}^{2} y_{2}^{2})-d(x_{1}^{2} y_{2}^{2}+y_{1}^{2} x_{2}^{2})\\ =&\left(x_{1} x_{2}+d y_{1} y_{2}\right)^{2}-d\left(x_{1} y_{2}+x_{2} y_{1}\right)^{2}\\ =&1 \end{aligned}

逐次迭代可得上述迭代公式。

暴力迭代法

y=1y=1 开始枚举验证,每次+1+1

矩阵迭代法

kk 个迭代解用矩阵表示如下:

[xkyk]=[x1dy1y1x1]k1[x1y1]\left[\begin{array}{l}x_{k} \\ y_{k}\end{array}\right]=\left[\begin{array}{ll}x_{1} & d y_{1} \\ y_{1} & x_{1}\end{array}\right]^{k-1}\left[\begin{array}{l}x_{1} \\ y_{1}\end{array}\right]

求出第一个特解后用矩阵快速幂求得第kk 个解。

第二类佩尔方程

形如:x2dy2=k,d>1x^2-dy^2=k,d>1

解有迭代公式:

x=px1+dqy1y=py1+qx1x=p x_{1}+d q y_{1}\\y=p y_{1}+q x_{1}

其中(p,q)(p,q) 是第二类佩尔方程的一个特解,(x1,y1)(x_1,y_1) 是第一类佩尔方程的最小特解。

推导:根据上述,有p2dq2=k,x12dy12=1p^2-dq^2=k,x_{1}^{2}-d y_{1}^{2}=1 ,则(x12dy12)(p2dy2)=k\left(x_{1}^{2}-d y_{1}^{2}\right)\left(p^{2}-d y^{2}\right)=k

展开,有

x12p2dx12q2dy12p2+d2y12q2=(x12p2+d2y12q2)d(x12q2+y12p2)=(x1p+dy1q)2d(x1q+py1)2=k\begin{aligned} &x_{1}^{2} p^{2}-d x_{1}^{2} q^{2}-d y_{1}^{2} p^{2}+d^{2} y_{1}^{2} q^{2}\\=&(x_{1}^{2} p^{2}+d^{2} y_{1}^{2} q^{2})-d(x_{1}^{2} q^{2}+y_{1}^{2} p^{2})\\ =&\left(x_{1} p+d y_{1} q\right)^{2}-d\left(x_{1} q+p y_{1}\right)^{2}\\ =&k \end{aligned}

逐次迭代可得上述迭代公式。

对于每一组特解(p,q)(p,q) ,第kk 个迭代解用矩阵表示如下:

[xkyk]=[pdqqp]k1[x1y1]\left[\begin{array}{l}x_{k} \\ y_{k}\end{array}\right]=\left[\begin{array}{ll}p & d q \\ q& p\end{array}\right]^{k-1}\left[\begin{array}{l}x_{1} \\ y_{1}\end{array}\right]

求出第一个特解后用矩阵快速幂求得第kk 个解。