OFFSET
1,1
COMMENTS
If a^5 + b^5 = m, then (ma)^5 + (mb)^5 = m^6 = (m^3)^2 is square. Therefore A003347(n)^3 are terms of this sequence.
When k is in this sequence, k * (n^5) (n = 2, 3, ... ) is also in this sequence.
EXAMPLE
8^2 = 2^5 + 2^5, so 8 is in the sequence.
6655^2 = 22^5 + 33^5, so 6655 is in the sequence.
MATHEMATICA
lst={}; Do[If[IntegerQ[(n^2-a^5)^(1/5)], AppendTo[lst, n]], {n, 9000}, {a, (n^2/2)^(1/5)}]; lst
PROG
(PARI)
upto(n) = {
my(res = List(), u = n^2, i5);
for(i = 1, sqrtnint(u, 5),
i5 = i^5;
for(j = i, sqrtnint(u - i5, 5),
c = i5 + j^5;
if(issquare(c, &sc),
listput(res, sc))));
Set(res)} \\ David A. Corneth, Jun 17 2025
CROSSREFS
KEYWORD
nonn
AUTHOR
XU Pingya, Sep 04 2017
STATUS
approved
