A107715 - OEIS

%I #18 Mar 06 2016 11:09:13

%S 2,3,11,13,23,31,101,103,113,131,211,223,233,311,313,331,1013,1021,

%T 1031,1033,1103,1123,1201,1213,1223,1231,1301,1303,1321,2003,2011,

%U 2111,2113,2131,2203,2213,2221,2311,2333,3001,3011,3023,3121,3203,3221,3301,3313

%N Primes whose decimal representation contains only digits from the set {0,1,2,3}.

%H Reinhard Zumkeller, <a href="/A107715/b107715.txt">Table of n, a(n) for n = 1..1000</a>

%t Select[Prime[Range[500]],Max[IntegerDigits[#]]<4&] (* _Harvey P. Dale_, May 09 2012 *)

%t Select[FromDigits/@Tuples[{0,1,2,3},4],PrimeQ] (* _Harvey P. Dale_, Mar 06 2016 *)

%o (Haskell)

%o a107715 n = a107715_list !! (n-1)

%o a107715_list = filter ((== 1) . a010051) a007090_list

%o -- _Reinhard Zumkeller_, Aug 11 2011

%o (Python) from gmpy2 import digits

%o from sympy import isprime

%o [int(digits(n,4)) for n in range(1000) if isprime(int(digits(n,4)))] # _Chai Wah Wu_, Jul 31 2014

%Y Cf. A036953 (primes containing digits from set {0, 1, 2}).

%Y Cf. A010051, A007090, A193890.

%K base,nonn

%O 1,1

%A _Rick L. Shepherd_, May 22 2005