|
|
A352312
|
|
Start of the first run of exactly n consecutive primes using only prime digits.
|
|
1
|
|
|
|
OFFSET
|
1,1
|
|
COMMENTS
|
|
|
LINKS
|
|
|
EXAMPLE
|
a(1) = 7 because it is the first prime using only prime digits and whose next prime 11 does not use only prime digits.
a(3) = 37253 because 37253, 37273, 37277 is the first run of 3 consecutive primes using only prime digits, then next prime 37307 has a digit 0.
|
|
PROG
|
(Python)
from sympy import nextprime, isprime
from itertools import count, islice, product
def onlypd(n): return set(str(n)) <= set("2357")
def agen():
adict, n = {1:7, 4:2}, 1
yield 7
for digits in count(2):
for p in product("2357", repeat=digits-1):
for end in "37":
t0 = t = int("".join(p) + end)
run = 0
while isprime(t):
run += 1
t = nextprime(t)
if not onlypd(t): break
if run not in adict:
adict[run] = t0
if run > n:
for r in range(n+1, run+1):
if r in adict:
yield adict[r]
n += 1
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,base,more
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|