OFFSET
1,2
COMMENTS
Note that A027206 and A057429 treat Gaussian primes of a similar form. The remaining case, (1+i)^m + 1, is a Gaussian prime for m=1,2,3,4 only.
Let z = (1+i)^m - 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^m - sin(m*Pi/4)*2^(1+m/2). z is real when m=1. z is imaginary when m=4*k+2, in which case, z has magnitude 2^(2*k+1) - (-1)^k. These pure imaginary numbers are Gaussian primes when 2^(2*k+1)-1 is a Mersenne prime and 2*k+1 == 3 (mod 4); that is, when m 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,more
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
