OFFSET
1,2
COMMENTS
It is remarkable that this sequence should consist entirely of integers.
LINKS
Paul D. Hanna, Table of n, a(n) for n = 1..301
FORMULA
a(n) ~ sqrt(1-c) * 3^(3*n-1) * n^(2*n - 3/2) / (sqrt(2*Pi) * c^n * (3-c)^(2*n-1) * exp(2*n)), where c = -LambertW(-3*exp(-3)) = 0.1785606278779211065968086697... = -A226750. - Vaclav Kotesovec, Aug 06 2018
B(x) = x + (1/(m+1)) * A((m+1)^3 * x^(m+1)) satisfies [x^n] exp(n*x*(1+n^2*x^m) - n*B(x) ) = 0 for n >= 1. - Seiichi Manyama, Apr 09 2026
EXAMPLE
O.g.f.: A(x) = x + 9*x^2 + 552*x^3 + 85842*x^4 + 24653700*x^5 + 11219022936*x^6 + 7393496092416*x^7 + 6649411839351120*x^8 + ...
ILLUSTRATION OF DEFINITION.
The table of coefficients of x^k/k! in exp( n^3*x - n*A(x) ) begins:
n=1: [1, 0, -18, -3312, -2059236, ...];
n=2: [1, 6, 0, -7056, -4281984, ...];
n=3: [1, 24, 522, 0, -6980580, ...];
n=4: [1, 60, 3528, 189792, 0, ...];
n=5: [1, 120, 14310, 1679040, 181358460, 0, ...]; ...
in which the coefficient of x^n in row n forms a diagonal of zeros.
RELATED SERIES.
exp(A(x)) = 1 + x + 19*x^2/2! + 3367*x^3/3! + 2074537*x^4/4! + 2969379361*x^5/5! + 8096147776171*x^6/6! + 37321188279552199*x^7/7! + ...
where [x^n] exp(n^3*x) / exp(A(x))^n = 0 for n >= 1.
PROG
(PARI) {a(n, k=0) = my(A=[0], m); for(i=1, n+1, m=#A; A=concat(A, 0); A[m+1] = Vec( exp(m^2*(m+k)*x - m*Ser(A)) )[m+1]/m ); polcoef( Ser(A), n)}
for(n=1, 15, print1(a(n, 0), ", "))
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Jul 29 2018
EXTENSIONS
Entry revised by Paul D. Hanna, Mar 01 2026
STATUS
approved
