OFFSET
0,3
COMMENTS
Sequences of the type b(m)+m*b(m-1), where b is a polygonal number:
LINKS
Bruno Berselli, Table of n, a(n) for n = 0..1000
Index entries for linear recurrences with constant coefficients, signature (4,-6,4,-1).
FORMULA
G.f.: x*(1+5*x+9*x^2)/(1-x)^4.
a(n) - a(-n) = A008531(n) for n>0.
MATHEMATICA
Table[n (5 n^2 - 8 n + 5)/2, {n, 0, 40}]
CoefficientList[Series[x (1 + 5 x + 9 x^2)/(1 - x)^4, {x, 0, 45}], x] (* Vincenzo Librandi, Aug 18 2013 *)
LinearRecurrence[{4, -6, 4, -1}, {0, 1, 9, 39}, 50] (* Harvey P. Dale, May 19 2017 *)
PROG
(Magma) [n*(5*n^2-8*n+5)/2: n in [0..40]];
(Magma) I:=[0, 1, 9, 39]; [n le 4 select I[n] else 4*Self(n-1)-6*Self(n-2)+4*Self(n-3)-Self(n-4): n in [1..45]]; // Vincenzo Librandi, Aug 18 2013
(PARI) a(n)=n*(5*n^2-8*n+5)/2 \\ Charles R Greathouse IV, Oct 07 2015
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Bruno Berselli, Jun 07 2013
STATUS
approved