A327716 Expansion of Product_{k>=1} B(x^k), where B(x) is the g.f. of A003823.
1, 1, 1, 1, 2, 3, 3, 3, 4, 6, 7, 9, 10, 12, 14, 17, 21, 23, 26, 32, 40, 45, 51, 58, 69, 80, 89, 102, 116, 135, 154, 177, 198, 224, 253, 288, 326, 361, 408, 459, 521, 583, 650, 723, 812, 909, 1009, 1122, 1244, 1393, 1547, 1716, 1898, 2101, 2326, 2575, 2845, 3132, 3456, 3809
Offset: 0
Keywords
Links
- Vaclav Kotesovec, Table of n, a(n) for n = 0..10000
- Eric Weisstein's World of Mathematics, Rogers-Ramanujan Continued Fraction.
Crossrefs
Programs
-
Mathematica
nmax = 60; CoefficientList[Series[Product[QPochhammer[x^(5*j - 3)] * QPochhammer[x^(5*j - 2)]/(QPochhammer[x^(5*j - 4)] * QPochhammer[x^(5*j - 1)]), {j, 1, nmax}], {x, 0, nmax}], x] (* Vaclav Kotesovec, Sep 23 2019 *) -
PARI
N=66; x='x+O('x^N); Vec(1/prod(k=1, N, (1-x^k)^sumdiv(k, d, kronecker(5, d))))
Formula
G.f.: Product_{i>=1} Product_{j>=1} (1-x^(i*(5*j-2))) * (1-x^(i*(5*j-3))) / ((1-x^(i*(5*j-1))) * (1-x^(i*(5*j-4)))).
G.f.: Product_{k>=1} (1-x^k)^(-A035187(k)).
a(n) ~ c * exp(Pi*sqrt(r*n)) / n^(3/4), where r = 4*log((1+sqrt(5))/2) / (3*sqrt(5)) = 0.2869392939760026925..., c = 0.203427046022096... - Vaclav Kotesovec, Sep 24 2019, updated May 09 2020
Comments