OFFSET
1,3
COMMENTS
a(n) (mod 3) yields period 6 sequence: [0,0,0,0,0,2] starting at n=3.
a(6*n+2) == 2 (mod 3) for n>=1.
FORMULA
E.g.f. satisfies:
(1) A(x - x^2*cosh(x)^2/2) = x.
(2) A(x) = x + Sum_{n>=1} d^(n-1)/dx^(n-1) cosh(x)^(2*n)*x^(2*n) / (2^n*n!).
(2) A(x) = x*exp( Sum_{n>=1} d^(n-1)/dx^(n-1) cosh(x)^(2*n)*x^(2*n-1) / (2^n*n!) ).
a(n) ~ n^(n-1) * sqrt(2/(1 + (1+2*s^2)*cosh(2*s) + 4*s*sinh(2*s))) / (exp(n) * r^(n-1/2)), where s = 0.568824148293580379787367453... is the root of the equation s*cosh(s)*(cosh(s) + s*sinh(s)) = 1, and r = s - s^2*(cosh(s))^2/2 = 0.3488028982491643456675... - Vaclav Kotesovec, Jan 12 2014
EXAMPLE
E.g.f: A(x) = x + x^2/2! + 3*x^3/3! + 27*x^4/4! + 285*x^5/5! + 3585*x^6/6! +...
Related expansions:
A(x)^2 = 2*x^2/2! + 6*x^3/3! + 30*x^4/4! + 330*x^5/5! + 4410*x^6/6! + 67830*x^7/7! +...
A(x) = x + cosh(x)^2*x^2/2 + d/dx cosh(x)^4*x^4/(4*2!) + d^2/dx^2 cosh(x)^6*x^6/(8*3!) + d^3/dx^3 cosh(x)^8*x^8/(16*4!) +...
log(A(x)/x) = 1 + cosh(x)^2*x/2 + d/dx cosh(x)^4*x^3/(4*2!) + d^2/dx^2 cosh(x)^6*x^5/(8*3!) + d^3/dx^3 cosh(x)^8*x^7/(16*4!) +...
MATHEMATICA
Rest[CoefficientList[InverseSeries[Series[x - (x^2*Cosh[x]^2)/2, {x, 0, 20}], x], x] * Range[0, 20]!] (* Vaclav Kotesovec, Jan 12 2014 *)
PROG
(PARI) {a(n)=n!*polcoeff(serreverse(x-x^2/2*cosh(x+x*O(x^n))^2), n)}
for(n=1, 30, print1(a(n), ", "))
(PARI) {Dx(n, F)=local(D=F); for(i=1, n, D=deriv(D)); D}
{a(n)=local(A=x); A=x+sum(m=1, n, Dx(m-1, cosh(x+x*O(x^n))^(2*m)*x^(2*m)/2^m/m!)); n!*polcoeff(A, n)}
(PARI) {Dx(n, F)=local(D=F); for(i=1, n, D=deriv(D)); D}
{a(n)=local(A=x+x^2+x*O(x^n)); A=x*exp(sum(m=1, n, Dx(m-1, cosh(x+x*O(x^n))^(2*m)*x^(2*m-1)/2^m/m!)+x*O(x^n))); n!*polcoeff(A, n)}
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Jul 13 2012
STATUS
approved