OFFSET
1,4
LINKS
Alois P. Heinz, Table of n, a(n) for n = 1..500
FORMULA
G.f.: B(x)-B(x)^2, where B(x) is g.f. for A005753.
a(n) ~ c * d^n / n^(5/2), where d = A246312 = 5.249032491228170579164952216..., c = 0.17807103914078424643862998... . - Vaclav Kotesovec, Aug 25 2014
MAPLE
b:= proc(n, i) option remember; `if`(n=0, 1, `if`(i<1, 0,
add(binomial(2*b(i-1$2), j)*b(n-i*j, i-1), j=0..n/i)))
end:
a:= n-> b(n-1$2)-add(b(j-1$2)*b(n-j-1$2), j=1..n-1):
seq(a(n), n=1..35); # Alois P. Heinz, Aug 01 2013
MATHEMATICA
b[n_, i_] := b[n, i] = If[n == 0, 1, If[i<1, 0, Sum[Binomial[2*b[i-1, i-1], j]*b[n - i*j, i-1], {j, 0, n/i}]]] // FullSimplify; a[n_] := b[n-1, n-1] - Sum[b[j-1, j-1]*b[n-j-1, n-j-1], {j, 1, n-1}]; Table[a[n], {n, 1, 35}] (* Jean-François Alcover, Feb 24 2015, after Alois P. Heinz *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Vladeta Jovovic, Feb 10 2005
STATUS
approved