OFFSET
1,1
COMMENTS
All terms are multiples of 12.
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000 (terms 1..1000 from Harvey P. Dale)
FORMULA
k/2 +/-1 and also (k^3)/2 +/- 1 are twin primes.
EXAMPLE
1044=521+523, 1044^3=1137893184=568946591+568946593.
MATHEMATICA
sttpQ[n_]:=Module[{x=(n-2)/2}, AllTrue[{x, x+2}, PrimeQ]]; Select[Range[ 12, 300000, 12], sttpQ[#]&&sttpQ[#^3]&] (* The program uses the AllTrue function from Mathematica version 10 *) (* Harvey P. Dale, Jul 25 2017 *)
PROG
(PARI) isok(n) = !(n % 2) && isprime(n/2-1) && isprime(n/2+1) && isprime(n^3/2-1) && isprime(n^3/2+1); \\ Michel Marcus, Oct 19 2013
CROSSREFS
KEYWORD
nonn
AUTHOR
Zak Seidov, Jun 20 2012
STATUS
approved