OFFSET
1,1
LINKS
J. K. Andersen, Consecutive Congruent Primes.
EXAMPLE
a(4) = 20809 because this number is the first in a sequence of 4 consecutive primes all of the form 8n + 1.
MATHEMATICA
NextPrime[ n_Integer ] := Module[ {k = n + 1}, While[ ! PrimeQ[ k ], k++ ]; Return[ k ] ]; PrevPrime[ n_Integer ] := Module[ {k = n - 1}, While[ ! PrimeQ[ k ], k-- ]; Return[ k ] ]; p = 0; Do[ a = Table[ -1, {n} ]; k = Max[ 1, p ]; While[ Union[ a ] != {1}, k = NextPrime[ k ]; a = Take[ AppendTo[ a, Mod[ k, 8 ] ], -n ] ]; p = NestList[ PrevPrime, k, n ]; Print[ p[ [ -2 ] ] ]; p = p[ [ -1 ] ], {n, 1, 9} ] a(10) > 123700000.
Module[{nn=75100000, pr8}, pr8=Table[If[Mod[p, 8]==1, 1, 0], {p, Prime[Range[nn]]}]; Prime[#]&/@Table[SequencePosition[ pr8, PadRight[ {}, n, 1], 1], {n, 10}]][[;; , 1, 1]] (* The program generates the first 10 terms of the sequence. *) (* Harvey P. Dale, Sep 29 2024 *)
CROSSREFS
KEYWORD
nonn,changed
AUTHOR
Robert G. Wilson v, Oct 10 2000
EXTENSIONS
More terms from Jens Kruse Andersen, May 28 2006
a(16)-a(18) from Giovanni Resta, Aug 04 2013
STATUS
approved