OFFSET
1,1
COMMENTS
Primes in A254636.
FORMULA
PROG
(Python)
import sympy
from sympy import isprime
TOP = 10000
a = [0]*TOP
for y in range(2, TOP//2):
for x in range(y, TOP//2):
k = x*y + x + y
if k>=TOP: break
a[k]+=1
print([n for n in range(TOP) if a[n]==0 and isprime(n)])
(PARI) v=[]; for(m=2, 500, for(k=m, 500, if(isprime(P=k*m+k+m), v=concat(v, P)))); v=vecsort(v, , 8); forprime(p=1, 2000, if(!vecsearch(v, p), print1(p, ", "))) \\ Derek Orr, Mar 21 2015
CROSSREFS
KEYWORD
nonn
AUTHOR
Alex Ratushnyak, Mar 14 2015
STATUS
approved