A318001 E.g.f. A(x) satisfies: cosh(A(x) - A(-x)) - sinh(A(x) + A(-x)) = 1.

1, 2, 6, 56, 600, 8432, 144816, 2892416, 66721920, 1732489472, 50144683776, 1604936139776, 56236356234240, 2137961925773312, 87642967518836736, 3863105286629851136, 182345733925971394560, 9130908475775186173952, 481864839159167717277696, 27108466364634568866922496, 1642481780780610712999034880

Offset

1,2

Comments

First negative term is a(27).

Links

Paul D. Hanna, Table of n, a(n) for n = 1..100

Formula

E.g.f. A(x) satisfies:

(1) A(-A(-x)) = x.

(2a) 1 = Sum_{n>=0} (-1)^n * ( A(x) - (-1)^n*A(-x) )^n/n!.

(2b) 1 = Sum_{n>=0} ( x + (-1)^n*A(A(x)) )^n/n!.

(3a) 1 = cosh(A(x) - A(-x)) - sinh(A(x) + A(-x)).

(3b) 1 = cosh(A(-x))*exp(-A(x)) - sinh(A(-x))*exp(A(x)).

(3c) 1 = cosh(x)*exp(-A(A(x))) + sinh(x)*exp(A(A(x))).

(4a) A(x) = log( 2*cosh(A(-x)) / (1 + sqrt(1 + 2*sinh(2*A(-x)))) ).

(4b) A(x) = log( (sqrt(1 + 2*sinh(2*A(-x))) - 1) / (2*sinh(A(-x))) ).

(5) A(A(x)) = log( 2*cosh(x) / (1 + sqrt(1 - 2*sinh(2*x))) ), which is the e.g.f. of A318000.

Example

E.g.f.: A(x) = x + 2*x^2/2! + 6*x^3/3! + 56*x^4/4! + 600*x^5/5! + 8432*x^6/6! + 144816*x^7/7! + 2892416*x^8/8! + 66721920*x^9/9! + 1732489472*x^10/10! + 50144683776*x^11/11! + 1604936139776*x^12/12! + 56236356234240*x^13/13! + 2137961925773312*x^14/14! + 87642967518836736*x^15/15! + ...
such that cosh(A(x) - A(-x)) - sinh(A(x) + A(-x)) = 1.
RELATED SERIES.
(1) exp(A(x)) = 1 + x + 3*x^2/2! + 13*x^3/3! + 105*x^4/4! + 1141*x^5/5! + 16083*x^6/6! + 276193*x^7/7! + 5561265*x^8/8! + 128834761*x^9/9! + 3365571363*x^10/10! + ...
which equals (sqrt(1 + 2*sinh(2*A(-x))) - 1) / (2*sinh(A(-x))).
(2) A(A(x)) = x + 4*x^2/2! + 24*x^3/3! + 256*x^4/4! + 3840*x^5/5! + 73024*x^6/6! + 1688064*x^7/7! + 45991936*x^8/8! + 1443102720*x^9/9! + ... + A318000(n)*x^n/n! + ...
which equals log( 2*cosh(x) / (1 + sqrt(1 - 2*sinh(2*x))) ).

Prog

(PARI) {a(n) = my(A=x+x^2 +x*O(x^n), S=x); for(i=1, n, S = (A - subst(A, x, -x))/2;
A = S + log(cosh(2*S) - 1 + sqrt(1 + (cosh(2*S) - 1)^2))/2;
A = (A - subst(serreverse(A), x, -x))/2 ); n!*polcoeff(A, n)}
for(n=1, 25, print1(a(n), ", "))

Crossrefs

Cf. A318000 (A(A(x))), A318006 (variant).

Sequence in context: A074023 A354315 A000146 * A211933 A167010 A014070

Adjacent sequences: A317998 A317999 A318000 * A318002 A318003 A318004

Keyword

sign

Author

Paul D. Hanna, Aug 20 2018

Status

approved