|
%I #52 Aug 05 2022 07:45:13
%S 2,3,5,7,13,17,23,29,31,37,43,47,53,59,67,71,73,79,83,97,113,131,137,
%T 139,167,173,179,197,223,229,233,239,271,283,293,311,313,317,331,337,
%U 347,353,359,367,373,379,383,397,431,433,439,443,467,479,523,547,571
%N Left- or right-truncatable primes.
%C Repeatedly removing a digit from either the left or right produces only primes. There are 149677 terms in this sequence, ending with 8939662423123592347173339993799.
%C The number of n-digit terms is A298048(n). - _Jon E. Schoenfield_, Jan 28 2022
%H T. D. Noe, <a href="/A137812/b137812.txt">Table of n, a(n) for n = 1..10000</a>
%H I. O. Angell and H. J. Godwin, <a href="http://dx.doi.org/10.1090/S0025-5718-1977-0427213-2">On Truncatable Primes</a>, Math. Comput. 31, 265-267, 1977.
%H T. D. Noe, <a href="/A137812/a137812.png">Plot of all terms</a>
%H Carlos Rivera, <a href="http://www.primepuzzles.net/puzzles/puzz_002.htm">Puzzle 2: Prime Strings</a>, The Prime Puzzles and Problems Connection.
%H Daniel Starodubtsev, <a href="/A137812/a137812_1.txt">Full sequence</a>
%H Eric Weisstein, <a href="http://mathworld.wolfram.com/TruncatablePrime.html">MathWorld: Truncatable Prime</a>
%H <a href="/index/Tri#tprime">Index entries for sequences related to truncatable primes</a>
%e 139 is here because (removing 9 from the right) 13 is prime and (removing 1 from the left) 3 is prime.
%t Clear[s]; s[0]={2,3,5,7}; n=1; While[s[n]={}; Do[k=s[n-1][[i]]; Do[p=j*10^n+k; If[PrimeQ[p], AppendTo[s[n],p]], {j,9}]; Do[p=10*k+j; If[PrimeQ[p], AppendTo[s[n],p]], {j,9}], {i,Length[s[n-1]]}]; s[n]=Union[s[n]]; Length[s[n]]>0, n++ ];t=s[0]; Do[t=Join[t,s[i]], {i,n}]; t
%o (Python)
%o from sympy import isprime
%o def agen():
%o primes = [2, 3, 5, 7]
%o while len(primes) > 0:
%o yield from primes
%o cands = set(int(d+str(p)) for p in primes for d in "123456789")
%o cands |= set(int(str(p)+d) for p in primes for d in "1379")
%o primes = sorted(c for c in cands if isprime(c))
%o afull = [an for an in agen()]
%o print(afull[:60]) # _Michael S. Branicky_, Aug 04 2022
%Y Cf. A024770 (right-truncatable primes), A024785 (left-truncatable primes), A077390 (left-and-right truncatable primes), A080608.
%Y Cf. A298048 (number of n-digit terms).
%K base,fini,nonn
%O 1,1
%A _T. D. Noe_, Feb 11 2008
|
