OFFSET
1,1
COMMENTS
Primes of the form A087681(k)/7, any index k.
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 1..1000
EXAMPLE
For p=5, both 7*5-6=29 and 7*5+6=41 are prime,
for p=11, both 7*11-6=71 and 7*11+6=83 are prime.
MATHEMATICA
lst={}; Do[p=Prime[n]; If[PrimeQ[7*p-6]&&PrimeQ[7*p+6], AppendTo[lst, p]], {n, 6!}]; lst
Select[Prime[Range[700]], And @@ PrimeQ/@{7 # + 6, 7 # - 6}&] (* Vincenzo Librandi, Apr 09 2013 *)
PROG
(Magma) [p: p in PrimesUpTo(3000) | IsPrime(7*p-6) and IsPrime(7*p+6)]; // Vincenzo Librandi, Apr 09 2013
(PARI) is(n)=isprime(n) && isprime(7*n-6) && isprime(7*n+6) \\ Charles R Greathouse IV, Mar 28 2017
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Vladimir Joseph Stephan Orlovsky, Aug 16 2009
EXTENSIONS
Examples rephrased by R. J. Mathar, Aug 20 2009
STATUS
approved