login
Number of 4-level rooted trees with n leaves.
17

%I #40 Feb 18 2023 14:45:11

%S 1,1,4,10,30,75,206,518,1344,3357,8429,20759,51044,123973,299848,

%T 719197,1716563,4070800,9607797,22555988,52718749,122655485,284207304,

%U 655894527,1508046031,3454808143,7887768997,17949709753,40719611684,92096461012,207697731344

%N Number of 4-level rooted trees with n leaves.

%H Alois P. Heinz, <a href="/A007713/b007713.txt">Table of n, a(n) for n = 0..1000</a>

%H P. J. Cameron, <a href="http://www.cs.uwaterloo.ca/journals/JIS/VOL3/groups.html">Sequences realized by oligomorphic permutation groups</a>, J. Integ. Seqs. Vol. 3 (2000), #00.1.5.

%H B. A. Huberman and T. Hogg, <a href="https://doi.org/10.1016/0167-2789(86)90308-1">Complexity and adaptation</a>, Evolution, games and learning (Los Alamos, N.M., 1985). Phys. D 22 (1986), no. 1-3, 376-384.

%H N. J. A. Sloane, <a href="/transforms.txt">Transforms</a>

%H <a href="/index/Ro#rooted">Index entries for sequences related to rooted trees</a>

%F Euler transform applied thrice to all-1's sequence.

%e From _Gus Wiseman_, Oct 11 2018: (Start)

%e Also the number of multiset partitions of multiset partitions of integer partitions of n. For example, the a(1) = 1 through a(4) = 30 multiset partitions are:

%e ((1)) ((2)) ((3)) ((4))

%e ((11)) ((12)) ((13))

%e ((1)(1)) ((111)) ((22))

%e ((1))((1)) ((1)(2)) ((112))

%e ((1)(11)) ((1111))

%e ((1))((2)) ((1)(3))

%e ((1))((11)) ((2)(2))

%e ((1)(1)(1)) ((1)(12))

%e ((1))((1)(1)) ((2)(11))

%e ((1))((1))((1)) ((1)(111))

%e ((11)(11))

%e ((1))((3))

%e ((2))((2))

%e ((1))((12))

%e ((1)(1)(2))

%e ((2))((11))

%e ((1))((111))

%e ((1)(1)(11))

%e ((11))((11))

%e ((1))((1)(2))

%e ((2))((1)(1))

%e ((1))((1)(11))

%e ((1)(1)(1)(1))

%e ((11))((1)(1))

%e ((1))((1))((2))

%e ((1))((1))((11))

%e ((1))((1)(1)(1))

%e ((1)(1))((1)(1))

%e ((1))((1))((1)(1))

%e ((1))((1))((1))((1))

%e (End)

%p with(numtheory): etr:= proc(p) local b; b:=proc(n) option remember; local d,j; if n=0 then 1 else add(add(d*p(d), d=divisors(j)) *b(n-j), j=1..n)/n fi end end: b0:= etr(1): b1:= etr(b0): a:= etr(b1): seq(a(n), n=0..30); # _Alois P. Heinz_, Sep 08 2008

%t i[ n_, m_ ] := 1 /; m==1 || n==0; i[ n_, m_ ] := (i[ n, m ]=1/n Sum[ i[ k, m ] Plus @@ ((# i[ #, m-1 ])& /@ Divisors[ n-k ]), {k, 0, n-1} ]) /; n>0 && m>1

%t etr[p_] := Module[{b}, b[n_] := b[n] = If[n == 0, 1, Sum[Sum[d*p[d], {d, Divisors[ j]}]*b[n-j], {j, 1, n}]/n]; b]; b0 = etr[Function[1]]; b1 = etr[b0]; a = etr[b1]; Table[a[n], {n, 1, 30}] (* _Jean-François Alcover_, Mar 05 2015, after _Alois P. Heinz_ *)

%Y Column k=4 of A290353.

%Y Main diagonal of A055886.

%Y Cf. A001970, A047968, A050342, A089259, A141268, A258466, A261049, A319066, A320328, A320330, A320331.

%K easy,nonn

%O 0,3

%A _N. J. A. Sloane_.