OFFSET
0,3
COMMENTS
Conjecture: given o.g.f. A(x), then sqrt( (A(x) - 1)/x ) is an integer series.
LINKS
Vaclav Kotesovec, Table of n, a(n) for n = 0..240
FORMULA
G.f.: 1/Q(0), where Q(k) = 1 -(2*k+1)^2*x/(1 -(2*k+1)^2*x/Q(k+1) ); (continued fraction). - Sergei N. Gladkovskii, Sep 17 2013
G.f.: Q(0), where Q(k) = 1 - x*(2*k+1)^2/( x*(2*k+1)^2 - 1/(1 - x*(2*k+1)^2/( x*(2*k+1)^2 - 1/Q(k+1)))); (continued fraction). - Sergei N. Gladkovskii, Oct 09 2013
a(n) ~ 2^(4*n+2) * n^(2*n-1/2) / (Pi^(2*n+1/2) * exp(2*n)). - Vaclav Kotesovec, Aug 25 2017
EXAMPLE
O.g.f.: A(x) = 1 + x + 2*x^2 + 13*x^3 + 206*x^4 + 5794*x^5 + 252068*x^6 +...
Let A(x) = 1 + x*B(x)^2, then B(x) appears to be an integer series:
B(x) = 1 + x + 6*x^2 + 97*x^3 + 2782*x^4 + 122670*x^5 + 7687932*x^6 + 649446621*x^7 + 71136143478*x^8 + 9806113041658*x^9 +...
the coefficients of B(x) are integer for at least the initial 500 terms.
MATHEMATICA
nmax = 20; CoefficientList[Series[1/Fold[(1 - #2/#1) &, 1, Reverse[(2*Range[nmax + 1] - 2*Floor[Range[nmax + 1]/2] - 1)^2*x]], {x, 0, nmax}], x] (* Vaclav Kotesovec, Aug 25 2017 *)
PROG
(PARI) /* Continued fraction: */
{a(n)=local(A=1+x, CF); CF=1+x; for(k=0, n, CF=1/(1-(2*((n-k)\2)+1)^2*x*CF+x*O(x^n))); A=CF; polcoeff(A, n)}
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Jul 30 2011
STATUS
approved