OFFSET
16,7
LINKS
FORMULA
a(n) = Sum_{k=1..floor(n/4)} Sum_{j=k..floor((n-k)/3)} Sum_{i=j..floor((n-j-k)/2)} [Omega(k) = Omega(j) = Omega(i) = Omega(n-i-j-k) = 2], where Omega is the number of prime factors with multiplicity (A001222) and [ ] is the (generalized) Iverson bracket.
a(n) = [x^n y^4] 1/Product_{j>=1} (1-y*x^A001358(j)). - Alois P. Heinz, May 21 2021
MATHEMATICA
Table[Sum[Sum[Sum[KroneckerDelta[PrimeOmega[k], PrimeOmega[j], PrimeOmega[i], PrimeOmega[n - i - j - k], 2], {i, j, Floor[(n - j - k)/2]}], {j, k, Floor[(n - k)/3]}], {k, Floor[n/4]}], {n, 16, 100}]
Table[Count[IntegerPartitions[n, {4}], _?(PrimeOmega[#]=={2, 2, 2, 2}&)], {n, 16, 95}] (* Harvey P. Dale, May 14 2022 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Wesley Ivan Hurt, Jan 19 2021
STATUS
approved
