<?xml version="1.0" encoding="UTF-8"?>
<rss version="2.0"
	xmlns:content="http://purl.org/rss/1.0/modules/content/"
	xmlns:wfw="http://wellformedweb.org/CommentAPI/"
	xmlns:dc="http://purl.org/dc/elements/1.1/"
	xmlns:atom="http://www.w3.org/2005/Atom"
	xmlns:sy="http://purl.org/rss/1.0/modules/syndication/"
	xmlns:slash="http://purl.org/rss/1.0/modules/slash/"
	>

<channel>
	<title>Michael Nielsen &#187; Polymath</title>
	<atom:link href="http://michaelnielsen.org/blog/category/polymath/feed/" rel="self" type="application/rss+xml" />
	<link>http://michaelnielsen.org/blog</link>
	<description></description>
	<lastBuildDate>Fri, 10 Feb 2012 19:10:31 +0000</lastBuildDate>
	<language>en</language>
	<sy:updatePeriod>hourly</sy:updatePeriod>
	<sy:updateFrequency>1</sy:updateFrequency>
	<generator>http://wordpress.org/?v=3.3</generator>
		<item>
		<title>Finding Primes: A Fun Subproblem</title>
		<link>http://michaelnielsen.org/blog/finding-primes-a-fun-subproblem/</link>
		<comments>http://michaelnielsen.org/blog/finding-primes-a-fun-subproblem/#comments</comments>
		<pubDate>Tue, 01 Sep 2009 21:02:15 +0000</pubDate>
		<dc:creator>Michael Nielsen</dc:creator>
				<category><![CDATA[Polymath]]></category>
		<category><![CDATA[Uncategorized]]></category>
		<category><![CDATA[polymath1]]></category>

		<guid isPermaLink="false">http://michaelnielsen.org/blog/?p=664</guid>
		<description><![CDATA[The ongoing open mathematics project finding primes aims to find a deterministic algorithm to efficiently generate -digit primes. The fastest known algorithm seems to be a method of Odlyzko which generates a -digit prime in time . The people working on the project have made some observations which come tantalizingly close to breaking that barrier. [...]]]></description>
			<content:encoded><![CDATA[<p>The ongoing open mathematics project <a href="http://michaelnielsen.org/polymath1/index.php?title=Finding_primes">finding   primes</a> aims to find a deterministic algorithm to efficiently generate <img src='/blog/latexrender/pictures/8ce4b16b22b58894aa86c421e8759df3.png' title='k' alt='k' align=absmiddle>-digit primes.  The fastest known algorithm seems to be a <a href="http://michaelnielsen.org/polymath1/index.php?title=Odlyzko's_method">method   of Odlyzko</a> which generates a <img src='/blog/latexrender/pictures/8ce4b16b22b58894aa86c421e8759df3.png' title='k' alt='k' align=absmiddle>-digit prime in time <img src='/blog/latexrender/pictures/fad1f0b21b11b544eccd1fb41b5ec4ee.png' title='O(10^{k/2})' alt='O(10^{k/2})' align=absmiddle>. The people working on the project have made some observations which come tantalizingly close to breaking that barrier.  The obstruction is a beautiful little problem that I thought many people might enjoy, and which may well be tractable.  If you&#8217;re interested in participating in the finding primes project, it might be a good entry point.  So if you&#8217;ve got a good idea to solve the problem, pitch in and help over at the <a href="http://polymathprojects.org">Polymath blog</a>.  But please be polite: read some <a href="http://polymathprojects.org/about/">background</a> <a href="http://michaelnielsen.org/polymath1/index.php?title=Finding_primes">first</a>, and take a look at some of the <a href="http://polymathprojects.org/2009/08/28/research-thread-iv-determinstic-way-to-find-primes/">research   threads</a> to get a feel for how things work, and what&#8217;s already known.</p>
<p>One small caveat: the argument that follows is not in any way mine, it&#8217;s all the work of other people!  A recent comment thread on this argument starts <a href="http://polymathprojects.org/2009/08/28/research-thread-iv-determinstic-way-to-find-primes/#comment-715">here</a>.</p>
<p>Let <img src='/blog/latexrender/pictures/e5cb16c20d9f01bbbfe8f299e28d1f4b.png' title='\pi(x)' alt='\pi(x)' align=absmiddle> denote the number of primes less than or equal to <img src='/blog/latexrender/pictures/9dd4e461268c8034f5c8564e155c67a6.png' title='x' alt='x' align=absmiddle>. <img src='/blog/latexrender/pictures/e5cb16c20d9f01bbbfe8f299e28d1f4b.png' title='\pi(x)' alt='\pi(x)' align=absmiddle> has a lot of structure, and there&#8217;s a surprising amount that can be said about it.  In particular, the people working on <a href="http://michaelnielsen.org/polymath1/index.php?title=Finding_primes">finding   primes</a> have figured out a clever way of computing the parity of <img src='/blog/latexrender/pictures/e5cb16c20d9f01bbbfe8f299e28d1f4b.png' title='\pi(x)' alt='\pi(x)' align=absmiddle> in time about <a href="http://polymathprojects.org/2009/08/28/research-thread-iv-determinstic-way-to-find-primes/#comment-729"><img src='/blog/latexrender/pictures/3cb473b589160a4f3fbcc93e185481a3.png' title='x^{5/11+o(1)}' alt='x^{5/11+o(1)}' align=absmiddle></a>.</p>
<p>Suppose you can find two <img src='/blog/latexrender/pictures/8ce4b16b22b58894aa86c421e8759df3.png' title='k' alt='k' align=absmiddle>-digit numbers <img src='/blog/latexrender/pictures/9dd4e461268c8034f5c8564e155c67a6.png' title='x' alt='x' align=absmiddle> and <img src='/blog/latexrender/pictures/415290769594460e2e485922904f345d.png' title='y' alt='y' align=absmiddle> such that the parity of <img src='/blog/latexrender/pictures/e5cb16c20d9f01bbbfe8f299e28d1f4b.png' title='\pi(x)' alt='\pi(x)' align=absmiddle> and <img src='/blog/latexrender/pictures/2b2d28ff25cab00621660d2de85d9ce8.png' title='\pi(y)' alt='\pi(y)' align=absmiddle> are different.  Set <img src='/blog/latexrender/pictures/61735b4fe89b403d2936b6a32648ad14.png' title='z = \lfloor (x+y)/2 \rfloor' alt='z = \lfloor (x+y)/2 \rfloor' align=absmiddle>, i.e., take <img src='/blog/latexrender/pictures/fbade9e36a3f36d3d676c1b808451dd7.png' title='z' alt='z' align=absmiddle> to be the midpoint between <img src='/blog/latexrender/pictures/9dd4e461268c8034f5c8564e155c67a6.png' title='x' alt='x' align=absmiddle> and <img src='/blog/latexrender/pictures/415290769594460e2e485922904f345d.png' title='y' alt='y' align=absmiddle>.  Compute the parity of <img src='/blog/latexrender/pictures/fbade9e36a3f36d3d676c1b808451dd7.png' title='z' alt='z' align=absmiddle>.  It must have either a different parity to <img src='/blog/latexrender/pictures/9dd4e461268c8034f5c8564e155c67a6.png' title='x' alt='x' align=absmiddle> or a different parity to <img src='/blog/latexrender/pictures/415290769594460e2e485922904f345d.png' title='y' alt='y' align=absmiddle>.  Repeating this procedure <img src='/blog/latexrender/pictures/18bfba88133de0dabea899170b5e687b.png' title='O(k)' alt='O(k)' align=absmiddle> times, we can use a binary search to find adjacent <img src='/blog/latexrender/pictures/8ce4b16b22b58894aa86c421e8759df3.png' title='k' alt='k' align=absmiddle>-digit numbers <img src='/blog/latexrender/pictures/70ff83728fb2701cbc52fac8df762bf1.png' title='p-1' alt='p-1' align=absmiddle> and <img src='/blog/latexrender/pictures/83878c91171338902e0fe0fb97a8c47a.png' title='p' alt='p' align=absmiddle> such that <img src='/blog/latexrender/pictures/eafecf0d375175caa1d62d26effd515c.png' title='\pi(p-1)' alt='\pi(p-1)' align=absmiddle> and <img src='/blog/latexrender/pictures/2474ab5473c6609c4ead44d551c965fa.png' title='\pi(p)' alt='\pi(p)' align=absmiddle> have different parity.  We conclude that <img src='/blog/latexrender/pictures/83878c91171338902e0fe0fb97a8c47a.png' title='p' alt='p' align=absmiddle> must be prime.  That takes time <img src='/blog/latexrender/pictures/9df4e0e237c9a6b2f98d5d1247bfc953.png' title='O(10^{5/11 k + o(1)})' alt='O(10^{5/11 k + o(1)})' align=absmiddle>, and so breaks the barrier set by Odlyzko&#8217;s method.</p>
<p>What&#8217;s the problem with this algorithm?  The problem is finding the initial <img src='/blog/latexrender/pictures/8ce4b16b22b58894aa86c421e8759df3.png' title='k' alt='k' align=absmiddle>-digit numbers <img src='/blog/latexrender/pictures/9dd4e461268c8034f5c8564e155c67a6.png' title='x' alt='x' align=absmiddle> and <img src='/blog/latexrender/pictures/415290769594460e2e485922904f345d.png' title='y' alt='y' align=absmiddle> such that <img src='/blog/latexrender/pictures/e5cb16c20d9f01bbbfe8f299e28d1f4b.png' title='\pi(x)' alt='\pi(x)' align=absmiddle> and <img src='/blog/latexrender/pictures/2b2d28ff25cab00621660d2de85d9ce8.png' title='\pi(y)' alt='\pi(y)' align=absmiddle> have different parity.  It would surprise me a great deal if this weren&#8217;t possible, but it&#8217;s not obvious (at least to me) how to do it quickly.  Is there a fast way of doing this?</p>
]]></content:encoded>
			<wfw:commentRss>http://michaelnielsen.org/blog/finding-primes-a-fun-subproblem/feed/</wfw:commentRss>
		<slash:comments>0</slash:comments>
		</item>
		<item>
		<title>Bertrand&#8217;s postulate</title>
		<link>http://michaelnielsen.org/blog/bertrands-postulate/</link>
		<comments>http://michaelnielsen.org/blog/bertrands-postulate/#comments</comments>
		<pubDate>Tue, 18 Aug 2009 13:53:23 +0000</pubDate>
		<dc:creator>Michael Nielsen</dc:creator>
				<category><![CDATA[Polymath]]></category>

		<guid isPermaLink="false">http://michaelnielsen.org/blog/?p=655</guid>
		<description><![CDATA[This post is motivated by the ongoing Polymath Project &#8220;finding primes&#8221;, an open source mathematics project whose aim (roughly speaking) is to find an efficient algorithm to deterministically generate a prime number of at least digits. The post isn&#8217;t a contribution to the ongoing research discussion, but rather describes background material I wanted to understand [...]]]></description>
			<content:encoded><![CDATA[<p>This post is motivated by the ongoing Polymath Project <a href="http://michaelnielsen.org/polymath1/index.php?title=Finding_primes">&#8220;finding   primes&#8221;</a>, an open source mathematics project whose aim (roughly speaking) is to find an efficient algorithm to deterministically generate a prime number of at least <img src='/blog/latexrender/pictures/8ce4b16b22b58894aa86c421e8759df3.png' title='k' alt='k' align=absmiddle> digits. The post isn&#8217;t a contribution to the ongoing <a href="http://polymathprojects.org/2009/08/13/research-thread-iii-determinstic-way-to-find-primes/">research   discussion</a>, but rather describes background material I wanted to understand better: a proof of a mathematical theorem known as Bertrand&#8217;s Postulate, a simple result from the 19th century about how frequently prime numbers occur.  It&#8217;s being cross-posted to the <a href="http://michaelnielsen.org/polymath1/index.php?title=Bertrand's_postulate">Polymath   wiki</a>.  The post requires some familiarity with elementary number theory.</p>
<p>Despite its name, Bertrand&#8217;s postulate is a theorem, not an example of what we&#8217;d ordinarily call a postulate in mathematics.  Here&#8217;s what it says:</p>
<p><strong>Theorem (Bertrand&#8217;s Postulate):</strong> <em>For every integer <img src='/blog/latexrender/pictures/d9111381cc3bce0d928f6865909ba5f0.png' title='n   \geq 2' alt='n   \geq 2' align=absmiddle>, there is a prime <img src='/blog/latexrender/pictures/83878c91171338902e0fe0fb97a8c47a.png' title='p' alt='p' align=absmiddle> satisfying <img src='/blog/latexrender/pictures/635729275f58d8d3d6003589d105f73c.png' title='n &lt; p &lt; 2n' alt='n &lt; p &lt; 2n' align=absmiddle>.</em></p>
<p>Bertrand&#8217;s postulate is relevant to the <a href="http://michaelnielsen.org/polymath1/index.php?title=Finding_primes">finding   primes</a> problem because it guarantees the existence of a <img src='/blog/latexrender/pictures/8ce4b16b22b58894aa86c421e8759df3.png' title='k' alt='k' align=absmiddle>-digit prime for any <img src='/blog/latexrender/pictures/8ce4b16b22b58894aa86c421e8759df3.png' title='k' alt='k' align=absmiddle>. Brute force search thus yields a <img src='/blog/latexrender/pictures/8ce4b16b22b58894aa86c421e8759df3.png' title='k' alt='k' align=absmiddle>-digit prime after about <img src='/blog/latexrender/pictures/8b9ac9b3a2caf53c8d97b9ab854976f5.png' title='O(10^k)' alt='O(10^k)' align=absmiddle> steps; this can be considered the &#8220;trivial bound&#8221; for the problem.</p>
<p>Bertrand&#8217;s postulate was apparently first proved by Chebyshev. For large <img src='/blog/latexrender/pictures/7b8b965ad4bca0e41ab51de7b31363a1.png' title='n' alt='n' align=absmiddle>, the claim follows as a consequence of the prime number theorem. We will give an elementary proof due to Erdos.</p>
<p>Our strategy is to analyse the prime factorization of the binomial coefficient <img src='/blog/latexrender/pictures/c92da943df73dc077dbee5514376346a.png' title='{2n \choose n}' alt='{2n \choose n}' align=absmiddle>.  What we&#8217;ll discover is that for primes <img src='/blog/latexrender/pictures/cdb071c4c9a56c35818230afe995b24e.png' title='p \leq n' alt='p \leq n' align=absmiddle>, the corresponding prime power in the prime factorization of <img src='/blog/latexrender/pictures/c92da943df73dc077dbee5514376346a.png' title='{2n \choose n}' alt='{2n \choose n}' align=absmiddle> is never very big.  In fact, we can guarantee that those powers are quite small &#8211; much smaller than what one might <em>a priori</em> believe could be the case.  What we&#8217;ll find as a consequence is that when we take the product of all the prime powers for <img src='/blog/latexrender/pictures/cdb071c4c9a56c35818230afe995b24e.png' title='p \leq n' alt='p \leq n' align=absmiddle>, the product can&#8217;t possibly be as big as <img src='/blog/latexrender/pictures/cbbfd7725796feb5ee07d20250bee011.png' title='{2n \choose   n}' alt='{2n \choose   n}' align=absmiddle>.  And that means that there must be primes in the range <img src='/blog/latexrender/pictures/635729275f58d8d3d6003589d105f73c.png' title='n &lt; p &lt; 2n' alt='n &lt; p &lt; 2n' align=absmiddle> which appear in the prime factorization of <img src='/blog/latexrender/pictures/c92da943df73dc077dbee5514376346a.png' title='{2n \choose n}' alt='{2n \choose n}' align=absmiddle>.</p>
<p>Rather than prove Bertrand's postulate directly, we'll first prove two results of independent interest.  The first is a bound on the product of all primes less than <img src='/blog/latexrender/pictures/7b8b965ad4bca0e41ab51de7b31363a1.png' title='n' alt='n' align=absmiddle>, due to Chebyshev.</p>
<p><strong>Lemma (Chebyshev bound):</strong> <em>For integers <img src='/blog/latexrender/pictures/96bee3237a70bf8577da50732cca3b6b.png' title='n \geq 2' alt='n \geq 2' align=absmiddle>,     <img src='/blog/latexrender/pictures/d0762fbf955ae375d425bd50f3887c30.png' title='\prod_{p \leq n} p \leq 4^n' alt='\prod_{p \leq n} p \leq 4^n' align=absmiddle>, where the product is over all     primes <img src='/blog/latexrender/pictures/83878c91171338902e0fe0fb97a8c47a.png' title='p' alt='p' align=absmiddle> less than or equal to <img src='/blog/latexrender/pictures/7b8b965ad4bca0e41ab51de7b31363a1.png' title='n' alt='n' align=absmiddle>.</em></p>
<p>The Chebyshev bound tells us that primes can&#8217;t occur too frequently (for example, with constant density) in the positive integers, for if they did, the product of primes on the left would rise much faster than <img src='/blog/latexrender/pictures/0e727e07f1a22829d8c383d1dbf6cdba.png' title='4^n' alt='4^n' align=absmiddle>.  Although it&#8217;s not needed for a proof of Bertrand&#8217;s postulate, insight can be gained by expressing this idea in a simple quantitative way.  Observe that if <img src='/blog/latexrender/pictures/bb11019dbf0bb779bd6788eb81561b14.png' title='\pi(n)' alt='\pi(n)' align=absmiddle> is the number of primes less than or equal to <img src='/blog/latexrender/pictures/7b8b965ad4bca0e41ab51de7b31363a1.png' title='n' alt='n' align=absmiddle>, then <img src='/blog/latexrender/pictures/77d1a8d811d593dfa64a1eea4b41e14c.png' title='\pi(n)! \leq \prod_{p \leq n} p' alt='\pi(n)! \leq \prod_{p \leq n} p' align=absmiddle>, simply because the <img src='/blog/latexrender/pictures/8ce4b16b22b58894aa86c421e8759df3.png' title='k' alt='k' align=absmiddle>th prime must be at least <img src='/blog/latexrender/pictures/8ce4b16b22b58894aa86c421e8759df3.png' title='k' alt='k' align=absmiddle>.  Combining this observation with Chebyshev&#8217;s bound gives <img src='/blog/latexrender/pictures/7950241cd048a8311d2bacab41f3cdb6.png' title='\pi(n)! \leq 4^n' alt='\pi(n)! \leq 4^n' align=absmiddle>. Taking logarithms of both sides, and applying <a href="http://en.wikipedia.org/wiki/Stirling's_approximation">Stirling&#8217;s   approximation</a>, we see that up to constant factors and small corrections, we have <img src='/blog/latexrender/pictures/b5ce8b532308bd4d2f0b1164c4b3420d.png' title='\pi(n) \ln \pi(n) \leq n' alt='\pi(n) \ln \pi(n) \leq n' align=absmiddle>.  This, in turn, implies that, up to constant factors and small corrections, <img src='/blog/latexrender/pictures/094dd9c3c4f2b7aeaafaf8b4c5e12eef.png' title='\pi(n) \leq n \ln \ln n / \ln n' alt='\pi(n) \leq n \ln \ln n / \ln n' align=absmiddle>.  Roughly speaking, the prime numbers are at most logarithmically dense in the positive integers.  In fact, the <a href="http://en.wikipedia.org/wiki/Prime_number_theorem">prime   number theorem</a> tells us that <img src='/blog/latexrender/pictures/d1680572e0df048dde82f5321b37d120.png' title='\pi(n) \approx n / \ln n' alt='\pi(n) \approx n / \ln n' align=absmiddle>, as <img src='/blog/latexrender/pictures/7b8b965ad4bca0e41ab51de7b31363a1.png' title='n' alt='n' align=absmiddle> becomes large.</p>
<p><strong>Proof of Chebyshev&#8217;s bound:</strong> We will induct on <img src='/blog/latexrender/pictures/7b8b965ad4bca0e41ab51de7b31363a1.png' title='n' alt='n' align=absmiddle>.  The case when <img src='/blog/latexrender/pictures/c303081f7a16f603112b0375bdc84883.png' title='n=2' alt='n=2' align=absmiddle> is obviously true.  Furthermore, if we assume the result when <img src='/blog/latexrender/pictures/7b8b965ad4bca0e41ab51de7b31363a1.png' title='n' alt='n' align=absmiddle> is odd, it follows immediately for <img src='/blog/latexrender/pictures/40b85027598d87611b1c8d5d11e46812.png' title='n+1' alt='n+1' align=absmiddle>, since the left-hand side is not affected by incrementing <img src='/blog/latexrender/pictures/7b8b965ad4bca0e41ab51de7b31363a1.png' title='n' alt='n' align=absmiddle>.  So we assume the result is true up to some even number, <img src='/blog/latexrender/pictures/b0c015bc7ddc9904ba6458a6dd623bbd.png' title='n = 2m' alt='n = 2m' align=absmiddle>, and try to prove the result for <img src='/blog/latexrender/pictures/be9e7408448baf2fd39d415d2da4ea10.png' title='2m+1' alt='2m+1' align=absmiddle>.  What we&#8217;ll aim to show is that:</p>
<p><img src='/blog/latexrender/pictures/b38a3ae161fd75c7f1feaa4257ca92e9.png' title='   (*) \prod_{m+1 &lt; p \leq 2m+1} p \leq 4^m ' alt='   (*) \prod_{m+1 &lt; p \leq 2m+1} p \leq 4^m ' align=absmiddle></p>
<p>If we can prove this, then multiplying by the inequality <img src='/blog/latexrender/pictures/1a524fd30e52fc586de249eba5352f97.png' title='\prod_{p   \leq m+1} p \leq 4^{m+1}' alt='\prod_{p   \leq m+1} p \leq 4^{m+1}' align=absmiddle> (which is true by the inductive hypothesis) will give us the desired result.  To prove (*), note that every prime <img src='/blog/latexrender/pictures/83878c91171338902e0fe0fb97a8c47a.png' title='p' alt='p' align=absmiddle> satisfying <img src='/blog/latexrender/pictures/0b52fe890cb4b66f5cf4b4b12d746540.png' title='m+1 &lt; p \leq 2m+1' alt='m+1 &lt; p \leq 2m+1' align=absmiddle> must divide <img src='/blog/latexrender/pictures/1173127beede094ae401fe4914ff094e.png' title='{2m+1   \choose m}' alt='{2m+1   \choose m}' align=absmiddle>, since <img src='/blog/latexrender/pictures/83878c91171338902e0fe0fb97a8c47a.png' title='p' alt='p' align=absmiddle> appears in the numerator, but not the denominator.  Furthermore, the uniqueness of prime factorization means that the product of all such primes must also divide <img src='/blog/latexrender/pictures/3d4d61c64ab40d6b2633237658d13dad.png' title='{2m+1 \choose   m}' alt='{2m+1 \choose   m}' align=absmiddle>, and so we have</p>
<p><img src='/blog/latexrender/pictures/e488ebcb12a63ac84eac865116e152a0.png' title='   \prod_{m+1 &lt; p \leq 2m+1} p \leq {2m+1 \choose m}. ' alt='   \prod_{m+1 &lt; p \leq 2m+1} p \leq {2m+1 \choose m}. ' align=absmiddle></p>
<p>The proof is completed by noting that <img src='/blog/latexrender/pictures/3dc28bd4d0d609b2741b5091d07f7ceb.png' title='{2m+1 \choose m}' alt='{2m+1 \choose m}' align=absmiddle> appears twice in the binomial expansion of <img src='/blog/latexrender/pictures/506627755f6629f96d89c56f4b3b40e7.png' title='(1+1)^{2m+1}' alt='(1+1)^{2m+1}' align=absmiddle>, and thus <img src='/blog/latexrender/pictures/446dbe0e76fdd26764355f2d7dfb243a.png' title='{2m+1 \choose   m} \leq 4^m' alt='{2m+1 \choose   m} \leq 4^m' align=absmiddle>.  <strong>QED</strong></p>
<p>There are two main ideas in our proof of Chebyshev&#8217;s bound, and it&#8217;s worth pausing to appreciate them.  The first idea is that there should be a relationship between products of primes and factorials, e.g.,</p>
<p><img src='/blog/latexrender/pictures/e488ebcb12a63ac84eac865116e152a0.png' title='   \prod_{m+1 &lt; p \leq 2m+1} p \leq {2m+1 \choose m}. ' alt='   \prod_{m+1 &lt; p \leq 2m+1} p \leq {2m+1 \choose m}. ' align=absmiddle></p>
<p>The second is that there should be a relationship between combinatorial coefficients like <img src='/blog/latexrender/pictures/3dc28bd4d0d609b2741b5091d07f7ceb.png' title='{2m+1 \choose m}' alt='{2m+1 \choose m}' align=absmiddle> and exponentials with fixed bases, e.g., <img src='/blog/latexrender/pictures/fd5322da7e17f8c7b04f38066965784b.png' title='{2m+1 \choose m} \leq 4^m' alt='{2m+1 \choose m} \leq 4^m' align=absmiddle>.  Once these two ideas are firmly in mind, the proof of the Chebyshev bound is obvious.</p>
<p>The second lemma used in our proof of Bertrand's postulate is a useful fact about the prime factorization of <img src='/blog/latexrender/pictures/c92da943df73dc077dbee5514376346a.png' title='{2n \choose n}' alt='{2n \choose n}' align=absmiddle>, telling us that <img src='/blog/latexrender/pictures/c92da943df73dc077dbee5514376346a.png' title='{2n \choose n}' alt='{2n \choose n}' align=absmiddle> has a rather special prime factorization which only involves low powers.  After the main statement of the lemma we will single out two special cases where the proof can be used to make statements that are a little stronger than the main statement.  These cases are included explicitly only because they are useful later. Note that although they complicate the statement of the lemma somewhat, they are trivial consequences of the proof of the main statement in the lemma.</p>
<p><strong>Lemma:</strong> <em>Let <img src='/blog/latexrender/pictures/83878c91171338902e0fe0fb97a8c47a.png' title='p' alt='p' align=absmiddle> be a prime, and define <img src='/blog/latexrender/pictures/dfafc1058f6edabb7f41b6e0b156dbe8.png' title='r(p,n)' alt='r(p,n)' align=absmiddle> to be   the largest number such that <img src='/blog/latexrender/pictures/88377c4d9984ba0f820c6fd2c6a92dab.png' title='p^{r(p,n)}' alt='p^{r(p,n)}' align=absmiddle> divides <img src='/blog/latexrender/pictures/c92da943df73dc077dbee5514376346a.png' title='{2n \choose n}' alt='{2n \choose n}' align=absmiddle>.   Then <img src='/blog/latexrender/pictures/a9fbed768a68122fc80e99b8917c3f62.png' title='p^{r(p,n)} \leq 2n' alt='p^{r(p,n)} \leq 2n' align=absmiddle>.  There are two special cases where we can   make stronger statements: (1) if <img src='/blog/latexrender/pictures/b1393237013f9d5bb3041ae6d70fc2a9.png' title='p &gt; \sqrt{2n}' alt='p &gt; \sqrt{2n}' align=absmiddle>, then <img src='/blog/latexrender/pictures/dfafc1058f6edabb7f41b6e0b156dbe8.png' title='r(p,n)' alt='r(p,n)' align=absmiddle> is   <img src='/blog/latexrender/pictures/c4ca4238a0b923820dcc509a6f75849b.png' title='1' alt='1' align=absmiddle> or <img src='/blog/latexrender/pictures/cfcd208495d565ef66e7dff9f98764da.png' title='0' alt='0' align=absmiddle>; (2) if <img src='/blog/latexrender/pictures/f2c0ed823881eadfa1347431940736f9.png' title='2n/3 &lt; p \leq n' alt='2n/3 &lt; p \leq n' align=absmiddle>, then <img src='/blog/latexrender/pictures/5481cc4cea0948af9b690474ea594be2.png' title='r(p,n) = 0' alt='r(p,n) = 0' align=absmiddle>.</em></p>
<p>This lemma contains the essential insight underlying Bertrand&#8217;s postulate, and that we explained earlier: because we can guarantee that the prime powers are low, <img src='/blog/latexrender/pictures/c92da943df73dc077dbee5514376346a.png' title='{2n \choose n}' alt='{2n \choose n}' align=absmiddle> must have many prime factors.  A careful accounting using the Chebyshev bound will show that we simply don&#8217;t have enough if one of those prime factors isn&#8217;t in the range <img src='/blog/latexrender/pictures/635729275f58d8d3d6003589d105f73c.png' title='n &lt; p &lt; 2n' alt='n &lt; p &lt; 2n' align=absmiddle>.</p>
<p><strong>Proof:</strong> To prove the lemma, first note that the prime power of <img src='/blog/latexrender/pictures/83878c91171338902e0fe0fb97a8c47a.png' title='p' alt='p' align=absmiddle> in the factorization of <img src='/blog/latexrender/pictures/388f554901ba5d77339eec8b26beebea.png' title='n!' alt='n!' align=absmiddle> is given by <img src='/blog/latexrender/pictures/2f316ef04cacd19d4149282a51dc49ef.png' title='\sum_{j=1}^{\infty} \lfloor n / p^j \rfloor' alt='\sum_{j=1}^{\infty} \lfloor n / p^j \rfloor' align=absmiddle>.  As a result, <img src='/blog/latexrender/pictures/7bccb23d063d614019457f0cee01943c.png' title='r(p,n) = \sum_{j=1}^\infty ( \lfloor 2n / p^j \rfloor &amp;#8211; 2 \lfloor n / p^j \rfloor)' alt='r(p,n) = \sum_{j=1}^\infty ( \lfloor 2n / p^j \rfloor &amp;#8211; 2 \lfloor n / p^j \rfloor)' align=absmiddle>.  All terms in this sum are either <img src='/blog/latexrender/pictures/cfcd208495d565ef66e7dff9f98764da.png' title='0' alt='0' align=absmiddle> or <img src='/blog/latexrender/pictures/c4ca4238a0b923820dcc509a6f75849b.png' title='1' alt='1' align=absmiddle>, depending on whether <img src='/blog/latexrender/pictures/c437c37fbebdea6db114e5367842b568.png' title='\lfloor 2n/p^j \rfloor' alt='\lfloor 2n/p^j \rfloor' align=absmiddle> is even or odd, and they are always <img src='/blog/latexrender/pictures/cfcd208495d565ef66e7dff9f98764da.png' title='0' alt='0' align=absmiddle> when <img src='/blog/latexrender/pictures/febd5a71aa0c2120073fd8ccb0e33e61.png' title='j &gt; \log_p (2n)' alt='j &gt; \log_p (2n)' align=absmiddle>, whence <img src='/blog/latexrender/pictures/4ddc2f9f8fb52d8c2f4be557c9ccb60b.png' title='r(p,n) \leq \log_p (2n)' alt='r(p,n) \leq \log_p (2n)' align=absmiddle>, which gives the main result.</p>
<p>In the special case of (1), <img src='/blog/latexrender/pictures/b1393237013f9d5bb3041ae6d70fc2a9.png' title='p &gt; \sqrt{2n}' alt='p &gt; \sqrt{2n}' align=absmiddle>, simply observe that all but the first term in the sum vanishes, and so <img src='/blog/latexrender/pictures/d3565f936ef584f453fe0e4da31ec147.png' title='r(p,n) = \lfloor 2n/ p \rfloor &amp;#8211; 2 \lfloor n/p \rfloor' alt='r(p,n) = \lfloor 2n/ p \rfloor &amp;#8211; 2 \lfloor n/p \rfloor' align=absmiddle>, which is either <img src='/blog/latexrender/pictures/cfcd208495d565ef66e7dff9f98764da.png' title='0' alt='0' align=absmiddle> or <img src='/blog/latexrender/pictures/c4ca4238a0b923820dcc509a6f75849b.png' title='1' alt='1' align=absmiddle>, according to whether <img src='/blog/latexrender/pictures/bf04c2e333a9b6abd0cd922d31d7b9bd.png' title='\lfloor 2n/p \rfloor' alt='\lfloor 2n/p \rfloor' align=absmiddle> is even or odd.</p>
<p>In the special case of (2), <img src='/blog/latexrender/pictures/f2c0ed823881eadfa1347431940736f9.png' title='2n/3 &lt; p \leq n' alt='2n/3 &lt; p \leq n' align=absmiddle>, observe again that all but the first term in the sum vanishes, and so <img src='/blog/latexrender/pictures/fb4d481f572221f04f3f9aecc70552ef.png' title='r(p,n) = \lfloor 2n/p \rfloor - 2 \lfloor n/p \rfloor' alt='r(p,n) = \lfloor 2n/p \rfloor - 2 \lfloor n/p \rfloor' align=absmiddle>.  But for any such value of <img src='/blog/latexrender/pictures/83878c91171338902e0fe0fb97a8c47a.png' title='p' alt='p' align=absmiddle> we see that <img src='/blog/latexrender/pictures/bf04c2e333a9b6abd0cd922d31d7b9bd.png' title='\lfloor 2n/p \rfloor' alt='\lfloor 2n/p \rfloor' align=absmiddle> is equal to <img src='/blog/latexrender/pictures/c81e728d9d4c2f636f067f89cc14862c.png' title='2' alt='2' align=absmiddle>, which is even, and so <img src='/blog/latexrender/pictures/5481cc4cea0948af9b690474ea594be2.png' title='r(p,n) = 0' alt='r(p,n) = 0' align=absmiddle>. <strong>QED</strong></p>
<p>The proof of Bertrand&#8217;s postulate is now straightforward.</p>
<p><strong>Proof of Bertrand&#8217;s postulate:</strong> For a contradiction, suppose there are no primes satisfying <img src='/blog/latexrender/pictures/635729275f58d8d3d6003589d105f73c.png' title='n &lt; p &lt; 2n' alt='n &lt; p &lt; 2n' align=absmiddle>.  In light of case (2) of the last lemma, we can assume that all the prime factors in <img src='/blog/latexrender/pictures/e25dc50d34d9862186ff1c39580591fa.png' title='{2n   \choose n}' alt='{2n   \choose n}' align=absmiddle> satisfy <img src='/blog/latexrender/pictures/09b54b5332d3a7f79eec802d95962bfd.png' title='p \leq 2n/3' alt='p \leq 2n/3' align=absmiddle>.  We split the prime factorizing of <img src='/blog/latexrender/pictures/c92da943df73dc077dbee5514376346a.png' title='{2n \choose n}' alt='{2n \choose n}' align=absmiddle> up into the cases where <img src='/blog/latexrender/pictures/7eb8c06079a5340592ac98b757c46e6a.png' title='p \leq \sqrt{2n}' alt='p \leq \sqrt{2n}' align=absmiddle> and where <img src='/blog/latexrender/pictures/c963f69cdb695898e3afa5341871df32.png' title='\sqrt{2n} &lt; p \leq 2n/3' alt='\sqrt{2n} &lt; p \leq 2n/3' align=absmiddle>.  Applying the last lemma we obtain:</p>
<p><img src='/blog/latexrender/pictures/c92e8b85eb957e35bcac5742ba9d1bf3.png' title='   {2n \choose n} \leq \prod_{p \leq \sqrt{2n}} (2n) \prod_{\sqrt{2n} &lt; p \leq 2n/3} p ' alt='   {2n \choose n} \leq \prod_{p \leq \sqrt{2n}} (2n) \prod_{\sqrt{2n} &lt; p \leq 2n/3} p ' align=absmiddle></p>
<p>We can bound the first product by noting that there are no more than <img src='/blog/latexrender/pictures/37d9e328b24ee5601a25c0923dc30db3.png' title='\sqrt{2n}' alt='\sqrt{2n}' align=absmiddle> primes of size at most <img src='/blog/latexrender/pictures/37d9e328b24ee5601a25c0923dc30db3.png' title='\sqrt{2n}' alt='\sqrt{2n}' align=absmiddle>, and bound the second product using Chebyshev's bound:</p>
<p><img src='/blog/latexrender/pictures/461ef0dc1b4e6241280611ea72d7547d.png' title='   {2n \choose n} \leq (2n)^{\sqrt 2n} 4^{2n/3}. ' alt='   {2n \choose n} \leq (2n)^{\sqrt 2n} 4^{2n/3}. ' align=absmiddle></p>
<p>Observing that <img src='/blog/latexrender/pictures/c92da943df73dc077dbee5514376346a.png' title='{2n \choose n}' alt='{2n \choose n}' align=absmiddle> is the largest of the <img src='/blog/latexrender/pictures/d2b2d9fec288403faf6e85ebf2c58972.png' title='2n+1' alt='2n+1' align=absmiddle> terms in the binomial expansion of <img src='/blog/latexrender/pictures/7137f9916b1410e79dec2e15b0e91727.png' title='(1+1)^{2n}=4^n' alt='(1+1)^{2n}=4^n' align=absmiddle>, we see that</p>
<p><img src='/blog/latexrender/pictures/d1b6366652b2585faaae37dff6093376.png' title='   \frac{4^n}{2n+1} \leq {2n \choose n} ' alt='   \frac{4^n}{2n+1} \leq {2n \choose n} ' align=absmiddle></p>
<p>Combinining the last two inequalities and rearranging a little gives</p>
<p><img src='/blog/latexrender/pictures/fef1bf31c7abf68a5a206744afb8391f.png' title='   4^{n/3} \leq (2n+1)(2n)^{\sqrt 2n}. ' alt='   4^{n/3} \leq (2n+1)(2n)^{\sqrt 2n}. ' align=absmiddle></p>
<p>The left-hand side rises much faster than the right, and so this inequality can only be true over some finite set of <img src='/blog/latexrender/pictures/7b8b965ad4bca0e41ab51de7b31363a1.png' title='n' alt='n' align=absmiddle>.  It should be straightforward to convince yourself that it fails for <img src='/blog/latexrender/pictures/14fdfd4e619edc4e7a0698e282d86c71.png' title='n &gt; 2048' alt='n &gt; 2048' align=absmiddle>, for example, and thus we must have <img src='/blog/latexrender/pictures/83d8fdac653b3d79d6ece509d4c43511.png' title='n \leq 2048' alt='n \leq 2048' align=absmiddle>.  In fact, we can do quite a bit better than <img src='/blog/latexrender/pictures/83d8fdac653b3d79d6ece509d4c43511.png' title='n \leq 2048' alt='n \leq 2048' align=absmiddle> with a bit more work, but the bound suffices for our purposes.  What we have now shown is that our initial assumption &#8211; no prime in the range <img src='/blog/latexrender/pictures/635729275f58d8d3d6003589d105f73c.png' title='n &lt; p &lt; 2n' alt='n &lt; p &lt; 2n' align=absmiddle> - can only possibly hold when <img src='/blog/latexrender/pictures/83d8fdac653b3d79d6ece509d4c43511.png' title='n \leq 2048' alt='n \leq 2048' align=absmiddle>.  But it's easily seen that the assumption is also false in that range, just by considering the sequence of primes <img src='/blog/latexrender/pictures/787812206553a33181020112c548c4ae.png' title='2,3,5,7,13,23,43,83,163,317,631,1259' alt='2,3,5,7,13,23,43,83,163,317,631,1259' align=absmiddle>, and <img src='/blog/latexrender/pictures/e71e5cd119bbc5797164fb0cd7fd94a4.png' title='2503' alt='2503' align=absmiddle>. Each prime in the sequence is less than twice its predecessor, and so there must be a prime between <img src='/blog/latexrender/pictures/7b8b965ad4bca0e41ab51de7b31363a1.png' title='n' alt='n' align=absmiddle> and <img src='/blog/latexrender/pictures/21e2c0c0472b331622877accbe29b91b.png' title='2n' alt='2n' align=absmiddle> for any <img src='/blog/latexrender/pictures/83d8fdac653b3d79d6ece509d4c43511.png' title='n \leq 2048' alt='n \leq 2048' align=absmiddle>. <strong>QED</strong></p>
]]></content:encoded>
			<wfw:commentRss>http://michaelnielsen.org/blog/bertrands-postulate/feed/</wfw:commentRss>
		<slash:comments>3</slash:comments>
		</item>
		<item>
		<title>The Polymath blog</title>
		<link>http://michaelnielsen.org/blog/the-polymath-blog/</link>
		<comments>http://michaelnielsen.org/blog/the-polymath-blog/#comments</comments>
		<pubDate>Wed, 29 Jul 2009 02:34:05 +0000</pubDate>
		<dc:creator>Michael Nielsen</dc:creator>
				<category><![CDATA[Polymath]]></category>

		<guid isPermaLink="false">http://michaelnielsen.org/blog/the-polymath-blog/</guid>
		<description><![CDATA[Earlier this year, Tim Gowers started a project in massively collaborative mathematics &#8211; an open approach to solving mathematical problems using blogs and wikis. The first iteration of this &#8220;Polymath Project&#8221; was very successful (see also Terry Tao’s recent mini-Polymath), and new iterations are now being planned. To help with that process, Terry Tao has [...]]]></description>
			<content:encoded><![CDATA[<p>Earlier this year, Tim Gowers started a project in <a href="http://gowers.wordpress.com/2009/01/27/is-massively-collaborative-mathematics-possible/">massively collaborative mathematics</a> &#8211; an open approach to solving mathematical problems using blogs and wikis.  The first iteration of this &#8220;Polymath Project&#8221; was <a href="http://michaelnielsen.org/blog/?p=584">very successful</a> (see also Terry Tao’s recent <a href="http://terrytao.wordpress.com/2009/07/22/imo-2009-q6-mini-polymath-project-impressions-reflections-analysis/#comment-40605">mini-Polymath</a>), and new iterations are now being planned.  To help with that process, <a href="http://terrytao.wordpress.com/2009/07/22/imo-2009-q6-mini-polymath-project-impressions-reflections-analysis/#comment-40605">Terry Tao</a> has set up a <a href="http://polymathprojects.org">Polymath blog</a>, and there is now a very lively discussion going on about possible problems, including the very interesting problem of finding an efficient deterministic algorithm to generate prime numbers above a specified size.  </p>
]]></content:encoded>
			<wfw:commentRss>http://michaelnielsen.org/blog/the-polymath-blog/feed/</wfw:commentRss>
		<slash:comments>2</slash:comments>
		</item>
	</channel>
</rss>

