cp's OEIS Frontend

A338018 Primes not in A228629.

%I A338018 #10 Oct 07 2020 11:25:04
%S A338018 3,5,11,13,17,19,29,31,37,41,43,47,53,59,71,73,79,89,97,101,103,113,
%T A338018 131,137,139,149,151,157,173,179,191,193,197,199,229,233,257,281,311,
%U A338018 313,317,331,337,347,349,353,359,373,379,397,419,431,433,439,443,449,457,479,491,499,503,541,547,557
%N A338018 Primes not in A228629.
%C A338018 Primes p such that there is no prime whose decimal digits are (in some order) the 9's complements of the decimal digits of p.
%C A338018 Leading zeros are not allowed.
%H A338018 Robert Israel, <a href="/A338018/b338018.txt">Table of n, a(n) for n = 1..10000</a>
%e A338018 a(4)=13 is a member because the 9's complements of the digits 1,3 are 8,2, and none of the integers 28 or 82 is prime.
%e A338018 a(26)=149 is a member because the 9's complements of its digits are 8,5,0, and none of the integers with those digits is prime.
%p A338018 R:= 3,5:
%p A338018 for d from 2 to 4 do
%p A338018   P:= select(isprime,[seq(i,i=10^(d-1)+1..10^d-1,2)]);
%p A338018   nP:= nops(P);
%p A338018   Pd:= map(sort@convert,P,base,10);
%p A338018   Ps:= convert(map(t -> ListTools:-Reverse([9$d]-t), Pd),set);
%p A338018   S:= remove(t -> member(Pd[t],Ps),[$1..nP]);
%p A338018   R:= R, op(P[S]);
%p A338018 od:
%p A338018 R; # _Robert Israel_, Oct 06 2020
%Y A338018 Cf. A228629. Contains all members of A030096 except 7.
%K A338018 nonn,base,look
%O A338018 1,1
%A A338018 _Robert Israel_, Oct 06 2020