OFFSET
0,6
FORMULA
E.g.f. of column k: exp(-k*x) / (1-x)^(k+1).
A(0,k) = A(1,k) = 1; A(n,k) = n*A(n-1,k) + k*(n-1)*A(n-2,k).
EXAMPLE
Square array begins:
1, 1, 1, 1, 1, 1, 1, ...
1, 1, 1, 1, 1, 1, 1, ...
2, 3, 4, 5, 6, 7, 8, ...
6, 11, 16, 21, 26, 31, 36, ...
24, 53, 88, 129, 176, 229, 288, ...
120, 309, 568, 897, 1296, 1765, 2304, ...
720, 2119, 4288, 7317, 11296, 16315, 22464, ...
PROG
(PARI) a(n, k) = n!*sum(j=0, n, (-k)^(n-j)*binomial(j+k, j)/(n-j)!);
CROSSREFS
KEYWORD
nonn,tabl
AUTHOR
Seiichi Manyama, Apr 24 2025
STATUS
approved
