OFFSET
1,1
LINKS
Michael De Vlieger, Table of n, a(n) for n = 1..10000
EXAMPLE
5 is a term because (5-1)/2 = 2 is prime and 2*5+1 = 11 is prime.
MATHEMATICA
Select[Range[1, 2000, 2], AllTrue[{(# - 1)/2, 2 # + 1}, PrimeQ] &] (* Michael De Vlieger, Apr 19 2024 *)
PROG
(Python)
from sympy import isprime
def a(n): return n%2 == 1 and isprime((n-1)>>1) and isprime(2*n+1)
print([n for n in range(2, 1900) if a(n)])
CROSSREFS
KEYWORD
nonn
AUTHOR
Alexandre Herrera, Apr 19 2024
STATUS
approved