OFFSET
1,1
MATHEMATICA
f[x_, y_]:=x^2+2^x+y^2+2^y; lst={}; Do[Do[p=f[x, y]; If[PrimeQ[p], AppendTo[lst, p]], {y, 4*5!}], {x, 4*5!}]; Take[Union[lst], 3*4! ]
PROG
(PARI) for(x=2, 100, for(y=1, x-1, p=x^2+2^x+y^2+2^y; if(isprime(p), print1(p", "))))
CROSSREFS
KEYWORD
nonn
AUTHOR
Vladimir Joseph Stephan Orlovsky, Jul 06 2009
EXTENSIONS
Program by Charles R Greathouse IV, Oct 12 2009
STATUS
approved