OFFSET
1,1
COMMENTS
The odd terms of A348604.
The numbers of terms not exceeding 10^k, for k = 4, 5, ..., are 1, 51, 360, 4117, 39803, 418663, 4099004, ... Apparently this sequence has an asymptotic density 0.0004...
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
EXAMPLE
8505 is a term since A160135(8505) = 8862 > 8505.
MATHEMATICA
esigma[n_] := Times @@ (Sum[First[#]^d, {d, Divisors[Last[#]]}] &) /@ FactorInteger[n]; q[n_] := DivisorSigma[1, n] - esigma[n] > n; Select[Range[1, 65000, 2], q]
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Oct 25 2021
STATUS
approved