OFFSET
1,1
COMMENTS
LINKS
Antti Karttunen, Table of n, a(n) for n = 1..20000
MAPLE
N:= 1000: # for terms <= N
S:= {seq(i, i=3..N, 4)}:
P:= select(isprime, [seq(i, i=3..N/3, 2)]):
for i from 1 to nops(P) do
p:= P[i];
for j from i to nops(P) do
q:= P[j];
if 2*p*q + q^2 > N then break fi;
for k from j to nops(P) do
r:= P[k];
v:= p*q + p*r + q*r;
if v > N then break fi;
S:= S minus {v};
od od od:
sort(convert(S, list)); # Robert Israel, Apr 17 2024
CROSSREFS
Cf. also A369250 (4m+3 primes missing from this sequence).
KEYWORD
nonn
AUTHOR
Antti Karttunen, Jan 20 2024
STATUS
approved