OFFSET
0,3
LINKS
Seiichi Manyama, Table of n, a(n) for n = 0..10000 (terms 0..1000 from Alois P. Heinz)
FORMULA
G.f.: 1/Product_{k>=1} (1-x^(k*(3*k-1)/2))*(1-x^(k*(3*k+1)/2)).
MATHEMATICA
nmax = 100; CoefficientList[Series[1/Product[(1-x^(k*(3*k-1)/2)) * (1-x^(k*(3*k+1)/2)), {k, 1, nmax}], {x, 0, nmax}], x] (* Vaclav Kotesovec, Dec 10 2017 *)
PROG
(PARI)
b(n) = (3*n^2 + 2*n + (n%2) * (2*n + 1)) / 8; \\ A001318
N=66; x='x+O('x^N);
Vec(1/prod(k=1, N, (1-x^b(k))) )
\\ Joerg Arndt, Oct 13 2014
CROSSREFS
KEYWORD
nonn
AUTHOR
Jon Perry, Jul 06 2004
STATUS
approved