OFFSET
0,2
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..1000
Index entries for linear recurrences with constant coefficients, signature (8,-28,56,-70,56,-28,8,-1).
FORMULA
a(n) = n*Stirling2(4 + n, 1 + n).
-a(-n-4) = (n+4)*abs(Stirling1(n+3, n)) for n >= 0.
-a(-n-4) = a(n) + 5*binomial(n+4, 5)*(n+2) for n >= 0.
From Colin Barker, Nov 21 2017: (Start)
G.f.: 15*x*(1 + 4*x + 2*x^2) / (1 - x)^8.
a(n) = (1/48)*(n*(2 + 3*n + n^2)^2*(12 + 7*n + n^2)).
a(n) = 8*a(n-1) - 28*a(n-2) + 56*a(n-3) - 70*a(n-4) + 56*a(n-5) - 28*a(n-6) + 8*a(n-7) - a(n-8) for n>7.
(End)
MATHEMATICA
LinearRecurrence[{8, -28, 56, -70, 56, -28, 8, -1}, {0, 15, 180, 1050, 4200, 13230, 35280, 83160}, 32]
Table[n*StirlingS2[4 + n, 1 + n], {n, 0, 50}] (* G. C. Greubel, Nov 20 2017 *)
PROG
(PARI) for(n=0, 30, print1(n*stirling(n+4, n+1, 2), ", ")) \\ G. C. Greubel, Nov 20 2017
(Magma) [0] cat [((n + 1)/2)*(n + 2)*Factorial(n+4)/(Factorial(4)*Factorial(n-1)): n in [1..30]]; // G. C. Greubel, Nov 20 2017
(PARI) concat(0, Vec(15*x*(1 + 4*x + 2*x^2) / (1 - x)^8 + O(x^40))) \\ Colin Barker, Nov 21 2017
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Peter Luschny, Oct 20 2017
STATUS
approved