OFFSET
0,3
LINKS
Seiichi Manyama, Table of n, a(n) for n = 0..387
K. Imatomi, M. Kaneko, E. Takeda, Multi-Poly-Bernoulli Numbers and Finite Multiple Zeta Values, J. Int. Seq. 17 (2014) # 14.4.5
M. Kaneko, Poly-Bernoulli numbers
Masanobu Kaneko, Poly-Bernoulli numbers, Journal de théorie des nombres de Bordeaux, 9 no. 1 (1997), Pages 221-228.
FORMULA
a(n) = numerator of Sum_{j=0..n} (-1)^(n+j) * j! * Stirling2(n, j) * (j+1)^(-k), for k = 4.
MAPLE
a:= (n, k)-> numer((-1)^n*add((-1)^m*m!*Stirling2(n, m)/(m+1)^k, m=0..n)):
seq(a(n, 4), n = 0..30);
MATHEMATICA
With[{k = 4}, Table[Numerator@ Sum[((-1)^(m + n))*m!*StirlingS2[n, m]*(m + 1)^(-k), {m, 0, n}], {n, 0, 17}]] (* Michael De Vlieger, Mar 18 2017 *)
PROG
(Magma)
A027647:= func< n, k | Numerator( (&+[(-1)^(j+n)*Factorial(j)*StirlingSecond(n, j)/(j+1)^k: j in [0..n]]) ) >;
[A027647(n, 4): n in [0..20]]; // G. C. Greubel, Aug 02 2022
(SageMath)
def A027647(n, k): return numerator( sum((-1)^(n+j)*factorial(j)*stirling_number2(n, j)/(j+1)^k for j in (0..n)) )
[A027647(n, 4) for n in (0..20)] # G. C. Greubel, Aug 02 2022
CROSSREFS
KEYWORD
sign,frac
AUTHOR
STATUS
approved