OFFSET
1,1
COMMENTS
For n > 0, k = (n + 1)(2n + 1)^2 is a term in this sequence, because k^2 = (n * (2n + 1)^2)^2 + (2n + 1)^5. Examples: 18, 75, 196, 405, 726, 1183.
When z^2 = x^2 + y^2 (i.e., z = A009003(n)), (z * y^4)^2 = (x * y^4)^2 + (y^2)^5. Thus z * y^4 is a term in this sequence. For example, 1200. More generally, for positive integer i, j and k, x^(5i - 5) * y^(5j - 1) * z^(5k - 5) is in this sequence.
When z^2 = x^2 + y^3 (i.e., z = A070745(n)), (z * y)^2 = (x * y)^2 + y^5. Thus z * y is in this sequence. E.g. 6, 18, 40, ... . More generally, for positive integer i, j and k, x^(5i - 5) * y^(5j - 4) * z^(5k - 4) is in this sequence.
When z^2 = x^2 + y^4 (i.e., z = A271576(n)), (z * y^3)^2 = (x * y^3)^2 + (y^2)^5. Thus z * y^3 is also in this sequence. E.g. 40, 405, 1107, ... . More generally, for positive integer i, j and k, x^(5i - 5) * y^(5j - 2) * z^(5k - 4) is in this sequence.
LINKS
Chai Wah Wu, Table of n, a(n) for n = 1..10000
EXAMPLE
6^2 = 2^2 + 2^5.
9^2 = 7^2 + 2^5.
MATHEMATICA
c[n_]: = Count[n^2 - Range[(n^2 - 1)^(1/5)]^5, _?(IntegerQ[Sqrt[#]] &)] > 0;
Select[Range[1200], c]
PROG
(PARI) isok(n) = for (k=1, n-1, if (ispower(n^2-k^2, 5), return (1)); ); return (0); \\ Michel Marcus, Oct 06 2017
CROSSREFS
KEYWORD
nonn
AUTHOR
XU Pingya, Oct 04 2017
STATUS
approved