%I #12 Jan 08 2015 17:36:30
%S 3,7,43,7757,31333
%N Primes p such that digits of p do not appear in p^7.
%C a(6)> 10^7.
%C a(6) > 10^9. - _Chai Wah Wu_, Jan 05 2015
%e 3 and 3^7 = 2187 have no digits in common, hence 3 is in the sequence.
%t Select[Prime[Range[1000000]], Intersection[IntegerDigits[#], IntegerDigits[#^7]]=={} &]
%o (Python)
%o from sympy import isprime
%o A253576_list = [n for n in range(1,10**6) if set(str(n)) & set(str(n**7)) == set() and isprime(n)]
%o # _Chai Wah Wu_, Jan 05 2015
%Y Cf. similar sequences listed in A253574.
%K nonn,base,more
%O 1,1
%A _Vincenzo Librandi_, Jan 04 2015
|