OFFSET
1,1
COMMENTS
These are the values of m such that A240081(m) 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
