文章目录

  • 不定方程
    • 二元一次不定方程
    • nnn 元一次不定方程
    • 毕达哥拉斯定理
    • 费马大定理
    • 佩尔方程
      • 第一类佩尔方程
        • 暴力迭代法
        • 矩阵迭代法
      • 第二类佩尔方程

不定方程

二元一次不定方程

转化为一元线性同余方程

nnn 元一次不定方程

转化为多元线性同余方程

毕达哥拉斯定理

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

本原的毕达哥拉斯三元组(x,y,z)且y为偶数⇔∃m,n(m>n,m,n互素,不同奇偶),x=m2−n2,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本原的毕达哥拉斯三元组(x,y,z)且y为偶数⇔∃m,n(m>n,m,n互素,不同奇偶),x=m2−n2,y=2mn,z=m2+n2

费马大定理

xn+yn=zn,n≥3,n∈Nx^n+y^n=z^n,n\geq 3,n\in Nxn+yn=zn,n≥3,n∈N 无非000 整数解

佩尔方程

第一类佩尔方程

形如:x2−dy2=1,d>1x^2-dy^2=1,d>1x2−dy2=1,d>1

ddd 是完全平方数⇒\Rightarrow⇒ 无解

解有迭代公式:
xn=xn−1x1+dyn−1y1yn=xn−1y1+yn−1x1x_{n}=x_{n-1} x_{1}+d y_{n-1} y_{1}\\ y_{n}=x_{n-1} y_{1}+y_{n-1} x_{1} xn​=xn−1​x1​+dyn−1​y1​yn​=xn−1​y1​+yn−1​x1​
推导:设特解(x1,y1),(x2,y2)(x_1,y_1),(x_2,y_2)(x1​,y1​),(x2​,y2​) ,则有x12−dy12=1,x22−dy22=1x_{1}^{2}-d y_{1}^{2}=1,x_{2}^{2}-d y_{2}^{2}=1x12​−dy12​=1,x22​−dy22​=1 则(x12−dy12)(x22−dy22)=1\left(x_{1}^{2}-d y_{1}^{2}\right)\left(x_{2}^{2}-d y_{2}^{2}\right)=1(x12​−dy12​)(x22​−dy22​)=1
展开,有x12x22−dx12y22−dy12x22+d2y12y22=(x12x22+d2y12y22)−d(x12y22+y12x22)=(x1x2+dy1y2)2−d(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} ===​x12​x22​−dx12​y22​−dy12​x22​+d2y12​y22​(x12​x22​+d2y12​y22​)−d(x12​y22​+y12​x22​)(x1​x2​+dy1​y2​)2−d(x1​y2​+x2​y1​)21​逐次迭代可得上述迭代公式。

暴力迭代法

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

矩阵迭代法

第kkk 个迭代解用矩阵表示如下:
[xkyk]=[x1dy1y1x1]k−1[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] [xk​yk​​]=[x1​y1​​dy1​x1​​]k−1[x1​y1​​]
求出第一个特解后用矩阵快速幂求得第kkk 个解。

第二类佩尔方程

形如:x2−dy2=k,d>1x^2-dy^2=k,d>1x2−dy2=k,d>1

解有迭代公式:
x=px1+dqy1y=py1+qx1x=p x_{1}+d q y_{1}\\y=p y_{1}+q x_{1} x=px1​+dqy1​y=py1​+qx1​
其中(p,q)(p,q)(p,q) 是第二类佩尔方程的一个特解,(x1,y1)(x_1,y_1)(x1​,y1​) 是第一类佩尔方程的最小特解。

推导:根据上述,有p2−dq2=k,x12−dy12=1p^2-dq^2=k,x_{1}^{2}-d y_{1}^{2}=1p2−dq2=k,x12​−dy12​=1则(x12−dy12)(p2−dy2)=k\left(x_{1}^{2}-d y_{1}^{2}\right)\left(p^{2}-d y^{2}\right)=k(x12​−dy12​)(p2−dy2)=k 展开,有
x12p2−dx12q2−dy12p2+d2y12q2=(x12p2+d2y12q2)−d(x12q2+y12p2)=(x1p+dy1q)2−d(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} ===​x12​p2−dx12​q2−dy12​p2+d2y12​q2(x12​p2+d2y12​q2)−d(x12​q2+y12​p2)(x1​p+dy1​q)2−d(x1​q+py1​)2k​
逐次迭代可得上述迭代公式。

对于每一组特解(p,q)(p,q)(p,q) ,第kkk 个迭代解用矩阵表示如下:
[xkyk]=[pdqqp]k−1[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] [xk​yk​​]=[pq​dqp​]k−1[x1​y1​​]
求出第一个特解后用矩阵快速幂求得第kkk 个解

【数论】【不定方程】n元一次不定方程、佩尔方程、毕达哥拉斯定理、费马大定理相关推荐

  1. 数论概论读书笔记 41.连分数与佩尔方程

    连分数与佩尔方程 通过反复地将小数部分翻到分母上并将整数部分分离,我们可以对任何一个数形成连分数: 连分数在视觉上滑向右下方,将它们写成分数却要花费很多的笔墨和空间,由于所有的分子都是1,故我们要做的 ...

  2. 数论读书笔记——线性丢番图方程——解不定方程

    线性丢番图方程: 背景:当我们需要求解特定方程的整数解的时候,那么就得到了一个丢番图方程,这些方程是根据丢番图(diophantus)命名的.他写下了一些方程并将解限定在有理数域上.方程ax+by=c ...

  3. 数论笔记 · 佩尔方程

    转载自:http://blog.sina.com.cn/s/blog_5d06e2390100ll92.html 佩尔方程实际上并不是佩尔提出的,而是费尔马提出,却被欧拉误记为佩尔提出,因此佩尔方程的 ...

  4. 数论笔记 · 佩尔方程(连分数)

    http://blog.sina.com.cn/s/blog_5d06e2390100ll92.html 佩尔方程实际上并不是佩尔提出的,而是费尔马提出,却被欧拉误记为佩尔提出,因此佩尔方程的名称沿用 ...

  5. 第三十四章 数论——高斯消元解线性方程组

    第三十四章 数论--高斯消元解线性方程组 一.高斯消元 1.线性方程组 2.高斯消元步骤 (1)数学知识铺垫 增广矩阵和阶梯矩阵 初等变换 (2)高斯消元步骤 二.代码模板 1.问题: 2.代码 一. ...

  6. 数论 —— 佩尔方程与连分数

    [概述] 连分数是一种特殊的繁分数,其形式为:  ,通常记为:,其中  和  称为连分数多项式,对于任意的 a 均为一次式,它们的比值称为第 n 个渐进值渐进分数. 佩尔(Pell)方程是一种不定二次 ...

  7. POJ 1320 Street Numbers 解佩尔方程

    传送门 Street Numbers Time Limit: 1000MS   Memory Limit: 10000K Total Submissions: 2529   Accepted: 140 ...

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

    title : 佩尔方程 date : 2021-10-31 tags : ACM,数学 author : Linno 佩尔方程 形如x2−dy2=1(d>1且d不为完全平方数)x^2-dy^2 ...

  9. 【枚举算法】佩尔方程

    佩尔方程是关于x.y的二次不定方程,表述为: x^2 - ny^2 = 1  (n为非平方正整数) 当x = 1或-1,y = 0时,满足方程.常把x.y中有一个零的解称为平凡解. 佩尔方程的非平凡解 ...

  10. 佩尔方程(超详细推导+例题讲解) 每日一遍,算法再见!

    这里写目录标题 佩尔方程 第一类佩尔方程 第一类佩尔方程例题讲解 第二类佩尔方程 佩尔方程 第一类佩尔方程 定义:形如x2−dy2=1x^2 - dy^2 = 1x2−dy2=1(d>1,且d不 ...

最新文章

  1. LeetCode 202. Happy Number--Python解法
  2. 009_Raphael绘制图形
  3. 非常全面的Android Bitmap知识点梳理
  4. PHP-代码审计-身份认证
  5. PAT 1007 Maximum Subsequence Sum
  6. 基于SignalR的站点有连接数限制问题及解决方案
  7. Log4j 与 Logback的ConversionPattern对比
  8. oracle支持ipv6,准备现有的网络以支持 IPv6
  9. Centos7安装vlc
  10. 前后端分离项目session跨域失效的解决方案
  11. 基于Silverlight4开发的相关工具
  12. 安装原生Win7-SP1重要补丁
  13. 双向链表、双向循环链表
  14. js initialize
  15. 实例分割最全综述(上):二阶段实例分割和一阶段实例分割
  16. 基于Android的医院挂号系统设计与实现(MySQL+Spring Boot+MyBatis+Android Studio+IDEA)
  17. OpenCV路在何方
  18. 【NLP自然语言处理】文本特征提取
  19. t检验中的t值和p值是什么关系_t检验和p值的关系
  20. cv.bitwise_and用法

热门文章

  1. 钟表维修管理系统技术解析(六) 录入保修单
  2. 相对湿度和绝对湿度的区别
  3. dell台式机进入安全模式_戴尔电脑如何进入安全模式?
  4. 双屏计算机主机是什么,电脑双屏显示有人知道关于一个主机两个显示器的问题吗 爱问知识人...
  5. windows系统c++/VS2019编译gRPC
  6. 计算机网络教程(第四版)
  7. 利用费马小定理进行素性测试
  8. iOS开发-class_ro_t和class_rw_t的区别
  9. (黎活明老师讲学)Android学习(一)---从网络获取图片
  10. [iuud8]ios不同版本下推送本地通知