OFFSET
1,1
COMMENTS
The terms in A007088 with a prime number of 1's are trivially contained in this sequence.
EXAMPLE
1011 is a member of this sequence because 1^1011 + 0^1011 + 1^1011 + 1^1011 = 3 is prime.
PROG
(Python)
import sympy
from sympy import isprime
def Pow(x):
..num = 0
..for i in str(x):
....num += int(i)**x
..if isprime(num):
....return True
x = 1
while x < 10**5:
..if Pow(x):
....print(x)
..x += 1
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Derek Orr, Mar 13 2014
EXTENSIONS
a(12)-a(37) from Giovanni Resta, Mar 14 2014
STATUS
approved