OFFSET
1,2
COMMENTS
Also numbers m such that A239472(m) > 0.
PROG
(Python)
import sympy
from sympy import isprime
def Lep(n):
for k in range(2*10**3):
num = k**n
for i in range(2, k):
num -= i**n
if num < 0:
return None
if isprime(num):
return k
n = 1
while n < 10**3:
if Lep(n) != None:
print(n)
n += 1
(PARI) a(n)=k=1; while((s=k^n-sum(i=2, k-1, i^n))>0, if(isprime(s), return(k)); k++)
for(n=1, 100, if(a(n), print1(n, ", "))) \\ Derek Orr, Mar 12 2015
CROSSREFS
KEYWORD
nonn
AUTHOR
Derek Orr, Mar 31 2014
STATUS
approved
