OFFSET
1,1
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 1..41118
MAPLE
A:={}: for n from 1 to 608 do if isprime(2*ithprime(n)+1)=true then A:=A union {ithprime(n)} else A:=A fi od: B:=convert(A, list): seq(B[2*j-1], j=1..nops(B)/2);
MATHEMATICA
Select[Prime[Range[1000]], PrimeQ[2 #+1]&][[1;; ;; 2]] (* Vincenzo Librandi, Dec 24 2025 *)
PROG
(Magma) [primes[i] : i in [1..#primes by 2]] where primes := [p : p in PrimesUpTo(5000) | IsPrime(2*p+1)]; // Vincenzo Librandi, Dec 24 2025
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
N. J. A. Sloane, Nov 20 2004
EXTENSIONS
More terms from Emeric Deutsch, Dec 09 2004
STATUS
approved
