OFFSET
1,2
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 1..200
FORMULA
a(n) = A093966(n, n).
From Vaclav Kotesovec, Nov 20 2012: (Start)
E.g.f.: x*(exp(x) - x)/(1-x)^2.
Recurrence: (n-1)*a(n) = n*(n+1)*a(n-1) - (n-1)*n*a(n-2) for n>2.
a(n) ~ n!*n*(e-1). (End)
MATHEMATICA
Rest[CoefficientList[Series[x(E^x-x)/(1-x)^2, {x, 0, 20}], x]* Range[0, 20]!] (* Vaclav Kotesovec, Nov 20 2012 *)
PROG
(PARI) my(x='x+O('x^66)); Vec(serlaplace(x*(exp(x)-x)/(1-x)^2)) \\ Joerg Arndt, May 11 2013
(Magma) [n le 2 select 4^(n-1) else n*((n+1)*Self(n-1) - (n-1)*Self(n-2))/(n-1): n in [1..30]]; // G. C. Greubel, Dec 29 2021
(Sage) [factorial(n)*( x*(exp(x) -x)/(1-x)^2 ).series(x, n+1).list()[n] for n in (1..30)] # G. C. Greubel, Dec 29 2021
CROSSREFS
KEYWORD
nonn
AUTHOR
Ralf Stephan, Apr 20 2004
EXTENSIONS
Name changed by Olivier Gérard, Aug 06 2016
STATUS
approved