%I #13 May 18 2024 19:36:23
%S 4,9,49,289,841,2209,2809,3481,4489,6241,24649,29929,66049,94249,
%T 97969,121801,124609,128881,167281,201601,218089,299209,310249,332929,
%U 434281,452929,458329,546121,619369,727609,863041,2093809
%N a(n) = prime^2 and digits of prime do not appear in a(n).
%H Harvey P. Dale, <a href="/A030088/b030088.txt">Table of n, a(n) for n = 1..442</a> (* All terms up to the 10 millionth prime. *)
%F a(n) = A030086(n)^2. - _Robert Israel_, Nov 02 2022
%p R:= NULL: count:= 0: p:= 0:
%p while count < 40 do
%p p:= nextprime(p);
%p ps:= p^2;
%p if convert(convert(p,base,10),set) intersect convert(convert(ps,base,10),set) = {} then
%p R:= R, ps; count:= count+1
%p fi
%p od:
%p R; # _Robert Israel_, Nov 02 2022
%t Select[Prime[Range[250]],Intersection[IntegerDigits[#],IntegerDigits[#^2]]=={}&]^2 (* _Harvey P. Dale_, May 18 2024 *)
%Y Cf. A030086.
%K nonn,base
%O 1,1
%A _Patrick De Geest_
|