OFFSET
1,2
COMMENTS
Also numbers n such that A239472(n) > 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