OFFSET
1,2
LINKS
Zak Seidov, Table of n, a(n) for n = 1..10000
EXAMPLE
3 belongs to the sequence because 4*3+1 and 4*3-1 are both primes.
MATHEMATICA
Select[Range[1023], And @@ PrimeQ[{-1, 1} + 4# ] &] (* Ray Chandler, Dec 06 2006 *)
PROG
(Magma) [n: n in [1..2000] | IsPrime(4*n+1) and IsPrime(4*n-1)] // Vincenzo Librandi, Nov 18 2010
(PARI) list(lim)=my(v=List(), p=2); forprime(q=3, 4*lim+1, if(q-p==2 && p%4==3, listput(v, q\4)); p=q); Vec(v) \\ Charles R Greathouse IV, Dec 03 2016
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
EXTENSIONS
More terms from Erich Friedman
STATUS
approved