OFFSET
0,3
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..506
MAPLE
b:= proc(n, i) option remember; `if`(n=0 or i=1, 1,
b(n, i-1) +b(n-i, min(i, n-i))/i)
end:
a:= n-> denom(b(n$2)):
seq(a(n), n=0..30);
MATHEMATICA
b[n_, i_] := b[n, i] = If[n==0 || i==1, 1, b[n, i-1]+b[n-i, Min[i, n-i]]/i];
a[n_] := Denominator[b[n, n]];
a /@ Range[0, 30] (* Jean-François Alcover, Apr 29 2020, after Alois P. Heinz *)
PROG
(PARI) a(n) = {my(s=0); forpart(p=n, s += 1/vecprod(Vec(p))); denominator(s); } \\ Michel Marcus, Apr 29 2020
CROSSREFS
KEYWORD
nonn,frac
AUTHOR
Alois P. Heinz, Dec 04 2018
STATUS
approved