login
Primes p such that not all decimal digits of p occur in p^3.
1

%I #8 Nov 15 2019 00:55:32

%S 2,3,7,13,17,19,23,31,37,41,43,47,53,79,83,89,103,107,113,127,131,139,

%T 149,157,163,167,173,179,181,193,197,199,211,223,227,239,241,257,263,

%U 269,271,277,283,293,307,313,317,349,353,359,367,373,379,383,397,409,419,421,431,433,439,443,457

%N Primes p such that not all decimal digits of p occur in p^3.

%C Complement of A030080 in the primes.

%C Includes 3*10^k-7 if k is in A102964.

%H Robert Israel, <a href="/A329487/b329487.txt">Table of n, a(n) for n = 1..10000</a>

%e a(4)=13 is in the sequence because the digit 3 is in 13 but not in 13^3=2197.

%p filter:= proc(p)

%p if not isprime(p) then return false fi;

%p not (convert(convert(p,base,10),set) subset convert(convert(p^3,base,10), set))

%p end proc:

%p select(filter, [2,seq(i,i=3..10000,2)]);

%Y Cf. A030080, A102964.

%K nonn,base

%O 1,1

%A _Robert Israel_, Nov 14 2019