找回密码
 欢迎注册
查看: 54750|回复: 26

[求助] 如何求解pell方程x^2-d*y^2=M

[复制链接]
发表于 2019-3-19 14:06:38 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?欢迎注册

×
x^2-d*y^2=M
其中M的绝对值大于d的平方根,
要是M的绝对值小于d的平方根,这个只要连分数就可以了!
但是大于怎么办?
比如
x^2-15y^2=61
我要的不是穷举法,我就是想知道如何求解


https://en.wikipedia.org/wiki/Pell%27s_equation
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2019-3-19 18:18:41 | 显示全部楼层
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2019-3-19 20:48:49 | 显示全部楼层
对于方程 `x^2-15y^2=61`,因为 `25^2\equiv 15 \pmod {61}`,故取二次无理分式展开 `\D\frac{25+\sqrt{15}}{61}=[0,2,8,\overline{1,6}]`,容易检验发现第三个完全商 `[8,\overline{1,6}]=5+\sqrt{15}`,其分母为1,故取第二个渐近分数 `p_2/q_2=[0,2]=1/2`,于是有 `(25q_2-61p_1)^2-15q_2^2=(-1)^261`.
因此得到最小解 `x_0=|25q_2-61p_1|=11,\;y=q_2=2`.
然后求标准形式 `x^2-15y^2=1` 的最小解,连分式展开 `\sqrt{15}=[3,\overline{1,6}]`,因为 `6=2\times 3`,故取渐近分数 `p_1/q_1=[3,1]=4`,得到最小解`x=p_1=4,\;y=q_1=1`,再根据上面链接中5楼的公式,可得通解。

点评

nyy
应该是25q2−61p2,而不是25q2−61p1,应该是这个  发表于 2025-7-13 21:47
确实,上面输入有误,应该是 `p_1`  发表于 2019-3-21 13:33
应该是25q2−61p2,而不是25q2−61p1  发表于 2019-3-21 13:28
我发现虽然你是力学毕业的,但是你很牛  发表于 2019-3-21 12:11
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2019-3-21 08:30:56 | 显示全部楼层
kastin 发表于 2019-3-19 20:48
对于方程 `x^2-15y^2=61`,因为 `25^2\equiv 15 \pmod {61}`,故取二次无理分式展开 `\D\frac{25+\sqrt{15} ...

我今天看华罗庚文集第274页,
感觉算看懂了,与你这个类似,
但是感觉他的办法似乎比你的好!

点评

但是他的办法似乎是穷举法!  发表于 2019-3-21 09:36
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2025-7-11 14:03:04 | 显示全部楼层
我不会
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
回复

使用道具 举报

发表于 2025-7-11 20:25:05 | 显示全部楼层
pell方程是genus=0的曲线,所有genus=0的曲线都有参数解,且参数解能够完整刻画曲线的所有解. 该pell方程$x^2-15y^2=61$的参数解是
$[x,y]=[\frac{14 U^2-90 U+210}{U^2-15},\frac{3 U^2-28 U+45}{U^2-15}] = [-14+\frac{30 n (3 m-14 n)}{m^2-15 n^2},-3+\frac{2 n (14 m-45 n)}{m^2-15 n^2}]$


递归公式是$a(n+4)-8 a(n+2)+a(n)=0$
  1. LinearRecurrence[{0,8,0,-1},{{11,2},{14,3},{74,19},{101,26}},40]
复制代码

  1. RecurrenceTable[{a[n+4]==8a[n+2]-a[n],a[1]=={11,2},a[2]=={14,3},a[3]=={74,19},a[4]=={101,26}},a,{n,40}]
复制代码

点评

nyy
看不懂  发表于 2025-7-13 21:47
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2025-7-12 13:14:11 | 显示全部楼层
wayne 发表于 2025-7-11 20:25
pell方程是genus=0的曲线,所有genus=0的曲线都有参数解,且参数解能够完整刻画曲线的所有解. 该pell方程$x^2 ...

你这个是什么意思呀?
看不懂
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2025-7-12 16:29:56 | 显示全部楼层
https://math.stackexchange.com/q ... n-to-pells-equation

拉马努金的圆周率原来也与pell方程的有关系
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2025-7-12 16:47:55 | 显示全部楼层
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2025-7-13 19:26:17 | 显示全部楼层
Lagrange-Matthews-Mollin算法

点评

nyy
搜索这个就能找到  发表于 2025-7-14 08:21
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
您需要登录后才可以回帖 登录 | 欢迎注册

本版积分规则

小黑屋|手机版|数学研发网 ( 苏ICP备07505100号 )

GMT+8, 2025-8-3 19:46 , Processed in 0.030789 second(s), 17 queries .

Powered by Discuz! X3.5

© 2001-2025 Discuz! Team.

快速回复 返回顶部 返回列表