1,1

NextPrime[n_Int]:=Module[{k}, k=n+1; While[ !PrimeQ[k], k++ ]; k]; PrevPrime[n_Int]:=Module[{k}, k=n-1; While[ !PrimeQ[k], k-- ]; k]; s=""; For[i=2, i< 10^5*5, p=Prime[i]; If[(Prime[i-1]+Prime[i+1])/2==p, r=p*(p-1)+1; a=PrevPrime[r]; b=NextPrime[r]; If[PrimeQ[r]&&r==(a+b)/2, (*Print[p, ":", a, ", ", b, "; ", r]*)s=s<>ToString[p]<>", "]]; i++ ]; Print[s]

Cf. A006562.

Sequence in context: A029566 A051651 A065324 this_sequence A086691 A112625 A118498

Adjacent sequences: A138261 A138262 A138263 this_sequence A138265 A138266 A138267

nonn

Vladimir Orlovsky (4vladimir(AT)gmail.com), May 05 2008