A030086 Primes p whose digits do not appear in p^2.
2, 3, 7, 17, 29, 47, 53, 59, 67, 79, 157, 173, 257, 307, 313, 349, 353, 359, 409, 449, 467, 547, 557, 577, 659, 673, 677, 739, 787, 853, 929, 1447, 1607, 1747, 2029, 2087, 2113, 2237, 3257, 3359, 3467, 3533, 3559, 3767, 3779, 4253, 4337, 4787, 5333, 5557, 5659
Offset: 1
Examples
2 and 2^2 = 4 have no digits in common, hence 2 is in the sequence.
Links
- Michael S. Branicky, Table of n, a(n) for n = 1..13322 (terms < 10^16; terms 502..676 from Harvey P. Dale and Vincenzo Librandi, terms 1..501 from Harvey P. Dale)
Crossrefs
Programs
-
Mathematica
Select[Prime[Range[700]],Intersection[IntegerDigits[#],IntegerDigits[ #^2]] == {}&] (* Harvey P. Dale, Oct 02 2014 *)
-
PARI
lista(nn) = {forprime (n=1, nn, if (#setintersect(Set(vecsort(digits(n^2))), Set(vecsort(digits(n)))) == 0, print1(n, ", ")); ); } \\ Michel Marcus, Jan 04 2015
Extensions
Offset changed from 0 to 1 from Vincenzo Librandi, Jan 04 2015
Comments