OFFSET
0,9
LINKS
FORMULA
a(n) = [x^n y^3] Product_{k>=1} 1/(1 - y*x^A246655(k)).
a(n) = Sum_{j=1..floor(n/3)} Sum_{i=j..floor((n-j)/2)} [omega(i) * omega(j) * omega(n-i-j) == 1], where omega(n) is the number of distinct prime factors of n and [==] is the Iverson bracket. - Wesley Ivan Hurt, Apr 25 2019
EXAMPLE
a(11) = 4 because we have [7, 2, 2], [5, 4, 2], [5, 3, 3] and [4, 4, 3].
MAPLE
f:= proc(n, k, pmax) option remember;
local t, p, j;
if n = 0 then return `if`(k=0, 1, 0) fi;
if k = 0 then return 0 fi;
if n > k*pmax then return 0 fi;
t:= 0:
for p in A246655 do
if p > pmax then return t fi;
t:= t + add(procname(n-j*p, k-j, min(p-1, n-j*p)), j=1..min(k, floor(n/p)))
od;
t
end proc:
seq(f(n, 3, n), n=0..80) # Robert Israel, Apr 25 2019
MATHEMATICA
Array[Count[IntegerPartitions[#, {3}], _?(AllTrue[#, PrimePowerQ] &)] &, 81, 0]
CROSSREFS
KEYWORD
nonn,look
AUTHOR
Ilya Gutkovskiy, Apr 24 2019
STATUS
approved