cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A338018 Primes not in A228629.

Original entry on oeis.org

3, 5, 11, 13, 17, 19, 29, 31, 37, 41, 43, 47, 53, 59, 71, 73, 79, 89, 97, 101, 103, 113, 131, 137, 139, 149, 151, 157, 173, 179, 191, 193, 197, 199, 229, 233, 257, 281, 311, 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
Offset: 1

Views

Author

Robert Israel, Oct 06 2020

Keywords

Comments

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.
Leading zeros are not allowed.

Examples

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

Crossrefs

Cf. A228629. Contains all members of A030096 except 7.

Programs

  • Maple
    R:= 3,5:
    for d from 2 to 4 do
      P:= select(isprime,[seq(i,i=10^(d-1)+1..10^d-1,2)]);
      nP:= nops(P);
      Pd:= map(sort@convert,P,base,10);
      Ps:= convert(map(t -> ListTools:-Reverse([9$d]-t), Pd),set);
      S:= remove(t -> member(Pd[t],Ps),[$1..nP]);
      R:= R, op(P[S]);
    od:
    R; # Robert Israel, Oct 06 2020