OFFSET
2,2
LINKS
R. J. Mathar, Topologically Distinct Sets of Non-intersecting Circles in the Plane, arXiv:1603.00077 [math.CO], 2016, row sums Table 7.
MAPLE
A280786 := proc(N)
if N < 2 then
0;
else
add(A280787(N, f), f=1..N-1) ;
end if;
end proc:
A280787 := proc(N, f)
option remember ;
local Npr, ct ;
if f = N then
return 0;
elif f = N-1 then
return 1;
elif f = 1 then
else
ct := 0 ;
for Npr from 1 to N-1 do
ct := ct+procname(Npr, 1)*A033185(N-Npr, f-1) ;
end do:
ct ;
end if;
end proc:
seq(A280786(n), n=2..30) ; # R. J. Mathar, Mar 06 2017
MATHEMATICA
a81[n_] := a81[n] = If[n <= 1, n, Sum[a81[n - j]*DivisorSum[j, #1*a81[#1] &], {j, n - 1}]/(n - 1)];
A027852[n_] := Module[{dh = 0, np}, For[np = 0, np <= n, np++, dh = a81[np]*a81[n - np] + dh]; If[EvenQ[n], dh = a81[n/2] + dh]; dh/2];
t[n_] := t[n] = Module[{d, j}, If[n == 1, 1, Sum[Sum[d*t[d], {d, Divisors[j]}]*t[n - j], {j, 1, n - 1}]/(n - 1)]];
b[1, 1, 1] = 1;
b[n_, i_, p_] := b[n, i, p] = If[p > n, 0, If[n == 0, 1, If[Min[i, p] < 1, 0, Sum[b[n - i*j, i - 1, p - j]*Binomial[t[i] + j - 1, j], {j, 0, Min[n/i, p]}]]]]; A033185[n_, k_] := b[n, n, k];
A280787[n_, f_] := A280787[n, f] = Module[{ct}, Which[f == n, Return[0], f == n - 1, Return[1], f == 1, Return[A280786[n - 1] + A280788[n - 2]], True, ct = 0; Do[ct += A280787[np, 1]*A033185[n - np, f - 1], {np, 1, n - 1}]]; ct];
Table[A280786[n], {n, 2, 30}] (* Jean-François Alcover, Nov 23 2017, after R. J. Mathar and Alois P. Heinz *)
CROSSREFS
KEYWORD
nonn
AUTHOR
N. J. A. Sloane, Jan 20 2017
STATUS
approved