- 注册时间
- 2009-2-12
- 最后登录
- 1970-1-1
- 威望
- 星
- 金币
- 枚
- 贡献
- 分
- 经验
- 点
- 鲜花
- 朵
- 魅力
- 点
- 上传
- 次
- 下载
- 次
- 积分
- 26215
- 在线时间
- 小时
|
楼主 |
发表于 2025-5-12 21:41:52
|
显示全部楼层
我来走一下流程,熟悉熟悉Maple和PARI
我发现没有一个数学工具可以一步到位,比如Maple有Weierstrass变换,但是没有后续的椭圆曲线相关的点的运算。 而PARI/Gp没有Weierstrass变换。
首先 对$\frac{a}{b+c}+\frac{c}{a+b}+\frac{b}{a+c}-n $通分,并令$c=1$,就是$f(a,b,c,n) = a^3-a^2 b n+a^2 b-a^2 c n+a^2 c-a b^2 n+a b^2-2 a b c n+3 a b c-a c^2 n+a c^2+b^3-b^2 c n+b^2 c-b c^2 n+b c^2+c^3$
用Maple发现
- gg := -a^2*b*n - a*b^2*n + a^3 + a^2*b - a^2*n + a*b^2 - 2*a*b*n + b^3 - b^2*n + a^2 + 3*a*b - a*n + b^2 - b*n + a + b + 1;
- Weierstrassform(gg, a, b, x, y)
复制代码
得到对应的Weierstrass方程是$-\frac{1}{48} (2 n+3) \left(8 n^3+36 n^2+6 n-93\right) x-\frac{1}{864} \left(4 n^2+12 n-3\right) \left(16 n^4+96 n^3+120 n^2-216 n-423\right)+x^3+y^2 = 0$,
变换关系是$x\to a^2 (-n-2)+a \left(b (n+2) n+n^2+n-3\right)+b^2 (-n-2)+b \left(n^2+n-3\right)+\frac{2 n^2}{3}+n-\frac{7}{4}$
$y\to a^3 (n+2)^2+a^2 \left(-b n (n+2)^2-\frac{1}{2} (2 n (n+1)-5) (n+2)\right)+a \left(b^2 (n+2)^2-\frac{1}{2} b (2 n-3) (n+2)^2+\frac{1}{2} (11-n (n (2 n+7)+1))\right)+b^2 \left(-\frac{n}{2}-1\right)+\frac{1}{2} b \left(n^2+n-3\right)+\frac{1}{2} (n+2)^2$
逆变换表达式就是${a\to \frac{4 n^2+36 n+12 x-24 y+69}{2 \left(4 n^3+20 n^2+12 n x+9 n+24 x-42\right)},b\to \frac{4 n^2+36 n+12 x+24 y+69}{2 \left(4 n^3+20 n^2+12 n x+9 n+24 x-42\right)}}$
然后该PARI出场了。
- n=4;
- a4=-1/48*(3 + 2*n)*(-93 + 6*n + 36*n^2 + 8*n^3);
- a6=(1/864)*(-3+12*n+4*n^2)*(-423-216*n+120*n^2+96*n^3+16*n^4);
- e=ellinit([0,0,0,a4,a6]);
- pts=ellratpoints(e,10000000)
- file=fileopen("E:/pari_data.txt", "a");foreach(pts, p, filewrite(file, p));fileclose(file)
复制代码
- [[103/36, 728/27], [103/36, -728/27], [5353/576, 28171/13824], [5353/576, -28171/13824], [-191/12, 65/2], [-191/12, -65/2], [-59/12, 49], [-59/12, -49], [97/12, 7/2], [97/12, -7/2], [109/12, 0], [121/12, 13/2], [121/12, -13/2], [277/12, 91], [277/12, -91], [2557/300, 91/125], [2557/300, -91/125], [29233/300, 237679/250], [29233/300, -237679/250], [7369/588, 13949/686], [7369/588, -13949/686], [109957/1452, 860496/1331], [109957/1452, -860496/1331], [26653/2028, 52577/2197], [26653/2028, -52577/2197], [244981/26508, 176995/103823], [244981/26508, -176995/103823], [41/6, 39/4], [41/6, -39/4], [-449023/71286, 263764935/5180116], [-449023/71286, -263764935/5180116]]
复制代码
通过变换,得到$\frac{a}{b+c}+\frac{c}{a+b}+\frac{b}{a+c}=4$的时候,a,b,c的值是
- {-1,4,11}
- {-5,9,11}
- {-11,-9,5}
- {-11,-4,1}
- {-8784,5165,9499}
- {-9071,5951,9841}
- {-9499,-5165,8784}
- {-9841,-5951,9071}
- {-375326521,679733219,883659076}
- {-934668779,-396650011,137430135}
复制代码
得到$\frac{a}{b+c}+\frac{c}{a+b}+\frac{b}{a+c}=14$的时候,a,b,c的值是
- {-7,8,13}
- {-13,-8,7}
- {-33,37,47}
- {-47,-37,33}
- {-2075,683,2117}
- {-2117,-683,2075}
- {-2944,2769,3031}
- {-3031,-2769,2944}
- {-70242391,79121224,87621389}
- {-77350969,96250941,258686551}
复制代码
|
|