%I #9 Apr 19 2014 19:16:33
%S 5,23,53,113,173,197,269,317,383,443,557,563,587,647,659,773,797,827,
%T 947,983,1097,1103,1187,1217,1229,1889,1913,1949,2039,2099,2153,2213,
%U 2339,2357,2399,2417,2447,2579,2693,2837,2879,2897,2903,2939,2969,3089,3203
%N Primes p such that (p^3 + 4)/3 is prime.
%H K. D. Bajpai, <a href="/A241099/b241099.txt">Table of n, a(n) for n = 1..8857</a>
%e 5 is prime and appears in the sequence because (5^3 + 4)/3 = 43 which is a prime.
%e 23 is prime and appears in the sequence because (23^3 + 4)/3 = 4057 which is a prime.
%p KD:= proc() local a,b;a:=ithprime(n); b:=(a^3+4)/3; if b=floor(b) and isprime(b) then RETURN (a); fi; end: seq(KD(), n=1..1000);
%t Select[Prime[Range[500]], PrimeQ[(#^3 + 4)/3] &]
%t n = 0; Do[If[PrimeQ[(Prime[k]^3 + 4)/3], n = n + 1; Print[n, " ", Prime[k]]], {k, 1, 200000}] (* b-file *)
%Y Cf. A109953 (primes p:(p^2+1)/3 is prime).
%Y Cf. A118915 (primes p:(p^2+5)/6 is prime).
%Y Cf. A118918 (primes p:(p^2+11)/12 is prime).
%K nonn
%O 1,1
%A _K. D. Bajpai_, Apr 15 2014
|