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