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
