OFFSET
1,1
COMMENTS
Includes all the perfect numbers (A000396). The nonperfect terms have an abundant proper divisor which is not pseudoperfect, i.e., a proper divisor which is a weird number (A006037).
The first term with one weird divisor is a(3) = 350, having the weird divisor 70.
The first term with 2 weird divisors is a(202) = 658312, having the 2 weird divisors 9272 and 10792.
The first term with 3 weird divisors is a(353) = 1574930, having the 3 weird divisors 70, 10430 and 10570.
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..2500
EXAMPLE
350 is a term since it is pseudoperfect: 1 + 5 + 14 + 35 + 50 + 70 + 175 = 350. All of its proper divisors, {1, 2, 5, 7, 10, 14, 25, 35, 50, 70, 175} are not pseudoperfect, and it is not primitive abundant, since its divisor 70 is abundant.
MATHEMATICA
pspQ[n_] := Module[{d = Most @ Divisors[n], x}, Plus @@d >= n && SeriesCoefficient[Series[Product[1 + x^d[[i]], {i, Length[d]}], {x, 0, n}], n] > 0]; primPspQ[n_] := pspQ[n] && AllTrue[Most @ Divisors[n], !pspQ[#] &]; primAbQ[n_] := DivisorSigma[1, n] > 2*n && AllTrue[Most @ Divisors[n], DivisorSigma[1, #] < 2*# &]; Select[Range[1000], primPspQ[#] && !primAbQ[#] &]
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, May 30 2020
STATUS
approved