OFFSET
0,11
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)], where omega is the number of distinct prime factors (A001221) and [ ] is the (generalized) Iverson bracket.
MATHEMATICA
Table[Sum[Sum[Sum[KroneckerDelta[PrimeNu[k], PrimeNu[j], PrimeNu[i], PrimeNu[n - i - j - k]], {i, j, Floor[(n - j - k)/2]}], {j, k, Floor[(n - k)/3]}], {k, Floor[n/4]}], {n, 0, 80}]
CROSSREFS
KEYWORD
nonn
AUTHOR
Wesley Ivan Hurt, Jan 19 2021
STATUS
approved