OFFSET
1,1
LINKS
Colin Barker, Table of n, a(n) for n = 1..1000
Index entries for linear recurrences with constant coefficients, signature (7,-21,35,-35,21,-7,1).
FORMULA
From Colin Barker, Sep 05 2019: (Start)
G.f.: 9*x*(1 + x)*(1 + 8*x + 22*x^2 + 8*x^3 + x^4) / (1 - x)^7.
a(n) = n^2*(2 + n^2)^2.
a(n) = 7*a(n-1) - 21*a(n-2) + 35*a(n-3) - 35*a(n-4) + 21*a(n-5) - 7*a(n-6) + a(n-7) for n>7.
(End)
E.g.f.: exp(x)*x*(9 + 63*x + 114*x^2 + 69*x^3 + 15*x^4 + x^5). - Conjectured by Stefano Spezia, Sep 05 2019 after Colin Barker
From Chai Wah Wu, Sep 10 2019: (Start)
Above conjectures are true. Proof: k + (k+1)^2 + (k+2)^3 = (k+1)*(k+3)^2 and thus is a perfect square if and only if k+1 = n^2 is a perfect square. This implies that (k+1)*(k+3)^2 = n^2*(n^2+2)^2.
(End)
EXAMPLE
9 = 0 + 1^2 + 2^3. 0,1,2 are consecutive numbers and 9 is a perfect square. Hence, 9 is a member of the sequence.
18225 = 24 + 25^2 + 26^3. 24,25,26 are consecutive numbers and 18225 is a perfect square. Hence 18225 is a member of the sequence.
PROG
(PARI) a(n) = n^2*(2 + n^2)^2 \\ David A. Corneth, Sep 11 2019
(PARI) Vec(9*x*(1 + x)*(1 + 8*x + 22*x^2 + 8*x^3 + x^4) / (1 - x)^7 + O(x^40)) \\ Colin Barker, Sep 11 2019
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Philip Mizzi, Sep 05 2019
STATUS
approved