OFFSET
1,2
COMMENTS
Agrees with A155138 except for insertion of zero after a(6) = 15.
LINKS
Index entries for linear recurrences with constant coefficients, signature (4,-6,4,-1).
FORMULA
G.f.: 3*x*(9-20*x+9*x^2+16*x^5+x^6-19*x^7+x^8+5*x^9)/(1-x)^4.
EXAMPLE
MATHEMATICA
CoefficientList[Series[3x (9-20x+9x^2+16x^5+x^6-19x^7+x^8+5x^9)/(1-x)^4, {x, 0, 40}], x] (* or *) LinearRecurrence[{4, -6, 4, -1}, {0, 27, 48, 57, 48, 15, 0, 48, 147, 288, 477}, 50] (* Harvey P. Dale, Sep 02 2021 *)
PROG
(Magma) [ Integers()!SquareRoot(a) : n in [ -30..1500] | IsSquare(a) where a is n^3+28*n^2 ];
(PARI) a(n)=if(n>6, n^3 - 3*n^2 - 25*n + 27, [0, 27, 48, 57, 48, 15, 0][n+1]) \\ Charles R Greathouse IV, Oct 18 2022
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Klaus Brockhaus, Jan 21 2009
STATUS
approved