OFFSET
1,1
COMMENTS
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
MATHEMATICA
b = {}; a = {}; Do[If[PrimeQ[((Prime[n + 2] + Prime[n + 1])/2 + (Prime[n + 1] + Prime[n])/2)/2], AppendTo[a, ((Prime[n + 2] + Prime[n + 1])/2 + (Prime[n + 1] + Prime[n])/2)/2]], {n, 1, 10000}]; Do[If[PrimeQ[(a[[k + 1]] + a[[k]])/2], AppendTo[b, (a[[k + 1]] + a[[k]])/2]], {k, 1, Length[a] - 1}]; b
PROG
(PARI) {m=250000; a=0; p=2; q=3; r=5; while(r<=m, if((p+r)/2==q, if(a>0, if(isprime(b=(a+q)/2), print1(b, ", "))); a=q); p=q; q=r; r=nextprime(r+1))} \\ Klaus Brockhaus, Jan 05 2007
CROSSREFS
KEYWORD
nonn
AUTHOR
Artur Jasinski, Dec 27 2006
EXTENSIONS
Edited and extended by Klaus Brockhaus, Jan 05 2007
STATUS
approved