OFFSET
1,1
COMMENTS
Conjecture: all terms are of the form 2*m^2 (twice a perfect square; A001105).
EXAMPLE
For k = 18: sigma(18) = 39 and gcd(18, 39) = 3, so 18 is a term; moreover 18 = 2*3^2.
For k = 5652: sigma(5652) = 14378 and gcd(5652, 14378) = 2, so 5652 is not a term.
PROG
(Python)
from sympy import divisor_sigma, gcd, isprime
def ok(k):
return isprime(k-1) and isprime(k+1) and gcd(k, divisor_sigma(k)) == 3
print([k for k in range(1, 20000000) if ok(k)])
CROSSREFS
KEYWORD
nonn
AUTHOR
Aied Sulaiman, Apr 24 2026
STATUS
approved
