OFFSET
0,1
LINKS
Robert Israel, Table of n, a(n) for n = 0..10000
MAPLE
P:= select(isprime, [seq(i, i=7..1000, 2)]): nP:= nops(P):
R:= 3, 5: s:= 8:
for i from 3 to 100 do
found:= false;
for j from 1 to nP do
if not isprime(s+P[j]) then
found:= true; R:= R, P[j]; s:= R[-2]+P[j]; P:= subsop(j=NULL, P); nP:= nP-1; break
fi;
od;
if not found then break fi;
od:
R; # Robert Israel, Jun 06 2025
PROG
(PARI) print1("3, 5, "); pr=vector(10000); v=vector(100); v[1]=3; v[2]=5; pr[3]=1; pr[5]=1; for(n=3, 100, forprime(p=3, 10000, if(!pr[p]&&!isprime(v[n-2]+v[n-1]+p), print1(p", "); v[n]=p; pr[p]=1; break)))
CROSSREFS
KEYWORD
nonn
AUTHOR
Amarnath Murthy, Sep 25 2002
EXTENSIONS
Corrected and extended by Ralf Stephan, Mar 27 2003
STATUS
approved
