OFFSET
0,3
LINKS
Vaclav Kotesovec, Table of n, a(n) for n = 0..1000
FORMULA
a(n) ~ -c * 2^(2*n - 1) / (sqrt(Pi) * n^(3/2)), where c = Product_{k>=1} (1 + sqrt(1 - 4*(1/4)^k))/2 = 0.4567034206737725013365271429022657551331606541289778092649... - Vaclav Kotesovec, May 06 2021
MATHEMATICA
nmax = 30; CoefficientList[Series[Product[(1+Sqrt[1-4*x^k])/2, {k, 1, nmax}], {x, 0, nmax}], x] (* Vaclav Kotesovec, May 06 2021 *)
PROG
(PARI) N=66; x='x+O('x^N); Vec(prod(k=1, N, (1+sqrt(1-4*x^k))/2))
(PARI) N=66; x='x+O('x^N); Vec(prod(i=1, N, 1-sum(j=1, N\i, binomial(2*j-2, j-1)*x^(i*j)/j)))
CROSSREFS
KEYWORD
sign
AUTHOR
Seiichi Manyama, Sep 22 2019
STATUS
approved