OFFSET
0,2
COMMENTS
Numerators are A133002.
LINKS
Daniel Suteu, Table of n, a(n) for n = 0..200
Hector Blandin and Rafael Diaz, Compositional Bernoulli numbers, arXiv:0708.0809 [math.CO], 2007-2008, p. 11, 1st table.
FORMULA
a(n) = denominator(f(n) * n!), where f(0) = 1, f(n) = -Sum_{k=0..n-1} f(k) / ((n-k+1)!)^2. - Daniel Suteu, Feb 23 2018
E.g.f. for fractions: x / (BesselI(0,2*sqrt(x)) - 1). - Ilya Gutkovskiy, Sep 01 2021
EXAMPLE
1, -1/4, 5/72, -1/48, 139/21600, -1/540, 859/2540160, 71/483840, -9769/36288000 (corrected by Daniel Suteu, Feb 24 2018).
MATHEMATICA
f[0] = 1; f[n_] := f[n] = -Sum[f[k]/((n - k + 1)!)^2, {k, 0, n - 1}]; a[n_] := Denominator[f[n]*n!]; Table[a[n], {n, 0, 19}] (* Jean-François Alcover, Feb 25 2018, after Daniel Suteu *)
CROSSREFS
KEYWORD
nonn,frac
AUTHOR
Jonathan Vos Post, Aug 09 2007
EXTENSIONS
Terms beyond a(8) from Daniel Suteu, Feb 24 2018
STATUS
approved