|
%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
|
