reviewed
approved
reviewed
approved
proposed
reviewed
editing
proposed
(PARI) isab(x) = sigma(x) > 2*x; \\ A005101
isok(k) = isab(k) && isab(k+1) && isab(k+2); \\ Michel Marcus, Nov 19 2022
approved
editing
Amiram Eldar, <a href="/A096536/b096536_1.txt">Table of n, a(n) for n = 1..22583</a> (terms below 10^12; terms 1..1000 from Donovan Johnson)
reviewed
approved
proposed
reviewed
editing
proposed
From Amiram Eldar, Oct 02 2022: (Start)
The least term that is == 1 (mod 3) is a(1292) = 55959128224, and the least term that is divisible by 3 is a(1590) = 68972878974.
The numbers of terms not exceeding 10^k, for k = 9, 10, ..., are 9, 226, 2298, 22583, ... . Apparently, the asymptotic density of this sequence exists and equals 2.2...*10^(-8). (End)
For 171078830 = 2*5*13*23*29*1973, sigma(n)/n = 2.09355, for 171078831 = 3^3*7*11*19*61*71, sigma(n)/n = 2.00396 and for 171078832 = 2^4*31*344917, sigma(n)/n = 2.00000579.