Homeに戻る  一覧に戻る 

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

[2025.10.13]A^4+B^4+C^4=3362*D^4の整点

■平方因子を持たない正整数nを固定したとき、不定方程式
       A^4+B^4+C^4=2*n^2*D^4 ----------(1)
を満たす自明でない整数の組(A,B,C,D) (ただし C!=0かつgcd(A,B,C,D)=1)を探す。

以下では、Elkiesの論文(参考文献[1])の方法およびTom Womackの文書(参考文献[5])を参考にして、(1)を満たす整数の組(A,B,C,D)を探す。
ここで、整数A,B,C,Dは0以上として良い。


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

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


■2次曲線(3a±)がsingularであるかどうかを議論する。
(3a±)の右辺の２次式の判別式D1は、
      D1=4*(u^2-2*u+2)^2-4*(-u^2+4*u-2)^2=4*u*(u-1)*(u-2)
であるので、２次曲線(3a±)がsingular(２本の直線または1点)になるのは、u=0,1,2のとき、また、そのときに限る。

(i)u=0のとき、(3a±)は、
       ±y^2=(x-1)^2
となり、左辺の復号は+のみ成立するとして良い。よって、
       y=±(x-1)
となるので、(2)より、
       x^4 - 2*x^3 + 3*x^2 - 2*x + 1=n^2*t^4 -----(*)
       (x^2 - x + 1)^2=(n*t^2)^2
       x^2 - x + 1=±(n*t^2)
n > 0 かつ x^2-x+1 > 0なので、右辺の復号±のうち、+のみが成立する。よって、
       x^2 - x + 1=n*t^2　--------(**)
となる。
２次曲線(**)が有理点を持つかどうかの議論が必要である。

(ii)u=1のとき、(3a±)は、
       ±(-y^2)=(x+1)^2
となるので、左辺の復号のうち、-のみが成立するとして良い。よって、
       y=±(x+1)
となる。(i)と同様にして、
       x^2 + x + 1=n*t^2　--------(***)
となる。
(***)のxを-xで置き換えると、(**)に一致するので、(***)が有理点を持つことと、(**)が有理点を持つことは同値である。
(i)と同じく、２次曲線(**)が有理点を持つかどうかの議論が必要である。

(iii)u=2のとき、(3a±)は、
       ±y^2=(x+1)^2
となるので、左辺の復号のうち、+のみが成立するとして良い。よって、
       y=±(x+1)
となるので、(2)より、
       x^4 + 2*x^3 + 3*x^2 + 2*x + 1=n^2*t^4
となり、これは、(*)のxを-xで置き換えたものであるので、(*)または(**)に帰着できる。

以下では、u!=0,1,2のとき、つまり、(3a)がnon-singularである場合のみを考察する。


■2次曲線(3b±)がsingularであるかどうかを議論する。
(3b±)の右辺の判別式D2は、
      D2=(-u^2+4*u-2)^2-4*(u^2-2*u+2)^2=-(u^2+2)*(3*u^2-8*u+8)=(u^2+2)*(3*(u-4/3)^2+8/3)
であり、有理数の根を持たない。
よって、任意の有理数uに対して、(3b±)はnon-singularである。

さらに、判別式D2は常に負であり、実数の根を持たない。
よって、(3b±)の右辺は常に正であり、復号±のどちらか１つのみが成立する。
(3b±)の左辺より、u^2 > 2のとき、(3b+)のみ、u^2 < 2のとき、(3b-)のみが成立する。


■２次曲線(3a±),(3b±)は、u,xの変換τ(u,x)=((u-2)/(u-1),-x)に対する対称性を持つ。
ここで、τ(τ(u,v))=τ((u-2)/(u-1),-x))=(u,x)である。つまり、τはinvolutionである。
(3a)をinvolution τで変換すると、不変であり、(3b+),(3b-)をτで変換すると、それぞれ(3b-),(3b+)になる。
よって、有理数uとτ(u)=(u-2)/(u-1)のどちらか一方(例えば、高さが大ききくない方)について、(3a),(3b±)が有理点を持つかどうか調べれば良い。

[pari/gpによる計算]
(02:59) gp > YY2(u,x)-YY2((u-2)/(u-1),-x)
%39 = 0
(07:07) gp > TT2(n,u,x)+TT2(n,(u-2)/(u-1),-x)
%40 = 0
[2025.11.08追記]τと異なるinvoclution σ(u,x)=((2*u-2)/(u-2),-x)および(στ)(u,x)=(2/u,-x)を見つけた。これにより、(3a-),(3b±)の共通有理点は、調べなくてもs良い。

3362=2*41^2であるので、以下では、n=41の場合について、議論する。


■n=41のとき、２次曲線(**)は、有理点を持たないことが確認できる。

{MAGMAでの計算]
> P2< x,y,z > := ProjectiveSpace(Rationals(), 2);
> N:=41;
> C := Conic(P2,-N*y^2+x^2-x*z+z^2);
> HasRationalPoint(C);
false
>

以下のように、証明も可能である。
２次曲線(**)を同次化すると、
       x^2 - x*z + z^2=n*t^2　--------(***)
となる。
以下の２つは同値である。
・平面上の２次曲線(**)は有理点(x,t)を持つ。
・射影平面上の２次曲線(***)はz!=0である有理点[x:z:t]を持つ。

n=41のとき、整数環Z/41Z上で、(***)の整数解はx=z=0(mod 41),t(mod 41)は任意,に限る。
整数x,zはどちらも素数41で割り切れるので、(***)の左辺は41^2で割り切れる。よって、(***)の右辺も41^2で割り切れるので、整数tは素数41で割り切れる。
x,y,tはいずれも41で割り切れるので、(***)はより小さい整数解(x/41,z/41,t/41)を持つ。
無限降下法により、(***)はz!=0である整数解を持たない。
よって、２次曲線(***)はz!=0である有理点を持たない。つまり、２次曲線(**)は有理点を持たない。

■有理数u(u!=0,1,2)の高さが小さいものから、順に調べる。
さらに、(3a±),(3b±)を共に満たす有理数の組(x,y,t)が存在すれば、2次曲線(3a+)と(3a-)の少なくとも一方が有理点を持ち、かつ、(3b+)と(3b-)の少なくとも一方が有理点を持つ。

例えば、有理数uの高さが200以下の範囲で、２次曲線(3a+)と２つの２次曲線の和集合(3b±)が共に有理点を持つようなuを選択すると、以下のように144個のuが抽出される。
これらのuについて、(3a+),(3b±)を共に満たす有理数の組(x,y,t)を見つければ良い。

[MAGMAによる計算]
> PP(41,1,200);
** u= -1 ; tau(u)= 3/2 ; -7*x^2 + y^2 + 10*x*z - 7*z^2
  (1 : -2 : 1)  C2b (3/4 : -1/4 : 1)
** u= 1/5 ; tau(u)= 9/4 ; -31*x^2 + 49*y^2 + 82*x*z - 31*z^2
  (11/25 : -24/175 : 1)  C1b (-5880/4223 : -8437/29561 : 1)
** u= 1/17 ; tau(u)= 33/16 ; -511*x^2 + 577*y^2 + 1090*x*z - 511*z^2
  (347/237 : -32/237 : 1)  C1b (10929/8632 : -1559/8632 : 1)
** u= -1/49 ; tau(u)= 99/50 ; -4999*x^2 + 4801*y^2 + 9802*x*z - 4999*z^2
  (-353/543 : 910/543 : 1)  C2b (13955/35901 : 4921/35901 : 1)
** u= 1/53 ; tau(u)= 105/52 ; -5407*x^2 + 5617*y^2 + 11026*x*z - 5407*z^2
  (-523/244295 : -240208/244295 : 1)  C1b (16240/27893 : -155059/1143613 : 1)
** u= 3/2 ; tau(u)= -1 ; 7*x^2 - y^2 + 10*x*z + 7*z^2
  (-7/3 : -14/3 : 1)  C1a (13/21 : -11/21 : 1)
** u= 4/17 ; tau(u)= 30/13 ; -322*x^2 + 562*y^2 + 916*x*z - 322*z^2
  (9/22 : -1/22 : 1)  C1b (-565/548 : -129/548 : 1)
** u= -4/45 ; tau(u)= 94/49 ; -4786*x^2 + 4034*y^2 + 8852*x*z - 4786*z^2
  (395/1376 : -1113/1376 : 1)  C2b (-562715/78501 : -99583/78501 : 1)
** u= -4/49 ; tau(u)= 102/53 ; -5602*x^2 + 4786*y^2 + 10420*x*z - 5602*z^2
  (88/81 : -35/81 : 1)  C2b (-45761/5087 : -7933/5087 : 1)
** u= 4/173 ; tau(u)= 342/169 ; -57106*x^2 + 59842*y^2 + 116980*x*z - 57106*z^2
  (45732/59503 : 7423/59503 : 1)  C1b (1403257/183156 : 204353/183156 : 1)
** u= 7/9 ; tau(u)= 11/2 ; 41*x^2 + 113*y^2 + 170*x*z + 41*z^2
  (-137/187 : -138/187 : 1)  C1b (61139/47837 : -11681/47837 : 1)
** u= 8/9 ; tau(u)= 10 ; 62*x^2 + 98*y^2 + 164*x*z + 62*z^2
  (-7/5 : 24/35 : 1)  C1b (0 : -1/7 : 1)
** u= 8/25 ; tau(u)= 42/17 ; -514*x^2 + 1186*y^2 + 1828*x*z - 514*z^2
  (-151/14 : 115/14 : 1)  C1b (-126120/1223 : -17339/1223 : 1)
** u= -8/45 ; tau(u)= 98/53 ; -5554*x^2 + 3986*y^2 + 9668*x*z - 5554*z^2
  (172/241 : -147/241 : 1)  C2b (111440/153 : 19151/153 : 1)
** u= 8/81 ; tau(u)= 154/73 ; -10594*x^2 + 13058*y^2 + 23780*x*z - 10594*z^2
  (1973/7193 : -4392/7193 : 1)  C1b (-252817/353256 : -77189/353256 : 1)
** u= -8/101 ; tau(u)= 210/109 ; -23698*x^2 + 20338*y^2 + 44164*x*z - 23698*z^2
  (401/2841 : 2668/2841 : 1)  C2b (10152968/433235 : -1615479/433235 : 1)
** u= 9/4 ; tau(u)= 1/5 ; 31*x^2 - 49*y^2 + 82*x*z + 31*z^2
  (1 : -12/7 : 1)  C1a (0 : -1/7 : 1)
** u= 10 ; tau(u)= 8/9 ; -62*x^2 - 98*y^2 + 164*x*z - 62*z^2
  (7/5 : -24/35 : 1)  C1a (-5880/4223 : -8437/29561 : 1)
** u= 11/2 ; tau(u)= 7/9 ; -41*x^2 - 113*y^2 + 170*x*z - 41*z^2
  (17/43 : 18/43 : 1)  C1a (484/215 : 13/43 : 1)
** u= -11/149 ; tau(u)= 309/160 ; -51079*x^2 + 44281*y^2 + 95602*x*z - 51079*z^2
  (-16633/5217 : 23192/5217 : 1)  C2b (21917447/134753 : -3544743/134753 : 1)
** u= 12/109 ; tau(u)= 206/97 ; -18674*x^2 + 23618*y^2 + 42580*x*z - 18674*z^2
  (-67/12586 : 78815/88102 : 1)  C1b (-30332/24551 : -48417/171857 : 1)
** u= -12/185 ; tau(u)= 382/197 ; -77474*x^2 + 68306*y^2 + 146068*x*z - 77474*z^2
  (311/402 : -1121/2814 : 1)  C2b (182585/1948 : 8414757/559076 : 1)
** u= 13/53 ; tau(u)= 93/40 ; -3031*x^2 + 5449*y^2 + 8818*x*z - 3031*z^2
  (-1707/4423 : -4972/4423 : 1)  C1b (64656/7135 : -8789/7135 : 1)
** u= 15/17 ; tau(u)= 19/2 ; 217*x^2 + 353*y^2 + 586*x*z + 217*z^2
  (-2563/5209 : -1202/5209 : 1)  C1b (-484/361 : -69/361 : 1)
** u= 15/113 ; tau(u)= 211/98 ; -18983*x^2 + 25313*y^2 + 44746*x*z - 18983*z^2
  (78759773/1076725509 : 850981054/1076725509 : 1)  C1b (907747/368117 : 120717/368117 : 1)
** u= -16/97 ; tau(u)= 210/113 ; -25282*x^2 + 18562*y^2 + 44356*x*z - 25282*z^2
  (-1115/271 : -1586/271 : 1)  C2b (-292923/8533 : 50849/8533 : 1)
** u= 17/81 ; tau(u)= 145/64 ; -7903*x^2 + 12833*y^2 + 21314*x*z - 7903*z^2
  (79/1393 : -144/199 : 1)  C1b (436480/222199 : -2391413/9110159 : 1)
** u= 19/2 ; tau(u)= 15/17 ; -217*x^2 - 353*y^2 + 586*x*z - 217*z^2
  (5/3 : -2/3 : 1)  C1a (-3460/8471 : 1461/8471 : 1)
** u= 19/149 ; tau(u)= 279/130 ; -33439*x^2 + 44041*y^2 + 78202*x*z - 33439*z^2
  (7299/4079 : 466/4079 : 1)  C1b (17471/64225 : -8669/64225 : 1)
** u= -19/193 ; tau(u)= 405/212 ; -89527*x^2 + 74137*y^2 + 164386*x*z - 89527*z^2
  (269/163 : -1044/1141 : 1)  C2b (-90224/13185 : 112841/92295 : 1)
** u= -20/81 ; tau(u)= 182/101 ; -20002*x^2 + 12722*y^2 + 33524*x*z - 20002*z^2
  (-106/677 : -963/677 : 1)  C2b (895023/107623 : 149737/107623 : 1)
** u= 20/109 ; tau(u)= 198/89 ; -15442*x^2 + 23362*y^2 + 39604*x*z - 15442*z^2
  (-899/770 : -207/110 : 1)  C1b (-227215/49164 : -36019/49164 : 1)
** u= 21/37 ; tau(u)= 53/16 ; -71*x^2 + 2297*y^2 + 3250*x*z - 71*z^2
  (9/431 : 16/431 : 1)  C1b (32384/16693 : 4743/16693 : 1)
** u= 21/61 ; tau(u)= 101/40 ; -2759*x^2 + 7001*y^2 + 10642*x*z - 2759*z^2
  (-987/26695 : -17924/26695 : 1)  C1b (-4800680/310193 : 665823/310193 : 1)
** u= -28/153 ; tau(u)= 334/181 ; -64738*x^2 + 46034*y^2 + 112340*x*z - 64738*z^2
  (18401/5756 : -16257/5756 : 1)  C2b (37573/47109 : 6881/47109 : 1)
** u= 29/49 ; tau(u)= 69/20 ; 41*x^2 + 3961*y^2 + 5602*x*z + 41*z^2
  (-95/12687 : 196/12687 : 1)  C1b (-23933/12693 : -3547/12693 : 1)
** u= 30/13 ; tau(u)= 4/17 ; 322*x^2 - 562*y^2 + 916*x*z + 322*z^2
  (-324/107 : 101/107 : 1)  C1a (-3101/524 : 417/524 : 1)
** u= -32/53 ; tau(u)= 138/85 ; -13426*x^2 + 4594*y^2 + 20068*x*z - 13426*z^2
  (-6585/28567 : -8248/4081 : 1)  C2b (-2550880/334331 : 642951/334331 : 1)
** u= -32/117 ; tau(u)= 266/149 ; -43378*x^2 + 26354*y^2 + 71780*x*z - 43378*z^2
  (4453/4972 : 3609/4972 : 1)  C2b (-17636703/234536 : -3240329/234536 : 1)
** u= 33/16 ; tau(u)= 1/17 ; 511*x^2 - 577*y^2 + 1090*x*z + 511*z^2
  (87/47 : 128/47 : 1)  C1a (-240816/6653 : -36089/6653 : 1)
** u= -37/181 ; tau(u)= 399/218 ; -93679*x^2 + 64153*y^2 + 160570*x*z - 93679*z^2
  (-13631/8029 : 25286/8029 : 1)  C2b (-23385180/5062633 : -4657375/5062633 : 1)
** u= 40/113 ; tau(u)= 186/73 ; -9058*x^2 + 23938*y^2 + 36196*x*z - 9058*z^2
  (1085/8289 : 3584/8289 : 1)  C1b (167672/14257 : 22363/14257 : 1)
** u= 42/17 ; tau(u)= 8/25 ; 514*x^2 - 1186*y^2 + 1828*x*z + 514*z^2
  (-157/721 : 248/721 : 1)  C1a (306785/4327 : -42231/4327 : 1)
** u= -43/153 ; tau(u)= 349/196 ; -74983*x^2 + 44969*y^2 + 123650*x*z - 74983*z^2
  (-46909/149 : -60732/149 : 1)  C2b (9266448/237469 : 1670851/237469 : 1)
** u= 44/169 ; tau(u)= 294/125 ; -29314*x^2 + 55186*y^2 + 88372*x*z - 29314*z^2
  (-3998/253 : 3185/253 : 1)  C1b (-17148/13823 : -145217/566743 : 1)
** u= 48/193 ; tau(u)= 338/145 ; -39746*x^2 + 72194*y^2 + 116548*x*z - 39746*z^2
  (15393/1559 : -9646/1559 : 1)  C1b (7625875/3031328 : -1007673/3031328 : 1)
** u= -51/157 ; tau(u)= 365/208 ; -83927*x^2 + 46697*y^2 + 135826*x*z - 83927*z^2
  (375/4969 : 43832/34783 : 1)  C2b (4413779/3742165 : -4766859/26195155 : 1)
** u= -52/73 ; tau(u)= 198/125 ; -28546*x^2 + 7954*y^2 + 41908*x*z - 28546*z^2
  (602/843 : -1085/843 : 1)  C2b (-20/321 : -3481/13161 : 1)
** u= 52/137 ; tau(u)= 222/85 ; -11746*x^2 + 34834*y^2 + 51988*x*z - 11746*z^2
  (391/4038 : -1787/4038 : 1)  C1b (582631/34604 : -77681/34604 : 1)
** u= 53/16 ; tau(u)= 21/37 ; 71*x^2 - 2297*y^2 + 3250*x*z + 71*z^2
  (-43/4101 : -520/4101 : 1)  C1a (4414352/46499 : 580023/46499 : 1)
** u= -55/49 ; tau(u)= 153/104 ; -18607*x^2 + 1777*y^2 + 26434*x*z - 18607*z^2
  (-1/181 : 588/181 : 1)  C2b (-395328/67435 : -189697/67435 : 1)
** u= -56/89 ; tau(u)= 234/145 ; -38914*x^2 + 12706*y^2 + 57892*x*z - 38914*z^2
  (24903/68491 : 92168/68491 : 1)  C2b (515895/53792 : 113473/53792 : 1)
** u= 56/145 ; tau(u)= 234/89 ; -12706*x^2 + 38914*y^2 + 57892*x*z - 12706*z^2
  (-4351/433281 : 253196/433281 : 1)  C1b (223593/135833 : 31621/135833 : 1)
** u= 57/109 ; tau(u)= 161/52 ; -2159*x^2 + 20513*y^2 + 29170*x*z - 2159*z^2
  (-21933/71317 : 53012/71317 : 1)  C1b (2023/3641 : 531/3641 : 1)
** u= 61/113 ; tau(u)= 165/52 ; -1687*x^2 + 21817*y^2 + 30946*x*z - 1687*z^2
  (-38355/358771 : -24572/51253 : 1)  C1b (-65575/168047 : -24159/168047 : 1)
** u= -61/117 ; tau(u)= 295/178 ; -59647*x^2 + 23657*y^2 + 90746*x*z - 59647*z^2
  (163435/1673839 : -2466186/1673839 : 1)  C2b (21868/11365 : 148559/465965 : 1)
** u= -64/49 ; tau(u)= 162/113 ; -21442*x^2 + 706*y^2 + 30340*x*z - 21442*z^2
  (-19/86 : 553/86 : 1)  C2b (934495/34727 : 658985/34727 : 1)
** u= 64/153 ; tau(u)= 242/89 ; -11746*x^2 + 42722*y^2 + 62660*x*z - 11746*z^2
  (3797/25331 : -6270/25331 : 1)  C1b (2019/1232 : -11839/50512 : 1)
** u= -64/157 ; tau(u)= 378/221 ; -93586*x^2 + 45202*y^2 + 146980*x*z - 93586*z^2
  (-122/49 : -235/49 : 1)  C2b (-45483/216101 : 49297/216101 : 1)
** u= -64/169 ; tau(u)= 402/233 ; -104482*x^2 + 53026*y^2 + 165700*x*z - 104482*z^2
  (81143/30651 : -4940/1803 : 1)  C2b (9040925/650456 : -1687725/650456 : 1)
** u= -67/61 ; tau(u)= 189/128 ; -28279*x^2 + 2953*y^2 + 40210*x*z - 28279*z^2
  (-271/511 : 2256/511 : 1)  C2b (954776/14113 : 385049/14113 : 1)
** u= -67/81 ; tau(u)= 229/148 ; -39319*x^2 + 8633*y^2 + 56930*x*z - 39319*z^2
  (1073/11243 : 22392/11243 : 1)  C2b (709585/385336 : 148525/385336 : 1)
** u= 67/101 ; tau(u)= 135/34 ; 2177*x^2 + 15913*y^2 + 22714*x*z + 2177*z^2
  (-428387/713321 : 83478/101903 : 1)  C1b (17747/26773 : -4421/26773 : 1)
** u= 69/20 ; tau(u)= 29/49 ; -41*x^2 - 3961*y^2 + 5602*x*z - 41*z^2
  (107/2203 : 532/2203 : 1)  C1a (52800/1937 : 6937/1937 : 1)
** u= 76/113 ; tau(u)= 150/37 ; 3038*x^2 + 19762*y^2 + 28276*x*z + 3038*z^2
  (-899/8232 : 31/1176 : 1)  C1b (60/721 : -3953/29561 : 1)
** u= -76/153 ; tau(u)= 382/229 ; -99106*x^2 + 41042*y^2 + 151700*x*z - 99106*z^2
  (113/134 : 135/134 : 1)  C2b (288300/293027 : 51395/293027 : 1)
** u= -76/193 ; tau(u)= 462/269 ; -138946*x^2 + 68722*y^2 + 219220*x*z - 138946*z^2
  (-261804/3787507 : -5683805/3787507 : 1)  C2b (-2404821/42515 : 95935/8503 : 1)
** u= 76/197 ; tau(u)= 318/121 ; -23506*x^2 + 71842*y^2 + 106900*x*z - 23506*z^2
  (-81883/626328 : 454817/626328 : 1)  C1b (-2664725/615316 : -385755/615316 : 1)
** u= 77/109 ; tau(u)= 141/32 ; 3881*x^2 + 17833*y^2 + 25810*x*z + 3881*z^2
  (-14779/36579 : 21064/36579 : 1)  C1b (-353656/34453 : 46517/34453 : 1)
** u= -79/153 ; tau(u)= 385/232 ; -101407*x^2 + 40577*y^2 + 154466*x*z - 101407*z^2
  (-775/18251 : 29796/18251 : 1)  C2b (2013129/658295 : -357043/658295 : 1)
** u= -80/149 ; tau(u)= 378/229 ; -98482*x^2 + 38002*y^2 + 149284*x*z - 98482*z^2
  (56954/352381 : -501357/352381 : 1)  C2b (1583149/1458800 : 275759/1458800 : 1)
** u= -84/65 ; tau(u)= 214/149 ; -37346*x^2 + 1394*y^2 + 52852*x*z - 37346*z^2
  (-4968/6581 : -55327/6581 : 1)  C2b (8969/964 : 231837/39524 : 1)
** u= -84/101 ; tau(u)= 286/185 ; -61394*x^2 + 13346*y^2 + 88852*x*z - 61394*z^2
  (659/280 : -1063/280 : 1)  C2b (185669/71476 : -41529/71476 : 1)
** u= 88/117 ; tau(u)= 146/29 ; 6062*x^2 + 19634*y^2 + 29060*x*z + 6062*z^2
  (-40897/67007 : -46404/67007 : 1)  C1b (-7891593/386951 : 1050619/386951 : 1)
** u= -88/145 ; tau(u)= 378/233 ; -100834*x^2 + 34306*y^2 + 150628*x*z - 100834*z^2
  (-885/9682 : -17761/9682 : 1)  C2b (-4183600/1292103 : -1189159/1292103 : 1)
** u= 89/121 ; tau(u)= 153/32 ; 5873*x^2 + 21361*y^2 + 31330*x*z + 5873*z^2
  (-4161/6553 : 4840/6553 : 1)  C1b (-62755/6816 : 339095/279456 : 1)
** u= 91/149 ; tau(u)= 207/58 ; 1553*x^2 + 36121*y^2 + 51130*x*z + 1553*z^2
  (-659/3361 : 1622/3361 : 1)  C1b (-16924/4257 : -93329/174537 : 1)
** u= -91/181 ; tau(u)= 453/272 ; -139687*x^2 + 57241*y^2 + 213490*x*z - 139687*z^2
  (1166807/2368673 : 2589616/2368673 : 1)  C2b (82815/10664 : 16265/10664 : 1)
** u= 92/157 ; tau(u)= 222/65 ; 14*x^2 + 40834*y^2 + 57748*x*z + 14*z^2
  (-3276/113 : -721/113 : 1)  C1b (704260/387007 : -105543/387007 : 1)
** u= 93/40 ; tau(u)= 13/53 ; 3031*x^2 - 5449*y^2 + 8818*x*z + 3031*z^2
  (-203/1185 : -644/1185 : 1)  C1a (-4289/3296 : 621/3296 : 1)
** u= 94/49 ; tau(u)= -4/45 ; 4786*x^2 - 4034*y^2 + 8852*x*z + 4786*z^2
  (-863/1202 : -567/1202 : 1)  C1a (-68267/23187 : 9659/23187 : 1)
** u= -95/149 ; tau(u)= 393/244 ; -110047*x^2 + 35377*y^2 + 163474*x*z - 110047*z^2
  (11575/13669 : -16336/13669 : 1)  C2b (229720/945593 : 192043/945593 : 1)
** u= 98/53 ; tau(u)= -8/45 ; 5554*x^2 - 3986*y^2 + 9668*x*z + 5554*z^2
  (-3877/2269 : -2604/2269 : 1)  C1a (1477/723 : 341/723 : 1)
** u= 99/50 ; tau(u)= -1/49 ; 4999*x^2 - 4801*y^2 + 9802*x*z + 4999*z^2
  (-4759/5571 : -1330/5571 : 1)  C1a (-136169/32095 : 19399/32095 : 1)
** u= 100/153 ; tau(u)= 206/53 ; 4382*x^2 + 36818*y^2 + 52436*x*z + 4382*z^2
  (-4189/18056 : -8175/18056 : 1)  C1b (-1175/2364 : 13781/96924 : 1)
** u= 101/40 ; tau(u)= 21/61 ; 2759*x^2 - 7001*y^2 + 10642*x*z + 2759*z^2
  (-4345/27411 : 11068/27411 : 1)  C1a (-714935/452992 : -100827/452992 : 1)
** u= -101/117 ; tau(u)= 335/218 ; -84847*x^2 + 17177*y^2 + 122426*x*z - 84847*z^2
  (26171/37907 : 58398/37907 : 1)  C2b (14859121/1525159 : 4082023/1525159 : 1)
** u= 101/121 ; tau(u)= 141/20 ; 9401*x^2 + 19081*y^2 + 30082*x*z + 9401*z^2
  (-33991/86547 : 1144/5091 : 1)  C1b (5617880/832073 : 822357/832073 : 1)
** u= 102/53 ; tau(u)= -4/49 ; 5602*x^2 - 4786*y^2 + 10420*x*z + 5602*z^2
  (-2432/4917 : -3031/4917 : 1)  C1a (-185719/22172 : 28483/22172 : 1)
** u= -104/125 ; tau(u)= 354/229 ; -94066*x^2 + 20434*y^2 + 136132*x*z - 94066*z^2
  (-1577/3783 : 10820/3783 : 1)  C2b (-1219297/173063 : -383289/173063 : 1)
** u= 105/52 ; tau(u)= 1/53 ; 5407*x^2 - 5617*y^2 + 11026*x*z + 5407*z^2
  (-1583/473 : 1076/473 : 1)  C1a (-1897/783 : -10529/32103 : 1)
** u= -112/117 ; tau(u)= 346/229 ; -92338*x^2 + 14834*y^2 + 132260*x*z - 92338*z^2
  (-1187/1987 : 7374/1987 : 1)  C2b (308072/103445 : 959/1217 : 1)
** u= 116/121 ; tau(u)= 126/5 ; 13406*x^2 + 15826*y^2 + 29332*x*z + 13406*z^2
  (-2603/3656 : 759/3656 : 1)  C1b (6023/412 : 38353/16892 : 1)
** u= -116/149 ; tau(u)= 414/265 ; -126994*x^2 + 30946*y^2 + 184852*x*z - 126994*z^2
  (-1510/2383 : 7359/2383 : 1)  C2b (203055/130316 : -39901/130316 : 1)
** u= -116/173 ; tau(u)= 462/289 ; -153586*x^2 + 46402*y^2 + 226900*x*z - 153586*z^2
  (36268/70407 : 90967/70407 : 1)  C2b (-19901/337492 : -85779/337492 : 1)
** u= -119/89 ; tau(u)= 297/208 ; -72367*x^2 + 1681*y^2 + 102370*x*z - 72367*z^2
  (51/509 : -3112/509 : 1)  C2b (16336/6737 : 442229/276217 : 1)
** u= 119/145 ; tau(u)= 171/26 ; 12809*x^2 + 27889*y^2 + 43402*x*z + 12809*z^2
  (-1399/613 : 85662/102371 : 1)  C1b (4275052/9617 : 98418917/1606039 : 1)
** u= -120/101 ; tau(u)= 322/221 ; -83282*x^2 + 6002*y^2 + 118084*x*z - 83282*z^2
  (821/4465 : 14624/4465 : 1)  C2b (-171688/46507 : 101457/46507 : 1)
** u= 120/169 ; tau(u)= 218/49 ; 9598*x^2 + 42722*y^2 + 61924*x*z + 9598*z^2
  (-2683/814 : -1183/814 : 1)  C1b (-1241/2384 : -13689/97744 : 1)
** u= 126/5 ; tau(u)= 116/121 ; -13406*x^2 - 15826*y^2 + 29332*x*z - 13406*z^2
  (393/256 : 11/256 : 1)  C1a (-11395/428 : 71483/17548 : 1)
** u= 127/153 ; tau(u)= 179/26 ; 14777*x^2 + 30689*y^2 + 48170*x*z + 14777*z^2
  (-22699/37949 : 20262/37949 : 1)  C1b (43039/38377 : -9113/38377 : 1)
** u= 128/145 ; tau(u)= 162/17 ; 15806*x^2 + 25666*y^2 + 42628*x*z + 15806*z^2
  (-3593/1595 : 36/1595 : 1)  C1b (1520/7111 : 45533/291551 : 1)
** u= -128/153 ; tau(u)= 434/281 ; -141538*x^2 + 30434*y^2 + 204740*x*z - 141538*z^2
  (-3337/503 : -8016/503 : 1)  C2b (-14874939/93208 : -4279093/93208 : 1)
** u= -132/97 ; tau(u)= 326/229 ; -87458*x^2 + 1394*y^2 + 123700*x*z - 87458*z^2
  (27/122 : 829/122 : 1)  C2b (-2500/967 : 138885/39647 : 1)
** u= -132/181 ; tau(u)= 494/313 ; -178514*x^2 + 48098*y^2 + 261460*x*z - 178514*z^2
  (-1261/12634 : -26171/12634 : 1)  C2b (441284/1557007 : 333801/1557007 : 1)
** u= 135/34 ; tau(u)= 67/101 ; -2177*x^2 - 15913*y^2 + 22714*x*z - 2177*z^2
  (175/57 : 98/57 : 1)  C1a (958332/1737047 : 250909/1737047 : 1)
** u= 138/85 ; tau(u)= -32/53 ; 13426*x^2 - 4594*y^2 + 20068*x*z + 13426*z^2
  (-17961/8711 : 21934/8711 : 1)  C1a (-4792/5243 : -947/5243 : 1)
** u= 141/20 ; tau(u)= 101/121 ; -9401*x^2 - 19081*y^2 + 30082*x*z - 9401*z^2
  (829/1965 : -572/1965 : 1)  C1a (-1511696/262795 : 17207/20215 : 1)
** u= 141/32 ; tau(u)= 77/109 ; -3881*x^2 - 17833*y^2 + 25810*x*z - 3881*z^2
  (683/163 : 232/163 : 1)  C1a (47913/20465 : -1307/4093 : 1)
** u= 145/64 ; tau(u)= 17/81 ; 7903*x^2 - 12833*y^2 + 21314*x*z + 7903*z^2
  (3257/12745 : 13248/12745 : 1)  C1a (1520/7111 : 45533/291551 : 1)
** u= 146/29 ; tau(u)= 88/117 ; -6062*x^2 - 19634*y^2 + 29060*x*z - 6062*z^2
  (5761/3263 : 3780/3263 : 1)  C1a (357024/101683 : 47053/101683 : 1)
** u= 147/157 ; tau(u)= 167/10 ; 21409*x^2 + 27689*y^2 + 49498*x*z + 21409*z^2
  (-33283/20869 : 6958/20869 : 1)  C1b (-10257700/959209 : 1458891/959209 : 1)
** u= 150/37 ; tau(u)= 76/113 ; -3038*x^2 - 19762*y^2 + 28276*x*z - 3038*z^2
  (49/24 : 35/24 : 1)  C1a (16805/7931 : 96377/325171 : 1)
** u= -152/197 ; tau(u)= 546/349 ; -220498*x^2 + 54514*y^2 + 321220*x*z - 220498*z^2
  (36083/46527 : 64264/46527 : 1)  C2b (141269/390432 : 83347/390432 : 1)
** u= 153/32 ; tau(u)= 89/121 ; -5873*x^2 - 21361*y^2 + 31330*x*z - 5873*z^2
  (20507/37859 : 25080/37859 : 1)  C1a (2019/1232 : -11839/50512 : 1)
** u= 153/104 ; tau(u)= -55/49 ; 18607*x^2 - 1777*y^2 + 26434*x*z + 18607*z^2
  (-6089/8763 : 868/381 : 1)  C1a (9941848/711897 : 4453103/711897 : 1)
** u= 154/73 ; tau(u)= 8/81 ; 10594*x^2 - 13058*y^2 + 23780*x*z + 10594*z^2
  (-1223/458 : 603/458 : 1)  C1a (164912/502291 : -88579/502291 : 1)
** u= 161/52 ; tau(u)= 57/109 ; 2159*x^2 - 20513*y^2 + 29170*x*z + 2159*z^2
  (-717/10423 : -928/10423 : 1)  C1a (2027656/334297 : -273819/334297 : 1)
** u= 161/169 ; tau(u)= 177/8 ; 25793*x^2 + 31201*y^2 + 57250*x*z + 25793*z^2
  (-47169/42661 : -18668/42661 : 1)  C1b (2279/513 : -15649/21033 : 1)
** u= 162/17 ; tau(u)= 128/145 ; -15806*x^2 - 25666*y^2 + 42628*x*z - 15806*z^2
  (823/985 : -576/985 : 1)  C1a (436480/222199 : -2391413/9110159 : 1)
** u= 162/113 ; tau(u)= -64/49 ; 21442*x^2 - 706*y^2 + 30340*x*z + 21442*z^2
  (667/1153 : -9324/1153 : 1)  C1a (-31643/126280 : 15391/25256 : 1)
** u= 165/52 ; tau(u)= 61/113 ; 1687*x^2 - 21817*y^2 + 30946*x*z + 1687*z^2
  (-79315/2559901 : 468196/2559901 : 1)  C1a (-487440/45293 : -64061/45293 : 1)
** u= 167/10 ; tau(u)= 147/157 ; -21409*x^2 - 27689*y^2 + 49498*x*z - 21409*z^2
  (2949/4651 : -1034/4651 : 1)  C1a (431341/365380 : 63201/365380 : 1)
** u= 168/173 ; tau(u)= 178/5 ; 28174*x^2 + 31634*y^2 + 59908*x*z + 28174*z^2
  (-3079/4197 : 584/4197 : 1)  C1b (265912/6533 : -40887/6533 : 1)
** u= 171/26 ; tau(u)= 119/145 ; -12809*x^2 - 27889*y^2 + 43402*x*z - 12809*z^2
  (9/7 : 1034/1169 : 1)  C1a (-15500/207 : 358033/34569 : 1)
** u= 173/181 ; tau(u)= 189/8 ; 29801*x^2 + 35593*y^2 + 65650*x*z + 29801*z^2
  (-26497/20221 : -7620/20221 : 1)  C1b (-464693/176904 : 4793/13608 : 1)
** u= -176/197 ; tau(u)= 570/373 ; -247282*x^2 + 46642*y^2 + 355876*x*z - 247282*z^2
  (21039/1523209 : 3472558/1523209 : 1)  C2b (436576/98845 : 114297/98845 : 1)
** u= 177/8 ; tau(u)= 161/169 ; -25793*x^2 - 31201*y^2 + 57250*x*z - 25793*z^2
  (5209/7661 : -1508/7661 : 1)  C1a (-513/2279 : 15649/93439 : 1)
** u= -177/193 ; tau(u)= 563/370 ; -242471*x^2 + 43169*y^2 + 348298*x*z - 242471*z^2
  (191/1303 : 19466/9121 : 1)  C2b (64868/96877 : -152811/678139 : 1)
** u= 177/197 ; tau(u)= 217/20 ; 30529*x^2 + 46289*y^2 + 78418*x*z + 30529*z^2
  (-1791/859 : -64/859 : 1)  C1b (63749/20560 : 437331/842960 : 1)
** u= 178/5 ; tau(u)= 168/173 ; -28174*x^2 - 31634*y^2 + 59908*x*z - 28174*z^2
  (2281/3162 : 347/3162 : 1)  C1a (105968/103385 : -16377/103385 : 1)
** u= 179/26 ; tau(u)= 127/153 ; -14777*x^2 - 30689*y^2 + 48170*x*z - 14777*z^2
  (1639/4541 : 678/4541 : 1)  C1a (-15782268/1153985 : 447535/230797 : 1)
** u= 182/101 ; tau(u)= -20/81 ; 20002*x^2 - 12722*y^2 + 33524*x*z + 20002*z^2
  (-3118/1795 : 2367/1795 : 1)  C1a (-11985308/2441183 : 6643/8447 : 1)
** u= 186/73 ; tau(u)= 40/113 ; 9058*x^2 - 23938*y^2 + 36196*x*z + 9058*z^2
  (2843/6609 : -6928/6609 : 1)  C1a (-4640/7999 : -1111/7999 : 1)
** u= 187/197 ; tau(u)= 207/10 ; 34769*x^2 + 42649*y^2 + 77818*x*z + 34769*z^2
  (-3967/3709 : -1674/3709 : 1)  C1b (-1051488549/428334749 : -140459269/428334749 : 1)
** u= -188/137 ; tau(u)= 462/325 ; -175906*x^2 + 2194*y^2 + 248788*x*z - 175906*z^2
  (-4538/19963 : 209465/19963 : 1)  C2b (73547/50963 : 61173/50963 : 1)
** u= -188/185 ; tau(u)= 558/373 ; -242914*x^2 + 33106*y^2 + 346708*x*z - 242914*z^2
  (29101/89360 : -193833/89360 : 1)  C2b (-2463612/88615 : -902899/88615 : 1)
** u= 189/8 ; tau(u)= 173/181 ; -29801*x^2 - 35593*y^2 + 65650*x*z - 29801*z^2
  (15399/22699 : 3844/22699 : 1)  C1a (652725/79984 : 93065/79984 : 1)
** u= 189/128 ; tau(u)= -67/61 ; 28279*x^2 - 2953*y^2 + 40210*x*z + 28279*z^2
  (1487/8383 : 29392/8383 : 1)  C1a (-364056/59795 : -26471/11959 : 1)
** u= -196/153 ; tau(u)= 502/349 ; -205186*x^2 + 8402*y^2 + 290420*x*z - 205186*z^2
  (-1384/1073 : 11235/1073 : 1)  C2b (592877/10949 : 379939/10949 : 1)
** u= -196/181 ; tau(u)= 558/377 ; -245842*x^2 + 27106*y^2 + 349780*x*z - 245842*z^2
  (43931/5704 : 120687/5704 : 1)  C2b (2807940/748957 : -929425/748957 : 1)
** u= 198/89 ; tau(u)= 20/109 ; 15442*x^2 - 23362*y^2 + 39604*x*z + 15442*z^2
  (-8/21 : -1/3 : 1)  C1a (1661052/154445 : 248333/154445 : 1)
** u= 198/125 ; tau(u)= -52/73 ; 28546*x^2 - 7954*y^2 + 41908*x*z + 28546*z^2
  (-15559/420776 : 775755/420776 : 1)  C1a (48140/827 : 506203/33907 : 1)
** u= -199/193 ; tau(u)= 585/392 ; -267727*x^2 + 34897*y^2 + 381826*x*z - 267727*z^2
  (729/1427 : 2884/1427 : 1)  C2b (5950927/6069465 : 1697711/6069465 : 1)
144
>


■以下では、u=198/125のときに、(3a+),(3b±)を満たす有理数の組(x,y,t)を求める。
u=198/125のとき、(3a),(3b±)より
     7954*y^2 = 28546*x^2 + 41908*x + 28546 ------- (4a)
     ±15816529*t^2=1016269*x^2 + 1384481*x + 1016269 ------- (4b±)
となる。

ここで、２次曲線(4a)の有理点の１つを探す。
例えば、(x0,y0)=(-15559/420776, 775755/420776)は(3a)の有理点なので、任意の有理数kに対して、直線
     y=k*(x+15559/420776) + 775755/420776 ---------------(6)

と２次曲線(4a)の交点は高々2個の有理点であり、(x0,y0)と異なる交点(x(k),y(k))は、以下のようになる。
     x(k)=(61878143*k^2 + 6170355270*k - 8594866697)/(-1673426152*k^2 + 6005735848) ------- (7)
     y(k)=(3085177635*k^2 - 8372793090*k + 11072351115)/(-1673426152*k^2 + 6005735848) ----- - (8)

さらに、(7)を(3b±)に代入して、
     ±t^2=(-7015712318801663789*k^4 + 35045844339021901740*k^3 - 97498048856953916558*k^2 + 146426811189378455340*k - 104374351035488335469)/(-114814558534098611264*k^4 + 824112745263861945472*k^2 - 1478823386114043443264) -------- (9)
を得る。
両辺に適当な有理数の平方数を掛けて、
     ±(t*(3977*k^2 - 14273)*17251816)^2=287644205070868215349*k^4 - 1436879617899897971340*k^3 + 3997420003135110578878*k^2 - 6003499258764516668940*k + 4279348392455021754229 ---------- (10±)
つまり、
     ±□=287644205070868215349*k^4 - 1436879617899897971340*k^3 + 3997420003135110578878*k^2 - 6003499258764516668940*k + 4279348392455021754229 ---------- (11±)
を満たす有理数kを求めれば良い。


■(11)の右辺の4次多項式は、実根を持たない(負にならない)ので、(11+)のみ、考慮すれば良い。
よって、楕円曲線

     E+: Y^2=287644205070868215349*X^4 - 1436879617899897971340*X^3 + 3997420003135110578878*X^2 - 6003499258764516668940*X + 4279348392455021754229 --------- (12)

の有理点を求める。

楕円曲線E+をsyzygy mapで写すと。楕円曲線
     E0: y^2=x^3 - 131533381367244164215248081536630341171740672*x - 496946952622325121564854492450656958760656582555086171408919691264
が得られる。
E0のminimal standard modelは、楕円曲線
     E2: y^2=x^3 + x^2 - 3237629437155959990*x - 1919096286007227296033058600
であることが分かる。

[pari/gpによるsyzygyの計算]
(08:08) gp > v=V(287644205070868215349*x^4 - 1436879617899897971340*x^3 + 3997420003135110578878*x^2 - 6003499258764516668940*x + 4279348392455021754229)
%1 = [287644205070868215349, -1436879617899897971340, 3997420003135110578878, -6003499258764516668940, 4279348392455021754229]
(08:09) gp > e0=E0(v)
v=[287644205070868215349, -1436879617899897971340, 3997420003135110578878, -6003499258764516668940, 4279348392455021754229]
I=4871606717305339415379558575430753376731136
J=18405442689715745243142758979653961435579873427966154496626655232
%2 = [0, 0, 0, -131533381367244164215248081536630341171740672, -496946952622325121564854492450656958760656582555086171408919691264, 0, -263066762734488328430496163073260682343481344, -1987787810489300486259417969802627835042626330220344685635678765056, -17301030413900894278369859986709427802384131133245294602033481221399655787353162419011584, 6313602305627719882331907913758256376243552256, 429362167065688905032034281477367612369207287327594452097306613252096, 38957323767280780602286963641277929776911445974775009083860804370163649169527768164395707769955957793908669755309328911621423104000000, 12345879053944502914377148189239145836992482212544/1911082634971811704047467743911970767517640625, Vecsmall([1]), [Vecsmall([128, 1])], [0, 0, 0, 0, 0, 0, 0, 0]]
(08:09) gp > e2=E2(v)
v=[287644205070868215349, -1436879617899897971340, 3997420003135110578878, -6003499258764516668940, 4279348392455021754229]
I=4871606717305339415379558575430753376731136
J=18405442689715745243142758979653961435579873427966154496626655232
rr=[2524656, 2124629306112, 0, 0]
%3 = [0, 1, 0, -3237629437155959990, -1919096286007227296033058600, 4, -6475258874311919980, -7676385144028909184132234400, -10482244380015203422289683962613034500, 155406212983486079536, 1658099190177807105871646153216, 580981870673987692214356057322685073501028465608089000000, 12345879053944502914377148189239145836992482212544/1911082634971811704047467743911970767517640625, Vecsmall([1]), [Vecsmall([128, 1])], [0, 0, 0, 0, 0, [64096668788529680160, 73728, [2, 5; 3, 1; 5, 1; 11, 1; 13, 1; 17, 1; 41, 2; 73, 1; 97, 1; 131, 1; 35227, 1], [[5, 3, 0, 2], [1, 10, 0, 6], [1, 10, 0, 6], [1, 6, 0, 2], [1, 6, 0, 2], [1, 6, 0, 2], [2, -8, 0, 4], [1, 6, 0, 2], [1, 8, 0, 4], [1, 6, 0, 2], [1, 6, 0, 2]]], 0, [[2, 41]~]]]


■楕円曲線E+をある楕円曲線の2-descentと見なして、MAGMAで4-descentを求めると、成功する。
E+は、楕円曲線E0の2-descentの１つになっているので、そのminimal standard modelである楕円曲線E2の有理点を求める。
さらに、E2の有理点から、E0の有理点を求める。

[MAGMA 4-descentによる計算]
> SetClassGroupBounds("GRH");
> P := PolynomialRing(Rationals());
> f :=  287644205070868215349*x^4 - 1436879617899897971340*x^3 + 3997420003135110578878*x^2 - 6003499258764516668940*x\
 + 4279348392455021754229   ;
> C := HyperellipticCurve((f));
> time fd := FourDescent(C : RemoveTorsion);
Time: 22.078
> #fd;
2
> RP4(fd,10^6);
J=1
Elliptic Curve defined by y^2 = x^3 + x^2 - 3237629437155959990*x - 1919096286007227296033058600 over Rational Field
rootno=1
(-174942651055478264675987/244780606357156 : -20891124644686233492499046337781371/3829706654220268054696 : 1)
height 36.9962224017345060609216575337
true (-174942651055478264675987/244780606357156 : 20891124644686233492499046337781371/3829706654220268054696 : 1)
(-174942651055478264675987/244780606357156 : -20891124644686233492499046337781371/3829706654220268054696 : 1)
height 36.9962224017345060609216575337
true (-174942651055478264675987/244780606357156 : 20891124644686233492499046337781371/3829706654220268054696 : 1)
(-174942651055478264675987/244780606357156 : -20891124644686233492499046337781371/3829706654220268054696 : 1)
height 36.9962224017345060609216575337
true (-174942651055478264675987/244780606357156 : 20891124644686233492499046337781371/3829706654220268054696 : 1)
(-174942651055478264675987/244780606357156 : -20891124644686233492499046337781371/3829706654220268054696 : 1)
height 36.9962224017345060609216575337
true (-174942651055478264675987/244780606357156 : 20891124644686233492499046337781371/3829706654220268054696 : 1)
J=2
Elliptic Curve defined by y^2 = x^3 + x^2 - 3237629437155959990*x - 1919096286007227296033058600 over Rational Field
rootno=1
(-14294068250419325295136330646500671045398507/10710488752832501604605089844433796 :
-5513451734757779144704564173941245101175055697275239743153358521/1108444452616653816973134639524599540568006219348344 :
1)
height 82.3227669260234809866404587323
true (-14294068250419325295136330646500671045398507/10710488752832501604605089844433796 :
5513451734757779144704564173941245101175055697275239743153358521/1108444452616653816973134639524599540568006219348344 :
1)
(-14294068250419325295136330646500671045398507/10710488752832501604605089844433796 :
-5513451734757779144704564173941245101175055697275239743153358521/1108444452616653816973134639524599540568006219348344 :
1)
height 82.3227669260234809866404587323
true (-14294068250419325295136330646500671045398507/10710488752832501604605089844433796 :
5513451734757779144704564173941245101175055697275239743153358521/1108444452616653816973134639524599540568006219348344 :
1)
(-12809308108101805039622169849007356722051087806079663963/10662986919353356508233584268494411612549069124 :
16938726381485529637262909781345059397306307345660154988348340726522206651166214127/11010785814195173067531823505807795\
34492121078438995214271032277893432 : 1)
height 109.928798827137085827881942845
true (-12809308108101805039622169849007356722051087806079663963/10662986919353356508233584268494411612549069124 :
16938726381485529637262909781345059397306307345660154988348340726522206651166214127/11010785814195173067531823505807795\
34492121078438995214271032277893432 : 1)
(-14294068250419325295136330646500671045398507/10710488752832501604605089844433796 :
-5513451734757779144704564173941245101175055697275239743153358521/1108444452616653816973134639524599540568006219348344 :
1)
height 82.3227669260234809866404587323
true (-14294068250419325295136330646500671045398507/10710488752832501604605089844433796 :
5513451734757779144704564173941245101175055697275239743153358521/1108444452616653816973134639524599540568006219348344 :
1)
(-12809308108101805039622169849007356722051087806079663963/10662986919353356508233584268494411612549069124 :
16938726381485529637262909781345059397306307345660154988348340726522206651166214127/11010785814195173067531823505807795\
34492121078438995214271032277893432 : 1)
height 109.928798827137085827881942845
true (-12809308108101805039622169849007356722051087806079663963/10662986919353356508233584268494411612549069124 :
16938726381485529637262909781345059397306307345660154988348340726522206651166214127/11010785814195173067531823505807795\
34492121078438995214271032277893432 : 1)
realtime=3.255
これらのE2の有理点から、一次独立なものを選択すると、
       P1(-174942651055478264675987/244780606357156, 20891124644686233492499046337781371/3829706654220268054696)
       P2(-14294068250419325295136330646500671045398507/10710488752832501604605089844433796, 5513451734757779144704564173941245101175055697275239743153358521/1108444452616653816973134639524599540568006219348344)
の2個である。よって、rank(E2)は2以上である。


[pari/gpによる計算]
(13:35) gp > v=V(287644205070868215349*x^4 - 1436879617899897971340*x^3 + 3997420003135110578878*x^2 - 6003499258764516668940*x + 4279348392455021754229)
%4 = [287644205070868215349, -1436879617899897971340, 3997420003135110578878, -6003499258764516668940, 4279348392455021754229]
(13:58) gp > e=E2(v)
v=[287644205070868215349, -1436879617899897971340, 3997420003135110578878, -6003499258764516668940, 4279348392455021754229]
I=4871606717305339415379558575430753376731136
J=18405442689715745243142758979653961435579873427966154496626655232
rr=[2524656, 2124629306112, 0, 0]
%5 = [0, 1, 0, -3237629437155959990, -1919096286007227296033058600, 4, -6475258874311919980, -7676385144028909184132234400, -10482244380015203422289683962613034500, 155406212983486079536, 1658099190177807105871646153216, 580981870673987692214356057322685073501028465608089000000, 12345879053944502914377148189239145836992482212544/1911082634971811704047467743911970767517640625, Vecsmall([1]), [Vecsmall([128, 1])], [0, 0, 0, 0, 0, [64096668788529680160, 73728, [2, 5; 3, 1; 5, 1; 11, 1; 13, 1; 17, 1; 41, 2; 73, 1; 97, 1; 131, 1; 35227, 1], [[5, 3, 0, 2], [1, 10, 0, 6], [1, 10, 0, 6], [1, 6, 0, 2], [1, 6, 0, 2], [1, 6, 0, 2], [2, -8, 0, 4], [1, 6, 0, 2], [1, 8, 0, 4], [1, 6, 0, 2], [1, 6, 0, 2]]], 0, [[2, 41]~]]]
(13:58) gp > P1=r(e,-174942651055478264675987/244780606357156)
%6 = [-174942651055478264675987/244780606357156, 20891124644686233492499046337781371/3829706654220268054696]
(14:00) gp > P2=r(e,-14294068250419325295136330646500671045398507/10710488752832501604605089844433796)
%7 = [-14294068250419325295136330646500671045398507/10710488752832501604605089844433796, 5513451734757779144704564173941245101175055697275239743153358521/1108444452616653816973134639524599540568006219348344]
(14:01) gp > P3=r(e,-12809308108101805039622169849007356722051087806079663963/10662986919353356508233584268494411612549069124)
%8 = [-12809308108101805039622169849007356722051087806079663963/10662986919353356508233584268494411612549069124, 16938726381485529637262909781345059397306307345660154988348340726522206651166214127/1101078581419517306753182350580779534492121078438995214271032277893432]
(14:01) gp > matdet(ellheightmatrix(e,[P1,P2]))
time = 2 ms.
%9 = 2139.9395575113731884578223884306279106
(14:01) gp > matdet(ellheightmatrix(e,[P1,P2,P3]))
time = 1 ms.
%10 = 3.622301079776301326 E-33

さらに、P1,P2を有理変換[2524656, 2124629306112, 0, 0]-1で写すことにより、(楕円曲線E+の有理点をsyzygy mapで写した)楕円曲線E0の有理点を求めると、
       Q1(-278766212361028884962560363157615040/61195151589289, 42022169170381781318268794555009598687378195785697792/478713331777533506837)
       Q2(-22777197225615531144367809302749081158104851178197540800/2677622188208125401151272461108449, 11090221586976991451411531078086682941672735310414994171896001971993851427077934592/138555556577081727121641829940574942571000777418543)
となる。ここで、E0の有理点Q1,Q2は一次独立である。

[pari/gpによる計算]
(14:02) gp > Q1=chpi(v,P1[1])
v=[287644205070868215349, -1436879617899897971340, 3997420003135110578878, -6003499258764516668940, 4279348392455021754229]
I=4871606717305339415379558575430753376731136
J=18405442689715745243142758979653961435579873427966154496626655232
rr=[2524656, 2124629306112, 0, 0]
v=[287644205070868215349, -1436879617899897971340, 3997420003135110578878, -6003499258764516668940, 4279348392455021754229]
I=4871606717305339415379558575430753376731136
J=18405442689715745243142758979653961435579873427966154496626655232
time = 1 ms.
%11 = [-278766212361028884962560363157615040/61195151589289, 42022169170381781318268794555009598687378195785697792/478713331777533506837]
(14:13) gp > Q2=chpi(v,P2[1])
v=[287644205070868215349, -1436879617899897971340, 3997420003135110578878, -6003499258764516668940, 4279348392455021754229]
I=4871606717305339415379558575430753376731136
J=18405442689715745243142758979653961435579873427966154496626655232
rr=[2524656, 2124629306112, 0, 0]
v=[287644205070868215349, -1436879617899897971340, 3997420003135110578878, -6003499258764516668940, 4279348392455021754229]
I=4871606717305339415379558575430753376731136
J=18405442689715745243142758979653961435579873427966154496626655232
%12 = [-22777197225615531144367809302749081158104851178197540800/2677622188208125401151272461108449, 11090221586976991451411531078086682941672735310414994171896001971993851427077934592/138555556577081727121641829940574942571000777418543]

■楕円曲線E0の有理点Q1,Q2から、E0の有理点をいくつか計算する。
さらに、syzygy mapの逆変換によって、楕円曲線E+の有理点P(x,y)を求めると、そのx座標は以下のようになる。
 30407075/26279287,
 945967865/1248057599,
 1588437129/13721456335,
 -1219384020107/2053661177915,
  7238621986098507/15401438000287805,
 7238621986098507/15401438000287805,
 2523484525026218281/616919842550050055,
 2523484525026218281/616919842550050055,
 7918680427205562475/1555531906528806611,
 7918680427205562475/1555531906528806611,
 206736470097243783595/455301048052810772797,
 206736470097243783595/455301048052810772797,
 67894850173916916719315/284272617120491254805807,
 -11406552876127364645246855/8312938824148780426752383,
 290157027339362326908805193/203403933781733781606373135,
 24836379138590189117837568873/37450259193857864779920643975,
 59119753168574500344511852042541599146030143212531/80531942916272290909196354054205076510511037257045,
 5899052009973968657607739374986290502199527513150091/4872598158246038860927066737600301097060213712856325,
 -7156549106013460647943623161255911730566650508310945/9904307575498282160550911204810848846693742352348679,
 306339257986856753681671914442787946614938121166815315/2070347241147865363291315221807210585802457450699932501,
...

[pari/gpによる計算]
(14:14) gp > L=ss2x(v,Q1,Q2,10);
v=[287644205070868215349, -1436879617899897971340, 3997420003135110578878, -6003499258764516668940, 4279348392455021754229]
I=4871606717305339415379558575430753376731136
J=18405442689715745243142758979653961435579873427966154496626655232
time = 1,605 ms.
(14:14) gp > #L
%14 = 440
(14:15) gp > L[1..20]
%15 = [2558149/3980783, 41021807/13655501, 14788729/49974203, 330413341/86871163, -31145871143/3862611419, -3019955969993/834385389601, 13717287202379/8893169459297, 126066305416193/73804006159651, 431808066142349/326200619036507, 482060239105477/246713455686239, 10394002144961011/192432661134137, -56530686689987759/35997071277915763, 42041898049222025570390197943/16054003360228530162598916101, -17343999215485048407387590807/618879358155435961929403117999, 8293517906945049688708880839333/9703099224128674332461216478319, 52368010678183032840070558459789/10827070089912940353261785022263, 34309115383062652319268813713938321/4130474606045023602688640458990547, 969414219124685069574749128934218271/876930469304553073908384414500020553, -3607802204449280119048946937332963149/5725928851961196570604549935514845593, -4686413274039099144252404080250569589/28653439382230384730327228889914091777]


■kをE+の有理点のx座標とする。kを(7),(8),(9)に代入すると、方程式系(3+a),(3b+)を満たす有理数の組(x,y,t)が求まる。
つまり、(2)を満たす有理数の組(x,y,t)が求まる。
(2)の両辺に適当な整数を掛けることにより、(1)の整点(A,B,C,D)をいくつか求めることができる。
ここで、0 <= A <= B <= C を満たすように、A,B,Cを交換して、Dの小さい順に(1)の等式を並べ替えると、以下のようになる。

[pari/gpによる計算]
(14:26) gp > sss4(41,u0,L[1..20])
35401855^4+40865628^4+53562031^4=3362*7822733^4
76298339723306940^4+144376098024837517^4+392097054273222611^4=3362*51745745604910607^4
22630934278564908444565^4+218426009168410193424516^4+383430470008039234123883^4=3362*51630869931062938919159^4
66362246478987836342552992091201664685210465639474020^4+96903091845971137820040019904874448516482466179139001^4+124994312898188111523490811243725150204657332635449303^4=3362*17983883123274685049847136904533151182677188556511541^4
37682872887674456007962886831727581567709687419316284566119455689^4+87890763112391594397658900321621501270731784581369209175727203300^4+123973647924119793711813764065522766229845905720563581825933541127^4=3362*17253275064778568487130075273315604954065636718785561424329636709^4
691048672272960450860044325926062322311036853674019263752393076169781180^4+1113195841691642261609432943293771082922611362472659021789669936692090693^4+1451927091024265743100085140482652585315995467075892655672274233243677261^4=3362*207292626492004064390780316253678173463653077794812959325479503602344027^4
time = 1 ms.
%18 = 6

[2025.10.17追記]
u=198/125のときとu=τ(198/125)=-52/73のときは、(1)の同じ整点が得られるので、高さの小さいu=-52/73の方を代表として扱う。
■さらに、u=52/137, 17/81, 44/169のときについても、(1)の整点を見つけることができたので、一緒にまとめると、以下のようになる。

[参考文献]

Last Update: 2025.11.25
H.Nakao

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