OFFSET
1,1
LINKS
Eric Weisstein's World of Mathematics, Quiteprime
PROG
(Python)
from math import isqrt
from sympy import primerange
def ok(n):
if n < 2: return False
for p in primerange(2, isqrt(n)+1):
isqrtp = isqrt(p)
if abs(2*(n%p)-p) > p + 1 - isqrt(p) - (isqrtp**2 < p): return False
return True
print(list(filter(ok, range(2, 350)))) # Michael S. Branicky, Aug 11 2021
CROSSREFS
KEYWORD
nonn,changed
AUTHOR
Eric W. Weisstein, Jim Ferry (jferry(AT)uiuc.edu)
EXTENSIONS
Title improved by Sean A. Irvine, Aug 11 2021
STATUS
approved