OFFSET
0,7
LINKS
FORMULA
a(n) = [x^n y^2] Product_{k>=1} 1/(1 - y*x^A246655(k)).
EXAMPLE
a(10) = 3 because we have [8, 2], [7, 3] and [5, 5].
MAPLE
# note that this requires A246655 to be pre-computed
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:
map(f, [$0..100]); # Robert Israel, Apr 29 2019
MATHEMATICA
Array[Count[IntegerPartitions[#, {2}], _?(AllTrue[#, PrimePowerQ] &)] &, 101, 0]
CROSSREFS
KEYWORD
nonn,look
AUTHOR
Ilya Gutkovskiy, Apr 24 2019
STATUS
approved