OFFSET
1,1
COMMENTS
The nonunitary version of A005835.
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
EXAMPLE
36 is in the sequence since its nonunitary divisors are 2, 3, 6, 12, 18 and 36 = 6 + 12 + 18.
MATHEMATICA
nudiv[n_] := Module[{d = Divisors[n]}, Select[d, GCD[#, n/#] > 1 &]]; s = {}; Do[d = nudiv[n]; If[Total[d] < n, Continue[]]; c = SeriesCoefficient[Series[Product[1 + x^d[[i]], {i, Length[d]}], {x, 0, n}], n]; If[c > 0, AppendTo[s, n]], {n, 1, 700}]; s
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Sep 30 2019
STATUS
approved