OFFSET
1,3
COMMENTS
From Christopher J. Smyth, Jan 26 2018: (Start)
The sequence is defined by the recurrence formula below. This recurrence is very similar to that of the sequence b(n) = A001147(n-1), which satisfies b(1)=1 and, for n >= 2, b(n) = Sum_{i=1..floor((n-1)/2)} binomial(n, i) * b(i) * b(n-i) + B, where B = 0 (n odd), = (1/2)*binomial(n, n/2)*b(n/2)^2 (n even) [see formula of Walsh on A001147 page]. Removal of the factor 1/2 from the definition of B gives, for n >= 3, the formula below for a(n).
This sequence seems to have been defined in the mistaken belief that it had applications. In fact the applications stated on earlier versions of this page actually belonged to A001147 -- see my comment on the A001147 page.
(End)
LINKS
Michael De Vlieger, Table of n, a(n) for n = 1..403
FORMULA
a(1) = 1, a(2) = 1 and a(n) = Sum_{i=1..floor(n/2)} binomial(n, i) * a(i) * a(n-i) for n >= 3.
EXAMPLE
a(3) = 3*a(1)*a(2) = 3, a(4) = 4*a(1)*a(3) + 6*a(2)^2 = 18.
MATHEMATICA
Fold[Append[#1, Sum[Binomial[#2, i] #1[[i]] #1[[#2 - i]], {i, Floor[#2/2]}]] &, {1, 1}, Range[3, 21]] (* Michael De Vlieger, Dec 13 2017 *)
PROG
(MATLAB) function m = a(n); if n==1 m = 1; elseif n==2 m = 1; else m = 0; for i=1:floor(n/2); f1 = binomial(n, i); f2 = a(i); f3 = a(n-i); m = m + f1*f2*f3; end; end;
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Niko Brummer (niko.brummer(AT)gmail.com), Aug 08 2005
STATUS
approved