OFFSET
0,3
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..500
Loïc Foissy, Algebraic structures on typed decorated rooted trees, arXiv:1811.07572 [math.RA], 2018.
FORMULA
a(n) ~ c * d^n / n^(3/2), where d = 24.654188324989308481297681284616076174..., c = 0.044677658292795889027452779147731609565... . - Vaclav Kotesovec, Aug 26 2014
G.f. A(x) satisfies: A(x) = x*exp(9*Sum_{k>=1} A(x^k)/k). - Ilya Gutkovskiy, Mar 20 2018
MAPLE
with(numtheory):
a:= proc(n) option remember; `if`(n<2, n, (add(add(d*
a(d), d=divisors(j))*a(n-j)*9, j=1..n-1))/(n-1))
end:
seq(a(n), n=0..25);
MATHEMATICA
a[n_] := a[n] = If[n<2, n, Sum[Sum[d*a[d], {d, Divisors[j]}]*a[n-j]*9, {j, 1, n-1}]/(n-1)];
Table[a[n], {n, 0, 25}] (* Jean-François Alcover, Feb 23 2019, from Maple *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Alois P. Heinz, Aug 19 2014
STATUS
approved