TAOCP 7.2.1.4 Exercise 45
Let $p(n)$ denote the partition function.
Section 7.2.1.4: Generating All Partitions
Exercise 45. [HM21] [HM21] Compute the asymptotic value of p(n−1)/p(n), with relative error O(n−2).
Verified: no
Solve time: 21m59s
Solution
Let $p(n)$ denote the partition function. From the Hardy–Ramanujan asymptotic formula (first Rademacher term suffices at the required precision),
$$ p(n)=\frac{1}{4n\sqrt{3}}\exp!\left(a\sqrt{n}\right)\left(1+O!\left(e^{-c\sqrt{n}}\right)\right), \qquad a=\pi\sqrt{\frac{2}{3}}, $$
so the error is exponentially small and does not affect any expansion in powers of $n^{-1/2}$.
Hence
$$ \frac{p(n-1)}{p(n)}
\frac{n}{n-1}\exp!\bigl(a(\sqrt{n-1}-\sqrt{n})\bigr)\left(1+O!\left(e^{-c\sqrt{n}}\right)\right). $$
Set
$$ \sqrt{n-1}=\sqrt{n}\left(1-\frac{1}{n}\right)^{1/2}. $$
Using the binomial expansion,
$$ \left(1-\frac{1}{n}\right)^{1/2}
1-\frac{1}{2n}-\frac{1}{8n^{2}}-\frac{1}{16n^{3}}+O(n^{-4}), $$
so
$$ \sqrt{n-1}
\sqrt{n} -\frac{1}{2\sqrt{n}} -\frac{1}{8n^{3/2}} -\frac{1}{16n^{5/2}} +O(n^{-7/2}). $$
Therefore
$$ \Delta:=\sqrt{n-1}-\sqrt{n}
-\frac{1}{2\sqrt{n}} -\frac{1}{8n^{3/2}} -\frac{1}{16n^{5/2}} +O(n^{-7/2}). $$
Define
$$ t_1=-\frac{a}{2\sqrt{n}},\quad t_2=-\frac{a}{8n^{3/2}},\quad t_3=-\frac{a}{16n^{5/2}}. $$
Then
$$ a\Delta=t_1+t_2+t_3+O(n^{-7/2}). $$
Expansion of the exponential
First,
$$ \exp(t_1)
1-\frac{a}{2\sqrt{n}}+\frac{a^2}{8n}-\frac{a^3}{48n^{3/2}}+\frac{a^4}{384n^2}+O(n^{-5/2}). $$
Next,
$$ \exp(t_2)=1-\frac{a}{8n^{3/2}}+O(n^{-3}), $$
and $t_3$ contributes only beyond $n^{-5/2}$ in the final product.
Multiplying,
$$ \exp(a\Delta)
1-\frac{a}{2\sqrt{n}}+\frac{a^2}{8n} -\left(\frac{a^3}{48}+\frac{a}{8}\right)n^{-3/2} +\left(\frac{a^4}{384}+\frac{a^2}{16}\right)n^{-2} +O(n^{-5/2}). $$
Prefactor expansion
$$ \frac{n}{n-1}=\frac{1}{1-1/n}=1+\frac{1}{n}+\frac{1}{n^2}+O(n^{-3}). $$
Final multiplication
Let
$$ S= 1+s_1 n^{-1/2}+s_2 n^{-1}+s_3 n^{-3/2}+s_4 n^{-2}, $$
where
$$ s_1=-\frac{a}{2},\quad s_2=\frac{a^2}{8},\quad s_3=-\left(\frac{a^3}{48}+\frac{a}{8}\right),\quad s_4=\frac{a^4}{384}+\frac{a^2}{16}. $$
Then
$$ \frac{p(n-1)}{p(n)}
\left(1+\frac{1}{n}+\frac{1}{n^2}\right)S+O(n^{-5/2}). $$
Collecting terms up to $n^{-2}$ gives
$$ \frac{p(n-1)}{p(n)}
1 +s_1 n^{-1/2} +(s_2+1)n^{-1} +(s_3+2s_1)n^{-3/2} +(s_4+s_2+1)n^{-2} +O(n^{-5/2}). $$
Substituting $a=\pi\sqrt{\frac{2}{3}}$ yields the explicit expansion
$$ \frac{p(n-1)}{p(n)}
1 -\frac{\pi}{\sqrt{6n}} +\left(\frac{\pi^2}{12}+1\right)\frac{1}{n} -\left(\frac{\pi^3}{36\sqrt{6}}+\frac{\pi}{\sqrt{6}}\right)\frac{1}{n^{3/2}} +\left(\frac{\pi^4}{288}+\frac{\pi^2}{24}+1\right)\frac{1}{n^2} +O(n^{-5/2}). $$
Thus the asymptotic value of $p(n-1)/p(n)$, with relative error $O(n^{-2})$, is
$$ \boxed{ 1 -\frac{\pi}{\sqrt{6n}} +\left(\frac{\pi^2}{12}+1\right)\frac{1}{n} -\left(\frac{\pi^3}{36\sqrt{6}}+\frac{\pi}{\sqrt{6}}\right)\frac{1}{n^{3/2}} +\left(\frac{\pi^4}{288}+\frac{\pi^2}{24}+1\right)\frac{1}{n^2} }. $$
This completes the solution. ∎