login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A104460 Number of hierarchical orderings for n unlabeled elements with 2 possible classes for levels l>=2. 13
1, 4, 13, 46, 154, 533, 1802, 6137, 20729, 69971, 235193, 789000, 2639004, 8807811, 29327841, 97456878, 323206002, 1069923013, 3535612108, 11664423298, 38422208659, 126374059558, 415069188175, 1361443135562, 4459861400156, 14591869576268, 47686017637926 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,2

COMMENTS

Consider a hierarchical ordering of n unlabeled elements into groups as defined in A034691. In addition assume that each level l with l >= 2 can fall into one of two classes A and B. Let | denote a separator among different groups and let : denote a separator between levels. Furthermore, let * denote an unlabeled element which is written as "a" if it falls into class A and as "b" if it falls into class B. As an example with n=4 one can have *|*:ab. In this example one has two groups, where the second group has tree elements, one on level l=1 and two on level l=2. One of the two elements on l=2 belongs to class A, the other to class B.

LINKS

Alois P. Heinz, Table of n, a(n) for n = 1..600

N. J. A. Sloane and Thomas Wieder, The Number of Hierarchical Orderings, Order 21 (2004), 83-89.

FORMULA

G.f.: 1 + Sum_{n>=1} a(n) * x^n = 1 / Product_{n>=1} (1-x^n)^(3^(n-1)).

A104460 is the Euler transform of powers of 3 [1, 3, 9, 27, 81, ...].

a(n) ~ exp(2*sqrt(n/3) - 1/6 + c/3) * 3^(n - 1/4) / (2*sqrt(Pi)*n^(3/4)), where c = Sum_{j>=2} 1/(j * (3^(j-1) - 1)). - Vaclav Kotesovec, Apr 12 2021

EXAMPLE

For n=3 there are 13 orderings:

*|*|*; *|**; *|*:a; *|*:b; ***; **|a; *:aa; *:a:a; **|b; *:bb; *:b:b; *:a:b; *:b:a.

MAPLE

etr:= proc(p) local b; b:=proc(n) option remember; `if`(n=0, 1, add(add(d*p(d), d=numtheory[divisors](j)) *b(n-j), j=1..n)/n) end end: a:= etr(n-> 3^(n-1)): seq(a(n), n=1..30); # Alois P. Heinz, Sep 08 2008

MATHEMATICA

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]; a = etr[Function[3^(#-1)]]; Table[a[n], {n, 1, 30}] (* Jean-Fran├žois Alcover, Mar 05 2015, after Alois P. Heinz *)

nmax = 30; Rest[CoefficientList[Series[Product[1/(1 - x^k)^(3^(k-1)), {k, 1, nmax}], {x, 0, nmax}], x]] (* Vaclav Kotesovec, Apr 12 2021 *)

CROSSREFS

Cf. A034691.

Sequence in context: A096353 A034553 A320743 * A095128 A149433 A047154

Adjacent sequences:  A104457 A104458 A104459 * A104461 A104462 A104463

KEYWORD

nonn

AUTHOR

Thomas Wieder, Mar 09 2005

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified January 28 14:26 EST 2022. Contains 350656 sequences. (Running on oeis4.)