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