OFFSET
1,1
COMMENTS
a(n) < 10 iff prime(n) is in A179336;
a(n) = prime(n) iff prime(n) is in A033274. [Corrected by M. F. Hasler, Aug 27 2012]
LINKS
Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
PROG
(PARI) A179335(n)={my(p=prime(n), m=0, M); for(d=1, n, M=10^d; n=p; until(n<=M || !n\=10, isprime(n%M) & (!m || m>n%M) & m=n%M); m & return(m))} \\ M. F. Hasler, Aug 27 2012
(Python)
from sympy import isprime, prime
def a(n):
s = str(prime(n))
ss = set(int(s[i:i+1+l]) for i in range(len(s)) for l in range(len(s)))
return min(t for t in ss if isprime(t))
print([a(n) for n in range(1, 94)]) # Michael S. Branicky, Jun 29 2022
CROSSREFS
KEYWORD
base,nonn
AUTHOR
Reinhard Zumkeller, Jul 11 2010
STATUS
approved