login
This site is supported by donations to The OEIS Foundation.
Logo

Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A027697 Odious primes: primes with odd number of 1's in binary expansion. 21
2, 7, 11, 13, 19, 31, 37, 41, 47, 59, 61, 67, 73, 79, 97, 103, 107, 109, 127, 131, 137, 151, 157, 167, 173, 179, 181, 191, 193, 199, 211, 223, 227, 229, 233, 239, 241, 251, 271, 283, 307, 313, 331, 367, 379, 397, 409, 419, 421, 431, 433, 439, 443, 457, 463, 487, 491, 499, 521, 541, 557, 563 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,1

COMMENTS

Comment from Vladimir Shevelev: Conjecture: a(n) < A027699(n) except for n = 2. I verified this conjecture up to n=5*10^7. Moreover, I conjecture that A027699(n) - a(n) tends to infinity.

Item 2 of Green is about this sequence: "This is an expository article to accompany my two lectures at the CDM conference. I have used this an excuse to make public two sets of notes I had lying around and also to put together a short reader's guide to some recent joint work with T. Tao. Contents: 1. An exposition, without much detail, of the work of Goldston, Pintz and Yildirim on gaps between primes; 2. A detailed discussion of the work of Mauduit and Rivat establishing that 50 percent of the primes have odd digit sum when written in base 2; 3. A reader's guide to recent work of T. Tao and the author on linear equations in primes. The sections can be read independently." - Jonathan Vos Post, Oct 04 2007

LINKS

T. D. Noe, Table of n, a(n) for n=1..10000

Ben Green, Three topics in additive prime number theory, Oct 03, 2007, pp. 12-27.

MAPLE

a:=proc(n) local nn: nn:= convert(ithprime(n), base, 2): if `mod`(sum(nn[j], j =1..nops(nn)), 2)=1 then ithprime(n) else end if end proc: seq(a(n), n=1..103); - Emeric Deutsch, Oct 24 2007

MATHEMATICA

Clear[BinSumOddQ]; BinSumOddQ[a_]:=Module[{i, s=0}, s=0; For[i=1, i<=Length[IntegerDigits[a, 2]], s+=Extract[IntegerDigits[a, 2], i]; i++ ]; OddQ[s]]; lst={}; Do[p=Prime[n]; If[BinSumOddQ[p], AppendTo[lst, p]], {n, 4!}]; lst [From Vladimir Joseph Stephan Orlovsky, Apr 06 2009]

PROG

(PARI) f(p)={v=binary(p); s=0; for(k=1, #v, if(v[k]==1, s++)); return(s%2)};

forprime(p=2, 563, if(f(p), print1(p, ", "))) - [W. Bomfim webonfim(AT)bol.com Jan 14, 2011]

CROSSREFS

Cf. A027699, A066148, A066149.

Cf. A000069 (odious numbers), A092246 (odd odious numbers)

Sequence in context: A161681 A020583 A140557 * A146315 A038892 A019346

Adjacent sequences:  A027694 A027695 A027696 * A027698 A027699 A027700

KEYWORD

nonn,easy,base

AUTHOR

N. J. A. Sloane.

EXTENSIONS

More terms from Scott Lindhurst (ScottL(AT)alumni.princeton.edu)

STATUS

approved

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

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

Last modified May 24 12:34 EDT 2013. Contains 225620 sequences.