|
|
A309867
|
|
Expansion of Product_{k>0} (1+sqrt(1-4*x^k))/2.
|
|
2
|
|
|
1, -1, -2, -2, -5, -9, -36, -104, -365, -1219, -4213, -14617, -51570, -183084, -656536, -2370066, -8613590, -31478538, -115632718, -426676244, -1580878746, -5878933054, -21936060630, -82100980070, -308146839623, -1159545407027, -4373730398473, -16533813947503
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,3
|
|
LINKS
|
|
|
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
|
|
|
STATUS
|
approved
|
|
|
|