A362066 - OEIS

%I #26 Apr 18 2023 10:24:38

%S 17,12491,14723,42437,57089,58193,61051,63131,63347,64553,64567,64577,

%T 64591,64601,64661,64679,64951,65071,65173,65293,65881,66863,69931,

%U 79817,99551,129083,165103,263071,284833,1407647,1515259,4303027,6440999,14968819,95517973,527737957,1893230839,1950929941,1964567161

%N Primes associated with the indices in A362060.

%C Is this the same as A114924, or are there base-10 expressions of pi(p) which become p after striking 2 or more digits? - _R. J. Mathar_, Apr 18 2023

%H Chai Wah Wu, <a href="/A362066/b362066.txt">Table of n, a(n) for n = 1..54</a> (n=1..50 from Jean-Marc Rebert)

%o (Python)

%o from sympy import sieve

%o def okA362060(n):

%o     p = sieve[n]

%o     while n and p:

%o         if n%10 == p%10:

%o             n //= 10

%o         p //= 10

%o     return n == 0

%o print([sieve[k] for k in range(1, 10**6) if okA362060(k)]) # _Michael S. Branicky_, Apr 07 2023

%o (Python)

%o from sympy import prime, nextprime

%o from itertools import count, islice

%o def A362066_gen(startvalue=1): # generator of terms >= startvalue

%o     p = prime(max(startvalue,1))

%o     for k in count(max(startvalue,1)):

%o         c = iter(str(p))

%o         if all(map(lambda b:any(map(lambda a:a==b,c)),str(k))):

%o             yield p

%o         p = nextprime(p)

%o A362066_list = list(islice(A362066_gen(),20)) # _Chai Wah Wu_, Apr 07 2023

%Y Cf. A362060.

%K nonn,base

%O 1,1

%A _Jean-Marc Rebert_, Apr 07 2023