OFFSET
1,1
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 1..1000
EXAMPLE
3 is a term because 3 is prime and 6*3^2+6*3-1 is 71 which is prime. 13 is a term because 13 is prime and 6*13^2+6*13-1 is 1091 which is prime. - Soumil Mandal, Apr 14 2016
MATHEMATICA
Select[Prime[Range[300]], PrimeQ[6 #^2 + 6 # - 1] &]
PROG
(Magma) [p: p in PrimesUpTo(2000) | IsPrime(6*p^2+6*p-1)];
(PARI) lista(nn) = forprime(p=2, nn, if(isprime(6*p^2+6*p-1), print1(p, ", "))); \\ Altug Alkan, Apr 12 2016
(Python)
from gmpy2 import is_prime
for p in range(3, 10**5, 2):
if(not is_prime(p)):continue
elif(is_prime(6*p**2+6*p-1)):print(p)
# Soumil Mandal, Apr 14 2016
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Vincenzo Librandi, Apr 12 2016
STATUS
approved