OFFSET
1,1
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
EXAMPLE
165375 is in the sequence since both 165375 and 165376 are nonunitary abundant: nusigma(165375) = 179280 > 165375, and nusigma(165376) = 183600 > 165376 (nusigma is the sum of nonunitary divisors, A048146).
MATHEMATICA
f[p_, e_] := (p^(e + 1) - 1)/(p - 1); nuabQ[n_] := Times @@ (f @@@ FactorInteger[n]) - Times @@ (1 + Power @@@ FactorInteger[n]) > n; s = {}; q1 = False; Do[q2 = nuabQ[n]; If[q1 && q2, AppendTo[s, n - 1]]; q1 = q2, {n, 2, 10^7}]; s
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Sep 30 2019
STATUS
approved