OFFSET
0,2
COMMENTS
a(n) is also the number of acyclic orientations of the complete bipartite graph K_{4,n}. - Vincent Pilaud, Sep 16 2020
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 0..500
K. Imatomi, M. Kaneko, and E. Takeda, Multi-Poly-Bernoulli Numbers and Finite Multiple Zeta Values, J. Int. Seq. 17 (2014) # 14.4.5.
Ken Kamano, Sums of Products of Bernoulli Numbers, Including Poly-Bernoulli Numbers, J. Int. Seq. 13 (2010), 10.5.2.
Ken Kamano, Sums of Products of Poly-Bernoulli Numbers of Negative Index, Journal of Integer Sequences, Vol. 15 (2012), #12.1.3.
Masanobu Kaneko, Poly-Bernoulli numbers, Journal de théorie des nombres de Bordeaux, 9 no. 1 (1997), Pages 221-228.
Hiroyuki Komaki, Relations between Multi-Poly-Bernoulli numbers and Poly-Bernoulli numbers of negative index, arXiv:1503.04933 [math.NT], 2015.
Index entries for linear recurrences with constant coefficients, signature (14,-71,154,-120).
FORMULA
a(n) = 24*5^n -36*4^n +14*3^n -2^n. - Vladeta Jovovic, Nov 14 2003
G.f.: (1+4*x)*(1-x)^2/((1-2*x)*(1-3*x)*(1-4*x)*(1-5*x)).
E.g.f.: 24*exp(5*x) - 36*exp(4*x) + 14*exp(3*x) - exp(2*x). - G. C. Greubel, Feb 07 2018
MAPLE
a:= (n, k) -> (-1)^n*sum((-1)^j*j!*Stirling2(n, j)/(j+1)^k, j=0..n);
seq(a(n, -4), n=0..30);
MATHEMATICA
Table[24*5^n -36*4^n +14*3^n -2^n, {n, 0, 30}] (* G. C. Greubel, Feb 07 2018 *)
LinearRecurrence[{14, -71, 154, -120}, {1, 16, 146, 1066}, 30] (* Harvey P. Dale, Nov 20 2019 *)
PROG
(Magma) [24*5^n-36*4^n+14*3^n-2^n: n in [0..30]]; // Vincenzo Librandi, Jul 17 2011
(PARI) Vec((1+4*x)*((1-x)^2)/((1-2*x)*(1-3*x)*(1-4*x)*(1-5*x)) + O(x^30)) \\ Michel Marcus, Feb 13 2015
(SageMath) [24*5^n -36*4^n +14*3^n -2^n for n in (0..30)] # G. C. Greubel, Aug 02 2022
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
STATUS
approved