OFFSET
0,5
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..1000 (first 201 terms from Reinhard Zumkeller)
FORMULA
G.f.: Product(k>=1, (1+x^(3*k))*(1+x^(3*k-2)) ). - Emeric Deutsch, Aug 26 2007
a(n) ~ exp(Pi*sqrt(2*n)/3) / (2^(19/12) * sqrt(3) * n^(3/4)). - Vaclav Kotesovec, Aug 24 2015
EXAMPLE
a(7)=3 because we have 7, 61 and 43.
MAPLE
g:=product((1+x^(3*k))*(1+x^(3*k-2)), k=1..30): gser:=series(g, x=0, 100): seq(coeff(gser, x, n), n=0..65); # Emeric Deutsch, Aug 26 2007
MATHEMATICA
nmax = 100; CoefficientList[Series[Product[((1+x^(3*k))*(1+x^(3*k-2))), {k, 1, nmax}], {x, 0, nmax}], x] (* Vaclav Kotesovec, Aug 24 2015 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Reinhard Zumkeller, Aug 22 2007
EXTENSIONS
Prepended a(0) = 1, Joerg Arndt, Feb 22 2015
STATUS
approved