OFFSET
1,2
COMMENTS
If a(n) = 1, then n is in A056994.
MATHEMATICA
lnk[n_]:=Module[{c=n^128, k}, k=If[EvenQ[c], 1, 2]; While[!PrimeQ[c+ k^128], k = k+2]; k]; Join[{1}, Array[lnk, 60, 2]] (* Harvey P. Dale, Mar 17 2015 *)
PROG
(Python)
import sympy
from sympy import isprime
def a(n):
for k in range(10**4):
if isprime(n**128+k**128):
return k
n = 1
while n < 100:
print(a(n))
n += 1
(PARI) a(n)=for(k=1, 10^4, if(ispseudoprime(n^128+k^128), return(k)));
n=1; while(n<100, print(a(n)); n+=1)
CROSSREFS
KEYWORD
nonn
AUTHOR
Derek Orr, May 17 2014
STATUS
approved