OFFSET
0,5
FORMULA
E.g.f.: 1 / (exp(x) + exp(y) - exp(x+y))^4.
A(n,k) = A(k,n).
A(n,k) = Sum_{j=0..min(n,k)} (j!)^2 * binomial(j+3,3) * Stirling2(n,j) * Stirling2(k,j).
EXAMPLE
Square array begins:
1, 0, 0, 0, 0, 0, ...
0, 4, 4, 4, 4, 4, ...
0, 4, 44, 124, 284, 604, ...
0, 4, 124, 1084, 5164, 19804, ...
0, 4, 284, 5164, 48044, 313804, ...
0, 4, 604, 19804, 313804, 3281404, ...
PROG
(PARI) a(n, k) = sum(j=0, min(n, k), j!^2*binomial(j+3, 3)*stirling(n, j, 2)*stirling(k, j, 2));
CROSSREFS
KEYWORD
nonn,tabl
AUTHOR
Seiichi Manyama, Apr 04 2025
STATUS
approved
