OFFSET
0,4
COMMENTS
Let G(x) = (-x)/LambertW(-x), then A(x)^A(x) = 1/G(x) where G(x)^G(x) = 1/exp(x).
a(n) = 0 (mod 3) when n = 6*k+5, k>=0, otherwise a(n) = 1 (mod 3).
LINKS
Paul D. Hanna, Table of n, a(n) for n = 0..200
FORMULA
E.g.f.: exp( -T(-T(x)) ), where T(x) = Sum_{n>=1} n^(n-1)*x^n/n!.
a(n) ~ n^(n-1)/(1+LambertW(1)). - Vaclav Kotesovec, Aug 26 2016
EXAMPLE
E.g.f.: A(x) = 1 + x + x^2/2! + 7*x^3/3! + 37*x^4/4! + 441*x^5/5! + 4771*x^6/6! + 79213*x^7/7! + 1320649*x^8/8! + 28318321*x^9/9! + 636978151*x^10/10! + 16863972621*x^11/11! + 475580960317*x^12/12! +...
such that A(x)^A(x) = LambertW(-x)/(-x), where
LambertW(-x)/(-x) = 1 + x + 3*x^2/2! + 16*x^3/3! + 125*x^4/4! + 1296*x^5/5! + 16807*x^6/6! + 262144*x^7/7! + 4782969*x^8/8! +...+ (n+1)^(n-1)*x^n/n! +...
The logarithm of the e.g.f. A(x) begins
log(A(x)) = x + 6*x^3/3! + 12*x^4/4! + 320*x^5/5! + 2190*x^6/6! + 51492*x^7/7! + 685496*x^8/8! + 17286768*x^9 +...+ A097174(n)*x^n/n! +...
which equals -T(-T(x)), where
T(x) = x + 2*x^2/2! + 9*x^3/3! + 64*x^4/4! + 625*x^5/5! + 7776*x^6/6! + 117649*x^7/7! + 2097152*x^8/8! +...+ n^(n-1)*x^n/n! +...
MATHEMATICA
CoefficientList[Series[E^LambertW[-LambertW[-x]], {x, 0, 20}], x] * Range[0, 20]! (* Vaclav Kotesovec, Aug 26 2016 *)
PROG
(PARI) {a(n) = my(A=1+x, W); W=serreverse(x*exp(-x +x^2*O(x^n)))/x;
for(i=0, n, A = W^(1/A) ); n!*polcoeff(A, n)}
for(n=0, 40, print1(a(n), ", "))
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Aug 24 2016
STATUS
approved