OFFSET
1,1
COMMENTS
LINKS
Robert Israel, Table of n, a(n) for n = 1..1000
EXAMPLE
a(1)=18 because 1 + 30*k for k=18..23 are 541, 571, 601, 631, 661, 691 all primes: A000040(k) for k={100, 105, 110, 115, 121, 125}.
MAPLE
filter:= t -> andmap(isprime, [seq(30*(t+k)+1, k=0..5)]):
select(filter, [seq(seq(77*k + i, i=[18, 39, 53, 60, 74]), k=0..10000)]); # Robert Israel, Apr 04 2017
MATHEMATICA
Select[Range[18, 1000000, 7], PrimeQ[1 + 30*#] && PrimeQ[1 + 30*(# + 1)] && PrimeQ[1 + 30*(# + 2)] && PrimeQ[1 + 30*(# + 3)] && PrimeQ[1 + 30*(# + 4)] && PrimeQ[1 + 30*(# + 5)] &]
Select[Range[4, 10^6, 7], AllTrue[30(#+Range[0, 5])+1, PrimeQ]&] (* Harvey P. Dale, Dec 03 2023 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Zak Seidov, Mar 31 2017
STATUS
approved