login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A030087 Primes such that digits of p do not appear in p^3. 2

%I #27 May 06 2022 13:13:51

%S 2,3,7,43,47,53,157,223,263,487,577,587,823,4657,5657,6653,7177,8287,

%T 9343,26777,36293,46477,58787,72727,75707,176777,363313,530353,566653,

%U 959953,1771787,2525557,2555353,2626277,3656363,4414447,7110707,8448343,20700077,54475457,71117177,72722977,135135113,393321293,457887457,505053053,672722627

%N Primes such that digits of p do not appear in p^3.

%C Primes of sequence A029785. - _Michel Marcus_, Jan 04 2015

%H Michael S. Branicky, <a href="/A030087/b030087.txt">Table of n, a(n) for n = 1..86</a> (terms < 10^19; terms 74..86 via A029785)

%e 2 and 2^3=8 have no digits in common, hence 2 is in the sequence.

%t Select[Prime[Range[1500000]], Intersection[IntegerDigits[#], IntegerDigits[#^3]]=={} &] (* _Vincenzo Librandi_, Jan 04 2015 *)

%o (PARI) lista(nn) = {forprime (n=1, nn, if (#setintersect(Set(vecsort(digits(n^3))), Set(vecsort(digits(n)))) == 0, print1(n, ", ")); ); } \\ _Michel Marcus_, Jan 04 2015

%o (Python)

%o from sympy import isprime

%o A030087_list = [n for n in range(1,10**6) if set(str(n)) & set(str(n**3)) == set() and isprime(n)]

%o # _Chai Wah Wu_, Jan 05 2015

%Y Cf. A029785 (digits of n are not present in n^3), A030086 (similar, with p^2), A253574 (similar, with p^4).

%K nonn,base

%O 1,1

%A _Patrick De Geest_, Dec 11 1999

%E Changed offset from 0 to 1 and more terms from _Vincenzo Librandi_, Jan 04 2015

%E a(40)-a(47) from _Chai Wah Wu_, Jan 05 2015

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

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 28 08:22 EDT 2024. Contains 371236 sequences. (Running on oeis4.)