OFFSET
1,1
MATHEMATICA
t = Table[Mod[Prime[nn], 3], {nn, 10^8}]; (* A039701 *)
u = Map[SequencePosition[t, Flatten[Join[{ConstantArray[2, #], 1}]], 1] &, Range[21]]
Prime[Table[First[Flatten[u[[n]]]], {n, 1, 21}]]
(* Peter J. C. Moses, Dec 11 2025 *)
PROG
(Python)
from sympy import nextprime, prime
from itertools import count, islice
def agen(NN=10000): # generator of terms for n <= NN
n, pk, residues = 1, 1, "X"*NN
for k in count(1):
pk = nextprime(pk)
resk = pk%3
if resk == 1:
while residues.endswith("2"*n) and n <= NN:
yield prime(k-n)
n += 1
residues = residues[1:] + str(resk)
print(list(islice(agen(), 12))) # Michael S. Branicky, Jan 11 2026
CROSSREFS
KEYWORD
nonn
AUTHOR
Clark Kimberling, Jan 08 2026
EXTENSIONS
a(22)-a(27) from Michael S. Branicky, Jan 11 2026
a(28) from Jinyuan Wang, Jan 16 2026
STATUS
approved
