cp's OEIS Frontend

A057876 Primes p with the following property: let d_1, d_2, ... be the distinct digits occurring in the decimal expansion of p. Then for each d_i, dropping all the digits d_i from p produces a prime number. Leading 0's are not allowed.

%I A057876 #11 Oct 27 2019 00:16:50
%S A057876 23,37,53,73,113,131,137,151,173,179,197,211,311,317,431,617,719,1531,
%T A057876 1831,1997,2113,2131,2237,2273,2297,2311,2797,3137,3371,4337,4373,
%U A057876 4733,4919,6173,7297,7331,7573,7873,8191,8311,8831,8837,12239,16673,19531
%N A057876 Primes p with the following property: let d_1, d_2, ... be the distinct digits occurring in the decimal expansion of p. Then for each d_i, dropping all the digits d_i from p produces a prime number.  Leading 0's are not allowed.
%H A057876 Robert Israel, <a href="/A057876/b057876.txt">Table of n, a(n) for n = 1..1000</a>
%e A057876 1531 gives primes 53, 131 and 151 after dropping digits 1, 5 and 3.
%e A057876 A larger example 1210778071 gives primes 12177871, 2077807, 110778071, 1210801 and 121077071 after dropping digits 0, 1, 2, 7 and 8.
%p A057876 filter:= proc(n) local L,d,Lp;
%p A057876   if not isprime(n) then return false fi;
%p A057876   L:= convert(n,base,10);
%p A057876   for d in convert(L,set) do
%p A057876     Lp:= subs(d=NULL,L);
%p A057876     if Lp=[] or Lp[-1] = 0 then return false fi;
%p A057876     if not isprime(add(Lp[i]*10^(i-1),i=1..nops(Lp))) then return false fi;
%p A057876   od;
%p A057876   true
%p A057876 end proc:
%p A057876 select(filter, [seq(i,i=13..20000,2)]); # _Robert Israel_, Jul 13 2018
%Y A057876 Cf. A057877-A057883, A051362, A034302-A034305.
%K A057876 nonn,base
%O A057876 1,1
%A A057876 _Patrick De Geest_, Oct 15 2000
%E A057876 Name edited by _Robert Israel_, Jul 13 2018