login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A166929 Positive integers m such that m^4 = a^2 + b^2 and a + b = c^2 for some coprime integers a, b, c. 4
1, 13, 1525, 2165017, 42422452969, 7658246457672229, 15512114571284835412957, 452005526897888844293504165425, 126314830357375266295717376544111167953, 368440923990671763222767414151367493861848396861 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,2
COMMENTS
Values of m in coprime solutions to 2m^4 = c^4 + d^2 (so that a, b = (c^2 +- d)/2).
Corresponding values of c are given in A167437.
Terms with positive a,b,c are given in A166930.
This is a generalized Somos-4 sequence. - Michael Somos, Jan 29 2023
LINKS
Keith Conrad, Fermat's Big Problem, 2008, slides 20-30.
FORMULA
a/m^2 = (-8*u^2 + 32*v + u^4 + 4*v^2 - 4*u^3 + 16*u*v)/(8 + u^2 + 4*u)^2 and b/m^2 = (4*u^3 - 8*u^2 - 4*v*u^2 - 16*u*v + 64)/(8 + u^2 + 4*u)^2 where (u,v) is a rational point on the elliptic curve v^2 = u^3 + 8*u.
a(n) = a(1-n) = (1764*a(n-1)*a(n-3) - 3107*a(n-2)^2)/a(n-4) = A360187(2*n-1) for all n in Z. - Michael Somos, Jan 29 2023
PROG
(PARI) {a(n) = my(A); if(n<1, n=1-n); A = vector(max(4, n+2)); A[1] = 13; A[2] = 1; A[3] = 1; A[4] = 13; for(k=5, n+2, A[k] = (1764*A[k-1]*A[k-3] - 3107*A[k-2]^2)/A[k-4]); A[n+2]}; /* Michael Somos, Jan 29 2023 */
CROSSREFS
Cf. A360187.
Sequence in context: A185193 A066457 A203515 * A079917 A028450 A201177
KEYWORD
nonn
AUTHOR
Max Alekseyev, Oct 23 2009
EXTENSIONS
Edited by Max Alekseyev, Nov 03 2009
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 29 06:44 EDT 2024. Contains 371265 sequences. (Running on oeis4.)