OFFSET
1,1
COMMENTS
These are the integers of the form (3*A014574(i))^(1/3), any index i. - R. J. Mathar, Dec 14 2008
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
EXAMPLE
6 is a term since (6^3)/3 = 72 and (71, 73) are twin primes.
30 is a term since (30^3)/3 = 9000 and (8999, 9001) are twin primes.
MATHEMATICA
lst1={}; lst2={}; Do[ p1=Prime[n]; p2=Prime[n+1]; If[p2-p1==2, e=(3*(p1+1))^(1/3); i=Floor[e]; If[e==i, AppendTo[lst1, (p1+1)]; AppendTo[lst2, i]]], {n, 2*10!}]; Print[lst1]; Print[lst2]
fQ[n_] := PrimeQ[n^3/3 - 1] && PrimeQ[n^3/3 + 1]; lst = {}; Do[If[fQ@n, AppendTo[lst, n]], {n, 3, 10^4, 3}]; lst
PROG
(Magma) [k:k in [3..9000 by 3]| IsPrime(k^3 div 3 -1) and IsPrime(k^3 div 3 +1)]; // Marius A. Burtea, Jan 01 2020
CROSSREFS
KEYWORD
nonn
AUTHOR
Vladimir Joseph Stephan Orlovsky, Dec 12 2008
EXTENSIONS
Edited and extended by Robert G. Wilson v, Dec 14 2008
Corrected divisor in definition. - R. J. Mathar, Dec 20 2008
STATUS
approved