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, pp]], {n, 16^3}]; Print[a];
Sequence in context: A028385 A190947 A203910 * A139777 A015327 A079230
Adjacent sequences: A138759 A138760 A138761 * A138763 A138764 A138765
nonn
Vladimir Orlovsky (4vladimir(AT)gmail.com), May 15 2008