1,1
a={}; Do[p1=Prime[n]; p2=Prime[n+1]; pp=p1^3+p2^2; If[PrimeQ[pp-1]&&PrimeQ[pp+1], AppendTo[a, p1]], {n, 16^3}]; Print[a];
Sequence in context: A002439 A132395 A064016 * A132937 A068655 A088385
Adjacent sequences: A138732 A138733 A138734 * A138736 A138737 A138738
nonn
Vladimir Orlovsky (4vladimir(AT)gmail.com), May 15 2008