A253576 Primes p such that digits of p do not appear in p^7.
3, 7, 43, 7757, 31333
Offset: 1
Examples
3 and 3^7 = 2187 have no digits in common, hence 3 is in the sequence.
Programs
-
Mathematica
Select[Prime[Range[1000000]], Intersection[IntegerDigits[#], IntegerDigits[#^7]]=={} &]
-
Python
from sympy import isprime A253576_list = [n for n in range(1,10**6) if set(str(n)) & set(str(n**7)) == set() and isprime(n)] # Chai Wah Wu, Jan 05 2015
Comments