OFFSET
0,3
LINKS
Robert Israel, Table of n, a(n) for n = 0..414
FORMULA
(n-1) a(n) + (n^2-2n) a(n-1) = Sum_{j=1..n-1} C(n,j) a(j) a(n-j). - Robert Israel, Jul 13 2015
For n > 1, a(n) = n!*Sum_{k=1..n-1} a(n-k)*(-1)^(k+1)/(k*(n-k)!). - Tani Akinari, Nov 25 2017
a(n) ~ n! * (-1)^n / (n * (log(n))^2) * (1 - 2*(gamma + 1) / log(n) + (3*gamma^2 + 6*gamma + 3 - Pi^2/2) / (log(n))^2), where gamma is the Euler-Mascheroni constant (A001620). - Vaclav Kotesovec, Nov 29 2017
MAPLE
G:=x/(1-ln(1+x)): Gser:=series(G, x=0, 28): 0, seq(n!*coeff(Gser, x^n), n=1..25); # yields the signed sequence
A:= proc(n) option remember;
(-(n^2-2*n)*procname(n-1) + add(binomial(n, j)*procname(j)*procname(n-j), j=1..n-1))/(n-1) end proc:
A(0):= 0: A(1):= 1:
map(A, [$0..50]); # Robert Israel, Jul 13 2015
MATHEMATICA
Range[0, 25]! CoefficientList[Series[x/(1 - Log[1 + x]), {x, 0, 25}], x] (* Vincenzo Librandi, Jul 06 2015 *)
PROG
(Maxima)
a[n]:=if n<2 then n else n!*sum(a[n-k]*(-1)^(k+1)/(k*
(n-k)!), k, 1, n-1); makelist(a[n], n, 0, 100); /* Tani Akinari, Nov 25 2017 */
(PARI) x='x+O('x^30); concat([0], Vec(serlaplace(x/(1 - log(1+x))))) \\ G. C. Greubel, Aug 30 2018
(Magma) m:=30; R<x>:=PowerSeriesRing(Rationals(), m); b:=Coefficients(R!(x/(1 - Log(1+x)))); [0] cat [Factorial(n)*b[n]: n in [1..m-1]]; // G. C. Greubel, Aug 30 2018
CROSSREFS
KEYWORD
sign
AUTHOR
Roger L. Bagula, Jun 26 2005
STATUS
approved