A053530 Expansion of e.g.f.: exp(-x - x^2/2 + x*exp(x)).

1, 0, 1, 3, 7, 35, 171, 847, 5041, 32643, 223705, 1659581, 13182159, 110802133, 984241363, 9212696235, 90477239521, 929604133343, 9969157068273, 111329454692485, 1291932988047775, 15550838026589061, 193833398512358011, 2498039016973836491

Offset: 0

Author

N. J. A. Sloane, Jan 16 2000

Keywords

• nonn

Comments

The number of simple labeled graphs on n nodes whose connected components are stars. - Geoffrey Critzer, Dec 10 2011

Equivalently, the number of minimal edge covers of the complete graph K_n. - Andrew Howroyd, Aug 04 2017

References

• R. P. Stanley, Enumerative Combinatorics, Cambridge, Vol. 2, 1999; see Problem 5.15(b).

Links

• Vincenzo Librandi, Table of n, a(n) for n = 0..200
• Vaclav Kotesovec, Asymptotic solution of the equations using the Lambert W-function
• Vladimir Kruchinin and D. V. Kruchinin, Composita and their properties, arXiv:1103.2582 [math.CO], 2011-2013.
• Eric Weisstein's World of Mathematics, Complete Graph
• Eric Weisstein's World of Mathematics, Minimal Edge Cover

Crossrefs

Cf. A000248, A210655.

Programs

Magma

m:=30; R:=PowerSeriesRing(Rationals(), m); b:=Coefficients(R!( Exp(-x -x^2/2 +x*Exp(x)) )); [Factorial(n-1)*b[n]: n in [1..m]]; // G. C. Greubel, May 15 2019

Mathematica

nn = 30; a = x Exp[x]; Range[0, nn]! CoefficientList[Series[Exp[a - x^2/2! - x], {x, 0, nn}], x] (* Geoffrey Critzer, Dec 10 2011 *)
CoefficientList[Series[Exp[-x - x^2/2 + x Exp[x]], {x, 0, 30}], x] Range[0, 30]! (* Eric W. Weisstein, Aug 10 2017 *)
Table[n! Sum[1/k! (Binomial[k, n-k] 2^(k-n) (-1)^k + Sum[Binomial[k, j] Sum[j^(i-j)/(i-j)! Binomial[k-j, n-i-k+j] 2^(i-j+k-n) (-1)^(k-j), {i, j, n-k+j}], {j, k}]), {k, n}], {n, 30}] (* Eric W. Weisstein, Aug 10 2017 *)

Maxima

a(n):=n!*sum((binomial(k,n-k)*2^(k-n)*(-1)^k +sum(binomial(k,j) *sum(j^(i-j)/(i-j)!*binomial(k-j,n-i-k+j)*(1/2)^(n-i-k+j)*(-1)^(k-j),i,j,n-k+j),j,1,k))/k!,k,1,n); /* Vladimir Kruchinin, Sep 10 2010 */

PARI

x='x+O('x^30); Vec(serlaplace(exp(-x-1/2*x^2+x*exp(x)))) \\ Altug Alkan, Aug 10 2017

Sage

m = 30; T = taylor(exp(-x -x^2/2 +x*exp(x)), x, 0, m); [factorial(n)*T.coefficient(x, n) for n in (0..m)] # G. C. Greubel, May 15 2019

Formula

a(n) = n!*Sum_{k=1..n} (1/k!)*( binomial(k, n-k)*2^(k-n)*(-1)^k + Sum_{j=1..k} binomial(k,j)* (Sum_{i=j..n-k+j} (j^(i-j)/(i-j)! * binomial(k-j,n-i-k+j)*(1/2)^(n-i-k+j)*(-1)^(k-j)) ) ), n>0. - Vladimir Kruchinin, Sep 10 2010

From Vaclav Kotesovec, Aug 06 2014: (Start)

a(n) ~ n^n / (exp(r^2/2 + n*r/(1+r)) * r^n * sqrt(r^2*(1+r)/n + 2+r-1/(1+r))), where r is the root of the equation r*(exp(r)*(1+r)-1-r) = n.

(a(n)/n!)^(1/n) ~ exp(1/(2*LambertW(sqrt(n)/2)))/(2*LambertW(sqrt(n)/2)).

(End)