OFFSET
0,6
COMMENTS
a(n) = A000586(n) for n <= 42.
LINKS
Charles R Greathouse IV, Table of n, a(n) for n = 0..10000
FORMULA
G.f.: Product_{k>=1} (1 + x^A109611(k)). - Andrew Howroyd, Dec 28 2017
MATHEMATICA
terms = 81;
gf = Times @@ (1 + x^SequencePosition[ PrimeOmega[ Range[terms]], {1, _, 1|2}][[All, 1]]) + O[x]^terms;
CoefficientList[gf, x] (* Jean-François Alcover, Jul 02 2018 *)
PROG
(PARI) P=1+O(x^1001); forprime(p=2, 1e3, if(bigomega(p+2)<3, P*=1+x^p)); Vec(P) \\ Charles R Greathouse IV, May 13 2013
CROSSREFS
KEYWORD
nonn
AUTHOR
Reinhard Zumkeller, Aug 26 2005
STATUS
approved