OFFSET
1,1
COMMENTS
a(n) = n+1 iff n is in A174156.
LINKS
Jens Kruse Andersen, Table of n, a(n) for n = 1..10000
EXAMPLE
34^32 + 35^32 = 35884563485651241417769982593434627193100364196481 is not prime. 34^2 + 36^32 = 2^32*(17^32+18^32) is not prime. 34^32 + 37^32 = 162384303092765940334383766635859112907663593431937 is prime. Thus a(34) = 37.
PROG
(Python)
import sympy
from sympy import isprime
def a(n):
..for k in range(n+1, 10**4):
....if isprime(k**32+n**32):
......return k
n = 1
while n < 100:
..print(a(n), end=', ')
..n += 1
(PARI) a(n)=for(k=n+1, 10^4, if(isprime(k^32+n^32), return(k)))
n=1; while(n<100, print1(a(n), ", "); n++)
CROSSREFS
KEYWORD
nonn
AUTHOR
Derek Orr, Jul 08 2014
STATUS
approved