Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
A038038
Number of labeled dyslexic rooted compound windmills with n nodes.
2
1, 2, 9, 64, 635, 8142, 128394, 2405784, 52237296, 1289735280, 35680244490, 1093124356560, 36734913680244, 1343439671120904, 53111942243599680, 2256996523715952480, 102589277119862240160, 4966450613410762025280, 255110144003301741908832
OFFSET
1,2
FORMULA
Divides by n and shifts left under "DIJ" (bracelet, indistinct, labeled) transform.
E.g.f: series reversion of 2*x/(2 + x + x^2/2 - log(1-x)). - Andrew Howroyd, Sep 19 2018
a(n) ~ sqrt(s) * (2 - s^2)^(n + 1/2) * n^n / (sqrt(2 - 2*s + s^2) * n * 2^n * (1-s)^(n - 1/2) * exp(n)), where s = 0.7579492001963653206343844374776312472163... is the root of the equation 4 - 6*s - s^2 + s^3 - 2*(1-s)*log(1-s) = 0. - Vaclav Kotesovec, Apr 21 2020
MATHEMATICA
m = 20;
CoefficientList[InverseSeries[2*x/(2 + x + x^2/2 - Log[1 - x]) + O[x]^m], x]*Range[0, m - 1]! // Rest (* Jean-François Alcover, Sep 08 2019 *)
PROG
(PARI) Vec(serlaplace(serreverse(2*x/(2 + x + x^2/2 - log(1-x + O(x^20)))))) \\ Andrew Howroyd, Sep 19 2018
CROSSREFS
Cf. A032289.
Sequence in context: A000169 A112319 A232552 * A368790 A048801 A152915
KEYWORD
nonn,eigen
EXTENSIONS
Terms a(17) and beyond from Andrew Howroyd, Sep 19 2018
STATUS
approved