OFFSET
1,1
COMMENTS
The unitary version of A064771.
It appears that most of the terms are divisible by 3. Terms that are not divisible by 3 are 3770, 5530, 7210, 7630, ... - Michel Marcus, Dec 15 2017
The least odd term is 442365. - Amiram Eldar, Jun 10 2020
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
EXAMPLE
150 is in the sequence since its aliquot unitary divisors are 1, 2, 3, 6, 25, 50, 75 and there is only one subset whose sum is 150: {25, 50, 75}.
MATHEMATICA
ud[n_] := Block[{d = Divisors[n]}, Select[d, GCD[#, n/#] == 1 &]];
a = {}; n = 0; While[Length[a] < 100, n++; d = Most[ud[n]];
c = SeriesCoefficient[Series[Product[1 + x^d[[i]], {i, Length[d]}], {x, 0, n}], n]; If[c ==1, AppendTo[a, n]]]; a
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Nov 28 2017
STATUS
approved