OFFSET
0,5
COMMENTS
FORMULA
Let D be the inverse Pascal operator: [D(f)]_m = Sum_{i=0..m} (-1)^(m - i)*binomial(m, i)*f(i). Then A(n, k) = [D^4(i^n*5^i)]_k = Sum_{j=0..n} Stirling2(n, j)*5^j*(k)_j.
For fixed k, the e.g.f. of A(n, k) in n is Sum_{n>=0} A(n, k)*u^n/n! = (1 + 5*(exp(u) - 1))^k.
A(n + 1, k) = k*(A(n, k) + 4*A(n, k - 1)) for k >= 1, with A(0, k) = 1 for k >= 0 and A(n, 0) = 0 for n > 0.
More generally, A_b(n, k) = Sum_{j=0..n} Stirling2(n, j)*b^j*(k)_j = Sum_{j=0..k} (-1)^(k - j)*binomial(k, j)*(b - 1)^(k - j)*b^j*j^n. For fixed k, the e.g.f. in n is Sum_{n>=0} A_b(n, k)*x^n/n! = (1 + b(e^x - 1))^k.
EXAMPLE
k | 0 1 2 3 4 5
n |--------------------------------------
0 | 1 1 1 1 1 1
1 | 0 5 10 15 20 25
2 | 0 5 60 165 320 525
3 | 0 5 160 1215 3920 9025
4 | 0 5 360 5565 35120 123525
5 | 0 5 760 21015 229520 1320025
A(2,2) = Stirling2(2, 1)*5*(2)_1 + Stirling2(2, 2)*25*(2)_2 = 1*5*2 + 1*25*2 = 60.
CROSSREFS
KEYWORD
nonn,tabl
AUTHOR
Dalton Heilig, May 05 2026
STATUS
approved
