OFFSET
0,3
COMMENTS
Binomial transform of A084900. Third binomial transform of heptagonal numbers A000566. Fourth binomial transform of (0,1,5,0,0,0,...).
Coefficients in the hypergeometric series identity 1 - 13*x/(x + 12) + 108*x*(x - 1)/((x + 12)*(x + 16)) - 736*x*(x - 1)*(x - 2)/((x + 12)*(x + 16)*(x + 20)) + ... = 0, valid in the half-plane Re(x) > 0. Cf. A276289 and A077616. - Peter Bala, May 30 2019
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..1000
Index entries for linear recurrences with constant coefficients, signature (12,-48,64).
FORMULA
G.f.: x*(1+x)/(1-4*x)^3.
E.g.f.: x*(2 + 5*x)*exp(4*x)/2. - G. C. Greubel, Jun 06 2019
a(n) = 12*a(n-1)-48*a(n-2)+64*a(n-3). - Wesley Ivan Hurt, May 28 2021
MATHEMATICA
Table[2^(2*n-5)*n*(5*n+3), {n, 0, 30}] (* G. C. Greubel, Jun 06 2019 *)
LinearRecurrence[{12, -48, 64}, {0, 1, 13}, 30] (* or *) CoefficientList[ Series[-((x (1+x))/(-1+4 x)^3), {x, 0, 30}], x] (* Harvey P. Dale, Jul 14 2021 *)
PROG
(PARI) vector(30, n, n--; 2^(2*n-5)*n*(5*n+3)) \\ G. C. Greubel, Jun 06 2019
(Magma) [2^(2*n-5)*n*(5*n+3): n in [0..30]]; // G. C. Greubel, Jun 06 2019
(Sage) [2^(2*n-5)*n*(5*n+3) for n in (0..30)] # G. C. Greubel, Jun 06 2019
(GAP) List([0..30], n-> 2^(2*n-5)*n*(5*n+3)) # G. C. Greubel, Jun 06 2019
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Paul Barry, Jun 10 2003
STATUS
approved