OFFSET
1,1
COMMENTS
All the terms are nonsquarefree (since squarefree numbers do not have nonunitary divisors).
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..400
EXAMPLE
840 is a term since its aliquot unitary divisors are {1, 3, 5, 7, 8, 15, 21, 24, 35, 40, 56, 105, 120, 168, 280} and 1 + 5 + 7 + 8 + 15 + 35 + 40 + 56 + 105 + 120 + 168 + 280 = 840, and its nonunitary divisors are {2, 4, 6, 10, 12, 14, 20, 28, 30, 42, 60, 70, 84, 140, 210, 420} and 70 + 140 + 210 + 420 = 840.
MATHEMATICA
pspQ[n_] := Module[{d = Divisors[n], ud, nd, x}, ud = Select[d, CoprimeQ[#, n/#] &]; nd = Complement[d, ud]; ud = Most[ud]; Plus @@ ud >= n && Plus @@ nd >= n && SeriesCoefficient[Series[Product[1 + x^ud[[i]], {i, Length[ud]}], {x, 0, n}], n] > 0 && SeriesCoefficient[Series[Product[1 + x^nd[[i]], {i, Length[nd]}], {x, 0, n}], n] > 0]; Select[Range[10^4], pspQ]
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, May 25 2020
STATUS
approved