This site is supported by donations to The OEIS Foundation.

Annual Appeal: Today, Nov 11 2014, is the 4th anniversary of the launch of the new OEIS web site. 70,000 sequences have been added in these four years, all edited by volunteers. Please make a donation (tax deductible in the US) to help keep the OEIS running.

 Hints (Greetings from The On-Line Encyclopedia of Integer Sequences!)
 A073019 Numbers n such that the absolute value of the real part of (1+2*I)^n is prime. 0
 2, 3, 4, 5, 7, 11, 13, 17, 19, 89, 163, 167, 233, 251, 349, 401, 757, 1423, 3049, 3643, 11987, 19801, 26393, 32303, 35363, 35863, 36607, 36683 (list; graph; refs; listen; history; text; internal format)
 OFFSET 1,1 COMMENTS In the first 28 terms only n = 4 is composite; does this continue? - Charles R Greathouse IV, Nov 25 2014 LINKS FORMULA {n: |A006495(n)| in A000040}. MAPLE A006495 := proc(n) add((-4)^l*binomial(n, 2*l), l=0..floor(n/2)) ; end proc: isA073019 := proc(n) isprime(abs(A006495(n))) ; end proc: for n from 1 do if isA073019(n) then print(n) ; end if; end do: # R. J. Mathar, Sep 06 2011 MATHEMATICA x=1; For[n=1, n<=14000, n++, {x=(1+2*I)*x, If[PrimeQ[Re[x]], Print[n]]}] PROG (PARI) is(n)=ispseudoprime(abs(real((1+2*I)^n))) \\ Charles R Greathouse IV, Nov 21 2014 CROSSREFS Cf. A006495. Sequence in context: A211781 A046022 A175787 * A174291 A007885 A192586 Adjacent sequences:  A073016 A073017 A073018 * A073020 A073021 A073022 KEYWORD nonn AUTHOR Benoit Cloitre, Aug 03 2002 EXTENSIONS Edited by John W. Layman, Aug 13 2002 a(21) from Robert G. Wilson v, Aug 22 2002 a(22)-a(28) from Charles R Greathouse IV, Nov 25 2014 STATUS approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Transforms | Superseeker | Recent | More pages
The OEIS Community | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement .

Last modified December 21 10:09 EST 2014. Contains 252304 sequences.