丢番图逼近之二 -- 用连分数的方法求解Pell方程
引言
上一篇的连分数主要也只是一个科普内容,还有一些内容不便于放在上篇里边,于是加在本文中.
由于本文并不像上一篇那样会面向学高数的同学,因此本文不会那么通俗,而是会比较严肃一点.
Pell方程
形如
最早碰到这个问题还是高中阶段,兜兜转转又要面对它了,而且还有可能出现在我初等数论的考试中😫.当初面对的还只是
但是!过去软糯的我已经死啦!现在是更软糯的我!拥有代数数论和连分数知识的我已经无往不利了!
利用代数数论的知识能够得到Pell方程解的结构,而利用连分数这个利器,求解这类一般的Pell方程已经成一个消磨时间的打趣罢了.
Pell方程解的结构
Dirichlet单位定理
在代数数论中,有如下事实:
定理1(Dirichlet单位定理): 设
为 次数域, 到 中有 个实嵌入,有 对复嵌入,则: 中单位群 ,其中 是单位根群, 是秩为 的自由Abel群.
其中单位群
并且由定理条件可知,存在元素
因此,若
回到Pell方程
由于
的情况
这种情况一般而言比较简单,因为
定理2: 当
时,设 是 的基本单位,令 . (1) 如果
:
a. 不存在有理整数解.
b. 的有理整数解为 . (2) 如果
:
a. 的有理整数解为 .
b. 的有理整数解为 .
当然对于另一种情况也是类似的考虑,只不过需要多转一个弯而已,比起其他代数里面的”一个弯”,这个简直可以说是友好至极了🤠.
的情况
在这种情况下,
从而
因此便有了和前一种情况的第一处不同,此处应该令
引理1: 当
时,设 是 的基本单位,令 . (1) 如果
:
a. 不存在有理整数解.
b. 的有理整数解为 . (2) 如果
:
a. 的有理整数解为 .
b. 的有理整数解为 .
可以发现的是,该引理和定理2是高度相似的,也只有
因此,对于
定理3: 当
时,设 是 的基本单位,令 . (1) 如果
:
a. 不存在有理整数解.
. 如果 ,则: 的有理整数解为
. 如果 ,则: 的有理整数解为 (2) 如果
:
. 如果 ,则: 的有理整数解为
. 如果 ,则: 的有理整数解为
. 如果 ,则: 的有理整数解为
. 如果 ,则: 的有理整数解为
这个定理看上去比较吓人,但实际上也确实比较吓人了.其中
基本单位的大小估计
上面只得到了Pell方程解的结构,虽然已经可以通过手动改变
事实上,对该问题基本单位的上界的估计也是有一些结果的.而最后的估计式由华罗庚给出:
其中
取
当
无平方因子,则 的基本单位为 . 当
无平方因子,则 的基本单位为 .
利用连分数求基本单位
与
在上一篇中介绍有
定理4(Lagrange定理): 循环连分数和二次无理数是一一对应的.
在Fatfish提供的证明二中,引入了两个新的符号
但是看这推导式,不如不看.总之在最后计算基本单位的时候,我们也只是用到了连分数展开的第一个循环而已,因此手撕就差不多够了,即使是61这个数(Fermat自己说是比较人性化的选了个比较小的数).
例如:
通过这种暴力手撕的方式,我们也能够得到
基本单位的确定
以下两个定理,能够说明Pell方程的有理整数解和连分数的渐近分数之间的关系:
定理5: 记
,其中 ,则 一定是 的渐近分数. 定理6: 设
的渐近分数为 ,则 .即 .
那么我们确定基本单位的方法基本上已经很明确了:
例如: 对于
最后来看一下Fermat所选取的两个比较人性化的数61和109到底有何特别之处,让我们先对
这里展开的循环节足足11位!而还有糕手!109更是重量级!它的循环节有15位!
因此这两个数对应的Pell方程的超大的基本单位,就能从这两个连分数展开中窥见一斑.因此Fermat这个糟老头子信不了一点.🤬
从19号开始写到25号写完本文,而且最后连分数里边还有一些内容被我很快的跳过去了,比如Legender判别准则来判断渐近分数,Liouville超越数这种十分逼近有理数的数,以及太超前的辛钦定理.这些内容之后可以再回过头来看一看,还是很神奇的.
参考资料
[1] 冯克勤.代数数论[M].哈尔滨:哈尔滨工业大学出版社,2018
参考第三章2.2节-实二次域的基本单位,Pell方程
[2] FatFish.连分数入门[Z]:超理论坛.https://chaoli.club/index.php/2756?see_lz=1
又是FatFish的帖子,主要参考3.2-Lagrange定理以及3.5-Pell方程