OFFSET
1,2
LINKS
Paul D. Hanna, Table of n, a(n) for n = 1..201
FORMULA
a(n) ~ sqrt(1-c) * 5^(5*n - 2) * n^(4*n - 5/2) / (sqrt(2*Pi) * (5-c)^(4*n-2) * c^n * exp(4*n)), where c = -LambertW(-5*exp(-5)) = 0.0348857682557236963012408686771... - Vaclav Kotesovec, May 18 2026
EXAMPLE
O.g.f.: A(x) = x + 98*x^2 + 214344*x^3 + 2422736288*x^4 + 86404408885000*x^5 + 7522353279784847520*x^6 + ...
ILLUSTRATION OF DEFINITION.
The table of coefficients of x^k/k! in exp( n^5*x - n^2*A(x) ) begins
n = 1: [1, 0, -196, -1286064, ...];
n = 2: [1, 28, 0, -5188160, ...];
n = 3: [1, 234, 52992, 0, ...];
n = 4: [1, 1008, 1012928, 994132224, 0, ...];
n = 5: [1, 3100, 9605100, 29713278400, 90217316417200, 0, ...]; ...
in which a diagonal, the coefficient of x^n in row n, is all zeros.
RELATED SERIES.
exp(A(x)) = 1 + x + 197*x^2/2! + 1286653*x^3/3! + 58150931593*x^4/4! + 10368822328678841*x^5/5! + 5416156760994815588941*x^6/6! + ...
where [x^n] exp(n^5*x) / exp(A(x))^(n^2) = 0 for n >= 1.
PROG
(PARI) {a(n, p=5, q=2) = my(A=[0], m); for(i=1, n+1, m=#A; A=concat(A, 0); A[m+1] = Vec( exp(m^p*x - m^q*Ser(A)) )[m+1]/(m^q) ); polcoef( Ser(A), n)}
for(n=1, 15, print1(a(n, 5, 2), ", "))
CROSSREFS
Below, sequences with the parameters (p,q) indicated have an o.g.f. F(x) satisfying [x^n] exp( n^p*x - n^q*F(x) ) = 0 for n >= 1.
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Mar 02 2026
STATUS
approved
