OFFSET
1,2
COMMENTS
Note that A027206 and A057429 treat Gaussian primes of a similar form. The remaining case, (1+i)^n + 1, is a Gaussian prime for n=1,2,3,4 only.
Let z = (1+i)^n - i. If z is not pure real or pure imaginary, then z is a Gaussian prime if the product of z and its conjugate is a rational prime. That product is 1 + 2^n - sin(n*Pi/4)*2^(1+n/2). z is real when n=1. z is imaginary when n=4k+2, in which case, z has magnitude 2^(2k+1) - (-1)^k. These pure imaginary numbers are Gaussian primes when 2^(2k+1)-1 is a Mersenne prime and 2k+1 = 3 (mod 4); that is, when n is twice an odd number in A112634. - T. D. Noe, Mar 07 2011
MATHEMATICA
fQ[n_] := PrimeQ[(1 + I)^n - I, GaussianIntegers -> True]; Select[ Range[0, 30000], fQ]
CROSSREFS
KEYWORD
nonn
AUTHOR
T. D. Noe, Jan 31 2005
EXTENSIONS
a(25)-a(29) from Robert G. Wilson v, Mar 02 2011.
0 prepended by T. D. Noe, Mar 07 2011
STATUS
approved