OFFSET
1,1
COMMENTS
n is odd iff a(n) is even.
LINKS
Lars-Erik Svahn, Table of n, a(n) for n = 1..10000
Lars-Erik Svahn, numbertheory.4th
Akshaa Vatwani, Bounded gaps between Gaussian primes, Journal of Number Theory, Volume 171, February 2017, Pages 449-473.
Eric Weisstein's World of Mathematics, Gaussian Prime.
EXAMPLE
a(1)=2 since (1 + 1)^2 + (1 + 1)^2 is not prime, but 1^2 + 2^2 = 5 and (1 + 1)^2 + (2 + 1)^2 = 13 are prime.
MATHEMATICA
Rest@ FoldList[Module[{k = 1}, While[Times @@ Boole@ Map[PrimeQ, {#2^2 + k^2, (#2 + 1)^2 + (k + 1)^2}] < 1, k++]; k] &, 1, Range@ 76] (* Michael De Vlieger, Mar 25 2017 *)
PROG
(ANS-Forth)
s" numbertheory.4th" included
: Gauss_twin \ n -- a(n)
locals| n | 0
begin 1+ dup dup * n dup * + isprime
over 1+ dup * n 1+ dup * + isprime and
until ;
(PARI) a(n) = my(k=0); while (! (isprime(n^2+k^2) && isprime((n+1)^2+(k+1)^2)), k++); k; \\ Michel Marcus, Mar 25 2017
CROSSREFS
KEYWORD
nonn
AUTHOR
Lars-Erik Svahn, Mar 25 2017
STATUS
approved