|
|
A003579
|
|
Dowling numbers: e.g.f. exp(x + (exp(b*x) - 1)/b), with b=7.
|
|
13
|
|
|
1, 2, 11, 99, 1066, 13283, 190933, 3117900, 56729565, 1132679479, 24564972756, 574431351673, 14394977015245, 384489778509034, 10894501505088695, 326149933663962479, 10280153573323314858
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,2
|
|
LINKS
|
|
|
FORMULA
|
E.g.f.: exp(x + (exp(7*x) - 1)/7).
a(n) = exp(-1/7) * Sum_{k>=0} (7*k + 1)^n / (7^k * k!). - Ilya Gutkovskiy, Apr 16 2020
a(n) ~ 7^(n + 1/7) * n^(n + 1/7) * exp(n/LambertW(7*n) - n - 1/7) / (sqrt(1 + LambertW(7*n)) * LambertW(7*n)^(n + 1/7)). - Vaclav Kotesovec, Jun 26 2022
|
|
MATHEMATICA
|
With[{m=20, b=7}, 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] * 7^k * BellB[k, 1/7], {k, 0, n}], {n, 0, 20}] (* Vaclav Kotesovec, Apr 17 2020 *)
|
|
PROG
|
(PARI) Vec(serlaplace( exp(z + 1/7 * exp(7 * z) - 1/7) ) ) \\ Joerg Arndt, Feb 24 2019
(Magma) m:=20; c:=7; 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=7; 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
(GAP) b:=7;; a:=[1, 2];; for n in [3..20] do a[n]:=2*a[n-1]+Sum([0..n-3], i->Binomial(n-2, i)*b^(n-2-i)*a[i+1]); od; Print(a); # Muniru A Asiru, Apr 10 2019
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|