A111579 Triangle A(r,c) read by rows, which contains the row sums of the triangle T(n,k)= T(n-1,k-1)+((c-1)*k+1)*T(n-1,k) in column c.
1, 1, 1, 1, 2, 1, 1, 4, 2, 1, 1, 8, 5, 2, 1, 1, 16, 15, 6, 2, 1, 1, 32, 52, 24, 7, 2, 1, 1, 64, 203, 116, 35, 8, 2, 1, 1, 128, 877, 648, 214, 48, 9, 2, 1, 1, 256, 4140, 4088, 1523, 352, 63, 10, 2, 1, 1, 512, 21147, 28640, 12349, 3008, 536, 80, 11, 2, 1
Offset: 0
Programs
-
Maple
T := proc(n,k,c) if k < 0 or k > n then 0 ; elif n <= 1 then 1; else procname(n-1,k-1,c)+((c-1)*k+1)*procname(n-1,k,c) ; fi; end: A111579 := proc(r,c) local n; if c = 0 then 1 ; else n := r-c ; add( T(n,k,c),k=0..n) ; end if; end: seq(seq(A111579(r,c),c=0..r),r=0..10) ; # R. J. Mathar, Oct 30 2009
-
Mathematica
T[n_, k_, c_] := T[n, k, c] = If[k < 0 || k > n, 0, If[n <= 1, 1, T[n-1, k-1, c] + ((c-1)*k+1)*T[n-1, k, c]]]; A111579[r_, c_] := Module[{n}, If[c == 0, 1, n = r - c; Sum[T[n, k, c], {k, 0, n}]]]; Table[A111579[r, c], {r, 0, 10}, {c, 0, r}] // Flatten (* Jean-François Alcover, Aug 01 2023, after R. J. Mathar *)
Formula
Extensions
Edited by R. J. Mathar, Oct 30 2009
Comments