A253575 Primes p such that digits of p do not appear in p^6.
2, 3, 13, 44449
Offset: 1
Examples
2 and 2^6 = 64 have no digits in common, hence 2 is in the sequence.
Programs
-
Maple
select(t -> isprime(t) and convert(convert(t,base,10),set) intersect convert(convert(t^6,base,10),set) = {}, {2, seq(i,i=3..10^5,2)}); # Robert Israel, May 25 2025
-
Mathematica
Select[Prime[Range[1000000]], Intersection[IntegerDigits[#], IntegerDigits[#^6]]=={} &]
-
Python
from sympy import isprime A253575_list = [n for n in range(1,10**6) if set(str(n)) & set(str(n**6)) == set() and isprime(n)] # Chai Wah Wu, Jan 05 2015
Comments