OFFSET
1,1
COMMENTS
Note that n must be odd. Terms of primitive Pythagorean triples: (n^2, (n^4-1)/2, (n^4+1)/2).
LINKS
Charles R Greathouse IV, Table of n, a(n) for n = 1..10000
EXAMPLE
a(1)=41 because (3^4 + 1)/2 = 82/2 = 41 is prime.
MATHEMATICA
Select[(Range[200]^4+1)/2, PrimeQ] (* Harvey P. Dale, Mar 09 2013 *)
PROG
(Magma) [ a: n in [0..2500] | IsPrime(a) where a is ((n^4+1) div 2) ]; // Vincenzo Librandi, Apr 15 2011
(PARI) list(lim)=my(v=List(), t); forstep(n=3, sqrtnint(lim\1*2-1, 4), 2, if(isprime(t=(n^4+1)/2), listput(v, t))); Vec(v) \\ Charles R Greathouse IV, Feb 14 2017
CROSSREFS
KEYWORD
nonn
AUTHOR
Hugo Pfoertner, Jun 19 2004
EXTENSIONS
Name edited by Zak Seidov, Apr 14 2011
STATUS
approved