A030086 - OEIS

%I #35 Nov 09 2022 19:24:25

%S 2,3,7,17,29,47,53,59,67,79,157,173,257,307,313,349,353,359,409,449,

%T 467,547,557,577,659,673,677,739,787,853,929,1447,1607,1747,2029,2087,

%U 2113,2237,3257,3359,3467,3533,3559,3767,3779,4253,4337,4787,5333,5557,5659

%N Primes p whose digits do not appear in p^2.

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

%H Michael S. Branicky, <a href="/A030086/b030086.txt">Table of n, a(n) for n = 1..13322</a> (terms < 10^16; terms 502..676 from Harvey P. Dale and Vincenzo Librandi, terms 1..501 from Harvey P. Dale)

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

%t Select[Prime[Range[700]],Intersection[IntegerDigits[#],IntegerDigits[ #^2]] == {}&] (* _Harvey P. Dale_, Oct 02 2014 *)

%o (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

%Y Cf. A029783 (numbers n whose digits are not present in n^2).

%Y Cf. similar sequences listed in A253574.

%K nonn,base

%O 1,1

%A _Patrick De Geest_, Dec 11 1999

%E Offset changed from 0 to 1 from _Vincenzo Librandi_, Jan 04 2015