如何求解pell方程x^2-d*y^2=M
x^2-d*y^2=M其中M的绝对值大于d的平方根，
要是M的绝对值小于d的平方根，这个只要连分数就可以了！
但是大于怎么办？
比如
x^2-15y^2=61
我要的不是穷举法，我就是想知道如何求解
https://en.wikipedia.org/wiki/Pell%27s_equation 见https://bbs.emath.ac.cn/thread-15533-1-1.html
以及http://mathworld.wolfram.com/PellEquation.html最后部分（和相应的参考文献） 对于方程 `x^2-15y^2=61`，因为 `25^2\equiv 15 \pmod {61}`，故取二次无理分式展开 `\D\frac{25+\sqrt{15}}{61}=`，容易检验发现第三个完全商 `=5+\sqrt{15}`，其分母为1，故取第二个渐近分数 `p_2/q_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}=`，因为 `6=2\times 3`，故取渐近分数 `p_1/q_1==4`，得到最小解`x=p_1=4,\;y=q_1=1`，再根据上面链接中5楼的公式，可得通解。 kastin 发表于 2019-3-19 20:48
对于方程 `x^2-15y^2=61`，因为 `25^2\equiv 15 \pmod {61}`，故取二次无理分式展开 `\D\frac{25+\sqrt{15} ...
我今天看华罗庚文集第274页,
感觉算看懂了,与你这个类似,
但是感觉他的办法似乎比你的好! 我不会 pell方程是genus=0的曲线,所有genus=0的曲线都有参数解,且参数解能够完整刻画曲线的所有解. 该pell方程$x^2-15y^2=61$的参数解是
$=[\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$
LinearRecurrence[{0,8,0,-1},{{11,2},{14,3},{74,19},{101,26}},40]
RecurrenceTable[{a==8a-a,a=={11,2},a=={14,3},a=={74,19},a=={101,26}},a,{n,40}] wayne 发表于 2025-7-11 20:25
pell方程是genus=0的曲线,所有genus=0的曲线都有参数解,且参数解能够完整刻画曲线的所有解. 该pell方程$x^2 ...
你这个是什么意思呀？
看不懂 https://math.stackexchange.com/questions/1045127/how-to-find-a-fundamental-solution-to-pells-equation
拉马努金的圆周率原来也与pell方程的有关系 https://math.stackexchange.com/questions/1496083/general-method-for-determining-if-ax2-bx-c-is-square Lagrange-Matthews-Mollin算法
页:
[1]
2