Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).
%I #12 Jul 14 2012 13:52:00
%S 0,3,4,5,8,10,16,26,29,34,73,113,122,157,178,241,353,457,997,1042,
%T 3041,4562,6434,8506,10141,19378,19882,22426,27529
%N Numbers n such that (1+i)^n - i is a Gaussian prime.
%C 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.
%C 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
%t fQ[n_] := PrimeQ[(1 + I)^n - I, GaussianIntegers -> True]; Select[ Range[0, 30000], fQ]
%Y Cf. A027206 ((1+i)^n + i is a Gaussian prime), A057429 ((1+i)^n - 1 is a Gaussian prime).
%K nonn
%O 1,2
%A _T. D. Noe_, Jan 31 2005
%E a(25)-a(29) from _Robert G. Wilson v_, Mar 02 2011.
%E 0 prepended by _T. D. Noe_, Mar 07 2011