OFFSET
1,1
COMMENTS
LINKS
Charles R Greathouse IV, Table of n, a(n) for n = 1..10000
FORMULA
a(n) == 1 (mod 60).
EXAMPLE
61 is a term because 61 and 73 are both terms in A005383.
MATHEMATICA
Select[Range@ 160000, Mod[#, 60] == 1 && PrimeQ[#] && PrimeQ[# + 12] && PrimeQ[(# + 1)/2] && PrimeQ[(# + 13)/2] &] (* Indranil Ghosh, Mar 17 2017 *)
Select[Range[1, 152000, 60], AllTrue[{#, #+12, (#+1)/2, (#+13)/2}, PrimeQ]&] (* Harvey P. Dale, Oct 14 2021 *)
PROG
(PARI) is(n)=n%60==1 && isprime(n) && isprime(n+12) && isprime((n+1)/2) && isprime((n+13)/2) \\ Charles R Greathouse IV, Mar 17 2017
(Python)
from sympy import isprime
[i for i in range(1, 160001) if i%60 == 1 and isprime(i) and isprime(i + 12) and isprime((i + 1)/2) and isprime((i + 13)/2)] # Indranil Ghosh, Mar 17 2017
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved