OFFSET
1,1
COMMENTS
Subsequence of A066649.
LINKS
Charles R Greathouse IV, Table of n, a(n) for n = 1..10000
EXAMPLE
2 is in the sequence because 2 is prime, 2=1^2+1^3 and 1^3-1^2=0;
31 is in the sequence because 31 is prime, 31=2^2+3^3 and 3^3-2^2>0;
43 is in the sequence because 43 is prime, 43=4^2+3^3 and 3^3-4^2>0;
PROG
(PARI) list(lim)=my(v=List(), t, B); lim\=1; for(b=1, sqrtn(lim+.5, 3), B=b^3; for(a=1, min(lim-B, sqrtint(B)), if(isprime(t=B+a^2), listput(v, t)))); vecsort(Vec(v), , 8) \\ Charles R Greathouse IV, Aug 28 2011
CROSSREFS
KEYWORD
nonn
AUTHOR
Juri-Stepan Gerasimov, Aug 27 2011
EXTENSIONS
Corrected by D. S. McNeil, Aug 27 2011
STATUS
approved