%I #21 Jan 01 2024 12:24:02
%S 20,30,50,70,101,103,107,109,110,130,170,190,200,203,209,230,290,300,
%T 301,307,310,370,401,403,407,410,430,470,500,503,509,530,590,601,607,
%U 610,670,700,701,703,709,710,730,790,803,809,830,890,907,970,1001,1003,1007,1009,1010,1013,1027,1030
%N Naught-y numbers (A011540) that after removing all zeros become zeroless primes (A038618).
%H Harvey P. Dale, <a href="/A256227/b256227.txt">Table of n, a(n) for n = 1..1000</a>
%p N:= 4: # to produce all terms with <= N digits
%p ZLO:= proc(d) # produce set of d-digit odd zeroless numbers
%p option remember;
%p if d = 1 then {1,3,5,7,9}
%p else
%p map(t -> seq(t+x*10^(d-1),x=1..9), ZLO(d-1))
%p fi
%p end proc:
%p addzeros:= proc(x,d) # d-digit numbers formed by inserting 0's into x
%p local L,n,R;
%p L:= convert(x,base,10);
%p n:= nops(L);
%p R:= map(t -> [op(t),d], combinat[choose](d-1,n-1));
%p seq(add(L[i]*10^(r[i]-1),i=1..n), r = R);
%p end proc:
%p Z[1]:= {2,3,5,7}:
%p for i from 2 to N-1 do Z[i]:= select(isprime,ZLO(i)) od:
%p `union`(seq(seq(map(addzeros,Z[i],d), i=1..d-1),d=2..N));
%p # if using Maple 11 or earlier, uncomment the next line
%p # sort(convert(%,list)); # _Robert Israel_, Mar 19 2015
%t ss={};Do[id=IntegerDigits[p];If[Min[id]<1&&PrimeQ[FromDigits[Delete[id,Position[id,0]]]],ss={ss,p}],{p,20,2000}];Flatten[ss]
%t Select[Range[1200],DigitCount[#,10,0]>0&&PrimeQ[FromDigits[DeleteCases[ IntegerDigits[ #],0]]]&] (* _Harvey P. Dale_, Jan 01 2024 *)
%o (PARI) is(n)=my(d=digits(n),e=select(x->x,d)); #e<#d && isprime(fromdigits(e)) \\ _Charles R Greathouse IV_, Mar 19 2015
%Y A256186 is the intersection of this sequence with A000040.
%Y Cf. A011540, A038618, A056709.
%K nonn,base
%O 1,1
%A _Charles R Greathouse IV_ and _Zak Seidov_, Mar 19 2015
|