找回密码
 欢迎注册
查看: 52608|回复: 19

[讨论] 解一个椭圆曲线方程

[复制链接]
发表于 2009-1-17 11:07:40 | 显示全部楼层 |阅读模式

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

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

×
Determine all pairs integer $(x,y)$ satisfying the equation :\[ y(x+y)=x^3-7x^2+11x-3 \]
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2009-1-17 13:15:33 | 显示全部楼层
{x,y} {1,-2} {1,1} {2,-1} {6,-9} {6,3}
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2009-1-17 13:58:14 | 显示全部楼层
原帖由 mathematica 于 2009-1-17 13:15 发表 {x,y} {1,-2} {1,1} {2,-1} {6,-9} {6,3}
能否列出解题过程?谢谢

点评

nyy
把标题改成解椭圆曲线方程的整数解  发表于 2023-5-8 12:20
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2009-1-17 14:50:28 | 显示全部楼层
$(2y+x)^2=4x^3-27x^2+44x-12=(x-2)(4x^2-19x+6)$ 其中$(x-2,4x^2-19x+6)|16$ 由此可以知道必然有$x-2=t^2$或$x-2=2t^2$,然后继续讨论吧

点评

nyy
你为什么不用pari/gp?  发表于 2023-5-9 14:42
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2023-5-6 11:42:41 | 显示全部楼层
northwolves 发表于 2009-1-17 13:58
能否列出解题过程?谢谢

我很笨的,只会穷举法搞!这个好像是椭圆曲线。只有有限个整点
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2023-5-8 19:44:19 | 显示全部楼层
mathe 发表于 2009-1-17 14:50
$(2y+x)^2=4x^3-27x^2+44x-12=(x-2)(4x^2-19x+6)$
其中$(x-2,4x^2-19x+6)|16$
由此可知必有$x-2=t^ ...

令`Y=2y+x,X=x-2`, 方程化为\[
Y^2=X(4X^2-3X-16)\tag1
\]由于`\gcd(X,4X^2-3X-16)|16`, 故\[\begin{split}
X&=s^2, &4X^2-3X-16&=t^2\\
\text{or  } X&=2s^2, &4X^2-3X-16&=2t^2
\end{split}\tag{2}\]
丢番图方程\[
4X^2-3X-16=t^2\tag3
\]可配方化为\[
(8X-3)^2-(4t)^2=265\\
\]易解得`X=4,t=±6`或者`X=17,t=±33`, 取`X`为平方数只有`X=4`, 对应`(x,y)=(6,3),(6,-9)`。
丢番图方程\[
4X^2-3X-16=2t^2\tag4
\]可配方化为\[
(8X-3)^2-2(4t)^2=265\\
\]这个扩展佩尔方程无解。(扩展佩尔方程`x^2-2y^2=z`有解的充要条件是`z`不含有`±3\pmod8`的素因子)

点评

@mathe 括号前多的那个系数4已被提前删除了,漏解与之无关。  发表于 2023-5-10 08:23
前面我的方程有个系数算错了  发表于 2023-5-9 11:27
咋漏解了?  发表于 2023-5-8 22:36
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2023-5-8 22:43:24 | 显示全部楼层
漏解是由于符号原因,还有以下两种情况:
\[\begin{split}
X&=-s^2, &4X^2-3X-16&=-t^2\\
\text{or  } X&=-2s^2, &4X^2-3X-16&=-2t^2
\end{split}\tag{2'}
\]对应于\[
(8X-3)^2+(4t)^2=265\\
(8X-3)^2+2(4t)^2=265
\]
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2023-5-9 09:15:49 | 显示全部楼层
  1. E:=EllipticCurve([1,-7,0,11,-3]);
  2. IntegralPoints(E);
复制代码


算出结果,但是结果不稳定,有两种

第一种结果
  1. [ (1 : 1 : 1), (2 : -1 : 1), (6 : 3 : 1) ]
  2. [ <(1 : 1 : 1), 1>, <(2 : -1 : 1), 1>, <(6 : 3 : 1), 1> ]
复制代码


第二种结果
  1. [ (1 : -2 : 1), (2 : -1 : 1), (6 : -9 : 1) ]
  2. [ <(1 : -2 : 1), 1>, <(2 : -1 : 1), 1>, <(6 : -9 : 1), 1> ]
复制代码


在线计算器
http://magma.maths.usyd.edu.au/calc/

点评

nyy
一次只能搞出一种结果,不知道为什么  发表于 2023-5-9 09:16
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2023-5-9 11:57:21 | 显示全部楼层
nyy 发表于 2023-5-9 09:15
算出结果,但是结果不稳定,有两种

第一种结果
  1. F:=EllipticCurve([1,-7,0,11,-3]);
  2. F;
  3. MW, f := MordellWeilGroup(F);
  4. MW;
  5. a := f(MW.1);
  6. b := f(MW.2);
  7. a;
  8. b;
复制代码


输出结果,谁能看懂这个结果?是这个群是由两个生成元生成的吗?
  1. Elliptic Curve defined by y^2 + x*y = x^3 - 7*x^2 + 11*x - 3 over Rational Field
  2. Abelian Group isomorphic to Z/2 + Z
  3. Defined on 2 generators
  4. Relations:
  5.     2*MW.1 = 0
  6. (2 : -1 : 1)
  7. (1 : 1 : 1)
复制代码
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2023-5-10 09:50:54 | 显示全部楼层
nyy 发表于 2023-5-9 11:57
输出结果,谁能看懂这个结果?是这个群是由两个生成元生成的吗?
  1. SetClassGroupBounds("GRH");
  2. E:=EllipticCurve([0,888889]);
  3. IntegralPoints(E);
  4. MW, f := MordellWeilGroup(E);
  5. MW;
  6. a := f(MW.1);
  7. b := f(MW.2);
  8. a;
  9. b;
复制代码


输出结果:
  1. [ (10950 : 1145833 : 1) ]
  2. [ <(10950 : 1145833 : 1), 1> ]
  3. Abelian Group isomorphic to Z + Z
  4. Defined on 2 generators (free)
  5. (10950 : 1145833 : 1)
  6. (283368/961 : 153436021/29791 : 1)
复制代码
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
您需要登录后才可以回帖 登录 | 欢迎注册

本版积分规则

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

GMT+8, 2024-12-23 01:26 , Processed in 0.026781 second(s), 19 queries .

Powered by Discuz! X3.5

© 2001-2024 Discuz! Team.

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