OFFSET
1,3
COMMENTS
a(n) modulo 8 has period 4: [1, 1, 5, 5].
a(n) = 5 (mod 10) for n > 2.
First negative term is a(21).
LINKS
Paul D. Hanna, Table of n, a(n) for n = 1..100
EXAMPLE
E.g.f.: A(x) = x + x^3/3! + 5*x^5/5! + 85*x^7/7! + 2825*x^9/9! + 151625*x^11/11! + 12098125*x^13/13! + 1339476125*x^15/15! + 196410020625*x^17/17! + 37062144900625*x^19/19! + 8772471210303125*x^21/21! + 2519410212081953125*x^23/23! + 854580849916226265625*x^25/25! +...
such that
A(A(x)) = x + 2*x^3 + 20*x^5 + 520*x^7 + 26000*x^9 + 2132000*x^11 + 260104000*x^13 + 44217680000*x^15 + ... + [Product_{k=1..n-1} ((2*k-1)^2 + 1)] * x^(2*n-1)/(2*n-1)! + ...
which equals sinh( arcsin(x) ).
RELATED SERIES.
Given e.g.f. F(x) of A280790 such that sin( F( sinh( F(x) ) ) ) = x,
then A(x) = sinh( F(x) ) where
F(x) = x + 4*x^5/5! + 2320*x^9/9! + 9857600*x^13/13! + 159122080000*x^17/17! + 7098806416000000*x^21/21! + 686863244097538560000*x^25/25! + ... + A280790(n)*x^(4*n-3)/(4*n-3)! + ...
also sin( F( A(x) ) ) = x.
PROG
(PARI) {a(n) = my(A=[1]); for(i=1, n, A = concat(A, [0, 0]); A[#A] = polcoeff(sinh(asin(x +x*O(x^#A))) - subst(x*Ser(A), x, x*Ser(A)), #A)/2 ); (2*n-1)!*A[2*n-1]}
for(n=1, 20, print1(a(n), ", "))
CROSSREFS
KEYWORD
sign
AUTHOR
Paul D. Hanna, Sep 05 2018
STATUS
approved