OFFSET
1,1
COMMENTS
All the terms are nondeficient numbers (A023196).
All the perfect numbers (A000396) are terms.
This sequence is infinite: if k = 2^(p-1)*(2^p-1) is an even perfect number and q > 2^p-1 is a prime, then k*q is a term.
Since the total sum of divisors of any term is even, none of the terms are squares or twice squares.
Are there odd terms in this sequence? There are none below 10^10.
The numbers of terms not exceeding 10^k, for k = 1, 2, ..., are 1, 6, 63, 605, 6164, 61291, 614045, 6139193, 61382607, 613861703, ... . Apparently, the asymptotic density of this sequence exists and equals 0.06138... .
LINKS
EXAMPLE
6 is a term since the sum of its deficient divisors, 1 + 2 + 3 is equal to 6, its only nondeficient divisor.
30 is a term since the sum of its deficient divisors, 1 + 2 + 3 + 5 + 10 + 15 = 36 is equal to the sum of its nondeficient divisors, 6 + 30 = 36.
MATHEMATICA
q[n_] := DivisorSum[n, If[DivisorSigma[-1, #] < 2, #, -#] &] == 0; Select[Range[1000], q]
PROG
(PARI) is(n) = sumdiv(n, d, if(sigma(d, -1) < 2, d, -d)) == 0;
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Sep 29 2022
STATUS
approved