<?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>Mathematical Food for Thought</title>
	<atom:link href="http://wangsblog.com/jeffrey/?feed=rss2" rel="self" type="application/rss+xml" />
	<link>http://wangsblog.com/jeffrey</link>
	<description>Serves a Daily Special and an All-You-Can-Eat Course in Problem Solving. Courtesy of me, Jeffrey Wang.</description>
	<lastBuildDate>Sun, 08 Jul 2007 00:22:04 +0000</lastBuildDate>
	<generator>http://wordpress.org/?v=2.9.2</generator>
	<language>en</language>
	<sy:updatePeriod>hourly</sy:updatePeriod>
	<sy:updateFrequency>1</sy:updateFrequency>
			<item>
		<title>What&#8217;s Your Function? Topic: Algebra. Level: AMC/AIME.</title>
		<link>http://wangsblog.com/jeffrey/?p=333</link>
		<comments>http://wangsblog.com/jeffrey/?p=333#comments</comments>
		<pubDate>Sat, 07 Jul 2007 23:56:15 +0000</pubDate>
		<dc:creator>paladin8</dc:creator>
				<category><![CDATA[Algebra]]></category>

		<guid isPermaLink="false">http://wangsblog.com/jeffrey/?p=333</guid>
		<description><![CDATA[Problem: Given two positive reals  and , show that there is a continuous function  that satisfies .
Solution: There are several special cases that are interesting to look at before we make a guess as to what type of function  will be. First we consider the case . This immediately gives
.
It should not [...]]]></description>
			<content:encoded><![CDATA[<p><strong>Problem</strong>: Given two positive reals <img src='http://www.wangsblog.com/jeffrey//pictures/9845198045ec71aa8304372c28b24630.gif' title=' \alpha ' alt=' \alpha ' align=absmiddle> and <img src='http://www.wangsblog.com/jeffrey//pictures/2fa6e1a5937d65ac9e67f692cebcddef.gif' title=' \beta ' alt=' \beta ' align=absmiddle>, show that there is a continuous function <img src='http://www.wangsblog.com/jeffrey//pictures/ce40937fdfbd06b8a15244e102a09356.gif' title=' f ' alt=' f ' align=absmiddle> that satisfies <img src='http://www.wangsblog.com/jeffrey//pictures/9d39b5f804c21d83f28230d3ef831935.gif' title=' f(x) = f(x+\alpha)+f(x+\beta) ' alt=' f(x) = f(x+\alpha)+f(x+\beta) ' align=absmiddle>.</p>
<p><strong>Solution</strong>: There are several special cases that are interesting to look at before we make a guess as to what type of function <img src='http://www.wangsblog.com/jeffrey//pictures/ce40937fdfbd06b8a15244e102a09356.gif' title=' f ' alt=' f ' align=absmiddle> will be. First we consider the case <img src='http://www.wangsblog.com/jeffrey//pictures/33d7a7bb70a50d0e88ff81f72e2fc183.gif' title=' \alpha = \beta = 1 ' alt=' \alpha = \beta = 1 ' align=absmiddle>. This immediately gives</p>
<p><img src='http://www.wangsblog.com/jeffrey//pictures/2f09355723738ab5036e3d0dfe37bdb4.gif' title=' f(x) = 2f(x+1) ' alt=' f(x) = 2f(x+1) ' align=absmiddle>.</p>
<p>It should not be too difficult to guess that <img src='http://www.wangsblog.com/jeffrey//pictures/a6260b26a112fcc8dffe89a182d44379.gif' title=' f(x) = 2^{-x} ' alt=' f(x) = 2^{-x} ' align=absmiddle> is a solution to this, as well as any constant multiple of it. Now try <img src='http://www.wangsblog.com/jeffrey//pictures/146e87d59e07356736448e22266bc6ce.gif' title=' \alpha = -1 ' alt=' \alpha = -1 ' align=absmiddle> and <img src='http://www.wangsblog.com/jeffrey//pictures/4a32f4a6b1c9f4bb4526458ffab34d16.gif' title=' \beta = -2 ' alt=' \beta = -2 ' align=absmiddle>, resulting in</p>
<p><img src='http://www.wangsblog.com/jeffrey//pictures/801875df4a642864874bcfb68a668de3.gif' title=' f(x) = f(x-1)+f(x-2) ' alt=' f(x) = f(x-1)+f(x-2) ' align=absmiddle>.</p>
<p>Looks a lot like Fibonacci, right? In fact, one possible function is just <img src='http://www.wangsblog.com/jeffrey//pictures/c8a3fbb1d467eb92567a6b040840fafd.gif' title=' f(x) = \phi^x ' alt=' f(x) = \phi^x ' align=absmiddle>. That&#8217;s pretty convenient.</p>
<p>Notice how both of these illuminating examples are exponential functions, which leads us to guess that our function will be exponential as well. So, following this track, we set</p>
<p><img src='http://www.wangsblog.com/jeffrey//pictures/2a7ab738136f42d38323542b8424dce2.gif' title=' f(x) = a^x ' alt=' f(x) = a^x ' align=absmiddle></p>
<p>so we simply need to solve</p>
<p><img src='http://www.wangsblog.com/jeffrey//pictures/344b11b3eebb997e8a589681f8258b05.gif' title=' a^x = a^{x+\alpha}+a^{x+\beta} ' alt=' a^x = a^{x+\alpha}+a^{x+\beta} ' align=absmiddle></p>
<p><img src='http://www.wangsblog.com/jeffrey//pictures/5c781d2ee37d28ef02d99f0fbcd3b7ec.gif' title=' a^x = a^x\left(a^{\alpha}+a^{\beta}\right) ' alt=' a^x = a^x\left(a^{\alpha}+a^{\beta}\right) ' align=absmiddle>.</p>
<p>Unfortunately, the equation <img src='http://www.wangsblog.com/jeffrey//pictures/d4e26c4ce4d23927d17c94714c9d5ae8.gif' title=' a^{\alpha}+a^{\beta} = 1 ' alt=' a^{\alpha}+a^{\beta} = 1 ' align=absmiddle> does not always have a solution (take <img src='http://www.wangsblog.com/jeffrey//pictures/e0f469259befd58058e9f8f726e4d159.gif' title=' \alpha = 1 ' alt=' \alpha = 1 ' align=absmiddle> and <img src='http://www.wangsblog.com/jeffrey//pictures/2612d4d6d9825fa08811f3b967188952.gif' title=' \beta = -1 ' alt=' \beta = -1 ' align=absmiddle> for example). But that&#8217;s ok and I&#8217;ll worry about it some other time. In any case we have found a function for whenever the equation <img src='http://www.wangsblog.com/jeffrey//pictures/d4e26c4ce4d23927d17c94714c9d5ae8.gif' title=' a^{\alpha}+a^{\beta} = 1 ' alt=' a^{\alpha}+a^{\beta} = 1 ' align=absmiddle> has a solution in the reals.</p>
]]></content:encoded>
			<wfw:commentRss>http://wangsblog.com/jeffrey/?feed=rss2&amp;p=333</wfw:commentRss>
		<slash:comments>5</slash:comments>
		</item>
		<item>
		<title>Addition At Its Finest. Topic: Calculus/S&amp;S.</title>
		<link>http://wangsblog.com/jeffrey/?p=332</link>
		<comments>http://wangsblog.com/jeffrey/?p=332#comments</comments>
		<pubDate>Sat, 30 Jun 2007 06:12:33 +0000</pubDate>
		<dc:creator>paladin8</dc:creator>
				<category><![CDATA[Calculus]]></category>
		<category><![CDATA[Sequences &#038; Series]]></category>

		<guid isPermaLink="false">http://wangsblog.com/jeffrey/?p=332</guid>
		<description><![CDATA[Problem: Evaluate  where  is a real number with .
Solution: Looking at that all too common denominator, we do a partial fraction decomposition in hopes of telescoping series. The summation becomes
.
Common Taylor series knowledge tells us that
,
which convenient fits the first part of the summation. As for the second part, we get

from the same [...]]]></description>
			<content:encoded><![CDATA[<p><strong>Problem</strong>: Evaluate <img src='http://www.wangsblog.com/jeffrey//pictures/681d434a5549ebfdded8c46598f1c172.gif' title=' \displaystyle \sum_{n=1}^{\infty} \frac{x^n}{n(n+1)} ' alt=' \displaystyle \sum_{n=1}^{\infty} \frac{x^n}{n(n+1)} ' align=absmiddle> where <img src='http://www.wangsblog.com/jeffrey//pictures/6722c218a6f30869ef6886dc4b050a37.gif' title=' x ' alt=' x ' align=absmiddle> is a real number with <img src='http://www.wangsblog.com/jeffrey//pictures/c21ca86348bd3cbcfe65d7684d8e3f2e.gif' title=' |x| &lt; 1 ' alt=' |x| &lt; 1 ' align=absmiddle>.</p>
<p>Solution: Looking at that all too common denominator, we do a partial fraction decomposition in hopes of telescoping series. The summation becomes</p>
<p><img src='http://www.wangsblog.com/jeffrey//pictures/7c18462c1c458d0f0e2466c144a7c905.gif' title=' \displaystyle \sum_{n=1}^{\infty} \left(\frac{x^n}{n}-\frac{x^n}{n+1}\right) ' alt=' \displaystyle \sum_{n=1}^{\infty} \left(\frac{x^n}{n}-\frac{x^n}{n+1}\right) ' align=absmiddle>.</p>
<p>Common Taylor series knowledge tells us that</p>
<p><img src='http://www.wangsblog.com/jeffrey//pictures/3c9b5ef2de0c398760e493cb16cb2b9b.gif' title=' \displaystyle \ln{(1-x)} = -\left(x+\frac{x^2}{2}+\frac{x^3}{3}+\cdots\right) = -\sum_{n=1}^{\infty} \frac{x^n}{n} ' alt=' \displaystyle \ln{(1-x)} = -\left(x+\frac{x^2}{2}+\frac{x^3}{3}+\cdots\right) = -\sum_{n=1}^{\infty} \frac{x^n}{n} ' align=absmiddle>,</p>
<p>which convenient fits the first part of the summation. As for the second part, we get</p>
<p><img src='http://www.wangsblog.com/jeffrey//pictures/58b5def4868ece420c256706d7fa3be4.gif' title=' \displaystyle \sum_{n=1}^{\infty} \frac{x^n}{n+1} = \frac{1}{x} \sum_{n=1}^{\infty} \frac{x^{n+1}}{n+1} = \frac{-\ln{(1-x)}-x}{x} ' alt=' \displaystyle \sum_{n=1}^{\infty} \frac{x^n}{n+1} = \frac{1}{x} \sum_{n=1}^{\infty} \frac{x^{n+1}}{n+1} = \frac{-\ln{(1-x)}-x}{x} ' align=absmiddle></p>
<p>from the same Taylor series. Combining the results, our answer is then</p>
<p><img src='http://www.wangsblog.com/jeffrey//pictures/ab5a302f3bca030c6d84db0bab22ee89.gif' title=' \displaystyle \sum_{n=1}^{\infty} \frac{x^n}{n(n+1)} = 1-\ln{(1-x)}+\frac{\ln{(1-x)}}{x} ' alt=' \displaystyle \sum_{n=1}^{\infty} \frac{x^n}{n(n+1)} = 1-\ln{(1-x)}+\frac{\ln{(1-x)}}{x} ' align=absmiddle>.</p>
<p>QED.</p>
<p>&#8212;&#8212;&#8212;&#8212;&#8212;&#8212;&#8211;</p>
<p>Comment: Even though the trick at the beginning didn&#8217;t actually get much to telescope, the idea certainly made it easier to recognize the Taylor series. Algebraic manipulations are nifty to carry around and can be applied in problems wherever you go.</p>
<p>&#8212;&#8212;&#8212;&#8212;&#8212;&#8212;&#8211;</p>
<p>Practice Problem: Show that <img src='http://www.wangsblog.com/jeffrey//pictures/9eaded4448a6d29039eeb505fbfb45e9.gif' title=' \displaystyle \int_0^{\frac{\pi}{2}} \ln{(\tan{x})} = 0 ' alt=' \displaystyle \int_0^{\frac{\pi}{2}} \ln{(\tan{x})} = 0 ' align=absmiddle>.</p>
]]></content:encoded>
			<wfw:commentRss>http://wangsblog.com/jeffrey/?feed=rss2&amp;p=332</wfw:commentRss>
		<slash:comments>4</slash:comments>
		</item>
		<item>
		<title>The Smaller The Better. Topic: Calculus.</title>
		<link>http://wangsblog.com/jeffrey/?p=331</link>
		<comments>http://wangsblog.com/jeffrey/?p=331#comments</comments>
		<pubDate>Tue, 19 Jun 2007 03:37:57 +0000</pubDate>
		<dc:creator>paladin8</dc:creator>
				<category><![CDATA[Calculus]]></category>

		<guid isPermaLink="false">http://wangsblog.com/jeffrey/?p=331</guid>
		<description><![CDATA[Problem: Given a complicated function , find an approximate local minimum.
Solution: The adjective complicated is only placed so that we assume there is no easy way to solve  to immediately give the solution. We seek an algorithm that will lead us to a local minimum (hopefully a global minimum as well).
We start at an [...]]]></description>
			<content:encoded><![CDATA[<p><strong>Problem</strong>: Given a <em>complicated</em> function <img src='http://www.wangsblog.com/jeffrey//pictures/c70846ba38580c22cd60162d2f57a10a.gif' title=' f: \mathbb{R}^n \rightarrow \mathbb{R} ' alt=' f: \mathbb{R}^n \rightarrow \mathbb{R} ' align=absmiddle>, find an approximate local minimum.</p>
<p><strong>Solution</strong>: The adjective <em>complicated</em> is only placed so that we assume there is no easy way to solve <img src='http://www.wangsblog.com/jeffrey//pictures/6802c8c0e23c9c8c49e95ff8088b8ae2.gif' title=' \bigtriangledown f = 0 ' alt=' \bigtriangledown f = 0 ' align=absmiddle> to immediately give the solution. We seek an algorithm that will lead us to a local minimum (hopefully a global minimum as well).</p>
<p>We start at an arbitrary point <img src='http://www.wangsblog.com/jeffrey//pictures/2d439ce9231f5ce9a8371ed060b2ebae.gif' title=' X_0 = (x_1, x_2, \ldots, x_n) ' alt=' X_0 = (x_1, x_2, \ldots, x_n) ' align=absmiddle>. Consider the following process (for <img src='http://www.wangsblog.com/jeffrey//pictures/df910baf4f24a5d3a3796c0fb3b1f4e9.gif' title=' k = 0, 1, 2, \ldots ' alt=' k = 0, 1, 2, \ldots ' align=absmiddle>), known as gradient descent:</p>
<p>1. Calculate (approximately) <img src='http://www.wangsblog.com/jeffrey//pictures/4385883096806996b3c211ddd0693311.gif' title=' \bigtriangledown f(X_k) ' alt=' \bigtriangledown f(X_k) ' align=absmiddle>.</p>
<p>2. Set <img src='http://www.wangsblog.com/jeffrey//pictures/2290c1ac4d4216ed444426740550283f.gif' title=' X_{k+1} = X_k &amp;#8211; \gamma_k \bigtriangledown f(X_k) ' alt=' X_{k+1} = X_k &amp;#8211; \gamma_k \bigtriangledown f(X_k) ' align=absmiddle>, where <img src='http://www.wangsblog.com/jeffrey//pictures/f0f79156083c7774f82be50a478502c4.gif' title=' \gamma_k ' alt=' \gamma_k ' align=absmiddle> is a constant that can be determined by a <a href="http://en.wikipedia.org/wiki/Line_search">linesearch</a>.</p>
<p>It is well-known that the direction of the gradient is the direction of maximum increase and the direction opposite the gradient is the direction of maximum decrease. Hence this algorithm is based on the idea that we always move in the direction that will decrease <img src='http://www.wangsblog.com/jeffrey//pictures/ce40937fdfbd06b8a15244e102a09356.gif' title=' f ' alt=' f ' align=absmiddle> the most. Sounds pretty good, right? Well, unfortunately gradient descent converges very slowly so it is only really useful for smaller optimization problems. Fortunately, there exist other algorithms but obviously they are more complex, such as the <a href="http://en.wikipedia.org/wiki/Nonlinear_conjugate_gradient">nonlinear conjugate gradient method</a> or <a href="http://en.wikipedia.org/wiki/Newton%27s_method_in_optimization">Newton&#8217;s method</a>, the latter of which involves the computation of the inverse of the Hessian matrix, which is a pain.</p>
]]></content:encoded>
			<wfw:commentRss>http://wangsblog.com/jeffrey/?feed=rss2&amp;p=331</wfw:commentRss>
		<slash:comments>0</slash:comments>
		</item>
		<item>
		<title>Colorful! Topic: Calculus.</title>
		<link>http://wangsblog.com/jeffrey/?p=330</link>
		<comments>http://wangsblog.com/jeffrey/?p=330#comments</comments>
		<pubDate>Thu, 07 Jun 2007 01:15:12 +0000</pubDate>
		<dc:creator>paladin8</dc:creator>
				<category><![CDATA[Calculus]]></category>

		<guid isPermaLink="false">http://wangsblog.com/jeffrey/?p=330</guid>
		<description><![CDATA[Theorem: (Green&#8217;s Theorem) Let  be a simply connected plane region whose boundary is a simple, closed, piecewise smooth curve  oriented counterclockwise. If  and  are continuous and have continuous first partial derivatives on some open set containing , then
.
&#8212;&#8212;&#8212;&#8212;&#8212;&#8212;&#8211;
Problem: Evaluate , where  is the boundary of the region enclosed by  [...]]]></description>
			<content:encoded><![CDATA[<p><strong>Theorem</strong>: (Green&#8217;s Theorem) Let <img src='http://www.wangsblog.com/jeffrey//pictures/75695b46abca7ce53dfa3b4e984a45ca.gif' title=' R ' alt=' R ' align=absmiddle> be a simply connected plane region whose boundary is a simple, closed, piecewise smooth curve <img src='http://www.wangsblog.com/jeffrey//pictures/699e902f5598ca623370f833cffb1a57.gif' title=' C ' alt=' C ' align=absmiddle> oriented counterclockwise. If <img src='http://www.wangsblog.com/jeffrey//pictures/c5a558ec2d4a1e8bb63c8e4c3e744d41.gif' title=' f(x, y) ' alt=' f(x, y) ' align=absmiddle> and <img src='http://www.wangsblog.com/jeffrey//pictures/dac1a8e648c0b3757f482b3134f16813.gif' title=' g(x, y) ' alt=' g(x, y) ' align=absmiddle> are continuous and have continuous first partial derivatives on some open set containing <img src='http://www.wangsblog.com/jeffrey//pictures/75695b46abca7ce53dfa3b4e984a45ca.gif' title=' R ' alt=' R ' align=absmiddle>, then</p>
<p><img src='http://www.wangsblog.com/jeffrey//pictures/a4dd96a899e6b188fad8352d3c2afbf7.gif' title=' \displaystyle \oint_C f(x, y) dx + g(x, y) dy = \int_R \int \left(\frac{\partial g}{\partial x}-\frac{\partial f}{\partial y}\right) dA ' alt=' \displaystyle \oint_C f(x, y) dx + g(x, y) dy = \int_R \int \left(\frac{\partial g}{\partial x}-\frac{\partial f}{\partial y}\right) dA ' align=absmiddle>.</p>
<p>&#8212;&#8212;&#8212;&#8212;&#8212;&#8212;&#8211;</p>
<p><strong>Problem</strong>: Evaluate <img src='http://www.wangsblog.com/jeffrey//pictures/a82445f7f88d3dc221e2667b24ac48d3.gif' title=' \displaystyle \oint_C x^2y dx + (y+xy^2) dy ' alt=' \displaystyle \oint_C x^2y dx + (y+xy^2) dy ' align=absmiddle>, where <img src='http://www.wangsblog.com/jeffrey//pictures/699e902f5598ca623370f833cffb1a57.gif' title=' C ' alt=' C ' align=absmiddle> is the boundary of the region enclosed by <img src='http://www.wangsblog.com/jeffrey//pictures/feb38954971364ac01eaf10660813d07.gif' title=' y = x^2 ' alt=' y = x^2 ' align=absmiddle> and <img src='http://www.wangsblog.com/jeffrey//pictures/ac0976bb0785fdd32136b778b6e25876.gif' title=' x = y^2 ' alt=' x = y^2 ' align=absmiddle>.</p>
<p><strong>Solution</strong>: First, verify that this region satisfies all of the requirements for Green&#8217;s Theorem &#8211; indeed, it does. So we may apply the theorem with <img src='http://www.wangsblog.com/jeffrey//pictures/d9238d809fc3ef6c257110639d2305c9.gif' title=' f(x, y) = x^2y ' alt=' f(x, y) = x^2y ' align=absmiddle> and <img src='http://www.wangsblog.com/jeffrey//pictures/76dc21e78238864feb12bca8b6837074.gif' title=' g(x, y) = y+xy^2 ' alt=' g(x, y) = y+xy^2 ' align=absmiddle>. From these, we have <img src='http://www.wangsblog.com/jeffrey//pictures/2434d4b7c274be22987cdb47779dea46.gif' title=' \frac{\partial g}{\partial x} = y^2 ' alt=' \frac{\partial g}{\partial x} = y^2 ' align=absmiddle> and <img src='http://www.wangsblog.com/jeffrey//pictures/b31f30944c96b9a0d9c5085e4428e6b1.gif' title=' \frac{\partial f}{\partial y} = x^2 ' alt=' \frac{\partial f}{\partial y} = x^2 ' align=absmiddle>. Then we obtain</p>
<p><img src='http://www.wangsblog.com/jeffrey//pictures/b51d79924039588a707d7e97b216028a.gif' title=' \displaystyle \oint_C x^2y dx + (y+xy^2) dy = \int_R \int (y^2-x^2) dA ' alt=' \displaystyle \oint_C x^2y dx + (y+xy^2) dy = \int_R \int (y^2-x^2) dA ' align=absmiddle>.<br />
But clearly this integral over the region <img src='http://www.wangsblog.com/jeffrey//pictures/75695b46abca7ce53dfa3b4e984a45ca.gif' title=' R ' alt=' R ' align=absmiddle> can be represented as <img src='http://www.wangsblog.com/jeffrey//pictures/25e0b38992df9249fda1b8ef90ea94d3.gif' title=' \displaystyle \int_0^1 \int_{x^2}^{\sqrt{x}} (y^2-x^2) dy dx ' alt=' \displaystyle \int_0^1 \int_{x^2}^{\sqrt{x}} (y^2-x^2) dy dx ' align=absmiddle>, so it remains a matter of calculation to get the answer. First, we evaluate the inner integral to get</p>
<p><img src='http://www.wangsblog.com/jeffrey//pictures/c3332ec8db9094ba053721e2ebdcadb5.gif' title=' \displaystyle \int_0^1 \int_{x^2}^{\sqrt{x}} (y^2-x^2) dy dx = \int_0^1 \left[\frac{y^3}{3}-x^2y\right]_{x^2}^{\sqrt{x}} dx = \int_0^1 \left(\frac{x^{3/2}}{3}-x^{5/2}-\frac{x^6}{3}+x^4\right)dx ' alt=' \displaystyle \int_0^1 \int_{x^2}^{\sqrt{x}} (y^2-x^2) dy dx = \int_0^1 \left[\frac{y^3}{3}-x^2y\right]_{x^2}^{\sqrt{x}} dx = \int_0^1 \left(\frac{x^{3/2}}{3}-x^{5/2}-\frac{x^6}{3}+x^4\right)dx ' align=absmiddle>.</p>
<p>Then finally we have</p>
<p><img src='http://www.wangsblog.com/jeffrey//pictures/d47fb50d217ebd922b1645f5847a9066.gif' title=' \displaystyle \int_0^1 \left(\frac{x^{3/2}}{3}-x^{5/2}-\frac{x^6}{3}+x^4\right)dx = \left[\frac{2x^{5/2}}{15}-\frac{2x^{7/2}}{7}-\frac{x^7}{21}+\frac{x^5}{5}\right]_0^1 = \frac{2}{15}-\frac{2}{7}-\frac{1}{21}+\frac{1}{5} = 0 ' alt=' \displaystyle \int_0^1 \left(\frac{x^{3/2}}{3}-x^{5/2}-\frac{x^6}{3}+x^4\right)dx = \left[\frac{2x^{5/2}}{15}-\frac{2x^{7/2}}{7}-\frac{x^7}{21}+\frac{x^5}{5}\right]_0^1 = \frac{2}{15}-\frac{2}{7}-\frac{1}{21}+\frac{1}{5} = 0 ' align=absmiddle>.</p>
<p>QED.</p>
<p>&#8212;&#8212;&#8212;&#8212;&#8212;&#8212;&#8211;</p>
<p>Comment: To me, Green&#8217;s Theorem is a very interesting result. It&#8217;s not at all obvious that a line integral along the boundary of a region is equivalent to an integral of some partial derivatives in the region itself. A simplified proof of the result can be obtained by proving that</p>
<p><img src='http://www.wangsblog.com/jeffrey//pictures/903d1a8937f622681a95027f9ee07065.gif' title=' \displaystyle \oint_C f(x, y) dx = -\int_R \int \frac{\partial f}{\partial y} dA ' alt=' \displaystyle \oint_C f(x, y) dx = -\int_R \int \frac{\partial f}{\partial y} dA ' align=absmiddle> and <img src='http://www.wangsblog.com/jeffrey//pictures/7a2d8c110b710829a1f538800f64979c.gif' title=' \displaystyle \oint_C g(x, y) dy = \int_R \int \frac{\partial g}{\partial x} dA ' alt=' \displaystyle \oint_C g(x, y) dy = \int_R \int \frac{\partial g}{\partial x} dA ' align=absmiddle>.</p>
<p>&#8212;&#8212;&#8212;&#8212;&#8212;&#8212;&#8211;</p>
<p>Practice Problem: Let <img src='http://www.wangsblog.com/jeffrey//pictures/75695b46abca7ce53dfa3b4e984a45ca.gif' title=' R ' alt=' R ' align=absmiddle> be a plane region with area <img src='http://www.wangsblog.com/jeffrey//pictures/951196ca354c5c72b8356494f97c3b5d.gif' title=' A ' alt=' A ' align=absmiddle> whose boundary is a piecewise smooth simple closed curve <img src='http://www.wangsblog.com/jeffrey//pictures/699e902f5598ca623370f833cffb1a57.gif' title=' C ' alt=' C ' align=absmiddle>. Show that the centroid <img src='http://www.wangsblog.com/jeffrey//pictures/f45a08812bddfe9b5913e7fcec7cc85c.gif' title=' (\overline{x}, \overline{y}) ' alt=' (\overline{x}, \overline{y}) ' align=absmiddle> of <img src='http://www.wangsblog.com/jeffrey//pictures/75695b46abca7ce53dfa3b4e984a45ca.gif' title=' R ' alt=' R ' align=absmiddle> is given by</p>
<p><img src='http://www.wangsblog.com/jeffrey//pictures/eedfa9adccd1596cd0bb146ca3ab08a5.gif' title=' \displaystyle \overline{x} = \frac{1}{2A} \oint_C x^2 dy ' alt=' \displaystyle \overline{x} = \frac{1}{2A} \oint_C x^2 dy ' align=absmiddle> and <img src='http://www.wangsblog.com/jeffrey//pictures/84e2d23913056041dd0977834ebda5d7.gif' title=' \displaystyle \overline{y} = -\frac{1}{2A} \oint_C y^2 dx ' alt=' \displaystyle \overline{y} = -\frac{1}{2A} \oint_C y^2 dx ' align=absmiddle>.</p>
]]></content:encoded>
			<wfw:commentRss>http://wangsblog.com/jeffrey/?feed=rss2&amp;p=330</wfw:commentRss>
		<slash:comments>0</slash:comments>
		</item>
		<item>
		<title>More Integrals&#8230; *whine*. Topic: Calculus.</title>
		<link>http://wangsblog.com/jeffrey/?p=329</link>
		<comments>http://wangsblog.com/jeffrey/?p=329#comments</comments>
		<pubDate>Tue, 05 Jun 2007 01:36:11 +0000</pubDate>
		<dc:creator>paladin8</dc:creator>
				<category><![CDATA[Calculus]]></category>

		<guid isPermaLink="false">http://wangsblog.com/jeffrey/?p=329</guid>
		<description><![CDATA[Definition: (Jacobian) If  is the transformation from the -plane to the -plane defined by the equations  and , then the Jacobian of  is denoted by  or by  and is defined by
,
i.e. the determinant of the matrix of the partial derivatives (also known as the Jacobian matrix). Naturally, this can be [...]]]></description>
			<content:encoded><![CDATA[<p><strong>Definition</strong>: (Jacobian) If <img src='http://www.wangsblog.com/jeffrey//pictures/39300282214603edc7b46a69b156bdcb.gif' title=' T ' alt=' T ' align=absmiddle> is the transformation from the <img src='http://www.wangsblog.com/jeffrey//pictures/7d63c29bfc8b07bb0c67b5eb7c726e14.gif' title=' uv ' alt=' uv ' align=absmiddle>-plane to the <img src='http://www.wangsblog.com/jeffrey//pictures/2341d8d15abec0908a302b4a02a24d0e.gif' title=' xy ' alt=' xy ' align=absmiddle>-plane defined by the equations <img src='http://www.wangsblog.com/jeffrey//pictures/a57eaf17862a4c93e3cbebd8483bcb14.gif' title=' x = x(u, v) ' alt=' x = x(u, v) ' align=absmiddle> and <img src='http://www.wangsblog.com/jeffrey//pictures/cdf3cc1ec328aca785bf90387a242558.gif' title=' y = y(u, v) ' alt=' y = y(u, v) ' align=absmiddle>, then the Jacobian of <img src='http://www.wangsblog.com/jeffrey//pictures/39300282214603edc7b46a69b156bdcb.gif' title=' T ' alt=' T ' align=absmiddle> is denoted by <img src='http://www.wangsblog.com/jeffrey//pictures/1c7c753f86050a7a53ed7f755beac57b.gif' title=' J(u, v) ' alt=' J(u, v) ' align=absmiddle> or by <img src='http://www.wangsblog.com/jeffrey//pictures/bc1f2c61234aaa86797b4a1959fdcde5.gif' title=' \partial(x, y)/\partial(u, v) ' alt=' \partial(x, y)/\partial(u, v) ' align=absmiddle> and is defined by</p>
<p><img src='http://www.wangsblog.com/jeffrey//pictures/ed33eb860caacef9fa7dfb62a51c8f02.gif' title=' J(u, v) = \frac{\partial(x, y)}{\partial(u, v)} = \frac{\partial x}{\partial u} \cdot \frac{\partial y}{\partial v} &amp;#8211; \frac{\partial y}{\partial u} \cdot \frac{\partial x}{\partial v} ' alt=' J(u, v) = \frac{\partial(x, y)}{\partial(u, v)} = \frac{\partial x}{\partial u} \cdot \frac{\partial y}{\partial v} &amp;#8211; \frac{\partial y}{\partial u} \cdot \frac{\partial x}{\partial v} ' align=absmiddle>,</p>
<p>i.e. the determinant of the matrix of the partial derivatives (also known as the Jacobian matrix). Naturally, this can be generalized to more variables.</p>
<p>&#8212;&#8212;&#8212;&#8212;&#8212;&#8212;&#8211;</p>
<p><strong>Theorem</strong>: If the transformation <img src='http://www.wangsblog.com/jeffrey//pictures/a57eaf17862a4c93e3cbebd8483bcb14.gif' title=' x = x(u, v) ' alt=' x = x(u, v) ' align=absmiddle>, <img src='http://www.wangsblog.com/jeffrey//pictures/cdf3cc1ec328aca785bf90387a242558.gif' title=' y = y(u, v) ' alt=' y = y(u, v) ' align=absmiddle> maps the region <img src='http://www.wangsblog.com/jeffrey//pictures/b9675b5d61a63b01419a46d9a48e3c7b.gif' title=' S ' alt=' S ' align=absmiddle> in the <img src='http://www.wangsblog.com/jeffrey//pictures/7d63c29bfc8b07bb0c67b5eb7c726e14.gif' title=' uv ' alt=' uv ' align=absmiddle>-plane into the region <img src='http://www.wangsblog.com/jeffrey//pictures/75695b46abca7ce53dfa3b4e984a45ca.gif' title=' R ' alt=' R ' align=absmiddle> in the <img src='http://www.wangsblog.com/jeffrey//pictures/2341d8d15abec0908a302b4a02a24d0e.gif' title=' xy ' alt=' xy ' align=absmiddle>-plane, and if the Jacobian <img src='http://www.wangsblog.com/jeffrey//pictures/bc1f2c61234aaa86797b4a1959fdcde5.gif' title=' \partial(x, y)/\partial(u, v) ' alt=' \partial(x, y)/\partial(u, v) ' align=absmiddle> is nonzero and does not change sign on <img src='http://www.wangsblog.com/jeffrey//pictures/b9675b5d61a63b01419a46d9a48e3c7b.gif' title=' S ' alt=' S ' align=absmiddle>, then (with appropriate restrictions on the transformation and the regions) it follows that</p>
<p><img src='http://www.wangsblog.com/jeffrey//pictures/4071a3800fa76f40b914f6cb502c9924.gif' title=' \displaystyle \int_R \int f(x, y) dA_{xy} = \int_S \int f(x(u, v), y(u, v)) \left|\frac{\partial(x, y)}{\partial(u, v)} \right| dA_{uv} ' alt=' \displaystyle \int_R \int f(x, y) dA_{xy} = \int_S \int f(x(u, v), y(u, v)) \left|\frac{\partial(x, y)}{\partial(u, v)} \right| dA_{uv} ' align=absmiddle>.</p>
<p>&#8212;&#8212;&#8212;&#8212;&#8212;&#8212;&#8211;</p>
<p><strong>Problem</strong>: Evaluate <img src='http://www.wangsblog.com/jeffrey//pictures/a4979020bc5a465e09e74d759e5e88f1.gif' title=' \displaystyle \int_R \int e^{(y-x)/(y+x)} dA ' alt=' \displaystyle \int_R \int e^{(y-x)/(y+x)} dA ' align=absmiddle>, where <img src='http://www.wangsblog.com/jeffrey//pictures/75695b46abca7ce53dfa3b4e984a45ca.gif' title=' R ' alt=' R ' align=absmiddle> is the region in the first quadrant enclosed by the trapezoid with vertices <img src='http://www.wangsblog.com/jeffrey//pictures/cb6aac7de3bc1a351c2b26b33682a4e2.gif' title=' (0, 1); (1, 0); (0, 4); (4, 0) ' alt=' (0, 1); (1, 0); (0, 4); (4, 0) ' align=absmiddle>.</p>
<p><strong>Solution</strong>: The bounding lines can be written as <img src='http://www.wangsblog.com/jeffrey//pictures/e689d4ddb0d1c08b60bdb5df049d4804.gif' title=' x = 0 ' alt=' x = 0 ' align=absmiddle>, <img src='http://www.wangsblog.com/jeffrey//pictures/479ed1e87b885147344b2ed44d2d48d6.gif' title=' y = 0 ' alt=' y = 0 ' align=absmiddle>, <img src='http://www.wangsblog.com/jeffrey//pictures/3d71c8f9e10c68365c080e53116fb97e.gif' title=' y = -x+1 ' alt=' y = -x+1 ' align=absmiddle>, and <img src='http://www.wangsblog.com/jeffrey//pictures/06c4345e76ac83f70d754dab561989f7.gif' title=' y = -x+4 ' alt=' y = -x+4 ' align=absmiddle>. Now consider the transformation <img src='http://www.wangsblog.com/jeffrey//pictures/3ffabe37850c180f36a379d7c422e87c.gif' title=' u = y+x ' alt=' u = y+x ' align=absmiddle> and <img src='http://www.wangsblog.com/jeffrey//pictures/ffb6d74c871b8d8d9cd89bd036289db3.gif' title=' v = y-x ' alt=' v = y-x ' align=absmiddle>. In the <img src='http://www.wangsblog.com/jeffrey//pictures/7d63c29bfc8b07bb0c67b5eb7c726e14.gif' title=' uv ' alt=' uv ' align=absmiddle>-plane, the bounding lines of the new region <img src='http://www.wangsblog.com/jeffrey//pictures/b9675b5d61a63b01419a46d9a48e3c7b.gif' title=' S ' alt=' S ' align=absmiddle> can now be written as <img src='http://www.wangsblog.com/jeffrey//pictures/c0ecc1013ea7e8e15df95f9dc44698e4.gif' title=' u = 1 ' alt=' u = 1 ' align=absmiddle>, <img src='http://www.wangsblog.com/jeffrey//pictures/0d2f499b82d74a019d32fcd6099029c7.gif' title=' u = 4 ' alt=' u = 4 ' align=absmiddle>, <img src='http://www.wangsblog.com/jeffrey//pictures/f96635fb637b1b979af799b35e2e1935.gif' title=' v = u ' alt=' v = u ' align=absmiddle>, and <img src='http://www.wangsblog.com/jeffrey//pictures/a3e32601e5dd6af4c32a46297f7279c0.gif' title=' v = -u ' alt=' v = -u ' align=absmiddle>.</p>
<p>We can write <img src='http://www.wangsblog.com/jeffrey//pictures/6722c218a6f30869ef6886dc4b050a37.gif' title=' x ' alt=' x ' align=absmiddle> and <img src='http://www.wangsblog.com/jeffrey//pictures/ee23c4f89cdc7b5be951059c2435fa2d.gif' title=' y ' alt=' y ' align=absmiddle> as functions of <img src='http://www.wangsblog.com/jeffrey//pictures/bc2c9f94d4206d09fef01fa877f5b516.gif' title=' u ' alt=' u ' align=absmiddle> and <img src='http://www.wangsblog.com/jeffrey//pictures/c4f4b9b0ab0a2eb771bf7decb3a53c8d.gif' title=' v ' alt=' v ' align=absmiddle>: simply <img src='http://www.wangsblog.com/jeffrey//pictures/dfcab720ff296ff7f4afc2f5a539598f.gif' title=' x = \frac{u-v}{2} ' alt=' x = \frac{u-v}{2} ' align=absmiddle> and <img src='http://www.wangsblog.com/jeffrey//pictures/d91b46ca4950786d6a76c410912f24ac.gif' title=' y = \frac{u+v}{2} ' alt=' y = \frac{u+v}{2} ' align=absmiddle>. So the Jacobian <img src='http://www.wangsblog.com/jeffrey//pictures/6fb2c2bfe1a772a002d51c893312d1ef.gif' title=' \displaystyle \frac{\partial(x, y)}{\partial(u, v)} = \frac{\partial x}{\partial u} \cdot \frac{\partial y}{\partial v} &amp;#8211; \frac{\partial y}{\partial u} \cdot \frac{\partial x}{\partial v} = \frac{1}{2} \cdot \frac{1}{2} &amp;#8211; \frac{1}{2} \cdot \left(-\frac{1}{2} \right) = \frac{1}{2} ' alt=' \displaystyle \frac{\partial(x, y)}{\partial(u, v)} = \frac{\partial x}{\partial u} \cdot \frac{\partial y}{\partial v} &amp;#8211; \frac{\partial y}{\partial u} \cdot \frac{\partial x}{\partial v} = \frac{1}{2} \cdot \frac{1}{2} &amp;#8211; \frac{1}{2} \cdot \left(-\frac{1}{2} \right) = \frac{1}{2} ' align=absmiddle>.</p>
<p>Then our original integral becomes <img src='http://www.wangsblog.com/jeffrey//pictures/dd93d421496d3974a4f5f7faefc3a8b6.gif' title=' \displaystyle \int_R \int e^{(y-x)/(y+x)} dA = \frac{1}{2} \int_S \int e^{v/u} dA ' alt=' \displaystyle \int_R \int e^{(y-x)/(y+x)} dA = \frac{1}{2} \int_S \int e^{v/u} dA ' align=absmiddle>. And this is equivalent to</p>
<p><img src='http://www.wangsblog.com/jeffrey//pictures/e5ffc613c9ce441f0774195e29fdc1e9.gif' title=' \displaystyle \frac{1}{2} \int_S \int e^{v/u} dA = \frac{1}{2} \int_1^4 \int_{-u}^u e^{v/u} dv du = \frac{1}{2} \int_1^4 \big[ u e^{v/u} \big]_{v=-u}^u du = \frac{1}{2} \int_1^4 u\left(e-\frac{1}{e}\right) du = \frac{15}{4}\left(e-\frac{1}{e}\right) ' alt=' \displaystyle \frac{1}{2} \int_S \int e^{v/u} dA = \frac{1}{2} \int_1^4 \int_{-u}^u e^{v/u} dv du = \frac{1}{2} \int_1^4 \big[ u e^{v/u} \big]_{v=-u}^u du = \frac{1}{2} \int_1^4 u\left(e-\frac{1}{e}\right) du = \frac{15}{4}\left(e-\frac{1}{e}\right) ' align=absmiddle>.</p>
<p>QED.</p>
<p>&#8212;&#8212;&#8212;&#8212;&#8212;&#8212;&#8211;</p>
<p>Comment: Note that the above theorem is probably very important in multivariable calculus, as it is the equivalent to <img src='http://www.wangsblog.com/jeffrey//pictures/bc2c9f94d4206d09fef01fa877f5b516.gif' title=' u ' alt=' u ' align=absmiddle>-substitution in one variable, which we all know is the ultimate integration technique. It functions in the same way, giving you a lot more flexibility on the function you are integrating and the region you are integrating on.</p>
<p>&#8212;&#8212;&#8212;&#8212;&#8212;&#8212;&#8211;</p>
<p>Practice Problem: Evaluate <img src='http://www.wangsblog.com/jeffrey//pictures/e22959fc847a929b05ef49b2781d4315.gif' title=' \displaystyle \int_R \int (x^2-y^2) dA ' alt=' \displaystyle \int_R \int (x^2-y^2) dA ' align=absmiddle>, where <img src='http://www.wangsblog.com/jeffrey//pictures/75695b46abca7ce53dfa3b4e984a45ca.gif' title=' R ' alt=' R ' align=absmiddle> is the rectangular region enclosed by the lines <img src='http://www.wangsblog.com/jeffrey//pictures/6c6ce383b972a67348be9b60c1179d29.gif' title=' y = -x ' alt=' y = -x ' align=absmiddle>, <img src='http://www.wangsblog.com/jeffrey//pictures/30802bd769e3df2073d039c723afabc1.gif' title=' y = 1-x ' alt=' y = 1-x ' align=absmiddle>, <img src='http://www.wangsblog.com/jeffrey//pictures/9421ffefd601c2ceb478d81421a07297.gif' title=' y = x ' alt=' y = x ' align=absmiddle>, <img src='http://www.wangsblog.com/jeffrey//pictures/c49ff99c0bb9e7210134496a6e8a476e.gif' title=' y = x+2 ' alt=' y = x+2 ' align=absmiddle>.</p>
]]></content:encoded>
			<wfw:commentRss>http://wangsblog.com/jeffrey/?feed=rss2&amp;p=329</wfw:commentRss>
		<slash:comments>0</slash:comments>
		</item>
		<item>
		<title>One By One, We&#8217;re Making It Fun. Topic: Calculus/S&amp;S.</title>
		<link>http://wangsblog.com/jeffrey/?p=328</link>
		<comments>http://wangsblog.com/jeffrey/?p=328#comments</comments>
		<pubDate>Mon, 28 May 2007 18:34:38 +0000</pubDate>
		<dc:creator>paladin8</dc:creator>
				<category><![CDATA[Calculus]]></category>
		<category><![CDATA[Sequences &#038; Series]]></category>

		<guid isPermaLink="false">http://wangsblog.com/jeffrey/?p=328</guid>
		<description><![CDATA[Theorem: (Stolz-Cesaro) Let  and  be sequences of real numbers such that  is positive, strictly increasing, and unbounded. If the limit

exists, then the following limit also exists and we have
.
&#8212;&#8212;&#8212;&#8212;&#8212;&#8212;&#8211;
Theorem: (Summation by Parts) If  and  are two sequences, then
.
&#8212;&#8212;&#8212;&#8212;&#8212;&#8212;&#8211;
Problem: Let  be a sequence of real numbers such that  converges. [...]]]></description>
			<content:encoded><![CDATA[<p><strong>Theorem</strong>: (Stolz-Cesaro) Let <img src='http://www.wangsblog.com/jeffrey//pictures/5b5aa2a8cb5f10be03d54bcf4c436297.gif' title=' \{a_n\} ' alt=' \{a_n\} ' align=absmiddle> and <img src='http://www.wangsblog.com/jeffrey//pictures/bd7a1d05c1ea245f609326f91eb515e9.gif' title=' \{b_n\} ' alt=' \{b_n\} ' align=absmiddle> be sequences of real numbers such that <img src='http://www.wangsblog.com/jeffrey//pictures/bd7a1d05c1ea245f609326f91eb515e9.gif' title=' \{b_n\} ' alt=' \{b_n\} ' align=absmiddle> is positive, strictly increasing, and unbounded. If the limit</p>
<p><img src='http://www.wangsblog.com/jeffrey//pictures/19313e8681049445852d48e8617b6f81.gif' title=' \displaystyle \lim_{n \rightarrow \infty} \frac{a_{n+1}-a_n}{b_{n+1}-b_n} = L ' alt=' \displaystyle \lim_{n \rightarrow \infty} \frac{a_{n+1}-a_n}{b_{n+1}-b_n} = L ' align=absmiddle></p>
<p>exists, then the following limit also exists and we have</p>
<p><img src='http://www.wangsblog.com/jeffrey//pictures/7563795de12e4a53aeab675cf568bfec.gif' title=' \displaystyle \lim_{n \rightarrow \infty} \frac{a_n}{b_n} = L ' alt=' \displaystyle \lim_{n \rightarrow \infty} \frac{a_n}{b_n} = L ' align=absmiddle>.</p>
<p>&#8212;&#8212;&#8212;&#8212;&#8212;&#8212;&#8211;</p>
<p><strong>Theorem</strong>: (Summation by Parts) If <img src='http://www.wangsblog.com/jeffrey//pictures/722562fb6c75b13cd0f93ace56f0a8c9.gif' title=' \{f_k\} ' alt=' \{f_k\} ' align=absmiddle> and <img src='http://www.wangsblog.com/jeffrey//pictures/73da3d781233d7f8ca0b826446459bb4.gif' title=' \{g_k\} ' alt=' \{g_k\} ' align=absmiddle> are two sequences, then</p>
<p><img src='http://www.wangsblog.com/jeffrey//pictures/5c7d5290a9be2723056ee20f8ba62f30.gif' title=' \displaystyle \sum_{k=m}^n f_k(g_{k+1}-g_k) = [f_{n+1}g_{n+1}-f_mg_m]-\sum_{k=m}^n g_{k+1}(f_{k+1}-f_k) ' alt=' \displaystyle \sum_{k=m}^n f_k(g_{k+1}-g_k) = [f_{n+1}g_{n+1}-f_mg_m]-\sum_{k=m}^n g_{k+1}(f_{k+1}-f_k) ' align=absmiddle>.</p>
<p>&#8212;&#8212;&#8212;&#8212;&#8212;&#8212;&#8211;</p>
<p><strong>Problem</strong>: Let <img src='http://www.wangsblog.com/jeffrey//pictures/5b5aa2a8cb5f10be03d54bcf4c436297.gif' title=' \{a_n\} ' alt=' \{a_n\} ' align=absmiddle> be a sequence of real numbers such that <img src='http://www.wangsblog.com/jeffrey//pictures/a584d5a20f0149140d7c561f477ebe13.gif' title=' \displaystyle \sum_{k=0}^{\infty} a_k ' alt=' \displaystyle \sum_{k=0}^{\infty} a_k ' align=absmiddle> converges. Show that <img src='http://www.wangsblog.com/jeffrey//pictures/df82fc5a3f00b04c085d3e9cc1973f3d.gif' title=' \displaystyle \lim_{n \rightarrow \infty} \frac{1}{n+1} \sum_{k=0}^n k \cdot a_k = 0 ' alt=' \displaystyle \lim_{n \rightarrow \infty} \frac{1}{n+1} \sum_{k=0}^n k \cdot a_k = 0 ' align=absmiddle>.</p>
<p><strong>Solution</strong>: Define the sequence <img src='http://www.wangsblog.com/jeffrey//pictures/bd7a1d05c1ea245f609326f91eb515e9.gif' title=' \{b_n\} ' alt=' \{b_n\} ' align=absmiddle> by <img src='http://www.wangsblog.com/jeffrey//pictures/ff45ee143caedabfe8a4c7319b06eb40.gif' title=' \displaystyle b_n = \sum_{k=0}^n a_k ' alt=' \displaystyle b_n = \sum_{k=0}^n a_k ' align=absmiddle> and let <img src='http://www.wangsblog.com/jeffrey//pictures/79afa3a2130ffa201e9c4c9cb343fd40.gif' title=' \displaystyle \lim_{n \rightarrow \infty} b_n = L ' alt=' \displaystyle \lim_{n \rightarrow \infty} b_n = L ' align=absmiddle>. Then, by summation by parts with <img src='http://www.wangsblog.com/jeffrey//pictures/e949d2667fc3604ee967ddeefe3b8893.gif' title=' \{f_n\} = \{n\} ' alt=' \{f_n\} = \{n\} ' align=absmiddle> and <img src='http://www.wangsblog.com/jeffrey//pictures/ea41fcf03a6a2a3d47585303641cbbaf.gif' title=' \{g_n\} = \{b_n\} ' alt=' \{g_n\} = \{b_n\} ' align=absmiddle>, we have</p>
<p><img src='http://www.wangsblog.com/jeffrey//pictures/570116b7a23943ff64f47568513fd166.gif' title=' \displaystyle \sum_{k=0}^n k \cdot a_k = \sum_{k=0}^n k \cdot (b_{k+1}-b_k) = (n+1)b_{n+1}-\sum_{k=0}^n b_{k+1} ' alt=' \displaystyle \sum_{k=0}^n k \cdot a_k = \sum_{k=0}^n k \cdot (b_{k+1}-b_k) = (n+1)b_{n+1}-\sum_{k=0}^n b_{k+1} ' align=absmiddle>.</p>
<p>The summation we wish to take the limit of is then</p>
<p><img src='http://www.wangsblog.com/jeffrey//pictures/9383fda849c455a3279f77ca04ed3fc3.gif' title=' \displaystyle \frac{1}{n+1} \sum_{k=0}^n k \cdot a_k = b_{n+1}-\frac{1}{n+1} \sum_{k=0}^n b_{k+1} ' alt=' \displaystyle \frac{1}{n+1} \sum_{k=0}^n k \cdot a_k = b_{n+1}-\frac{1}{n+1} \sum_{k=0}^n b_{k+1} ' align=absmiddle>.</p>
<p>But since, by Stolz-Cesaro,</p>
<p><img src='http://www.wangsblog.com/jeffrey//pictures/dce43f72b3672666c8b9bc77b298ebe6.gif' title=' \displaystyle \lim_{n \rightarrow \infty} \frac{1}{n+1} \sum_{k=0}^n b_{k+1} = \lim_{n \rightarrow \infty} b_{n+1} = L ' alt=' \displaystyle \lim_{n \rightarrow \infty} \frac{1}{n+1} \sum_{k=0}^n b_{k+1} = \lim_{n \rightarrow \infty} b_{n+1} = L ' align=absmiddle>,</p>
<p>we obtain</p>
<p><img src='http://www.wangsblog.com/jeffrey//pictures/e8bd3e9a8c533ce9d7b3c625218476fa.gif' title=' \displaystyle \lim_{n \rightarrow \infty} \frac{1}{n+1} \sum_{k=0}^n k \cdot a_k = \lim_{n \rightarrow \infty} \left(b_n-\frac{1}{n+1} \sum_{k=0}^n b_{k+1}\right) = L-L = 0 ' alt=' \displaystyle \lim_{n \rightarrow \infty} \frac{1}{n+1} \sum_{k=0}^n k \cdot a_k = \lim_{n \rightarrow \infty} \left(b_n-\frac{1}{n+1} \sum_{k=0}^n b_{k+1}\right) = L-L = 0 ' align=absmiddle>.</p>
<p>QED.</p>
<p>&#8212;&#8212;&#8212;&#8212;&#8212;&#8212;&#8211;</p>
<p>Comment: Summation by parts is a very useful technique to change sums around so that they are easier to evaluate. If you hadn&#8217;t noticed, it is the discrete analogue of integration by parts and is in fact very similar. Stolz-Cesaro is powerful as well and seems like a discrete analogue to L&#8217;Hopital (but I&#8217;m not sure about this one). Applying well-known calculus ideas to discrete things can often turn into neat results.</p>
<p>&#8212;&#8212;&#8212;&#8212;&#8212;&#8212;&#8211;</p>
<p>Practice Problem: If <img src='http://www.wangsblog.com/jeffrey//pictures/5b5aa2a8cb5f10be03d54bcf4c436297.gif' title=' \{a_n\} ' alt=' \{a_n\} ' align=absmiddle> is a decreasing sequence such that <img src='http://www.wangsblog.com/jeffrey//pictures/d92e409635193c3a32bb6a407edc3ad8.gif' title=' \displaystyle \lim_{n \rightarrow \infty} a_n = 0 ' alt=' \displaystyle \lim_{n \rightarrow \infty} a_n = 0 ' align=absmiddle>, show that</p>
<p><img src='http://www.wangsblog.com/jeffrey//pictures/488e7b1bdf1e9fdfd4e3999a27c848f3.gif' title=' \displaystyle \sum_{k=1}^{\infty} a_k \cdot \sin{(kx)} ' alt=' \displaystyle \sum_{k=1}^{\infty} a_k \cdot \sin{(kx)} ' align=absmiddle></p>
<p>converges for all <img src='http://www.wangsblog.com/jeffrey//pictures/6722c218a6f30869ef6886dc4b050a37.gif' title=' x ' alt=' x ' align=absmiddle>.</p>
]]></content:encoded>
			<wfw:commentRss>http://wangsblog.com/jeffrey/?feed=rss2&amp;p=328</wfw:commentRss>
		<slash:comments>0</slash:comments>
		</item>
		<item>
		<title>Bigger Means Better. Topic: Algebra/Inequalities/Sets. Level: Olympiad.</title>
		<link>http://wangsblog.com/jeffrey/?p=327</link>
		<comments>http://wangsblog.com/jeffrey/?p=327#comments</comments>
		<pubDate>Sun, 13 May 2007 01:55:01 +0000</pubDate>
		<dc:creator>paladin8</dc:creator>
				<category><![CDATA[Algebra]]></category>
		<category><![CDATA[Inequalities]]></category>
		<category><![CDATA[Olympiad]]></category>
		<category><![CDATA[Sets]]></category>

		<guid isPermaLink="false">http://wangsblog.com/jeffrey/?p=327</guid>
		<description><![CDATA[Definition: A set  is said to be convex if, for any  and , we have .
&#8212;&#8212;&#8212;&#8212;&#8212;&#8212;&#8211;
Definition: Let  be a real-valued function defined on a convex set . We say that  is convex if, for all  and , we have
.
&#8212;&#8212;&#8212;&#8212;&#8212;&#8212;&#8211;
Theorem: (Jensen&#8217;s Inequality) Let  be a real-valued, convex function defined on [...]]]></description>
			<content:encoded><![CDATA[<p><strong>Definition</strong>: A set <img src='http://www.wangsblog.com/jeffrey//pictures/b9675b5d61a63b01419a46d9a48e3c7b.gif' title=' S ' alt=' S ' align=absmiddle> is said to be <em>convex</em> if, for any <img src='http://www.wangsblog.com/jeffrey//pictures/cd1769d0f399535e1f64a3ae1f1cb8da.gif' title=' X, Y \in S ' alt=' X, Y \in S ' align=absmiddle> and <img src='http://www.wangsblog.com/jeffrey//pictures/aa92c2607bfb75f77a16af7c23fced7c.gif' title=' \lambda \in [0,1] ' alt=' \lambda \in [0,1] ' align=absmiddle>, we have <img src='http://www.wangsblog.com/jeffrey//pictures/d0ea4d391a4253e5f094cc02c671566b.gif' title=' \lambda X+(1-\lambda)Y \in S ' alt=' \lambda X+(1-\lambda)Y \in S ' align=absmiddle>.</p>
<p>&#8212;&#8212;&#8212;&#8212;&#8212;&#8212;&#8211;</p>
<p><strong>Definition</strong>: Let <img src='http://www.wangsblog.com/jeffrey//pictures/40bc56b5b134c6a6e2a2c6e54611877b.gif' title=' f: D \rightarrow \mathbb{R} ' alt=' f: D \rightarrow \mathbb{R} ' align=absmiddle> be a real-valued function defined on a convex set <img src='http://www.wangsblog.com/jeffrey//pictures/2f4396bab5869c1e0c9f8a7620bf2518.gif' title=' D ' alt=' D ' align=absmiddle>. We say that <img src='http://www.wangsblog.com/jeffrey//pictures/ce40937fdfbd06b8a15244e102a09356.gif' title=' f ' alt=' f ' align=absmiddle> is <em>convex</em> if, for all <img src='http://www.wangsblog.com/jeffrey//pictures/4368740990172e10de4635c8f7473023.gif' title=' X, Y \in D ' alt=' X, Y \in D ' align=absmiddle> and <img src='http://www.wangsblog.com/jeffrey//pictures/1383cb4f657c2663137af17a8d3fe721.gif' title=' \lambda \in [0, 1] ' alt=' \lambda \in [0, 1] ' align=absmiddle>, we have</p>
<p><img src='http://www.wangsblog.com/jeffrey//pictures/d47683d94269a4cb66b054f7dfa5064c.gif' title=' \lambda f(X) + (1-\lambda) f(Y) \ge f(\lambda X+(1-\lambda)Y) ' alt=' \lambda f(X) + (1-\lambda) f(Y) \ge f(\lambda X+(1-\lambda)Y) ' align=absmiddle>.</p>
<p>&#8212;&#8212;&#8212;&#8212;&#8212;&#8212;&#8211;</p>
<p><strong>Theorem</strong>: (Jensen&#8217;s Inequality) Let <img src='http://www.wangsblog.com/jeffrey//pictures/ce40937fdfbd06b8a15244e102a09356.gif' title=' f ' alt=' f ' align=absmiddle> be a real-valued, convex function defined on a convex set <img src='http://www.wangsblog.com/jeffrey//pictures/2f4396bab5869c1e0c9f8a7620bf2518.gif' title=' D ' alt=' D ' align=absmiddle>. Given <img src='http://www.wangsblog.com/jeffrey//pictures/cc988c2163e10f9a382009b9ab85604b.gif' title=' x_1, x_2, \ldots, x_n \in D ' alt=' x_1, x_2, \ldots, x_n \in D ' align=absmiddle> and nonnegative reals <img src='http://www.wangsblog.com/jeffrey//pictures/f64e21b203022438e3028ba7e3dc8e17.gif' title=' w_1, w_2, \ldots, w_n ' alt=' w_1, w_2, \ldots, w_n ' align=absmiddle> such that <img src='http://www.wangsblog.com/jeffrey//pictures/604ae197fbee13f065886396c967c78c.gif' title=' w_1+w_2+\cdots+w_n = 1 ' alt=' w_1+w_2+\cdots+w_n = 1 ' align=absmiddle>, we have</p>
<p><img src='http://www.wangsblog.com/jeffrey//pictures/83859f8d97eb98b6725b3678adc211ce.gif' title=' w_1f(x_1)+w_2f(x_2)+\cdots+w_nf(x_n) \ge f(w_1x_1+w_2x_2+\cdots+w_nx_n) ' alt=' w_1f(x_1)+w_2f(x_2)+\cdots+w_nf(x_n) \ge f(w_1x_1+w_2x_2+\cdots+w_nx_n) ' align=absmiddle></p>
<p>or, more concisely,</p>
<p><img src='http://www.wangsblog.com/jeffrey//pictures/97972dc6dc24bd5f3f8400b11b68a812.gif' title=' \displaystyle \sum_{i=1}^n w_if(x_i) \ge f\left(\sum_{i=1}^n w_ix_i\right) ' alt=' \displaystyle \sum_{i=1}^n w_if(x_i) \ge f\left(\sum_{i=1}^n w_ix_i\right) ' align=absmiddle>.</p>
<p><strong>Proof</strong>: We proceed by induction on <img src='http://www.wangsblog.com/jeffrey//pictures/bfbdd7d089006253c9a32f7c78c15270.gif' title=' n ' alt=' n ' align=absmiddle>.</p>
<p>Base Case &#8211; <img src='http://www.wangsblog.com/jeffrey//pictures/27c3e74b761d002ff70b040db97a607d.gif' title=' n = 1 ' alt=' n = 1 ' align=absmiddle>, we have <img src='http://www.wangsblog.com/jeffrey//pictures/fac72d054d5ef22f07d3ec5064916c12.gif' title=' f(x_1) \ge f(x_1) ' alt=' f(x_1) \ge f(x_1) ' align=absmiddle> which is trivially true. <img src='http://www.wangsblog.com/jeffrey//pictures/83e17fa6b2d07a6666380146b302bee1.gif' title=' n = 2 ' alt=' n = 2 ' align=absmiddle>, we have <img src='http://www.wangsblog.com/jeffrey//pictures/1262f0e572dfd56c4f417ad0fc88aa4f.gif' title=' \lambda_1 f(x_1)+\lambda_2 f(x_2) \ge f(\lambda_1 x_1+\lambda_2 x_2) ' alt=' \lambda_1 f(x_1)+\lambda_2 f(x_2) \ge f(\lambda_1 x_1+\lambda_2 x_2) ' align=absmiddle> with <img src='http://www.wangsblog.com/jeffrey//pictures/572f90e63366d9856513e3421074773a.gif' title=' \lambda_1+\lambda_2 = 1' alt=' \lambda_1+\lambda_2 = 1' align=absmiddle> which is true by the definition of convexity.</p>
<p>Induction Step &#8211; Suppose <img src='http://www.wangsblog.com/jeffrey//pictures/077664638d46872ef4cd6d5eadbebc99.gif' title=' \displaystyle \sum_{i=1}^k w_if(x_i) \ge f\left(\sum_{i=1}^k w_ix_i\right) ' alt=' \displaystyle \sum_{i=1}^k w_if(x_i) \ge f\left(\sum_{i=1}^k w_ix_i\right) ' align=absmiddle>. We wish to show</p>
<p><img src='http://www.wangsblog.com/jeffrey//pictures/8d7119d9bdce69c3c38d0b8544decd3e.gif' title=' \displaystyle \sum_{i=1}^{k+1} u_if(x_i) \ge f\left(\sum_{i=1}^{k+1} u_ix_i\right) ' alt=' \displaystyle \sum_{i=1}^{k+1} u_if(x_i) \ge f\left(\sum_{i=1}^{k+1} u_ix_i\right) ' align=absmiddle></p>
<p>for nonnegative reals <img src='http://www.wangsblog.com/jeffrey//pictures/a0af4f2bcb0316761db1d11c8b771710.gif' title=' u_1, u_2, \ldots, u_{k+1} ' alt=' u_1, u_2, \ldots, u_{k+1} ' align=absmiddle> that sum to <img src='http://www.wangsblog.com/jeffrey//pictures/35c0804c862a635e2fe8371dc43e25d0.gif' title=' 1 ' alt=' 1 ' align=absmiddle>. Well,</p>
<p><img src='http://www.wangsblog.com/jeffrey//pictures/9fe0c6bc5ab3dae011be1c523ad6e59a.gif' title=' \displaystyle f\left(\sum_{i=1}^{k+1} u_ix_i\right) \le u_{k+1}f(x_{k+1})+(1-u_{k+1})f\left(\frac{1}{1-u_{k+1}} \sum_{i=1}^k u_ix_i\right) ' alt=' \displaystyle f\left(\sum_{i=1}^{k+1} u_ix_i\right) \le u_{k+1}f(x_{k+1})+(1-u_{k+1})f\left(\frac{1}{1-u_{k+1}} \sum_{i=1}^k u_ix_i\right) ' align=absmiddle></p>
<p>by the definition of convexity. But since <img src='http://www.wangsblog.com/jeffrey//pictures/b659a111616f24d82382bcfd28221643.gif' title=' \displaystyle \frac{1}{1-u_{k+1}} \sum_{i=1}^k u_i = \sum_{i=1}^k \frac{u_i}{1-u_{k+1}} = 1 ' alt=' \displaystyle \frac{1}{1-u_{k+1}} \sum_{i=1}^k u_i = \sum_{i=1}^k \frac{u_i}{1-u_{k+1}} = 1 ' align=absmiddle>, by our inductive hypothesis (the <img src='http://www.wangsblog.com/jeffrey//pictures/365c0b3ff8a6fa58b7ae709949b55608.gif' title=' k ' alt=' k ' align=absmiddle> term case) we must have</p>
<p><img src='http://www.wangsblog.com/jeffrey//pictures/4b551dc9afe64ce09eafe11ccfe5d404.gif' title=' \displaystyle f\left(\frac{1}{1-u_{k+1}} \sum_{i=1}^k u_ix_i\right) \le \sum_{i=1}^k \frac{u_i}{1-u_{k+1}} f(x_i) ' alt=' \displaystyle f\left(\frac{1}{1-u_{k+1}} \sum_{i=1}^k u_ix_i\right) \le \sum_{i=1}^k \frac{u_i}{1-u_{k+1}} f(x_i) ' align=absmiddle>.</p>
<p>Combining this into the above inequality, we obtain</p>
<p><img src='http://www.wangsblog.com/jeffrey//pictures/5da4e3c361d58536528a95d7e8800557.gif' title=' \displaystyle f\left(\sum_{i=1}^{k+1} u_ix_i\right) \le u_{k+1}f(x_{k+1})+(1-u_{k+1})\sum_{i=1}^k \frac{u_i}{1-u_{k+1}} f(x_i) = \sum_{i=1}^{k+1} u_if(x_i) ' alt=' \displaystyle f\left(\sum_{i=1}^{k+1} u_ix_i\right) \le u_{k+1}f(x_{k+1})+(1-u_{k+1})\sum_{i=1}^k \frac{u_i}{1-u_{k+1}} f(x_i) = \sum_{i=1}^{k+1} u_if(x_i) ' align=absmiddle></p>
<p>as desired, completing the induction. QED.</p>
<p>&#8212;&#8212;&#8212;&#8212;&#8212;&#8212;&#8211;</p>
<p><strong>Definition</strong>: The <em>convex hull</em> of a set <img src='http://www.wangsblog.com/jeffrey//pictures/b9675b5d61a63b01419a46d9a48e3c7b.gif' title=' S ' alt=' S ' align=absmiddle> is defined to be the smallest convex set containing <img src='http://www.wangsblog.com/jeffrey//pictures/b9675b5d61a63b01419a46d9a48e3c7b.gif' title=' S ' alt=' S ' align=absmiddle>. It is the set of all points that can be written as <img src='http://www.wangsblog.com/jeffrey//pictures/abcbcaacfa2c0ecd73270e7334ef5054.gif' title=' \displaystyle \sum_{i=1}^n a_ix_i ' alt=' \displaystyle \sum_{i=1}^n a_ix_i ' align=absmiddle>, where <img src='http://www.wangsblog.com/jeffrey//pictures/bfbdd7d089006253c9a32f7c78c15270.gif' title=' n ' alt=' n ' align=absmiddle> is a natural number, <img src='http://www.wangsblog.com/jeffrey//pictures/8dedebd3e856a03087a40ed6ce31e026.gif' title=' a_1, a_2, \ldots, a_n ' alt=' a_1, a_2, \ldots, a_n ' align=absmiddle> are nonnegative reals that sum to <img src='http://www.wangsblog.com/jeffrey//pictures/35c0804c862a635e2fe8371dc43e25d0.gif' title=' 1 ' alt=' 1 ' align=absmiddle>, and <img src='http://www.wangsblog.com/jeffrey//pictures/e930d365a3380d7151442d25c4c00a01.gif' title=' x_1, x_2, \ldots, x_n \in S ' alt=' x_1, x_2, \ldots, x_n \in S ' align=absmiddle>.</p>
<p>&#8212;&#8212;&#8212;&#8212;&#8212;&#8212;&#8211;</p>
<p><strong>Definition</strong>: A <em>vertex</em> of a set <img src='http://www.wangsblog.com/jeffrey//pictures/2f4396bab5869c1e0c9f8a7620bf2518.gif' title=' D ' alt=' D ' align=absmiddle> is a point <img src='http://www.wangsblog.com/jeffrey//pictures/ba04fd8adfd221875589b2922ef3ffcb.gif' title=' v \in D ' alt=' v \in D ' align=absmiddle> such that for all <img src='http://www.wangsblog.com/jeffrey//pictures/d7aae3c674ec344e4c5581d3d580e0ef.gif' title=' x \in D ' alt=' x \in D ' align=absmiddle> not equal to <img src='http://www.wangsblog.com/jeffrey//pictures/c4f4b9b0ab0a2eb771bf7decb3a53c8d.gif' title=' v ' alt=' v ' align=absmiddle> and <img src='http://www.wangsblog.com/jeffrey//pictures/8cb1d4aecc516b0b09877cbe868ebbcc.gif' title=' \lambda &gt; 1 ' alt=' \lambda &gt; 1 ' align=absmiddle> we have <img src='http://www.wangsblog.com/jeffrey//pictures/9b21fb5a6acab1b7fa67e900c5ad8ae9.gif' title=' \lambda v+(1-\lambda)x \not\in D ' alt=' \lambda v+(1-\lambda)x \not\in D ' align=absmiddle>.</p>
<p>&#8212;&#8212;&#8212;&#8212;&#8212;&#8212;&#8211;</p>
<p><strong>Theorem</strong>: Let <img src='http://www.wangsblog.com/jeffrey//pictures/d6024010775ec5820e03a25806860b19.gif' title=' V_D = \{v_1, v_2, \ldots, v_k\} ' alt=' V_D = \{v_1, v_2, \ldots, v_k\} ' align=absmiddle> be the set of vertices of a compact convex set <img src='http://www.wangsblog.com/jeffrey//pictures/2f4396bab5869c1e0c9f8a7620bf2518.gif' title=' D ' alt=' D ' align=absmiddle>. The convex hull of <img src='http://www.wangsblog.com/jeffrey//pictures/74fd22ca4f9458f85c2c2ec87727a102.gif' title=' V_D ' alt=' V_D ' align=absmiddle> is <img src='http://www.wangsblog.com/jeffrey//pictures/2f4396bab5869c1e0c9f8a7620bf2518.gif' title=' D ' alt=' D ' align=absmiddle>.</p>
<p><strong>Example</strong>: The set of vertices of a convex polygon is, in fact, its vertices as traditionally defined in geometry. Any point inside the polygon can be written as a convex combination of its vertices.</p>
<p>&#8212;&#8212;&#8212;&#8212;&#8212;&#8212;&#8211;</p>
<p><strong>Theorem</strong>: If <img src='http://www.wangsblog.com/jeffrey//pictures/ce40937fdfbd06b8a15244e102a09356.gif' title=' f ' alt=' f ' align=absmiddle> is a real-valued, convex function defined on a compact convex set <img src='http://www.wangsblog.com/jeffrey//pictures/2f4396bab5869c1e0c9f8a7620bf2518.gif' title=' D ' alt=' D ' align=absmiddle>, then <img src='http://www.wangsblog.com/jeffrey//pictures/b24695d0cf2355be6999103dcc2e1428.gif' title=' \displaystyle \max_{x \in D} f(x) = \max_{x \in V_D} f(x) ' alt=' \displaystyle \max_{x \in D} f(x) = \max_{x \in V_D} f(x) ' align=absmiddle>, where <img src='http://www.wangsblog.com/jeffrey//pictures/74fd22ca4f9458f85c2c2ec87727a102.gif' title=' V_D ' alt=' V_D ' align=absmiddle> is the set of vertices of <img src='http://www.wangsblog.com/jeffrey//pictures/2f4396bab5869c1e0c9f8a7620bf2518.gif' title=' D ' alt=' D ' align=absmiddle>.</p>
<p><strong>Proof</strong>: We will show that <img src='http://www.wangsblog.com/jeffrey//pictures/5f7bc5ff9cd75e5d60f7f92e9174c6a8.gif' title=' \displaystyle f(x) \le \max_{x \in V_D} f(x) ' alt=' \displaystyle f(x) \le \max_{x \in V_D} f(x) ' align=absmiddle> for all <img src='http://www.wangsblog.com/jeffrey//pictures/d7aae3c674ec344e4c5581d3d580e0ef.gif' title=' x \in D ' alt=' x \in D ' align=absmiddle>.</p>
<p>Let <img src='http://www.wangsblog.com/jeffrey//pictures/1c0ba1cc5cc83b6bf97ebac7732f7708.gif' title=' \displaystyle x = \sum_{i=1}^k \lambda_i v_i ' alt=' \displaystyle x = \sum_{i=1}^k \lambda_i v_i ' align=absmiddle> where <img src='http://www.wangsblog.com/jeffrey//pictures/8806bbef73f38c9d551d5d9ffc9cc81b.gif' title=' \lambda_1, \lambda_2, \ldots, \lambda_k ' alt=' \lambda_1, \lambda_2, \ldots, \lambda_k ' align=absmiddle> are nonnegative reals that sum to <img src='http://www.wangsblog.com/jeffrey//pictures/35c0804c862a635e2fe8371dc43e25d0.gif' title=' 1 ' alt=' 1 ' align=absmiddle> and <img src='http://www.wangsblog.com/jeffrey//pictures/d6024010775ec5820e03a25806860b19.gif' title=' V_D = \{v_1, v_2, \ldots, v_k\} ' alt=' V_D = \{v_1, v_2, \ldots, v_k\} ' align=absmiddle>. This is possible by the preceding theorem. Then by Jensen&#8217;s Inequality we get</p>
<p><img src='http://www.wangsblog.com/jeffrey//pictures/17876dddad4f8235b841bd23221c7fae.gif' title=' \displaystyle \sum_{i=1}^k \lambda_i f(v_i) \ge f\left(\sum_{i=1}^k \lambda_i v_i\right) = f(x) ' alt=' \displaystyle \sum_{i=1}^k \lambda_i f(v_i) \ge f\left(\sum_{i=1}^k \lambda_i v_i\right) = f(x) ' align=absmiddle>.</p>
<p>And since <img src='http://www.wangsblog.com/jeffrey//pictures/5980803dc14e9d259c763d768d281ae9.gif' title=' \displaystyle \max_{x \in V_D} f(x) \ge \sum_{i=1}^k \lambda_i f(v_i) ' alt=' \displaystyle \max_{x \in V_D} f(x) \ge \sum_{i=1}^k \lambda_i f(v_i) ' align=absmiddle>, we have <img src='http://www.wangsblog.com/jeffrey//pictures/6613ed0017e718a69f07a24892965a28.gif' title=' \displaystyle \max_{x \in V_D} f(x) \ge f(x) ' alt=' \displaystyle \max_{x \in V_D} f(x) \ge f(x) ' align=absmiddle> for all <img src='http://www.wangsblog.com/jeffrey//pictures/d7aae3c674ec344e4c5581d3d580e0ef.gif' title=' x \in D ' alt=' x \in D ' align=absmiddle>. Furthermore, since <img src='http://www.wangsblog.com/jeffrey//pictures/74fd22ca4f9458f85c2c2ec87727a102.gif' title=' V_D ' alt=' V_D ' align=absmiddle> is a subset of <img src='http://www.wangsblog.com/jeffrey//pictures/2f4396bab5869c1e0c9f8a7620bf2518.gif' title=' D ' alt=' D ' align=absmiddle>, we know that this maximum is attained, thus</p>
<p><img src='http://www.wangsblog.com/jeffrey//pictures/b24695d0cf2355be6999103dcc2e1428.gif' title=' \displaystyle \max_{x \in D} f(x) = \max_{x \in V_D} f(x) ' alt=' \displaystyle \max_{x \in D} f(x) = \max_{x \in V_D} f(x) ' align=absmiddle></p>
<p>as desired. QED.</p>
<p>&#8212;&#8212;&#8212;&#8212;&#8212;&#8212;&#8211;</p>
<p>Comment: This is a very useful result, as it allows us to limit our search for the maximum of a function to its set of vertices, which is usually a considerably smaller set, though it may still be infinite (think sphere). In any case, this works well for constrained optimization problems in which the domain is limited to a polygon, the easiest application of this theorem.</p>
<p>&#8212;&#8212;&#8212;&#8212;&#8212;&#8212;&#8211;</p>
<p>Practice Problem: Let <img src='http://www.wangsblog.com/jeffrey//pictures/6722c218a6f30869ef6886dc4b050a37.gif' title=' x ' alt=' x ' align=absmiddle> and <img src='http://www.wangsblog.com/jeffrey//pictures/ee23c4f89cdc7b5be951059c2435fa2d.gif' title=' y ' alt=' y ' align=absmiddle> be real numbers satisfying <img src='http://www.wangsblog.com/jeffrey//pictures/089e257151c877fdfe6e89bef12ad7b4.gif' title=' |2x-y| \le 3 ' alt=' |2x-y| \le 3 ' align=absmiddle> and <img src='http://www.wangsblog.com/jeffrey//pictures/de5d6552927ffcf182d9f57f01cfaf62.gif' title=' |y-3x| \le 1 ' alt=' |y-3x| \le 1 ' align=absmiddle>. Find the maximum value of <img src='http://www.wangsblog.com/jeffrey//pictures/6b2235cd9fa36dad9a3b753aa0dd7f3e.gif' title=' f(x, y) = (x-y)^2 ' alt=' f(x, y) = (x-y)^2 ' align=absmiddle>.</p>
]]></content:encoded>
			<wfw:commentRss>http://wangsblog.com/jeffrey/?feed=rss2&amp;p=327</wfw:commentRss>
		<slash:comments>3</slash:comments>
		</item>
		<item>
		<title>A Square of Divisors. Topic: Number Theory. Level: AIME/Olympiad.</title>
		<link>http://wangsblog.com/jeffrey/?p=326</link>
		<comments>http://wangsblog.com/jeffrey/?p=326#comments</comments>
		<pubDate>Wed, 09 May 2007 04:06:49 +0000</pubDate>
		<dc:creator>paladin8</dc:creator>
				<category><![CDATA[AIME]]></category>
		<category><![CDATA[Number Theory]]></category>
		<category><![CDATA[Olympiad]]></category>

		<guid isPermaLink="false">http://wangsblog.com/jeffrey/?p=326</guid>
		<description><![CDATA[Problem: (1999 Canada &#8211; #3) Determine all positive integers  with the property that . Here  denotes the number of positive divisors of .
Solution: So, testing some small numbers yields  as a solution. We claim that this is the only such solution.
Clearly, since  is a square, we can use a variant of [...]]]></description>
			<content:encoded><![CDATA[<p><strong>Problem</strong>: (1999 Canada &#8211; #3) Determine all positive integers <img src='http://www.wangsblog.com/jeffrey//pictures/05e9d7891b6c7bceed4d4b9a03612612.gif' title=' n &gt; 1 ' alt=' n &gt; 1 ' align=absmiddle> with the property that <img src='http://www.wangsblog.com/jeffrey//pictures/f3532e861d30d7ea3145595368c91c86.gif' title='n = (d(n))^2' alt='n = (d(n))^2' align=absmiddle>. Here <img src='http://www.wangsblog.com/jeffrey//pictures/a3cbe72badc6a5e482ba0d642aa689fa.gif' title='d(n)' alt='d(n)' align=absmiddle> denotes the number of positive divisors of <img src='http://www.wangsblog.com/jeffrey//pictures/7b8b965ad4bca0e41ab51de7b31363a1.gif' title='n' alt='n' align=absmiddle>.</p>
<p><strong>Solution</strong>: So, testing some small numbers yields <img src='http://www.wangsblog.com/jeffrey//pictures/ee0c812da8c3234207661b58e97facf1.gif' title=' n = 9 ' alt=' n = 9 ' align=absmiddle> as a solution. We claim that this is the only such solution.</p>
<p>Clearly, since <img src='http://www.wangsblog.com/jeffrey//pictures/bfbdd7d089006253c9a32f7c78c15270.gif' title=' n ' alt=' n ' align=absmiddle> is a square, we can use a variant of the usual prime decomposition and say that <img src='http://www.wangsblog.com/jeffrey//pictures/2bd61140af6fff8fdbe1be5e50c16e92.gif' title=' n = p_1^{2e_1} p_2^{2e_2} \cdots p_k^{2e_k} ' alt=' n = p_1^{2e_1} p_2^{2e_2} \cdots p_k^{2e_k} ' align=absmiddle>.</p>
<p>Furthermore, again since <img src='http://www.wangsblog.com/jeffrey//pictures/bfbdd7d089006253c9a32f7c78c15270.gif' title=' n ' alt=' n ' align=absmiddle> is a square, we know</p>
<p><img src='http://www.wangsblog.com/jeffrey//pictures/f17e44cb09ffab1c7cabc8908d560cd3.gif' title=' d(n) = (2e_1+1)(2e_2+1) \cdots (2e_k+1) ' alt=' d(n) = (2e_1+1)(2e_2+1) \cdots (2e_k+1) ' align=absmiddle></p>
<p>is odd, so <img src='http://www.wangsblog.com/jeffrey//pictures/427317d04878a103a5eb6a06af097423.gif' title=' n = (d(n))^2 ' alt=' n = (d(n))^2 ' align=absmiddle> must be odd as well, i.e. <img src='http://www.wangsblog.com/jeffrey//pictures/8c0f68e1eb5091d1974d1cc03a0d1f9b.gif' title=' 2 ' alt=' 2 ' align=absmiddle> is not one of the <img src='http://www.wangsblog.com/jeffrey//pictures/7ff32fdf5361cfc3af72db5e27ea9c90.gif' title=' p_i ' alt=' p_i ' align=absmiddle>. Then we use the equation given to us to get</p>
<p><img src='http://www.wangsblog.com/jeffrey//pictures/e3917188d4b65062b91e37e417f49ff8.gif' title=' p_1^{2e_1} p_2^{2e_2} \cdots p_k^{2e_k} = (2e_1+1)^2(2e_2+1)^2 \cdots (2e_k+1)^2 ' alt=' p_1^{2e_1} p_2^{2e_2} \cdots p_k^{2e_k} = (2e_1+1)^2(2e_2+1)^2 \cdots (2e_k+1)^2 ' align=absmiddle></p>
<p><img src='http://www.wangsblog.com/jeffrey//pictures/984aa03b457c889c8145a9ab93153556.gif' title=' p_1^{e_1} p_2^{e_2} \cdots p_k^{e_k} = (2e_1+1)(2e_2+1) \cdots (2e_k+1) ' alt=' p_1^{e_1} p_2^{e_2} \cdots p_k^{e_k} = (2e_1+1)(2e_2+1) \cdots (2e_k+1) ' align=absmiddle>.</p>
<p>Note, however, that by Bernoulli&#8217;s Inequality (overkill, I know) we have for <img src='http://www.wangsblog.com/jeffrey//pictures/d8ab9d6414448e53b6ad36ec43a5c38d.gif' title=' i = 1, 2, \ldots, k ' alt=' i = 1, 2, \ldots, k ' align=absmiddle></p>
<p><img src='http://www.wangsblog.com/jeffrey//pictures/9692bd0afb8d162acfe32c5e6d083ddc.gif' title=' (1+(p_i-1))^{e_i} \ge 1+(p_i-1)e_i \ge 2e_i+1 ' alt=' (1+(p_i-1))^{e_i} \ge 1+(p_i-1)e_i \ge 2e_i+1 ' align=absmiddle></p>
<p>with equality iff <img src='http://www.wangsblog.com/jeffrey//pictures/47954eec5ba886bc4c36e856aaf1d305.gif' title=' p_i = 3 ' alt=' p_i = 3 ' align=absmiddle> and <img src='http://www.wangsblog.com/jeffrey//pictures/a0ad3cdbe321f6a45c58d08ac7749d25.gif' title=' e_i = 1 ' alt=' e_i = 1 ' align=absmiddle>. So</p>
<p><img src='http://www.wangsblog.com/jeffrey//pictures/16f7c7a4c833ee382fc2b641415b908d.gif' title=' p_1^{e_1} p_2^{e_2} \cdots p_k^{e_k} \ge (2e_1+1)(2e_2+1) \cdots (2e_k+1) ' alt=' p_1^{e_1} p_2^{e_2} \cdots p_k^{e_k} \ge (2e_1+1)(2e_2+1) \cdots (2e_k+1) ' align=absmiddle>.</p>
<p>Since we want equality, we must have <img src='http://www.wangsblog.com/jeffrey//pictures/47954eec5ba886bc4c36e856aaf1d305.gif' title=' p_i = 3 ' alt=' p_i = 3 ' align=absmiddle> and <img src='http://www.wangsblog.com/jeffrey//pictures/a0ad3cdbe321f6a45c58d08ac7749d25.gif' title=' e_i = 1 ' alt=' e_i = 1 ' align=absmiddle> for all <img src='http://www.wangsblog.com/jeffrey//pictures/72f8ab13f56f855e098e0ea6e73251c1.gif' title=' i ' alt=' i ' align=absmiddle>. But since the primes are supposed to be distinct we can have exactly one prime so that <img src='http://www.wangsblog.com/jeffrey//pictures/6362e1440fbc34177c937be90418f380.gif' title=' n = 3^2 = 9 ' alt=' n = 3^2 = 9 ' align=absmiddle> is the only solution. QED.</p>
<p>&#8212;&#8212;&#8212;&#8212;&#8212;&#8212;&#8211;</p>
<p>Comment: Another one of those problems that you kind of look at the result and you&#8217;re like huh, that&#8217;s interesting. But anyway, just throwing in some weak inequalities led to a pretty straightforward solution. As long as you know how to find the number of divisors of a positive integer it isn&#8217;t too much of a stretch to figure the rest out, though it make take some time to get in the right direction since the problem is quite open-ended.</p>
<p>&#8212;&#8212;&#8212;&#8212;&#8212;&#8212;&#8211;</p>
<p>Practice Problem: (1999 Canada &#8211; #4) Suppose <img src='http://www.wangsblog.com/jeffrey//pictures/4d27475ebbb3474b3be6bf801410018c.gif' title='a_1,a_2,\cdots,a_8' alt='a_1,a_2,\cdots,a_8' align=absmiddle> are eight distinct integers from <img src='http://www.wangsblog.com/jeffrey//pictures/d25ad55b8e95922b638273c6fcb06fc2.gif' title='\{1,2,\ldots,16,17\}' alt='\{1,2,\ldots,16,17\}' align=absmiddle>. Show that there is an integer <img src='http://www.wangsblog.com/jeffrey//pictures/99fb54c066b724ea0b7bf0cf745f8232.gif' title='k &gt; 0' alt='k &gt; 0' align=absmiddle> such that the equation <img src='http://www.wangsblog.com/jeffrey//pictures/2af4008eefa336d10d3d4ab891e665f0.gif' title='a_i &amp;#8211; a_j = k' alt='a_i &amp;#8211; a_j = k' align=absmiddle> has at least three different solutions. Also, find a specific set of <img src='http://www.wangsblog.com/jeffrey//pictures/8f14e45fceea167a5a36dedd4bea2543.gif' title='7' alt='7' align=absmiddle> distinct integers from <img src='http://www.wangsblog.com/jeffrey//pictures/d25ad55b8e95922b638273c6fcb06fc2.gif' title='\{1,2,\ldots,16,17\}' alt='\{1,2,\ldots,16,17\}' align=absmiddle> such that the equation <img src='http://www.wangsblog.com/jeffrey//pictures/2af4008eefa336d10d3d4ab891e665f0.gif' title='a_i &amp;#8211; a_j = k' alt='a_i &amp;#8211; a_j = k' align=absmiddle> does not have three distinct solutions for any <img src='http://www.wangsblog.com/jeffrey//pictures/99fb54c066b724ea0b7bf0cf745f8232.gif' title='k &gt; 0' alt='k &gt; 0' align=absmiddle>.</p>
]]></content:encoded>
			<wfw:commentRss>http://wangsblog.com/jeffrey/?feed=rss2&amp;p=326</wfw:commentRss>
		<slash:comments>6</slash:comments>
		</item>
		<item>
		<title>Ready For AP Calculus? Topic: Calculus/S&amp;S. Level: AIME/Olympiad.</title>
		<link>http://wangsblog.com/jeffrey/?p=325</link>
		<comments>http://wangsblog.com/jeffrey/?p=325#comments</comments>
		<pubDate>Sat, 05 May 2007 20:31:20 +0000</pubDate>
		<dc:creator>paladin8</dc:creator>
				<category><![CDATA[AIME]]></category>
		<category><![CDATA[Calculus]]></category>
		<category><![CDATA[Olympiad]]></category>
		<category><![CDATA[Sequences &#038; Series]]></category>

		<guid isPermaLink="false">http://wangsblog.com/jeffrey/?p=325</guid>
		<description><![CDATA[Problem: Evaluate .
Solution: Let . Then , a well-known Taylor series. So we want to integrate this:

by parts using  and . Substituting  in the last integral, we have
.
So . Thus

for some constant . Using our knowledge that , , and  by L&#8217;Hopital twice, we see that

Then setting  we obtain  and [...]]]></description>
			<content:encoded><![CDATA[<p><strong>Problem</strong>: Evaluate <img src='http://www.wangsblog.com/jeffrey//pictures/525d3d736e96221704e20ee6621a8ce6.gif' title=' \displaystyle \sum_{n=1}^{\infty} \frac{1}{n^2 \cdot 2^n} ' alt=' \displaystyle \sum_{n=1}^{\infty} \frac{1}{n^2 \cdot 2^n} ' align=absmiddle>.</p>
<p><strong>Solution</strong>: Let <img src='http://www.wangsblog.com/jeffrey//pictures/2a0704cadb6de5ed74c6af9867e68748.gif' title=' \displaystyle f(x) = \sum_{n=1}^{\infty} \frac{x^n}{n^2} ' alt=' \displaystyle f(x) = \sum_{n=1}^{\infty} \frac{x^n}{n^2} ' align=absmiddle>. Then <img src='http://www.wangsblog.com/jeffrey//pictures/0444bca95c3e559bb82531f2bb1d52e9.gif' title=' \displaystyle f^{\prime}(x) = \sum_{n=1}^{\infty} \frac{x^{n-1}}{n} = -\frac{\ln{(1-x)}}{x} ' alt=' \displaystyle f^{\prime}(x) = \sum_{n=1}^{\infty} \frac{x^{n-1}}{n} = -\frac{\ln{(1-x)}}{x} ' align=absmiddle>, a well-known Taylor series. So we want to integrate this:</p>
<p><img src='http://www.wangsblog.com/jeffrey//pictures/b6a1f96007040bc8f2529ea05ba6a897.gif' title=' \displaystyle \int \frac{\ln{(1-x)}}{x}dx = \ln{x} \ln{(1-x)}-\int \frac{\ln{x}}{x-1}dx ' alt=' \displaystyle \int \frac{\ln{(1-x)}}{x}dx = \ln{x} \ln{(1-x)}-\int \frac{\ln{x}}{x-1}dx ' align=absmiddle></p>
<p>by parts using <img src='http://www.wangsblog.com/jeffrey//pictures/be082745e5296802aecca02fb9859baa.gif' title=' u = \ln{(1-x)} ' alt=' u = \ln{(1-x)} ' align=absmiddle> and <img src='http://www.wangsblog.com/jeffrey//pictures/bb306b30a9149a4f846be77a4532b602.gif' title=' dv = dx/x ' alt=' dv = dx/x ' align=absmiddle>. Substituting <img src='http://www.wangsblog.com/jeffrey//pictures/f92067800e29dce7227e18a45a4fd783.gif' title=' t = 1-x ' alt=' t = 1-x ' align=absmiddle> in the last integral, we have</p>
<p><img src='http://www.wangsblog.com/jeffrey//pictures/f7e760115d08dbc796aa3288f6c928cb.gif' title=' \displaystyle \int \frac{\ln{(1-x)}}{x}dx = \ln{x} \ln{(1-x)}-\int \frac{\ln{(1-t)}}{t}dt ' alt=' \displaystyle \int \frac{\ln{(1-x)}}{x}dx = \ln{x} \ln{(1-x)}-\int \frac{\ln{(1-t)}}{t}dt ' align=absmiddle>.</p>
<p>So <img src='http://www.wangsblog.com/jeffrey//pictures/db89ba8482c05e8d835e983b2e462191.gif' title=' -f(x) = \ln{x} \ln{(1-x)}+f(t)+C = \ln{x} \ln{(1-x)}+f(1-x)+C ' alt=' -f(x) = \ln{x} \ln{(1-x)}+f(t)+C = \ln{x} \ln{(1-x)}+f(1-x)+C ' align=absmiddle>. Thus</p>
<p><img src='http://www.wangsblog.com/jeffrey//pictures/cb1d1e7e3c3d68ca474766cde718511c.gif' title=' f(x)+f(1-x) = C-\ln{x} \ln{(1-x)} ' alt=' f(x)+f(1-x) = C-\ln{x} \ln{(1-x)} ' align=absmiddle></p>
<p>for some constant <img src='http://www.wangsblog.com/jeffrey//pictures/699e902f5598ca623370f833cffb1a57.gif' title=' C ' alt=' C ' align=absmiddle>. Using our knowledge that <img src='http://www.wangsblog.com/jeffrey//pictures/fe0c4947ffdad5dd8b5e3609021cb2ec.gif' title=' f(0) = 0 ' alt=' f(0) = 0 ' align=absmiddle>, <img src='http://www.wangsblog.com/jeffrey//pictures/61769437df02a7504e91a7bde833815c.gif' title=' \displaystyle f(1) = \frac{\pi^2}{6} ' alt=' \displaystyle f(1) = \frac{\pi^2}{6} ' align=absmiddle>, and <img src='http://www.wangsblog.com/jeffrey//pictures/b6842ac344f6bbb0f47b4bee04479088.gif' title=' \displaystyle \lim_{x \rightarrow 1} \ln{x} \ln{(1-x)}= 0 ' alt=' \displaystyle \lim_{x \rightarrow 1} \ln{x} \ln{(1-x)}= 0 ' align=absmiddle> by L&#8217;Hopital twice, we see that</p>
<p><img src='http://www.wangsblog.com/jeffrey//pictures/f159e4c93e18a993aa996956a3ba5ef8.gif' title=' f(0)+f(1) = C \Rightarrow C = \frac{\pi^2}{6} ' alt=' f(0)+f(1) = C \Rightarrow C = \frac{\pi^2}{6} ' align=absmiddle></p>
<p>Then setting <img src='http://www.wangsblog.com/jeffrey//pictures/ea674b31fb5fc34dfee7a9dbf435b585.gif' title=' x = \frac{1}{2} ' alt=' x = \frac{1}{2} ' align=absmiddle> we obtain <img src='http://www.wangsblog.com/jeffrey//pictures/8261ef1838aa499fed128b9d0e714251.gif' title=' \displaystyle 2f\left(\frac{1}{2}\right) = \frac{\pi^2}{6}-\ln{\frac{1}{2}} \ln{\frac{1}{2}} ' alt=' \displaystyle 2f\left(\frac{1}{2}\right) = \frac{\pi^2}{6}-\ln{\frac{1}{2}} \ln{\frac{1}{2}} ' align=absmiddle> and <img src='http://www.wangsblog.com/jeffrey//pictures/78d505a0961d908d1f1bdbc70e5f0c41.gif' title=' \displaystyle \sum_{n=1}^{\infty} \frac{1}{n^2 \cdot 2^n} = f\left(\frac{1}{2}\right) = \frac{\pi^2}{12}-\frac{\ln^2{(2)}}{2} ' alt=' \displaystyle \sum_{n=1}^{\infty} \frac{1}{n^2 \cdot 2^n} = f\left(\frac{1}{2}\right) = \frac{\pi^2}{12}-\frac{\ln^2{(2)}}{2} ' align=absmiddle>. QED.</p>
<p>&#8212;&#8212;&#8212;&#8212;&#8212;&#8212;&#8211;</p>
<p>Comment: This was a pretty tough problem that required you to compound a lot of calculus knowledge all into a single problem &#8211; series, integration by parts, limits. Recognizing all the steps was the first part; following through with the right computations was another. Still, there weren&#8217;t really any super clever tricks, mostly just standard substitutions and approaches applied in a somewhat non-standard way. Makes for a very nice problem.</p>
<p>&#8212;&#8212;&#8212;&#8212;&#8212;&#8212;&#8211;</p>
<p>Practice Problem #1: Show that <img src='http://www.wangsblog.com/jeffrey//pictures/2fec09949acad8c8f08b1c075a6637b0.gif' title=' \displaystyle \lim_{x \rightarrow 1} \ln{x} \ln{(1-x)} = 0 ' alt=' \displaystyle \lim_{x \rightarrow 1} \ln{x} \ln{(1-x)} = 0 ' align=absmiddle>.</p>
<p>Practice Problem #2: Evaluate <img src='http://www.wangsblog.com/jeffrey//pictures/4c4c303221f73b401c71a8dfb135202d.gif' title=' \displaystyle \sum_{n=1}^{\infty} \frac{x^n}{n} ' alt=' \displaystyle \sum_{n=1}^{\infty} \frac{x^n}{n} ' align=absmiddle>. Can you also find <img src='http://www.wangsblog.com/jeffrey//pictures/c5f2e70e99471324833a4bc45e340525.gif' title=' \displaystyle \sum_{n=1}^{\infty} \frac{x^n}{n^3} ' alt=' \displaystyle \sum_{n=1}^{\infty} \frac{x^n}{n^3} ' align=absmiddle>?</p>
]]></content:encoded>
			<wfw:commentRss>http://wangsblog.com/jeffrey/?feed=rss2&amp;p=325</wfw:commentRss>
		<slash:comments>0</slash:comments>
		</item>
		<item>
		<title>Just A &#8220;Natural&#8221; Thing. Topic: Algebra. Level: AIME/Olympiad.</title>
		<link>http://wangsblog.com/jeffrey/?p=324</link>
		<comments>http://wangsblog.com/jeffrey/?p=324#comments</comments>
		<pubDate>Sat, 28 Apr 2007 21:46:44 +0000</pubDate>
		<dc:creator>paladin8</dc:creator>
				<category><![CDATA[AIME]]></category>
		<category><![CDATA[Algebra]]></category>
		<category><![CDATA[Olympiad]]></category>

		<guid isPermaLink="false">http://wangsblog.com/jeffrey/?p=324</guid>
		<description><![CDATA[Problem: Factor .
Solution: Consider the polynomial . We have
,
so we want to factor the second term when . Call it  so that . Consider the relation
.
Since  is a root of the LHS, we factor it out of the RHS as well to get
.
Dividing through by  and rearranging, we obtain the nice expression
.
Letting [...]]]></description>
			<content:encoded><![CDATA[<p><strong>Problem</strong>: Factor <img src='http://www.wangsblog.com/jeffrey//pictures/4df1f2c7288e3b1ed1e28b768abbc1dc.gif' title=' 7^{6(2k-1)}-7^{5(2k-1)}+7^{4(2k-1)}-7^{3(2k-1)}+7^{2(2k-1)}-7^{2k-1}+1 ' alt=' 7^{6(2k-1)}-7^{5(2k-1)}+7^{4(2k-1)}-7^{3(2k-1)}+7^{2(2k-1)}-7^{2k-1}+1 ' align=absmiddle>.</p>
<p><strong>Solution</strong>: Consider the polynomial <img src='http://www.wangsblog.com/jeffrey//pictures/fd99add89f6d3de5bc79b4b92f38499c.gif' title=' x^7+1 ' alt=' x^7+1 ' align=absmiddle>. We have</p>
<p><img src='http://www.wangsblog.com/jeffrey//pictures/0463aab95360ee712653bbec3ee020e0.gif' title=' x^7+1 = (x+1)(x^6-x^5+x^4-x^3+x^2-1) ' alt=' x^7+1 = (x+1)(x^6-x^5+x^4-x^3+x^2-1) ' align=absmiddle>,</p>
<p>so we want to factor the second term when <img src='http://www.wangsblog.com/jeffrey//pictures/49dc08eaeb33eabe7bcfb4e2f0670d6f.gif' title=' x = 7^{2k-1} ' alt=' x = 7^{2k-1} ' align=absmiddle>. Call it <img src='http://www.wangsblog.com/jeffrey//pictures/1343923278289008b04ad5f121b3ca76.gif' title=' P(x) ' alt=' P(x) ' align=absmiddle> so that <img src='http://www.wangsblog.com/jeffrey//pictures/3c71acfee22b1cb0770a222bb586cb78.gif' title=' P(x) = \frac{x^7+1}{x+1} ' alt=' P(x) = \frac{x^7+1}{x+1} ' align=absmiddle>. Consider the relation</p>
<p><img src='http://www.wangsblog.com/jeffrey//pictures/83adb225b8235b8b2745469f1b73c7bd.gif' title=' (x+1)^7-(x^7+1) = 7x^6+21x^5+35x^4+35x^3+21x^2+7x ' alt=' (x+1)^7-(x^7+1) = 7x^6+21x^5+35x^4+35x^3+21x^2+7x ' align=absmiddle>.</p>
<p>Since <img src='http://www.wangsblog.com/jeffrey//pictures/d888afbeeba19927b1b7a6eb748ba7ea.gif' title=' -1 ' alt=' -1 ' align=absmiddle> is a root of the LHS, we factor it out of the RHS as well to get</p>
<p><img src='http://www.wangsblog.com/jeffrey//pictures/5d231d43b46bc8380d80f83b32d9ab82.gif' title=' (x+1)^7-(x^7+1) = 7x(x+1)(x^4+2x^3+3x^2+2x+1) = 7x(x+1)(x^2+x+1)^2 ' alt=' (x+1)^7-(x^7+1) = 7x(x+1)(x^4+2x^3+3x^2+2x+1) = 7x(x+1)(x^2+x+1)^2 ' align=absmiddle>.</p>
<p>Dividing through by <img src='http://www.wangsblog.com/jeffrey//pictures/2fefaab3a83b2ade501ffd890af15c0d.gif' title=' x+1 ' alt=' x+1 ' align=absmiddle> and rearranging, we obtain the nice expression</p>
<p><img src='http://www.wangsblog.com/jeffrey//pictures/e9827ee44df296d6579c6518076c65b7.gif' title=' P(x) = (x+1)^6-7x(x^2+x+1)^2 ' alt=' P(x) = (x+1)^6-7x(x^2+x+1)^2 ' align=absmiddle>.</p>
<p>Letting <img src='http://www.wangsblog.com/jeffrey//pictures/49dc08eaeb33eabe7bcfb4e2f0670d6f.gif' title=' x = 7^{2k-1} ' alt=' x = 7^{2k-1} ' align=absmiddle>, this becomes</p>
<p><img src='http://www.wangsblog.com/jeffrey//pictures/18612c4a25d26a531eada3e227b8ad04.gif' title=' P(7^{2k-1}) = (7^{3(2k-1)}+3 \cdot 7^{2(2k-1)}+3 \cdot 7^{2k-1}+1)^2-7^{2k}(7^{2(2k-1)}+7^{2k-1}+1)^2 ' alt=' P(7^{2k-1}) = (7^{3(2k-1)}+3 \cdot 7^{2(2k-1)}+3 \cdot 7^{2k-1}+1)^2-7^{2k}(7^{2(2k-1)}+7^{2k-1}+1)^2 ' align=absmiddle></p>
<p>which is conveniently enough a difference of two squares. And we&#8217;ll leave it as this because the factors are not particularly nice or anything. QED.</p>
<p>&#8212;&#8212;&#8212;&#8212;&#8212;&#8212;&#8211;</p>
<p>Comment: This &#8220;natural&#8221; factorization was basically the one crucial step to solving the USAMO #5 this year and most people did not see it, unsurprisingly. Taking the difference <img src='http://www.wangsblog.com/jeffrey//pictures/ab8b18191075fccbb561f7b3f08c7991.gif' title=' (x+1)^7-(x^7+1) ' alt=' (x+1)^7-(x^7+1) ' align=absmiddle> was the trickiest/cleverest part, and there were definitely a limited number of approaches to this factorization. Oh well, at least it seems sort of cool after you know about it.</p>
<p>&#8212;&#8212;&#8212;&#8212;&#8212;&#8212;&#8211;</p>
<p>Practice Problem: (2007 USAMO &#8211; #5) Prove that for every nonnegative integer n, the number <img src='http://www.wangsblog.com/jeffrey//pictures/78f5038c64a5e15817784309632c784a.gif' title='7^{7^{n}}+1' alt='7^{7^{n}}+1' align=absmiddle> is the product of at least <img src='http://www.wangsblog.com/jeffrey//pictures/ba5c7c09f61d4fe49d6d0209cb3c4b83.gif' title='2n+3' alt='2n+3' align=absmiddle> (not necessarily distinct) primes.</p>
]]></content:encoded>
			<wfw:commentRss>http://wangsblog.com/jeffrey/?feed=rss2&amp;p=324</wfw:commentRss>
		<slash:comments>1</slash:comments>
		</item>
	</channel>
</rss>
