STATUS
editing
approved
editing
approved
proposed
approved
editing
proposed
b[n_] := b[n] = If[n <= 1, n, (Sum[Sum[d*b[d], {d, Divisors[j]}]*b[n - j], {j, 1, n - 1}])/(n - 1)];
a[n_] := a[n] = If[n > 2, 1, 0] + b[n] - (Sum[b[k]*b[n - k], {k, 0, n}] - If[Mod[n, 2] == 0, b[n/2], 0])/2;
Array[a, 40] (* Jean-François Alcover, Mar 21 2017, translated from Maple *)
approved
editing
editing
approved
with (numtheory):
b:= proc(n) option remember; local d, j; `if` (n<=1, n,
(add (add (d*b(d), d=divisors(j)) *b(n-j), j=1..n-1))/(n-1))
(add (b(k)*b(n-k), k=0..n) -`if`(irem(n, 2)=0, b(n/2), 0))/2
seq (a(n), n=1..40);
approved
editing
editing
approved
Ta(5) = 4: .o-o-o. .o-o-o. .o-o-o. .o-o-o.
The labeled version is A215851.