OFFSET
0,2
LINKS
Muniru A Asiru, Table of n, a(n) for n = 0..180
Moussa Benoumhani, On Whitney numbers of Dowling lattices, Discrete Math. 159 (1996), no. 1-3, 13-33.
FORMULA
E.g.f.: exp(x + (exp(9*x) - 1)/9).
G.f.: 1/W(0), where W(k) = 1 - x - x/(1 - 9*(k+1)*x/W(k+1) ); (continued fraction). - Sergei N. Gladkovskii, Nov 07 2014
a(n) = exp(-1/9) * Sum_{k>=0} (9*k + 1)^n / (9^k * k!). - Ilya Gutkovskiy, Apr 16 2020
a(n) ~ 9^(n + 1/9) * n^(n + 1/9) * exp(n/LambertW(9*n) - n - 1/9) / (sqrt(1 + LambertW(9*n)) * LambertW(9*n)^(n + 1/9)). - Vaclav Kotesovec, Jun 26 2022
EXAMPLE
G.f. = 1 + 2*x + 13*x^2 + 143*x^3 + 1852*x^4 + 27563*x^5 + ...
MAPLE
seq(coeff(series(factorial(n)*exp(z+(1/9)*exp(9*z)-(1/9)), z, n+1), z, n), n = 0 .. 20); # Muniru A Asiru, Feb 24 2019
MATHEMATICA
With[{m=20, b=9}, CoefficientList[Series[Exp[x +(Exp[b*x]-1)/b], {x, 0, m}], x]*Range[0, m]!] (* G. C. Greubel, Feb 24 2019 *)
Table[Sum[Binomial[n, k] * 9^k * BellB[k, 1/9], {k, 0, n}], {n, 0, 20}] (* Vaclav Kotesovec, Apr 17 2020 *)
PROG
(PARI) Vec(serlaplace(exp(z + (exp(9*z) - 1)/9))) \\ Michel Marcus, Nov 07 2014
(Magma) m:=20; c:=9; R<x>:=PowerSeriesRing(Rationals(), m); b:=Coefficients(R!( Exp(x +(Exp(c*x)-1)/c) )); [Factorial(n-1)*b[n]: n in [1..m-1]]; // G. C. Greubel, Feb 24 2019
(Sage) m = 20; b=9; T = taylor(exp(x +(exp(b*x)-1)/b), x, 0, m); [factorial(n)*T.coefficient(x, n) for n in (0..m)] # G. C. Greubel, Feb 24 2019
CROSSREFS
KEYWORD
nonn
AUTHOR
EXTENSIONS
Name clarified by Muniru A Asiru, Feb 24 2019
STATUS
approved