OFFSET
1,1
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
EXAMPLE
30 is a term since its proper unitary divisors, 1 < d < 30, are {2, 3, 5, 6, 10, 15}, and 5 + 10 + 15 = 30.
MATHEMATICA
q[n_] := Module[{d = Most @ Select[Divisors[n], CoprimeQ[#, n/#] &], x}, Plus @@ d >= n && SeriesCoefficient[Series[Product[1 + x^d[[i]], {i, 2, Length[d]}], {x, 0, n}], n] > 0]; Select[Range[1000], q]
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Mar 10 2021
STATUS
approved