Homeに戻る  一覧に戻る 

Integer Points on A^4+20011*B^4=C^4+20011*D^4

[2026.07.16]A^4+20011*B^4=C^4+20011*D^4の整点

■Diophantine Equation
       A^4+n*B^4=C^4+n*D^4 ----------(1)
を満たす自明でない整数の組(A,B,C,D) (ただし 0 < A < C, 0 < D < Bかつgcd(A,B,C,D)=1)を探す。

以下では、Elkiesの論文(参考文献[1])の方法によって、(1)を満たす整数の組(A,B,C,D)を探す。

■(1)およびD!=0より、A/D=y,B/D=x.C/D=tとすると、
       y^4+n*x^4 = t^4+n ----------(2)
       y^4-t^4= y n(1-x^4)
       (y^2+t^2)*(x^2-t^2) = n*(1+x^2)*(1-x^2)
つまり、(2)を満たす有理数の組(x,y,t)を見つければ良い。

ここで、ある有理数u > 0に対して、
       y^2+t^2 = (n/u)*(1+x^2) ----------(2a)
       y^2-t^2 = u*(1-x^2) ----------(2b)
よって、
       2*u*y^2 = (n+u^2)*x^2+(n-u^2) ----------(3a)
       2*u*t^2 = (n-u^2)*x^2+(n+u^2) ----------(3b)
を満たす有理数の組(x,y,t)が存在すれば、(x,y,t)が(2)を満たすことが分かる。

[pari/gpによる計算]
gp > YY2(n,u,x)
%1 = (1/(2*u)*n - 1/2*u)*x^2 + (1/(2*u)*n + 1/2*u)
gp > TT2(n,u,x)
%2 = (1/(2*u)*n + 1/2*u)*x^2 + (1/(2*u)*n - 1/2*u)
gp > YY2(n,u,x)^2+n*x^4-TT2(n,u,x)^2-n
%3 = 0

■２次曲線(3a),(3b)は、n=u^2のときはsoncularであり、それ以外のときはnon-singularである。
２次曲線(3a)の右辺の判別式は
    -4*(n+u^2)*(n-u^2)
となる。ここで、nが平方数でなければ、任意の有理数uについて、判別式は0にならないので、２次曲線(3a)は常にnon-singularである。
nが平方数であるならば、n=u^2のときに限り、２次曲線(3a)はsingularであり。それ以外のuにつては、non-singularである。

同様に、２次曲線(3b)の右辺の判別式は
    -4*(n-u^2)*(n+u^2)
となる。nが平方数でなければ、任意の有理数uについて、判別式は0にならないので、２次曲線(3b)は常にnon-singularである。
nが平方数であるならば、n=u^2のときに限り、２次曲線(3b)はsingularであり。それ以外のuにつては、non-singularである。

■以下では、n=20011とする。

Pari/GPで簡単なプログラム("a45s.gp")を作成して、max{A,B,C,D} <= 100000の範囲で、小さい整数解を探すと、1個の整数解
    48159^4+20011*9898^4=88181^4+20011*9096^4
が見つかった。
これらの整数解から、(3a),(3b)を満たす有理数uをいくつか求めると、
    u = 3401870/9497
となる。
[gp2cによる計算]
gp > sss(20011,100000)
1:48159^4+20011*9898^4=88181^4+20011*9096^4
1 solutions,
time = 3h, 12min, 27,494 ms.
%1 = 1
{Pari/GPによる計算]
gp > uu(20011,48159,9898,88181,9096)
[u, -1; 802*u + 531493, 1; 9497*u - 3401870, 1]
[u, -1; 802*u - 531493, 1; 9497*u - 3401870, 1]
[u, -1; 802*u + 531493, 1; 9497*u + 3401870, 1]
time = 1 ms.
%1 = [-531493/802, 3401870/9497]

■有理数uの高さが小さいものから、順に調べる。
例えば、有理数uの高さが1000以下の範囲で、２つの２次曲線(3a)と(3b)が共に有理点を持つようなuを選択すると、
このように7361個のuが抽出される。
[MAGMAによる計算]
> PP(20011,1,1000);
**u= 1/24 , ht= 1229 ; ( -3403/629 , -1695828/629 ), ( 1229/83 , 603622/83 )
**u= 1/36 , ht= 6767 ; ( 2465/6769 , 4323504/6769 ), ( 2467/6767 , 4322787/6767 )
**u= 1/60 , ht= 53987 ; ( -53987/16493 , -43738092/16493 ), ( 40159/69619 , 62272481/69619 )
**u= 1/64 , ht= 853 ; ( -467/12797 , 20494463/25594 ), ( -31/853 , 2732153/3412 )
**u= 1/96 , ht= 231167 ; ( 231167/12985 , 90766359/5194 ), ( 243655/47801 , -243349830/47801 )
**u= 1/140 , ht= 94163 ; ( -132611/32511 , 53866145/10837 ), ( -23043/94163 , 114734244/94163 )
**u= 1/144 , ht= 4347 ; ( 21737/18903 , -103732298/56709 ), ( 4347/3781 , -20746318/11343 )
**u= 1/171 , ht= 32537 ; ( -676631/50195 , -177497023/10039 ), ( 9015/32537 , -132488185/97611 )
**u= 1/172 , ht= 22645 ; ( -234205/534947 , 766079290/534947 ), ( 11239/22645 , -6632875/4529 )
**u= 1/244 , ht= 8729 ; ( -8729/8363 , 18888251/8363 ), ( -33167/10875 , -1880599/375 )
**u= 1/256 , ht= 103 ; ( 103/39 , -117511/26 ), ( 368169/215833 , -2732076489/863332 )
**u= 1/264 , ht= 21623 ; ( -21623/5321 , 36191288/5321 ), ( -86963/151291 , -40516123/21613 )
**u= 1/291 , ht= 1276541 ; ( -1276541/364699 , -2265364701/364699 ), ( -5685011/217051 , 9707636411/217051 )
**u= 1/315 , ht= 87433 ; ( 816513/178393 , 4451276663/535179 ), ( -87433/15887 , 157762477/15887 )
**u= 1/316 , ht= 21337 ; ( 2377/21337 , 38174612/21337 ), ( 7133/33203 , -60386192/33203 )
**u= 1/324 , ht= 61077 ; ( -21707/61077 , 1050367906/549693 ), ( -1150047/305809 , 2754783748/393183 )
**u= 1/372 , ht= 210551 ; ( -125/210551 , 406207795/210551 ), ( -1085849/113975 , -421278816/22795 )
**u= 1/460 , ht= 1696531 ; ( 3260397/198803 , -7007686210/198803 ), ( -1696531/745619 , -567950990/106517 )
**u= 1/472 , ht= 321247 ; ( -45025/321247 , -704942600/321247 ), ( 1846771/383925 , 163964947/15357 )
**u= 1/484 , ht= 6443 ; ( 393/6443 , -14204849/6443 ), ( 5897/96647 , 213077396/96647 )
**u= 1/491 , ht= 402547 ; ( 402547/215967 , -337509011/71989 ), ( -2365719/1515145 , -6226753629/1515145 )
**u= 1/520 , ht= 92347 ; ( 72688/92347 , -1072263941/369388 ), ( -132917/17493 , 43685018/2499 )
**u= 1/536 , ht= 77237 ; ( 77237/2091 , -59643762/697 ), ( -2963283/1806997 , -8037636237/1806997 )
**u= 1/616 , ht= 33739 ; ( 5389/53899 , -134477749/53899 ), ( 31627/33739 , 114808288/33739 )
**u= 1/619 , ht= 3270469 ; ( 3270469/31861 , -8139455767/31861 ), ( 11702941/143735 , -29126781083/143735 )
**u= 1/628 , ht= 140843 ; ( 231815/100249 , -633095030/100249 ), ( 140843/37195 , -73030478/7439 )
**u= 1/636 , ht= 292319 ; ( 114640945/735253 , -289198905573/735253 ), ( 292319/9617 , 737802124/9617 )
**u= 1/676 , ht= 8755 ; ( -32293/131323 , -4572207706/1707199 ), ( -2153/8755 , -304819301/113815 )
**u= 1/691 , ht= 347311 ; ( 347311/202201 , -1056716201/202201 ), ( 72495569/20723085 , -198255712771/20723085 )
**u= 1/739 , ht= 11711 ; ( 4823/102371 , 278676287/102371 ), ( 11711/2843 , 32769517/2843 )
**u= 1/756 , ht= 190613 ; ( -190613/76619 , -565009834/76619 ), ( 186720167/5230355 , -1541214538309/15691065 )
**u= 1/820 , ht= 473249 ; ( 67349/473249 , 1369209505/473249 ), ( 200289691/4962411 , 573876152132/4962411 )
**u= 1/859 , ht= 296779 ; ( 5108633/20805679 , 62807313307/20805679 ), ( 296779/435 , 870060823/435 )
**u= 1/900 , ht= 2079 ; ( -151769/96873 , 8104503436/1453095 ), ( -2079/1327 , 111018892/19905 )
**u= 1/931 , ht= 1593 ; ( -61133/7173 , 20873563/797 ), ( 1593/73 , -4867047/73 )
**u= 1/976 , ht= 1620945 ; ( -1620945/623107 , -10853509049/1246214 ), ( -61266223/396031 , -191458427477/396031 )
**u= 2/35 , ht= 386 ; ( 386/141 , -114639/94 ), ( 40834/8271 , -34867549/16542 )
**u= 2/47 , ht= 115520 ; ( 115520/5963 , -112180741/11926 ), ( 290009/111268 , -301241647/222536 )
**u= 2/83 , ht= 1030 ; ( 56123/24470 , -15781067/9788 ), ( -1030/947 , 1803215/1894 )
**u= 2/87 , ht= 5936 ; ( 245000/28631 , -325465917/57262 ), ( -589/5936 , -1124391/1696 )
**u= 2/119 , ht= 2521 ; ( -1100/2521 , 4244485/5042 ), ( 261763/140100 , -91631291/56040 )
...省略...
**u= 986/975 , ht= 679429 ; ( 347731238/7936861 , -345962451589/79368610 ), ( -679429/425402 , -797361549/4254020 )
**u= 997/7 , ht= 4759 ; ( 4759/503 , -3733/503 ), ( -2131/25799 , -301/25799 )
**u= 997/8 , ht= 86750 ; ( -86750/9747 , 514327/12996 ), ( 2488221/185434 , -118893093/741736 )
**u= 997/60 , ht= 917099 ; ( 1756237/423697 , 44058723/423697 ), ( 917099/3221 , -22659066/3221 )
**u= 997/143 , ht= 49015 ; ( 49015/837 , -618277/279 ), ( 178569/1156403 , -44275413/1156403 )
**u= 997/172 , ht= 3973689 ; ( 733687411/13052223 , -10153830041/4350741 ), ( -3973689/752731 , -168160659/752731 )
**u= 997/183 , ht= 16811519 ; ( -169607555/41206703 , -7474971689/41206703 ), ( -16811519/4292041 , -43767261/252473 )
**u= 997/192 , ht= 140071 ; ( -140071/10679 , 12324513/21358 ), ( -632849/82477 , 112130159/329908 )
**u= 997/207 , ht= 2388227 ; ( 9782217/1878313 , 1361267669/5634939 ), ( -2388227/578727 , 336177299/1736181 )
**u= 997/247 , ht= 5213707 ; ( 5213707/44491 , 259481191/44491 ), ( -160463549/23147853 , 8074927601/23147853 )
**u= 997/255 , ht= 3160961 ; ( 2456549/3160961 , -202534773/3160961 ), ( 9551101/16823701 , -978461725/16823701 )
**u= 997/260 , ht= 94641497 ; ( 94641497/37448907 , -1732562269/12482969 ), ( 202210161/3438341 , 10334361603/3438341 )
**u= 997/287 , ht= 443325 ; ( 5812493/2917965 , -23265445/194531 ), ( 443325/215447 , 26457945/215447 )
**u= 997/312 , ht= 921091 ; ( 921091/88610 , 41412669/70888 ), ( -399281185/91347031 , -22924807605/91347031 )
**u= 997/316 , ht= 2077233 ; ( -96545/2077233 , 39044100/692411 ), ( -2077233/96545 , -23426460/19309 )
**u= 997/352 , ht= 702747 ; ( -55508153/2234532 , -8803056323/5958752 ), ( -702747/106963 , 84513969/213926 )
**u= 997/359 , ht= 7643155 ; ( 10329245/1007451 , -207606125/335817 ), ( -1333281/7643155 , 93121965/1528631 )
**u= 997/367 , ht= 546440877 ; ( 621093715/5374879 , 37687572139/5374879 ), ( -59517227/546440877 ,
33352676803/546440877 )
**u= 997/368 , ht= 1862571 ; ( -14731861/23784951 , -1133579377/15856634 ), ( -1862571/646315 , 119828352/646315 )
**u= 997/372 , ht= 2969863 ; ( -2969863/305783 , 182386748/305783 ), ( 7851723871/117653221 , 68554649407/16807603 )
**u= 997/388 , ht= 2473333 ; ( 49312037/1256879 , -236737321/96683 ), ( 976235/2473333 , 165904310/2473333 )
**u= 997/400 , ht= 1173745 ; ( 1173745/206947 , -755022241/2069470 ), ( 3464525/2547261 , -1362315397/12736305 )
**u= 997/408 , ht= 3502696 ; ( 697951/3502696 , -914280043/14010784 ), ( -52394413/21962771 , -3635657764/21962771 )
**u= 997/472 , ht= 20365235 ; ( 2164678221/150305944 , -597300144931/601223776 ), ( -854587/20365235 , 82511999/1197955
)
**u= 997/495 , ht= 35995463 ; ( 35995463/1004643 , 7613223299/3013929 ), ( 724115817/76865587 , 153985519087/230596761 )
**u= 997/516 , ht= 34727891 ; ( -1568785/34727891 , -2501831715/34727891 ), ( -213293963/32497195 , -3105467306/6499439
)
**u= 997/532 , ht= 1857841 ; ( -6492203867/1531893057 , -162453802149/510631019 ), ( 1857841/382411 , 138605939/382411 )
**u= 997/559 , ht= 193171 ; ( 32756115/6749017 , -2504761595/6749017 ), ( 193171/155645 , -530881/4447 )
**u= 997/568 , ht= 5078803 ; ( -5078803/3553581 , 155992696/1184527 ), ( 87333954/29324939 , -3974807361/16757108 )
**u= 997/575 , ht= 3712463 ; ( 3289755593/4539758121 , -709822785079/7566263535 ), ( -1496439/3712463 ,
1520220867/18562315 )
**u= 997/588 , ht= 12895375 ; ( 37248929/39053425 , -5804092768/54674795 ), ( 12895375/769103 , 6946955630/5383721 )
**u= 997/592 , ht= 333207919 ; ( -720995589/4168891 , -55570174726/4168891 ), ( 333207919/155012219 ,
-56655219343/310024438 )
**u= 997/607 , ht= 653027 ; ( -577969/653027 , 68063953/653027 ), ( -622131517/146142313 , -49881328057/146142313 )
**u= 997/640 , ht= 456732 ; ( 456732/171367 , -625494881/2741872 ), ( -4677639/2902244 , -7058926143/46435904 )
**u= 997/663 , ht= 14836439 ; ( -21005201/17572525 , -89354901/702901 ), ( 8875535/14836439 , 1410185795/14836439 )
**u= 997/684 , ht= 522172595 ; ( 522172595/29405581 , 43328990890/29405581 ), ( -1711909629/17698105 , 5006443157/624639
)
**u= 997/759 , ht= 18503305 ; ( 574370155/72955967 , 50529139345/72955967 ), ( 3429587/18503305 , 328465119/3700661 )
**u= 997/807 , ht= 297019621 ; ( -2024999675/199248883 , 183108917153/199248883 ), ( 297019621/240200365 ,
34376749111/240200365 )
**u= 997/816 , ht= 2392583 ; ( -156258155/16903879 , -28444625685/33807758 ), ( 2392583/1646005 , 105122085/658402 )
**u= 997/823 , ht= 30783345 ; ( 30783345/12242801 , -3010665637/12242801 ), ( -177998711/1774117 , -16178046737/1774117
)
**u= 997/852 , ht= 14282707 ; ( 15493109/15513557 , -2027364342/15513557 ), ( -14282707/4977889 , 1398645433/4977889 )
**u= 997/892 , ht= 9229773 ; ( -81227147/2057241 , -2562476537/685747 ), ( -9229773/3239101 , -925500192/3239101 )
**u= 997/928 , ht= 3665056 ; ( -1407068/21563721 , 5561260817/57503256 ), ( 3661547/3665056 , 3999663857/29320448 )
**u= 997/940 , ht= 1805647 ; ( 30924069/66657049 , 7137051031/66657049 ), ( -1805647/1304157 , -216337777/1304157 )
**u= 997/984 , ht= 4437409 ; ( 129909985/33409556 , -53317405705/133638224 ), ( -4437409/2292535 , 99267941/458507 )
**u= 997/988 , ht= 93298507 ; ( 126939493/3011575 , -2528646941/602315 ), ( -93298507/21059035 , -9524135984/21059035 )
**u= 997/991 , ht= 40247525 ; ( -60626015/2222421 , 2016633915/740807 ), ( 17024347/40247525 , -871590761/8049505 )
7361
>
また、小さい整点から求めたuについても調べる。

■これらのuについて、(2),(3a),(3b)を満たす有理数解(x,y,t)を持たないものもあれば、有理数解(x,y,t)を持つものもある。
これらのuを順に調べれば良い。

ここからは、A^4+B^4=C^4+D^4と同様なので、最終的に得られた整点のみ記述する。
ここで、対応する整点が見つかった各有理数uについて、0 < A < B, 0 < C < D, 0 < A < Cを満たすように、
A,B,C,Dの符号を変更したり、A,B,C,Dを交換して、Dの小さい順に並び替えると、以下のようになる。

■u=3401870/9497のときの楕円曲線およびそのMordell-Weil群の基底は、以下の通りである。

変数 値
n, u 20011, 3401870/9497
(3a)
 y^2 = -488124891/3228980*x^2 + 668510909/3228980
(3b),
(x0,t0) 		t^2 = 668510909/3228980*x^2 - 488124891/3228980,
(-4548/4949, 48159/9898)
xkt(k)
 (14685401040*k^2 + 155504447820*k + 3040387614132)/(-15980222020*k^2 + 3308460488641)
一変数多項式の等式
xx^2+n*y^4=z^4+n*w^4 		(81073841058447585744400*x^4 - 2761749046308083386366080*x^3 - 156186525536754891301504376*x^2 - 571777888180570929106834464*x + 3475097839005331995329307841)^2+20011*(29370802080*x^2 + 311008895640*x + 6080775228264)^4
= (155504447820*x^2 + 12161550456528*x + 32194816866531)^4+20011*(-31960444040*x^2 + 6616920977282)^4
楕円曲線 En,u y^2 = 81073841058447585744400*x^4 - 2761749046308083386366080*x^3 - 156186525536754891301504376*x^2 - 571777888180570929106834464*x + 3475097839005331995329307841)^2+20011*(29370802080*x^2 + 311008895640*x + 6080775228264
楕円曲線 Estdn.u ( En,uのminimal standard model ) y^2 = x^3 + x^2 - 50004566528123296237180656993911520*x - 492831845561205090049555477799439776999549996781132
Estdn.u(Q)の基底 (2532765438480553990368763746/7775888761, 91421630249318460389985117263764565135040/685685646833741)
height 16.0326980313085706081149797115
En.u(Q)の基底 (0, 58949960466529)
En.uの有理点(x,y)のx座標 0,
-16053/3032,
-506731269022/12958703985,
223879820934094/3421354908729,
-189236734896984828/24182697107313025,
-20026629864887984305/756946939587939312,
2287213080046035824661/722903464199770844120,
-1790289524054760751625436165470569932990866/50721081126660309931935990779467021903389,
102158436166424744281737173937044128511849748/152069345226187023246554322942823320845913275,
125934525711764207290919313898985281104373681355/408633744665698977126948695748176514047398992,
...

[Pari/GPによる計算]
gp > uxt(20011, 3401870/9497, -4548/4949, 48159/9898 )
(3b+) TT2(x)=668510909/3228980*x^2 - 488124891/3228980
%416 = (14685401040*k^2 + 155504447820*k + 3040387614132)/(-15980222020*k^2 + 3308460488641)
gp > xkt(k)=(14685401040*k^2 + 155504447820*k + 3040387614132)/(-15980222020*k^2 + 3308460488641)
%417 = (k)->(14685401040*k^2+155504447820*k+3040387614132)/(-15980222020*k^2+3308460488641)
gp > PP(YY2(20011, u0, xkt(x)))
4 real roots.
time = 2 ms.
%418 = 81073841058447585744400*x^4 - 2761749046308083386366080*x^3 - 156186525536754891301504376*x^2 - 571777888180570929106834464*x + 3475097839005331995329307841
[MAGMAによる計算]
> fd:=C0(81073841058447585744400*x^4 - 2761749046308083386366080*x^3 - 156186525536754891301504376*x^2 - 5717778881805\
70929106834464*x + 3475097839005331995329307841);
Time: 3539.516
4
> fd;
[
    Curve over Rational Field defined by
    71*$.1*$.2 + 404*$.1*$.3 + 248*$.1*$.4 + 105*$.2^2 - 5*$.2*$.3 + 891*$.2*$.4 - 375*$.3^2 + 229*$.3*$.4 - 1292*$.4^2,
    -481202*$.1*$.2 - 150354*$.1*$.3 - 224182*$.1*$.4 + 58158*$.2^2 - 475491*$.2*$.3 - 140944*$.2*$.4 + 224488*$.3^2 -
        672316*$.3*$.4 - 391634*$.4^2,
    Curve over Rational Field defined by
    129*$.1*$.2 + 311*$.1*$.3 + 329*$.1*$.4 + 193*$.2^2 + 374*$.2*$.3 - 195*$.2*$.4 - 362*$.3^2 + 1220*$.3*$.4 +
        1792*$.4^2,
    87077*$.1^2 - 41541*$.1*$.2 + 322698*$.1*$.3 - 582446*$.1*$.4 + 95888*$.2^2 + 187264*$.2*$.3 + 59538*$.2*$.4 -
        905156*$.3^2 + 401062*$.3*$.4 - 445800*$.4^2,
    Curve over Rational Field defined by
    11*$.1^2 + 313*$.1*$.2 + 239*$.1*$.3 + 203*$.1*$.4 + 116*$.2^2 - 579*$.2*$.3 - 213*$.2*$.4 - 180*$.3^2 - 215*$.3*$.4
        - 1360*$.4^2,
    344431*$.1^2 - 122162*$.1*$.2 + 238583*$.1*$.3 - 116085*$.1*$.4 - 246957*$.2^2 + 360220*$.2*$.3 + 129079*$.2*$.4 -
        223442*$.3^2 + 441363*$.3*$.4 - 236200*$.4^2,
    Curve over Rational Field defined by
    10*$.1^2 + 857*$.1*$.2 + 322*$.1*$.3 + 557*$.1*$.4 + 137*$.2^2 + 51*$.2*$.3 - 95*$.2*$.4 - 156*$.3^2 - 18*$.3*$.4 +
        303*$.4^2,
    200414*$.1^2 + 164292*$.1*$.2 - 860818*$.1*$.3 - 20185*$.1*$.4 + 67483*$.2^2 + 713177*$.2*$.3 - 282450*$.2*$.4 -
        304980*$.3^2 + 343708*$.3*$.4 + 104441*$.4^2
]
> RP4(fd,10^6);
J=1
Elliptic Curve defined by y^2 = x^3 + x^2 - 50004566528123296237180656993911520*x -
492831845561205090049555477799439776999549996781132 over Rational Field
rootno=-1
(2532765438480553990368763746/7775888761 : 91421630249318460389985117263764565135040/685685646833741 : 1)
height 16.0326980313085706081149797115
true (2532765438480553990368763746/7775888761 : 91421630249318460389985117263764565135040/685685646833741 : 1)
(2532765438480553990368763746/7775888761 : 91421630249318460389985117263764565135040/685685646833741 : 1)
height 16.0326980313085706081149797115
true (2532765438480553990368763746/7775888761 : 91421630249318460389985117263764565135040/685685646833741 : 1)
(2532765438480553990368763746/7775888761 : 91421630249318460389985117263764565135040/685685646833741 : 1)
height 16.0326980313085706081149797115
true (2532765438480553990368763746/7775888761 : 91421630249318460389985117263764565135040/685685646833741 : 1)
(2532765438480553990368763746/7775888761 : 91421630249318460389985117263764565135040/685685646833741 : 1)
height 16.0326980313085706081149797115
true (2532765438480553990368763746/7775888761 : 91421630249318460389985117263764565135040/685685646833741 : 1)
(2532765438480553990368763746/7775888761 : 91421630249318460389985117263764565135040/685685646833741 : 1)
height 16.0326980313085706081149797115
true (2532765438480553990368763746/7775888761 : 91421630249318460389985117263764565135040/685685646833741 : 1)
(2532765438480553990368763746/7775888761 : 91421630249318460389985117263764565135040/685685646833741 : 1)
height 16.0326980313085706081149797115
true (2532765438480553990368763746/7775888761 : 91421630249318460389985117263764565135040/685685646833741 : 1)
J=2
Elliptic Curve defined by y^2 = x^3 + x^2 - 50004566528123296237180656993911520*x -
492831845561205090049555477799439776999549996781132 over Rational Field
rootno=-1
J=3
Elliptic Curve defined by y^2 = x^3 + x^2 - 50004566528123296237180656993911520*x -
492831845561205090049555477799439776999549996781132 over Rational Field
rootno=-1
J=4
Elliptic Curve defined by y^2 = x^3 + x^2 - 50004566528123296237180656993911520*x -
492831845561205090049555477799439776999549996781132 over Rational Field
rootno=-1
realtime=30.910
4
> RP4(fd,10^8);
J=1
Elliptic Curve defined by y^2 = x^3 + x^2 - 50004566528123296237180656993911520*x -
492831845561205090049555477799439776999549996781132 over Rational Field
rootno=-1
(2532765438480553990368763746/7775888761 : 91421630249318460389985117263764565135040/685685646833741 : 1)
height 16.0326980313085706081149797115
true (2532765438480553990368763746/7775888761 : 91421630249318460389985117263764565135040/685685646833741 : 1)
(2532765438480553990368763746/7775888761 : 91421630249318460389985117263764565135040/685685646833741 : 1)
height 16.0326980313085706081149797115
true (2532765438480553990368763746/7775888761 : 91421630249318460389985117263764565135040/685685646833741 : 1)
(2532765438480553990368763746/7775888761 : 91421630249318460389985117263764565135040/685685646833741 : 1)
height 16.0326980313085706081149797115
true (2532765438480553990368763746/7775888761 : 91421630249318460389985117263764565135040/685685646833741 : 1)
(2532765438480553990368763746/7775888761 : 91421630249318460389985117263764565135040/685685646833741 : 1)
height 16.0326980313085706081149797115
true (2532765438480553990368763746/7775888761 : 91421630249318460389985117263764565135040/685685646833741 : 1)
(2532765438480553990368763746/7775888761 : 91421630249318460389985117263764565135040/685685646833741 : 1)
height 16.0326980313085706081149797115
true (2532765438480553990368763746/7775888761 : 91421630249318460389985117263764565135040/685685646833741 : 1)
(2532765438480553990368763746/7775888761 : 91421630249318460389985117263764565135040/685685646833741 : 1)
height 16.0326980313085706081149797115
true (2532765438480553990368763746/7775888761 : 91421630249318460389985117263764565135040/685685646833741 : 1)
(166849303893134883148437788160947636707966299408185651739623209284422929706808417426/226614460524646931384238701527448\
7132552862857777172552561183441 : -681529400035945165194051703966839797468301413245126713217186594874618139009332499451\
26293959402016207453208468574863913018560/10787763797991049584669397588220231747613493485505683322241825977481949739872\
7122871136773617511 : 1)
height 144.294282281777135473034817403
true (166849303893134883148437788160947636707966299408185651739623209284422929706808417426/2266144605246469313842387015\
274487132552862857777172552561183441 : 68152940003594516519405170396683979746830141324512671321718659487461813900933249\
945126293959402016207453208468574863913018560/1078776379799104958466939758822023174761349348550568332224182597748194973\
98727122871136773617511 : 1)
(2532765438480553990368763746/7775888761 : 91421630249318460389985117263764565135040/685685646833741 : 1)
height 16.0326980313085706081149797115
true (2532765438480553990368763746/7775888761 : 91421630249318460389985117263764565135040/685685646833741 : 1)
J=2
Elliptic Curve defined by y^2 = x^3 + x^2 - 50004566528123296237180656993911520*x -
492831845561205090049555477799439776999549996781132 over Rational Field
rootno=-1
J=3
Elliptic Curve defined by y^2 = x^3 + x^2 - 50004566528123296237180656993911520*x -
492831845561205090049555477799439776999549996781132 over Rational Field
rootno=-1
J=4
Elliptic Curve defined by y^2 = x^3 + x^2 - 50004566528123296237180656993911520*x -
492831845561205090049555477799439776999549996781132 over Rational Field
rootno=-1
realtime=319.424
4
> RP4(fd[2..4],10^10);
J=1
Elliptic Curve defined by y^2 = x^3 + x^2 - 50004566528123296237180656993911520*x -
492831845561205090049555477799439776999549996781132 over Rational Field
rootno=-1
J=2
Elliptic Curve defined by y^2 = x^3 + x^2 - 50004566528123296237180656993911520*x -
492831845561205090049555477799439776999549996781132 over Rational Field
rootno=-1
J=3
Elliptic Curve defined by y^2 = x^3 + x^2 - 50004566528123296237180656993911520*x -
492831845561205090049555477799439776999549996781132 over Rational Field
rootno=-1
realtime=6699.051
3
>
[Pari/GPによる計算]
gp > v=V(81073841058447585744400*x^4 - 2761749046308083386366080*x^3 - 156186525536754891301504376*x^2 - 571777888180570929106834464*x + 3475097839005331995329307841)
%666 = [81073841058447585744400, -2761749046308083386366080, -156186525536754891301504376, -571777888180570929106834464, 3475097839005331995329307841]
gp > e=E2(v)
v=[81073841058447585744400, -2761749046308083386366080, -156186525536754891301504376, -571777888180570929106834464, 3475097839005331995329307841]
I=23037784005440038223750867376553214137510594455490816
J=800805778113716386787325354887091273742056329377976827363953824978997520310272
rr=[59388, 1175644848, 0, 0]
time = 3 ms.
%667 = [0, 1, 0, -50004566528123296237180656993911520, -492831845561205090049555477799439776999549996781132, 4, -100009133056246592474361313987823040, -1971327382244820360198221911197759107998199987124528, -2500456673665508726544847043820517454593533712801108734972549575834928, 2400219193349918219384671535707752976, 425806714564881183401500772719206651019581982972380224, 7897266579200114239146217647740081315390195190133026449066314040685865911481424587758954841342557288960000, 3375924806402913965471924711927732647121270807223866051030648984925593772962490285768541065154421443200481/1928043598437527890416557042905293289890184372591070910416580576339322732295269674745838584312147775625, Vecsmall([1]), [Vecsmall([128, 1])], [0, 0, 0, 0, 0, [28113332979347621548284611760, 4194304, [2, 4; 3, 1; 5, 1; 11, 1; 17, 1; 139, 1; 409, 1; 1069, 1; 9497, 1; 20011, 1; 54236099, 1], [[4, -8, 0, 4], [1, 12, 0, 8], [1, 8, 0, 4], [1, 8, 0, 4], [1, 8, 0, 4], [1, 8, 0, 4], [1, 8, 0, 4], [1, 8, 0, 4], [1, 8, 0, 4], [1, 6, 0, 2], [1, 8, 0, 4]]], 0, [[2]~]]]
gp > Q=chpi(v,2532765438480553990368763746/7775888761)
v=[81073841058447585744400, -2761749046308083386366080, -156186525536754891301504376, -571777888180570929106834464, 3475097839005331995329307841]
I=23037784005440038223750867376553214137510594455490816
J=800805778113716386787325354887091273742056329377976827363953824978997520310272
rr=[59388, 1175644848, 0, 0]
v=[81073841058447585744400, -2761749046308083386366080, -156186525536754891301504376, -571777888180570929106834464, 3475097839005331995329307841]
I=23037784005440038223750867376553214137510594455490816
J=800805778113716386787325354887091273742056329377976827363953824978997520310272
time = 5 ms.
%668 = [8932897916826372750030319714832995152/7775888761, 19148954226960385258811376759712961508868794761602682880/685685646833741]
gp > L=ss1x(v,Q,40);
v=[81073841058447585744400, -2761749046308083386366080, -156186525536754891301504376, -571777888180570929106834464, 3475097839005331995329307841]
I=23037784005440038223750867376553214137510594455490816
J=800805778113716386787325354887091273742056329377976827363953824978997520310272
time = 1,265 ms.
gp > L[1..10]
%670 = [0, -16053/3032, -506731269022/12958703985, 223879820934094/3421354908729, -189236734896984828/24182697107313025, -20026629864887984305/756946939587939312, 2287213080046035824661/722903464199770844120, -1790289524054760751625436165470569932990866/50721081126660309931935990779467021903389, 102158436166424744281737173937044128511849748/152069345226187023246554322942823320845913275, 125934525711764207290919313898985281104373681355/408633744665698977126948695748176514047398992]
gp > sss4t(20011, u0, L[1..40])
[u0,x0,t0]=[3401870/9497, -4548/4949, 48159/9898]
48159^4+20011*9898^4=88181^4+20011*9096^4
47604039799647984432639087705271^4+20011*72737567159283845853998286455544^4=713950061055274851217747449605019^4+20011*62242121810044958270783231103278^4
21357324022572449761716682854703243050216860514479900645249684733487391929183583022657501^4+20011*3487367946599870038117639593139506820112723155621459760881098426846045288827009553377918^4=29024833603646511220329349574616368518363083700040292800637525977869637478646877976728609^4+20011*3329156904202553375387329714391431735670906428518372768567164774762941252959673717705464^4
4632128158847525255194044182804884954732262716783709108110666469019291209901558571725326734497872543463807506381460908791344908900986361873654630538696337325898609815405079^4+20011*1294375062956701887872148008553020182696245423115677891863247704610076639439233893123793364112968116078742409476120903638852835943300480003762377258332458120709962381624058^4=12094214291285880777738192051089320747745979893930736888445173639796667504308102093202197001472041831287397198079062943731101026504964777451780572494140460311662461157594339^4+20011*1151939768366086435716487392773312594780132256049479892707490156725601077228259064388667775688718704147245132525253583555284570131885157841722842582378477864601720597812616^4
...省略...
3137789139953206196209420812716326453408050498397265030195571403507639609303612211831264804766337421239375675328140628975266152377570266971399503921306099040006948538113357440624190036592590603362015066166037089641147158243772219941315587266329341808565482698534752934158423332688974391219578844987057450497828547064501825767086989773729487808517748759419589930932709391104657635498927595907868726275688342585608779726971925251751005091949756538029715327258522362712848764217739685835097498452686089979370798949344526179035527290912057064081482957398254902450383303036159213610471678734601929917803546729285081664184934862643306000154335923068390969634991695965057638679270843750818203419640825289203937515960090180916316737590767420161699098574780132735947768450825432494979845681479483211045727679405354176198381457334907519969196487992997073029464154025338226373701847185973476132353940173693545130468563049384586826107382662910714704909483607407972378143672372581197563554984799034896022748446953024883003302715558322060158745630167706333056586973037546444696950715401193161315243007618583565511649543733876551185690402454204960969345260811406584730753148751773329309147810940300924912907537031998682544132621728318107976846620956106454901731364542539820412261903584847988834313233481953931808880439418143009625094777461801046069255193110349001339838358359400125657452355061773842293357644053585123642023505854974547064704298588241258614761504265315823760835750514916330269619417820198820685547788399283463240464438490538223664802769^4+20011*688641088807358184532991360463526527328311326574089892350083900111712090898962673236552831590607179197292269046621301201670054274602212851441212202661886487951125654124890355256850224339530155556729776498146159856004244979247867045528059952430006423046013635134851716092765474195617207893877722234033016680422250337911233106821943327135560820615066305985596614104262339360561517523656336492444178501671139702596758818550605732471078694671665453608698077887902721663629651717516780395773547727427006366307140792812337151107526514870919164302501571167765850047284506899495640761472154398304462094129111945423794706708356283589456171067135034402884766673511606479261875252867114677756203161495231733999073679830963669642494486820440854562378209672109261935231195850525183568666332698676143310217990058411578219086205026464656138828548027023355290647186058600015480464732996412444519377742186313775729939518588375782102708189657642599721687730950555583348566691201492189187557631695788063818430620894833829364198198891339381725599249588073652166919833666424262444462276154541079714462266601609940122675605848422611393380415053241056738206675769745098433922601957289111240270171567422979627737496638982318226502216061622812500842979533316151040166370830762987408454120774335995091959620982193829598413256917979321415999262723505750308039314726992611446371330012890728571359586951099143957815649671791996274783403828362607318359360370043254746902866434421644215832518040528124006672572082632017406277101683141836334793818579397269910170195704^4=6216709919033421826367909752420051570612135814591607485994501681584816873282668190191473932049449478307653075087674429126471748624037991682190197905227255129334871576645240956313470902950121079799106849025364945716199620080179890776130057924555421086674390051046630828064375062187254285917856691032048928267915487626018188698993269056707593150032975303638765786052498967498153754357691514962082677223520496241399064531533854900717606604177220522800508978915751837915147604195019689215572620863624607033110967729360095909347131731529936113444472360692198423178643351054691448186297622915482152309232132271849014817967354554993196230912368320379140118824737349396069091713707531196417650153519586310171651206318394175221291101621988965825382243383975875100182375891027968085495185115822839307034853222810244889740330073957504418087906234332007337854574692003158313267288853383579773823287299809076211722204142303207450251738415206820516030157270956402503096965313452496254533252606468283618041160526161069852206492251831427046346960600762549506996044348774005616664126891815093154161639997765160975037518931642191911089556052525289479921400185006404208979295704699457376377653052219578199397692734524015751704171892937822847814479514654599510070157144577783261681857953244940559774915522792333864050214572027620144942101027595345011742690335280221312180572501791144545326597137109137820754516354359169796019296389922879365135308184247840589025851331043283498066614135467138384820699991576935119938772946319337195226306279873467591333591859^4+20011*627551333027216589763198260067362751158045522170009439938610323979543240533478950662189211744291272667666822169314005721652062008545944616103054874057730676019371329235330424234653848821307430984406881106702592442757600421564621779664540979653361874933472473928829642017044370080056668631091311440521469653396955494614756114522632566634065498804722308016297053918659864552160563784179157937324684953524993739054510729930358432590100339084022583259133176701991552898466273161340849225208537189997555178291732799177033419916999100825387001043912475323508942689924705223974383092435133082255073954503372421914383293762705272022689137161074985348044709720836421930878422698124889873790069091430727829789566666087693502039801218639625879445045709942896428718651083680731661293861551261139905343042839424428876530079702657255572668997305624814511845207528191410443805286841858232844672925445404699660255399914744557707091453090067856247333911963346011008579827023805296325647951402878609268580229641075492491995274335026376935710126891619659779236301109630327013244877991836268236977395879218076697297420003488458748057847675161325564089432104833126412157031882770231927021758025392504451304303289265292855065292441094290481576586914687787230438461682391819437346778166905531744359571679351975142787767075888813373409215906646650407514871276212118242070692113916687403620267540835251478370741417371735699372792797055726084438116469820637174955928225579024831218061588681434966318360833308639275201354870058492699173207945039457028467993843998^4
time = 9 ms.
%671 = 11
gp > Puxt(20011,3401870/9497, -4548/4949, 48159/9898 )
%672 = [20268460264611896436100*x^4 - 690437261577020846591520*x^3 - 39046631384188722825376094*x^2 - 142944472045142732276708616*x + 3475097839005331995329307841/4, 14685401040*x^2 + 155504447820*x + 3040387614132, 77752223910*x^2 + 6080775228264*x + 32194816866531/2, -15980222020*x^2 + 3308460488641]

[2026.07.28追記]u=3401870/9497のときの計算ログを追加した。

[参考文献]

Last Update: 2026.07.28
H.Nakao

Homeに戻る[Homeに戻る]  一覧に戻る[一覧に戻る]