OFFSET
1,1
COMMENTS
All terms are multiples of 6.
PROG
(Python)
from sympy import isprime, divisors
print([x for x in range(2, 2000) if any(all(isprime(v) for v in (x//y-1, x//y+1, x*y-1, x*y+1)) for y in divisors(x)[1:-1])])
(PARI) isok(k) = fordiv(k, d, if ((d>1) && (d<k) && isprime(k*d+1) && isprime(k*d-1) && isprime(k/d+1) && isprime(k/d-1), return(1))); \\ Michel Marcus, Mar 16 2026
CROSSREFS
KEYWORD
nonn
AUTHOR
Alex Ratushnyak, Mar 15 2026
STATUS
approved
