login
This site is supported by donations to The OEIS Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A156668 Positive integers k such that k^2 = (m^5 + n^5)/(m + n) for some coprime integers m, n. 4
1, 11, 101, 13361, 1169341, 1612186411, 1624763543401, 20188985439712961, 240020196429554642201, 29891946989942513908518251, 3506790234728288196345900732301, 5190947078637547438603476743093680561 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,2

COMMENTS

This sequence probably contains no more than 5 primes.

LINKS

Table of n, a(n) for n=1..12.

Kevin Acres and David Broadhurst, Rational points on y^2 = x^3 + 10*x^2 + 5*x

Dave Rusin, Re: Help with diophantine equ., sci.math newsgroup [Broken link]

Dave Rusin, Re: Help with diophantine equ., sci.math newsgroup [Cached copy]

FORMULA

Numerators of rational numbers (81*x^4 + 540*x^3 - 8370*x^2 + 33900*x - 47975)/(9*x^2 - 150*x + 445)^2, where x ranges over abscissas of rational points on the elliptic curve y^2 = x^3 - 85/3*x + 1550/27.

EXAMPLE

13361 belongs to this sequence since 13361^2 = (35^5 + 123^5) / (35 + 123) with gcd(35, 123)=1.

PROG

(PARI) { a(k) = local(P=ellpow(ellinit([0, 10, 0, 5, 0]), [-1, 2], k), s, t); s=P[1]^2; t=abs(numerator(P[2]^4/s-80*s)); while(t%2==0, t=t/2); t } /* David Broadhurst */

CROSSREFS

Cf. A156669, A156670, A186082.

Sequence in context: A089971 A082620 A199304 * A103992 A185949 A001387

Adjacent sequences:  A156665 A156666 A156667 * A156669 A156670 A156671

KEYWORD

nonn

AUTHOR

Max Alekseyev, Feb 13 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 | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified April 19 05:14 EDT 2019. Contains 322237 sequences. (Running on oeis4.)