OFFSET
1,1
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
MAPLE
N:= 1000: # for terms <= N
P:= select(isprime, [2, seq(i, i=3..N, 2)]):
nP:= nops(P):
S:= {}:
for i from 1 to nP do
for j from 1 to nP do
x:= P[i]*P[j]+P[i]-P[j];
if x > N then break fi;
if isprime(x) then S:= S union {x} fi
od od:
sort(convert(convert(P, set) minus S, list)); # Robert Israel, Jun 05 2025
MATHEMATICA
x = Take[ Select[ Union[ Flatten[ Table[ Prime[p]*Prime[q] + Prime[p] - Prime[q], {p, 116}, {q, 116}]]], PrimeQ[ # ] &], 200]; y = Table[ Prime[i], {i, PrimePi[ p[[ -1]] ]}]; Take[ Complement[y, x], 54] (* Robert G. Wilson v *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Zak Seidov, Feb 21 2004
STATUS
approved
