OFFSET
1,1
COMMENTS
These are the values of n such that A240081(n) is nonzero.
PROG
(Python)
import sympy
from sympy import isprime
def Leq(n):
..for k in range(1000):
....num = k**n
....for i in range(2, k):
......num -= i**n
......if num < 1:
........return None
....if isprime(num-1):
......return k
n = 1
while n < 10**3:
..if Leq(n) != None:
....print(n)
..n += 1
(PARI) f(n)=for(k=1, 10^3, num=k^n; for(i=2, k-1, num-=i^n; if(num<1, return(0))); if(ispseudoprime(num-1), return(k))); n=1; while(n<10^3, if(f(n), print(n)); n+=1)
CROSSREFS
KEYWORD
nonn
AUTHOR
Derek Orr, Mar 31 2014
EXTENSIONS
a(34)-a(52) from Giovanni Resta, Apr 02 2014
STATUS
approved