OFFSET
1,2
COMMENTS
No more terms below 10^9. - Amiram Eldar, Aug 28 2019
MATHEMATICA
Select[ Range[ 10^7 ], GCD[ #, DivisorSigma[ 1, # ] - # - 1 ] >= #/3 && !PrimeQ[ # ] & ]
PROG
(PARI) f(n) = if(n>1, sigma(n)-n-1, 0); \\ A048050
isok(n) = !isprime(n) && gcd(n, f(n)) >= n/3; \\ Michel Marcus, Aug 28 2019 [Corrected by Sean A. Irvine, Nov 25 2023]
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Robert G. Wilson v, Jan 23 2002
EXTENSIONS
Name clarified by Sean A. Irvine, Nov 25 2023
STATUS
approved