|
%I #33 Mar 05 2022 14:18:08
%S 2,3,5,7,11,13,31,37,41,43,67,73,79,101,137,239,271,353,449,757,859,
%T 1933,4649,8779,9091,9901,21401,21649,25601,27961,52579,62003,123551,
%U 333667,513239,538987,909091,1676321,2071723,2906161,5882353,10838689,35121409,52986961,99990001,265371653,1056689261,1058313049,1360682471
%N Primes p such that the decimal representation of 1/p does not contain every digit 0-9.
%C Every repunit prime (A004022) is here. There are 113 terms of A046107, having periods of up to 256, that are here. The only known unique-period prime (A007615) not here is the one having period 92092. Is this sequence finite? - _T. D. Noe_, Mar 13 2011
%e 4649 is in the sequence because 1/4649 = 0.00021510002151000215.... contain
%e only the digits 0, 1, 2 and 5.
%t Join[{2, 3, 5}, Select[Prime[Range[4, 10000]], Length[Union[RealDigits[1/#][[1, 1]]]] < 10 &]]
%o (Python)
%o from sympy import n_order, nextprime
%o from itertools import islice
%o def A187614_gen(): # generator of terms
%o yield from (2,3,5)
%o p = 7
%o while True:
%o if len(set('0'+str(10**(n_order(10, p))//p))) < 10:
%o yield p
%o p = nextprime(p)
%o A187614_list = list(islice(A187614_gen(),20)) # _Chai Wah Wu_, Mar 03 2022
%Y Cf. A187372.
%Y Cf. A352023 (does not contain digit 9)
%K nonn,base
%O 1,1
%A _Michel Lagneau_, Mar 12 2011
%E Extended by _T. D. Noe_, Mar 12 2011
|
