OFFSET
0,1
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