|
%I
%S 11,19,41,61,89,101,109,149,181,191,199,401,409,419,449,461,491,499,
%T 601,619,641,661,691,809,811,881,911,919,941,991,1009,1019,1049,1061,
%U 1069,1091,1109,1181,1409,1481,1489,1499,1601,1609,1619,1669,1699,1801
%N Primes with only nonprime decimal digits.
%C A109066(n) = 0 iff prime(n) is in this sequence. [From _Reinhard Zumkeller_, Jul 11 2010, corrected by _M. F. Hasler_, Aug 27 2012]
%C Or, primes p such that A193238(p) = 0. - _M. F. Hasler_, Aug 27 2012
%C Intersection of A084984 and A000040; complement of A179336 (within the primes A000040). [_Reinhard Zumkeller_, Jul 19 2011, edited by _M. F. Hasler_, Aug 27 2012]
%H T. D. Noe, <a href="/A034844/b034844.txt">Table of n, a(n) for n = 1..1000</a>
%H Chris Caldwell, <a href="http://primes.utm.edu/lists/small/1000.txt">The First 1,000 Primes</a>
%F a(n) >> n^1.285. [_Charles R Greathouse IV_, Feb 20 2012]
%e E.g. 149 is a prime made of nonprime digits(1,4,9).
%e 991 is a prime without any prime digits.
%t Select[Prime[Range[279]], Intersection[IntegerDigits[#], {2, 3, 5, 7}] == {} &] (* Jayanta Basu, Apr 18 2013 *)
%o (Haskell)
%o a034844 n = a034844_list !! (n-1)
%o a034844_list = filter (not . any (`elem` "2357") . show ) a000040_list
%o -- _Reinhard Zumkeller_, Jul 19 2011
%o (MAGMA) [p: p in PrimesUpTo(2000) | forall{d: d in [2,3,5,7] | d notin Set(Intseq(p))}]; // Bruno Berselli, Jul 27 2011
%o (PARI) is_A034844(n)=isprime(n)&!apply(x->isprime(x),eval(Vec(Str(n)))) \\ - _M. F. Hasler_, Aug 27 2012
%o (PARI) is_A034844(n)=isprime(n)&!setintersect(Set(Vec(Str(n))),Vec("2357")) \\ - _M. F. Hasler_, Aug 27 2012
%Y Cf. A000040, A084984.
%Y Cf. A034844, A033274, A152313, A141409.
%K nonn,base
%O 1,1
%A _Erich Friedman_
%E Edited by _N. J. A. Sloane_, Feb 22 2009 at the suggestion of _R. J. Mathar_
|
