OFFSET
1,1
COMMENTS
The least term that is not divisible by 3 is a(89047132) = 134785275625.
The numbers of terms not exceeding 10^k, for k = 4, 5, ..., are 4, 60, 640, 6650, 66044, 660230, 6604594, 66073470, ... . Apparently, the asymptotic density of this sequence exists and equals 0.000660... .
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
EXAMPLE
4725 is a term since it is odd, and A162296(4725) = 9728 > 2*4725.
MATHEMATICA
q[n_] := Module[{f = FactorInteger[n], p, e}, p = f[[;; , 1]]; e = f[[;; , 2]]; Times @@ ((p^(e + 1) - 1)/(p - 1)) - Times @@ (p + 1) > 2*n]; Select[Range[3, 60000, 2], q]
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Oct 06 2022
STATUS
approved