Дано выражение для числа сочетаний C(n,k). Объясните, как оценить C(1000,500) асимптотически и сравните применение формулы Стирлинга и логарифмических приближений
Коротко: C(1000,500)=1000!500!2C(1000,500)=\dfrac{1000!}{500!^2}C(1000,500)=500!21000!. Оценка по Стирлингу даёт простую и точную асимптику: 1) Стирлинг (основная формула) n!∼2πn (ne)n.
n!\sim\sqrt{2\pi n}\,\Big(\frac{n}{e}\Big)^n . n!∼2πn(en)n.
Применяя к 1000!1000!1000! и 500!500!500! получаем (для n=500n=500n=500, т.е. C(2n,n)C(2n,n)C(2n,n)) C(1000,500)∼4500π⋅500.
C(1000,500)\sim\frac{4^{500}}{\sqrt{\pi\cdot 500}}. C(1000,500)∼π⋅5004500.
Численно π⋅500≈39.636⇒C(1000,500)≈2100039.636≈2.703×10299.
\sqrt{\pi\cdot 500}\approx 39.636\quad\Rightarrow\quad C(1000,500)\approx\frac{2^{1000}}{39.636}\approx 2.703\times 10^{299}. π⋅500≈39.636⇒C(1000,500)≈39.63621000≈2.703×10299. 2) Логарифмические приближения (энтропийный подход) Общее асимптотическое поведение для C(n,k)C(n,k)C(n,k) записывается через энтропию H(p)=−plnp−(1−p)ln(1−p),p=kn,
H(p)=-p\ln p-(1-p)\ln(1-p),\qquad p=\frac{k}{n}, H(p)=−plnp−(1−p)ln(1−p),p=nk,
и lnC(n,k)=nH(p)+O(lnn).
\ln C(n,k)=nH(p)+O(\ln n). lnC(n,k)=nH(p)+O(lnn).
Для k=n/2k=n/2k=n/2 имеем H(1/2)=ln2H(1/2)=\ln 2H(1/2)=ln2, поэтому в нулевом порядке lnC(1000,500)≈1000ln2≈693.147.
\ln C(1000,500)\approx 1000\ln 2\approx 693.147. lnC(1000,500)≈1000ln2≈693.147.
Стирлингово уточнение даёт следующий член порядка lnn\ln nlnn: lnC(1000,500)=1000ln2−12ln(π⋅500)+⋯
\ln C(1000,500)=1000\ln 2-\tfrac12\ln(\pi\cdot 500)+\cdots lnC(1000,500)=1000ln2−21ln(π⋅500)+⋯
что даёт lnC≈689.4675\ln C\approx 689.4675lnC≈689.4675 и приводит к той же численной величине ≈2.70×10299 \approx2.70\times10^{299}≈2.70×10299. 3) Сравнение методов и погрешности - «Чистый» логарифмический (энтропийный) подход даёт правильный экспоненциальный порядок enH(p)e^{nH(p)}enH(p), но пропускает полиномиальный множитель (фактор n\sqrt{n}n и константы), который в логарифме даёт вклад O(lnn)O(\ln n)O(lnn) и может смещать десятки логарифмических единиц (здесь ≈3.68 в натуральных логарифмах). - Стирлинговоe приближение в форме для факториалов даёт и экспоненциальный член, и точный (порядка n−1/2n^{-1/2}n−1/2) множитель; дополнительные члены асимптотики (1/(12n)1/(12n)1/(12n), и т.д.) дают очень малую относительную поправку: для n=500n=500n=500 следующий член даёт множитель exp(−1/(8n)+… )≈0.99975\exp(-1/(8n)+\dots)\approx 0.99975exp(−1/(8n)+…)≈0.99975 (ошибка порядка 10−410^{-4}10−4 в относительном масштабе). - Практически: для численных оценок и корректных значащих цифр следует работать с логарифмами и со Стирлинговым разложением (включая член 12ln(2πn)\tfrac12\ln(2\pi n)21ln(2πn)); для оценки лишь порядка величины достаточно первого энтропийного члена. Итого: асимптотика C(1000,500)∼4500π⋅500≈2.703×10299,
C(1000,500)\sim\frac{4^{500}}{\sqrt{\pi\cdot 500}}\approx 2.703\times10^{299}, C(1000,500)∼π⋅5004500≈2.703×10299,
а использование Стирлинга в логарифмической форме даёт более точный результат и контролируемую малую относительную погрешность по сравнению с только экспоненциальным (энтропийным) приближением.
1) Стирлинг (основная формула)
n!∼2πn (ne)n. n!\sim\sqrt{2\pi n}\,\Big(\frac{n}{e}\Big)^n .
n!∼2πn (en )n. Применяя к 1000!1000!1000! и 500!500!500! получаем (для n=500n=500n=500, т.е. C(2n,n)C(2n,n)C(2n,n))
C(1000,500)∼4500π⋅500. C(1000,500)\sim\frac{4^{500}}{\sqrt{\pi\cdot 500}}.
C(1000,500)∼π⋅500 4500 . Численно
π⋅500≈39.636⇒C(1000,500)≈2100039.636≈2.703×10299. \sqrt{\pi\cdot 500}\approx 39.636\quad\Rightarrow\quad
C(1000,500)\approx\frac{2^{1000}}{39.636}\approx 2.703\times 10^{299}.
π⋅500 ≈39.636⇒C(1000,500)≈39.63621000 ≈2.703×10299.
2) Логарифмические приближения (энтропийный подход)
Общее асимптотическое поведение для C(n,k)C(n,k)C(n,k) записывается через энтропию
H(p)=−plnp−(1−p)ln(1−p),p=kn, H(p)=-p\ln p-(1-p)\ln(1-p),\qquad p=\frac{k}{n},
H(p)=−plnp−(1−p)ln(1−p),p=nk , и
lnC(n,k)=nH(p)+O(lnn). \ln C(n,k)=nH(p)+O(\ln n).
lnC(n,k)=nH(p)+O(lnn). Для k=n/2k=n/2k=n/2 имеем H(1/2)=ln2H(1/2)=\ln 2H(1/2)=ln2, поэтому в нулевом порядке
lnC(1000,500)≈1000ln2≈693.147. \ln C(1000,500)\approx 1000\ln 2\approx 693.147.
lnC(1000,500)≈1000ln2≈693.147. Стирлингово уточнение даёт следующий член порядка lnn\ln nlnn:
lnC(1000,500)=1000ln2−12ln(π⋅500)+⋯ \ln C(1000,500)=1000\ln 2-\tfrac12\ln(\pi\cdot 500)+\cdots
lnC(1000,500)=1000ln2−21 ln(π⋅500)+⋯ что даёт lnC≈689.4675\ln C\approx 689.4675lnC≈689.4675 и приводит к той же численной величине ≈2.70×10299 \approx2.70\times10^{299}≈2.70×10299.
3) Сравнение методов и погрешности
- «Чистый» логарифмический (энтропийный) подход даёт правильный экспоненциальный порядок enH(p)e^{nH(p)}enH(p), но пропускает полиномиальный множитель (фактор n\sqrt{n}n и константы), который в логарифме даёт вклад O(lnn)O(\ln n)O(lnn) и может смещать десятки логарифмических единиц (здесь ≈3.68 в натуральных логарифмах).
- Стирлинговоe приближение в форме для факториалов даёт и экспоненциальный член, и точный (порядка n−1/2n^{-1/2}n−1/2) множитель; дополнительные члены асимптотики (1/(12n)1/(12n)1/(12n), и т.д.) дают очень малую относительную поправку: для n=500n=500n=500 следующий член даёт множитель exp(−1/(8n)+… )≈0.99975\exp(-1/(8n)+\dots)\approx 0.99975exp(−1/(8n)+…)≈0.99975 (ошибка порядка 10−410^{-4}10−4 в относительном масштабе).
- Практически: для численных оценок и корректных значащих цифр следует работать с логарифмами и со Стирлинговым разложением (включая член 12ln(2πn)\tfrac12\ln(2\pi n)21 ln(2πn)); для оценки лишь порядка величины достаточно первого энтропийного члена.
Итого: асимптотика
C(1000,500)∼4500π⋅500≈2.703×10299, C(1000,500)\sim\frac{4^{500}}{\sqrt{\pi\cdot 500}}\approx 2.703\times10^{299},
C(1000,500)∼π⋅500 4500 ≈2.703×10299, а использование Стирлинга в логарифмической форме даёт более точный результат и контролируемую малую относительную погрешность по сравнению с только экспоненциальным (энтропийным) приближением.