OFFSET
0,6
LINKS
Alois P. Heinz, Rows n = 0..140, flattened
FORMULA
A(n,k) = Sum_{t=1..n} t^k * A122843(n,t).
For fixed k, A(n,k) ~ n! * n * sum(t>=1, t^k*(t^2+t-1)/(t+2)!) = n! * n * ((Bell(k) - Bell(k+1) + sum(j=0..k, (-1)^j*(2^j*((2*k-j+1)/(j+1))-1) *Bell(k-j)*C(k,j)))*exp(1) - (-1)^k*(2^k-1)), where Bell(k) are Bell numbers A000110. - Vaclav Kotesovec, Sep 12 2013
EXAMPLE
A(3,2) = 32 = 9+5+5+5+5+3 = 3^2+4*(2^2+1^2)+3*1^2: (1,2,3), (1,3,2), (2,1,3), (2,3,1), (3,1,2), (3,2,1).
Square array A(n,k) begins:
: 0, 0, 0, 0, 0, 0, 0, ...
: 1, 1, 1, 1, 1, 1, 1, ...
: 3, 4, 6, 10, 18, 34, 66, ...
: 12, 18, 32, 66, 152, 378, 992, ...
: 60, 96, 186, 426, 1110, 3186, 9846, ...
: 360, 600, 1222, 2964, 8254, 25620, 86782, ...
: 2520, 4320, 9086, 22818, 66050, 214410, 765506, ...
MAPLE
A:= (n, k)-> add(`if`(n=t, 1, n!/(t+1)!*(t*(n-t+1)+1
-((t+1)*(n-t)+1)/(t+2)))*t^k, t=1..n):
seq(seq(A(n, d-n), n=0..d), d=0..12);
MATHEMATICA
A[n_, k_] := Sum[If[n == t, 1, n!/(t + 1)!*(t*(n - t + 1) + 1 - ((t + 1)*(n - t) + 1)/(t + 2))]* t^k, {t, 1, n}]; Table[Table[A[n, d - n], {n, 0, d}], {d, 0, 12}] // Flatten (* Jean-François Alcover, Dec 27 2013, translated from Maple *)
CROSSREFS
KEYWORD
nonn,tabl
AUTHOR
Alois P. Heinz, Sep 10 2013
STATUS
approved