|
%I #15 Jan 24 2018 03:04:48
%S 59,157,251,24421,161059,2248123,6436351,6967903,15813283,20511157,
%T 22188073,58863901,86938339,90518881,131872261,263374753,440711113,
%U 553387693,865523209,1095912823,1194390013,1524845983,1573037779,2521008913,2979767551,3327970189
%N Primes of the form p^3 + q^5 where p and q are primes.
%C p and q cannot both be odd. Thus p=2 or q=2.
%H Robert Israel, <a href="/A122702/b122702.txt">Table of n, a(n) for n = 1..10000</a>
%F {a(n)} = {p^3 + q^5 in A000040 where p and q are in A000040}.
%e a(1) = 3^3 + 2^5 = 59.
%e a(2) = 5^3 + 2^5 = 157.
%e a(3) = 2^3 + 3^5 = 251.
%e a(4) = 29^3 + 2^5 = 24421.
%p N:= 10^10: # to get all terms <= N
%p A:= select(isprime,{seq(2^3 + q^5, q = select(isprime, [2,seq(i,i=3..floor((N-8)^(1/5)),2)])),
%p seq(2^5 + q^3, q = select(isprime, [2,seq(i,i=3..floor((N-2^5)^(1/3)),2)]))}):
%p sort(convert(A,list)); # _Robert Israel_, Jan 23 2018
%t With[{nn=10^8},Union[Select[Join[Prime[Range[Floor[Power[nn, (5)^-1]]]]^5+ 8,Prime[Range[Floor[Power[nn, (3)^-1]]]]^3+32],PrimeQ]]] (* _Harvey P. Dale_, Nov 26 2011 *)
%Y Cf. A000040, A045700 (Primes of form p^2+q^3 where p and q are primes).
%K easy,nonn
%O 1,1
%A _Jonathan Vos Post_, Sep 22 2006
%E More terms from _Harvey P. Dale_, Nov 26 2011
|
