OFFSET
1,1
COMMENTS
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
EXAMPLE
60 is a term since it is nonsquarefree (it is divisible by 4 = 2^2) and it is equal to a sum of its aliquot unitary divisors: 1 + 3 + 4 + 5 + 12 + 15 + 20 = 60.
MATHEMATICA
pspQ[n_] := !SquareFreeQ[n] && Module[{d = Most @ Select[Divisors[n], CoprimeQ[#, n/#] &], x}, Plus @@ d >= n && SeriesCoefficient[Series[Product[1 + x^d[[i]], {i, Length[d]}], {x, 0, n}], n] > 0]; Select[Range[1000], pspQ]
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, May 25 2020
STATUS
approved