|
|
A357607
|
|
Odd numbers k such that A162296(k) > 2*k.
|
|
2
|
|
|
4725, 6615, 7875, 8505, 11025, 14175, 15435, 17325, 19845, 20475, 22275, 23625, 24255, 25515, 26775, 28665, 29925, 31185, 33075, 36225, 36855, 37125, 37485, 38115, 39375, 40425, 41895, 42525, 46305, 47775, 48195, 50715, 51975, 53235, 53865, 55125, 57915, 59535
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
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
|
|
|
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
|
|
|
STATUS
|
approved
|
|
|
|