login
Array read by antidiagonals of higher order Fubini numbers.
7

%I #25 May 12 2023 15:30:26

%S 1,1,3,1,4,13,1,5,23,75,1,6,36,175,541,1,7,52,342,1662,4683,1,8,71,

%T 594,4048,18937,47293,1,9,93,949,8444,57437,251729,545835,1,10,118,

%U 1425,15775,143783,950512,3824282,7087261,1,11,146,2040,27146,313920,2854261,17975438,65361237,102247563

%N Array read by antidiagonals of higher order Fubini numbers.

%H Alois P. Heinz, <a href="/A153278/b153278.txt">Antidiagonals n = 1..150, flattened</a>

%H Istvan Mezo, <a href="https://arxiv.org/abs/0812.4047">On powers of Stirling matrices</a>, arXiv:0812.4047 [math.CO], 2008.

%H K. A. Penson, P. Blasiak, G. Duchamp, A. Horzela, and A. I. Solomon, <a href="https://arxiv.org/abs/quant-ph/0312202">Hierarchical Dobinski-type relations via substitution and the moment problem</a>, arXiv:quant-ph/0312202, 2003; J.Phys. A: Math.Gen. 37 3475-3487 (2004).

%e The table on p.6 of Mezo begins:

%e ===========================================================

%e F_p,n|n=1|n=2|n=3.|.n=4.|..n=5.|....n=6.|.....n=7.|comment

%e ===========================================================

%e p=1..|.1.|.3.|.13.|..75.|..541.|...4683.|...47293.|.A000670

%e p=2..|.1.|.4.|.23.|.175.|.1662.|..18937.|..251729.|.A083355

%e p=3..|.1.|.5.|.36.|.342.|.4048.|..57437.|..950512.|.A099391

%e p=4..|.1.|.6.|.52.|.594.|.8444.|.143783.|.2854261.|.A363008

%e p=5..|.1.|.7.|.71.|.949.|15775.|.313920.|.7279795.|.A363009

%e ===========================================================

%p f:= proc(n) option remember; local k; if n<=1 then 1 else

%p add(binomial(n, k) *f(n-k), k=1..n) fi

%p end:

%p stirtr:= proc(a) proc(n) option remember;

%p add( a(k) *Stirling2(n,k), k=0..n)

%p end end:

%p F:= (p,n)-> (stirtr@@(p-1))(f)(n):

%p seq(seq(F(d-n, n), n=1..d-1), d=1..13); # _Alois P. Heinz_, Feb 02 2009

%t f[n_] := f[n] = If[n <= 1, 1, Sum[Binomial[n, k]*f[n-k], {k, 1, n}]];

%t stirtr[a_] := Module[{g}, g[n_] := g[n] = Sum[a[k]*StirlingS2[n, k], {k, 0, n}]; g];

%t F[p_, n_] := (Composition @@ Table[stirtr, {p-1}])[f][n];

%t Table[Table[F[d-n, n], {n, 1, d-1}], {d, 1, 13}] // Flatten (* _Jean-François Alcover_, Mar 30 2016, after _Alois P. Heinz_ *)

%Y Cf. A000670, A083355, A099391, A153277, A363008, A363009.

%K easy,nonn,tabl

%O 1,3

%A _Jonathan Vos Post_, Dec 22 2008

%E More terms from _Alois P. Heinz_, Feb 02 2009