%I #33 Jun 18 2022 13:40:08
%S 3,5,7,11,13,17,19,31,37,53,59,71,73,79,97,113,131,137,139,151,157,
%T 173,179,191,193,197,199,311,313,317,331,337,353,359,373,379,397,557,
%U 571,577,593,599,719,733,739,751,757,773,797,911,919,937,953,971,977,991
%N Primes whose digits are all odd.
%H Seiichi Manyama, <a href="/A030096/b030096.txt">Table of n, a(n) for n = 1..10000</a> (terms 1..1000 from T. D. Noe)
%t Select[Prime[Range[500]],And@@OddQ[IntegerDigits[#]]&] (* _Harvey P. Dale_, Jan 28 2013 *)
%o (Haskell)
%o a030096 n = a030096_list !! (n-1)
%o a030096_list = filter f a000040_list where
%o f x = odd d && (x < 10 || f x') where (x', d) = divMod x 10
%o -- _Reinhard Zumkeller_, Apr 07 2014, Jan 29 2013
%o (PARI) is(n)=isprime(n) && #setintersect([0,2,4,6,8],Set(digits(n)))==0 \\ _Charles R Greathouse IV_, Feb 07 2017
%o (Magma) [p: p in PrimesUpTo(1000) | forall{d: d in [0,2,4,6,8] | d notin Set(Intseq(p))}]; // _Vincenzo Librandi_, Apr 29 2019
%Y Cf. A076704 = odd-digit prime powers of prime numbers; A091296 = odd-digit semiprimes; A000040 = prime numbers; A001358 = semiprimes.
%Y Intersection of A000040 and A014261.
%K nonn,base,nice
%O 1,1
%A _Patrick De Geest_
%E Edited by _N. J. A. Sloane_ at the suggestion of _T. D. Noe_ and _Jonathan Vos Post_, Sep 15 2007
|