OFFSET
1,3
COMMENTS
Related identity: Sum_{n>=0} (q^n + p)^n * r^n/n! = Sum_{n>=0} exp(p*q^n*r) * q^(n^2) * r^n/n!; here, q = A(x), p = -x, r = 1.
Conjecture: for n > 0, a(6*n + k) == [2,0,2,1,2,2] (mod 3) at k = [0,1,2,3,4,5], respectively.
LINKS
Paul D. Hanna, Table of n, a(n) for n = 1..300
FORMULA
E.g.f. A(x) = Sum_{n>=0} a(n)*x^n/n! satisfies the following formulas.
(1) Sum_{n>=0} (A(x)^n - x)^n / n! = 1.
(2) Sum_{n>=0} exp(-x*A(x)^n) * A(x)^(n^2) / n! = 1.
a(n) ~ -c * n^(n-1) / (exp(n) * r^n), where c = sqrt( (-r*Sum_{k>=0} ((k*(-r + s^k)^(-1 + k))/k!)) / Sum_{k>=0} ((-1 + k)*k^2 * s^(-2 + k)*(-r + s^k)^(-2 + k)*(-r + (1 + k)*s^k))/k! ) = 0.4187940561612508319941857365856199965995537726..., and where r = -0.3491010753747466229482161022113556770139942642631... and s = -0.6426571454882319173136663034041471668334385049965... are roots of the equations Sum_{k>=0} (-r + s^k)^k/k! = 1 and Sum_{k>=0} (k^2*s^(-1 + k)*(-r + s^k)^(-1 + k))/k! = 0. - Vaclav Kotesovec, Apr 11 2024
The values r and A(r) = s given above also satisfy Sum_{n>=0} n*(n - r*A(r)^n) * exp(-r*A(r)^n) * A(r)^(n^2) / n! = 0. - Paul D. Hanna, Apr 12 2024
EXAMPLE
E.g.f.: A(x) = x - x^2/2! + 7*x^3/3! - 37*x^4/4! + 371*x^5/5! - 4741*x^6/6! + 72885*x^7/7! - 1380541*x^8/8! + 29678815*x^9/9! - 726074821*x^10/10! + ...
where e.g.f. A(x) satisfies the following sums.
(1) 1 = 1 + (A(x) - x) + (A(x)^2 - x)^2/2! + (A(x)^3 - x)^3/3! + (A(x)^4 - x)^4/4! + (A(x)^5 - x)^5/5! + ...
(2) 1 = exp(-x) + exp(-x*A(x))*A(x) + exp(-x*A(x)^2)*A(x)^4/2! + exp(-x*A(x)^3)*A(x)^9/3! + exp(-x*A(x)^4)*A(x)^16/4! + ...
SPECIFIC VALUES.
A(1/3) = 0.309336999832107073180903710282149168034207161078640395207...
A(1/4) = 0.232922937634173409470673241764259081533730452334005659588...
A(1/5) = 0.187560199855301209894398645611115284037479048219241021351...
A(-1/3) = -0.5146620783815103062311605400508155869729182062358910349...
A(-1/4) = -0.3123628005245983090140211998639545568283470783996606926...
A(-1/5) = -0.2335665203884038676850050992335539648367581317265287642...
PROG
(PARI) /* Sum_{n>=0} exp(-x*A(x)^n) * A(x)^(n^2) / n! = 1 */
{a(n) = my(A=[0, 1]); for(i=1, n, A = concat(A, 0);
A[#A] = -polcoeff( sum(m=0, sqrtint(#A+1), exp(-x*Ser(A)^m +x*O(x^#A)) * Ser(A)^(m^2)/m! ), #A-1); ); n!*A[n+1]}
for(n=1, 30, print1(a(n), ", "))
(PARI) /* Sum_{n>=0} (A(x)^n - x)^n / n! = 1 */
{a(n) = my(A=[0, 1]); for(i=1, n, A = concat(A, 0);
A[#A] = -polcoeff( sum(m=0, #A+1, (Ser(A)^m - x)^m/m! ), #A-1) ); n!*A[n+1]}
for(n=1, 30, print1(a(n), ", "))
CROSSREFS
KEYWORD
sign
AUTHOR
Paul D. Hanna, Apr 10 2024
STATUS
approved