OFFSET
0,2
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..40
FORMULA
a(n) = [x^n] [ sqrt((1+x)/(1-x)) ]^(4^n).
More generally, the following coefficient of x^n in the series:
[x^n] Sum_{n>=0} arctanh(q^n*x)^n/n! = [x^n] [ sqrt((1+x)/(1-x)) ]^(q^n) is an integer for any even integer q.
EXAMPLE
G.f. A(x) = 1 + 4*x + 128*x^2 + 43712*x^3 + 178978816*x^4 + ...
MATHEMATICA
Table[SeriesCoefficient[(Sqrt[(1 + x)/(1 - x)])^(4^n), {x, 0, n}], {n, 0, 25}] (* G. C. Greubel, Apr 15 2017 *)
PROG
(PARI) {a(n)=polcoeff(sum(k=0, n, atanh(4^k*x +x*O(x^n))^k/k!), n)}
(PARI) {a(n)=polcoeff(((1+x)/(1-x +x*O(x^n)))^(4^n/2), n)}
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Jul 02 2008
STATUS
approved