%I #23 Mar 23 2022 16:19:21
%S 1,2,6,21,161,1357,18526,480106,12551561,683916417,51678803961,
%T 4881674119706,997454379905326,213822353304561757,
%U 79799551268268089761,53139223644814624290821,36631192030206080565822006,54202648602164057575419038802
%N Numerator of x-coordinate of (2n)*P where P is the generator for rational points on the curve y^2 + y = x^3 - x.
%H Seiichi Manyama, <a href="/A028936/b028936.txt">Table of n, a(n) for n = 1..106</a>
%H B. Mazur, <a href="https://doi.org/10.1090/S0273-0979-1986-15430-3">Arithmetic on curves</a>, Bull. Amer. Math. Soc. 14 (1986), 207-259; see p. 225.
%F P=(0, 0), 2P=(1, 0); if kP = (a, b) then (k+1)P = (a' = (b^2 - a^3)/a^2, b' = -1 - b*a'/a).
%F a(n) = A028940(2n). - _Seiichi Manyama_, Nov 19 2016
%F 0 = a(n)*a(n+6) - 5*a(n+1)*a(n+5) + 4*a(n+2)*a(n+4) - 20*a(n+3)^2 for all n in Z. a(n) = A006720(n+1)*A006720(n+2). - _Michael Somos_, Apr 12 2020
%e 4P =(2, -3).
%e a(3) = 6 = 2*3 = A006720(4)*A006720(5). - _Michael Somos_, Apr 12 2020
%o (PARI) a(n)=numerator(ellmul(E,[0,0],2*n)[1]) \\ _Charles R Greathouse IV_, Mar 23 2022
%Y Cf. A028937 (denominator), A028938, A028939, A028940.
%Y Cf. A006720.
%K nonn,frac
%O 1,2
%A _N. J. A. Sloane_
|