OFFSET
0,19
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..20000
FORMULA
a(n) = [x^n y^3] 1/Product_{j>=1} (1-y*x^A001358(j)).
MAPLE
h:= proc(n) option remember; `if`(n=0, 0,
`if`(numtheory[bigomega](n)=2, n, h(n-1)))
end:
b:= proc(n, i) option remember; series(`if`(n=0, 1, `if`(i<1, 0,
`if`(i>n, 0, x*b(n-i, h(min(n-i, i))))+b(n, h(i-1)))), x, 4)
end:
a:= n-> coeff(b(n, h(n)), x, 3):
seq(a(n), n=0..80);
MATHEMATICA
Table[Count[IntegerPartitions[n, {3}], _?(PrimeOmega[#]=={2, 2, 2}&)], {n, 0, 80}] (* Harvey P. Dale, Oct 12 2023 *)
CROSSREFS
KEYWORD
nonn,look
AUTHOR
Alois P. Heinz, May 19 2021
STATUS
approved