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: A:= NULL;
for i from 0 to 100 do
found:= false;
for j from 1 to nP do
if not isprime(s+P[j]) then
found:= true; A:= A, min(numtheory:-factorset(s+P[j])); 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:
A; # Robert Israel, Jun 06 2025
PROG
(PARI) 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(factor(v[n-2]+v[n-1]+p)[1, 1]", "):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
