OFFSET
1,1
COMMENTS
A number p is a term of this sequence iff p = q^2 + r^3 - 1 where p, q and r are all prime.
There are 1064 such numbers < 10^7.
EXAMPLE
173 is in the sequence because 173 = 7^2 + 5^3 - 1, and 173, 7, and 5 are all prime.
MATHEMATICA
max = 10^4; Take[Select[Sort[Flatten[Table[Prime[p]^2 + Prime[q]^3 - 1, {p, PrimePi[Sqrt[max]]}, {q, PrimePi[max^(1/3)]}]]], PrimeQ[#] && # < max &] (* Alonso del Arte, Mar 22 2013 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Kevin L. Schwartz & Christian N. K. Anderson, Mar 22 2013
STATUS
approved