OFFSET
1,8
LINKS
Alois P. Heinz, Table of n, a(n) for n = 1..1000
Eric Weisstein's World of Mathematics, Caterpillar Graph
MAPLE
with(numtheory):
b:= proc(n) option remember; `if`(n<=1, n,
(add(add(d*b(d), d=divisors(j))*b(n-j), j=1..n-1))/(n-1))
end:
a:= n-> b(n) -(add(b(k) *b(n-k), k=0..n)-`if`(irem(n, 2)=0,
b(n/2), 0))/2 -ceil(2^(n-4) + 2^(iquo(n-2, 2)-1)):
seq(a(n), n=1..40); # Alois P. Heinz, May 18 2013
MATHEMATICA
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_] := b[n] - (Sum[b[k]*b[n-k], {k, 0, n}] - If[ Mod[n, 2] == 0, b[n/2], 0])/2 - Ceiling[2^(n-4) + 2^(Quotient[n-2, 2] - 1)]; Table[a[n], {n, 1, 40}] (* Jean-François Alcover, Feb 19 2016, after Alois P. Heinz *)
CROSSREFS
KEYWORD
nonn
AUTHOR
EXTENSIONS
a(14) and up from Eric W. Weisstein, Jul 17 2004.
STATUS
approved