login
The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A190222 Noncomposite numbers all of whose decimal digits are noncomposite numbers (1,2,3,5,7). 3

%I

%S 1,2,3,5,7,11,13,17,23,31,37,53,71,73,113,127,131,137,151,157,173,211,

%T 223,227,233,251,257,271,277,311,313,317,331,337,353,373,521,523,557,

%U 571,577,727,733,751,757,773,1117,1123,1151,1153,1171,1213,1217,1223

%N Noncomposite numbers all of whose decimal digits are noncomposite numbers (1,2,3,5,7).

%H G. C. Greubel, <a href="/A190222/b190222.txt">Table of n, a(n) for n = 1..1000</a>

%H Chris Caldwell, The Prime Glossary, <a href="http://primes.utm.edu/glossary/xpage/UnholeyPrime.html">Unholey prime</a>

%F a(n) >> n^k where k = log(10)/log(5) = 1.43067.... - _Charles R Greathouse IV_, Dec 28 2017

%t Join[{1}, Select[Range[2000], PrimeQ[#] && Intersection[{0, 4, 6, 8, 9}, IntegerDigits[#]] == {} &]] (* _T. D. Noe_, May 09 2011 *)

%o (MAGMA) [1] cat [n: n in PrimesUpTo(1223) | Set(Intseq(n)) subset [1, 2, 3, 5, 7]]; // _Arkadiusz Wesolowski_, Apr 16 2014

%o (PARI) is(k) = if(!isprime(k) && k != 1, return(0)); setminus(vecsort(digits(k), , 8), [1, 2, 3, 5, 7]) == [] \\ _Iain Fox_, Dec 28 2017

%o (PARI) is(n) = if(isprime(n), #setminus(Set(digits(k)), [1,2,3,5,7])==0, n==1) \\ _Charles R Greathouse IV_, Dec 28 2017

%Y Subsequence of A001742.

%K nonn,base

%O 1,2

%A _Jaroslav Krizek_, May 06 2011

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified September 29 06:29 EDT 2020. Contains 337425 sequences. (Running on oeis4.)