OFFSET
1,1
COMMENTS
Primes of sequence A029785. - Michel Marcus, Jan 04 2015
LINKS
Michael S. Branicky, Table of n, a(n) for n = 1..86 (terms < 10^19; terms 74..86 via A029785)
EXAMPLE
2 and 2^3=8 have no digits in common, hence 2 is in the sequence.
MATHEMATICA
Select[Prime[Range[1500000]], Intersection[IntegerDigits[#], IntegerDigits[#^3]]=={} &] (* Vincenzo Librandi, Jan 04 2015 *)
PROG
(PARI) lista(nn) = {forprime (n=1, nn, if (#setintersect(Set(vecsort(digits(n^3))), Set(vecsort(digits(n)))) == 0, print1(n, ", ")); ); } \\ Michel Marcus, Jan 04 2015
(Python)
from sympy import isprime
A030087_list = [n for n in range(1, 10**6) if set(str(n)) & set(str(n**3)) == set() and isprime(n)]
# Chai Wah Wu, Jan 05 2015
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Patrick De Geest, Dec 11 1999
EXTENSIONS
Changed offset from 0 to 1 and more terms from Vincenzo Librandi, Jan 04 2015
a(40)-a(47) from Chai Wah Wu, Jan 05 2015
STATUS
approved