OFFSET
1,1
COMMENTS
These primes are such that p=1 (mod 4) and 6p+1 is prime, but there are other primes with these properties (13, 17, ...) not in this sequence.
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000 (terms 1..1313 from M. F. Hasler)
MATHEMATICA
Select[Range[10^4], PrimeQ[#] && PowerMod[2, #, 6# + 1] == 1 &] (* Amiram Eldar, Nov 13 2019 *)
PROG
(PARI) forprime(p=1, 1e5, Mod(2, p*6+1)^p-1||print1(p", "))
CROSSREFS
KEYWORD
nonn
AUTHOR
M. F. Hasler, Mar 21 2011
STATUS
approved