OFFSET
1,1
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 1..10431
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], j=1..nops(B)/2);
MATHEMATICA
Select[Prime[Range[1000]], PrimeQ[2 #+1]&][[2;; ;; 2]] (* Vincenzo Librandi, Dec 24 2025 *)
PROG
(Magma) A := []; for n in [2..1000] do p := NthPrime(n); if IsPrime(2*p + 1) then
Append(~A, p); end if; end for; B := A; C := [ B[2*j-1] : j in [1..(#B div 2)] ]; C; // 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
