OFFSET
1,1
COMMENTS
The first odd term is a(3577) = 99225.
The number of terms not exceeding 10^k, for k = 1, 2, ..., are 0, 3, 31, 360, 3605, 36160, 360840, 3618980, 36144059, ... Apparently this sequence has an asymptotic density 0.0361...
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
EXAMPLE
36 is a term since A348271(36) = 41 > 36.
MATHEMATICA
f[p_, e_] := Module[{b = IntegerDigits[e, 2], m}, m = Length[b]; Product[If[b[[j]] > 0, 1 + p^(2^(m - j)), 1], {j, 1, m}]]; isigma[1] = 1; isigma[n_] := Times @@ f @@@ FactorInteger[n]; s[n_] := DivisorSigma[1, n] - isigma[n]; Select[Range[1500], s[#] > # &]
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Oct 09 2021
STATUS
approved