|
|
A343471
|
|
Start of the first run of n or more consecutive primes using only prime digits.
|
|
1
|
|
|
|
OFFSET
|
1,1
|
|
COMMENTS
|
|
|
LINKS
|
|
|
EXAMPLE
|
a(1) = 2 because it is the first prime using only prime digits.
a(2) = 2 because 2, 3 is the first pair of consecutive primes using only prime digits.
a(5) = 2575723 because 2575723, 2575733, 2575753, 2575757, 2575777 is the first run of 5 consecutive primes using only prime digits.
|
|
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 = {i:2 for i in range(1, 5)}
for i in range(1, 5): yield 2
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:
for r in range(max(adict)+1, run+1):
adict[r] = t0
yield t0
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,base,hard,more
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|