10
правок
Изменения
Нет описания правки
{{Задача
|definition =
По заданному числу <tex>n</tex> найти количество его различных разбиений на положительные слагаемые<ref name="amount" /> <tex> m_0 + m_1 + m_2 + \dots + m_k = n </tex> так , что при всех <tex> i: m_i \leqslant m_{i+1} </tex>.}}
__TOC__
Существует много алгоритмов решения данной задачи. Вот некоторые из них:
Пусть <tex>P(n, m, k)</tex> — количество разбиений числа <tex>n</tex> на <tex>m</tex> слагаемых, каждое из которых не превосходит <tex>k</tex>.
Имеет место следующее рекуррентное соотношение:
<tex dpi = "145">P(n, m, k) = P(n, m, k - 1) + P(n - k, m - 1, k)</tex> при <tex>m,k \leqslant n</tex>
</p>
Начальные условия:
<tex>P(0,0,i) = 1, \forall i \geqslant 0</tex>
<tex>P(i, j, 0) = 0, \forall i,j>0</tex>
Рассмотрим множество разбиений числа <tex>n</tex> на <tex>m</tex> слагаемых, каждое из которых не больше <tex>k</tex>. Разделим его на две группы — в первой будут всеразбиения которые не содержат в качестве старшего слагаемого <tex>k</tex>. Таких разбиений <tex>P(n, m, k - 1)</tex>. Во второй — все разбиения со старшим слагаемым <tex>k</tex>. Их столько же, сколько разбиений числа <tex>n - k</tex> на <tex>m - 1</tex> слагаемое, каждое из которых не больше <tex>k</tex>, то есть <tex>P(n - k, m - 1, k)</tex>. Количество всех разбиений числа равно <tex>\sum\limits_{i=0}^nP(n, i, n)</tex>. Асимптотика <tex>O(n^{3})</tex> получается из-за использования трехмерной динамики. == Алгоритм за O(N<sup>2</sup>) ===
Обозначим <tex>P(n,k)</tex> — количество разбиений числа <tex>n</tex> на слагаемые, каждое из которых не превосходит <tex>k</tex>. Оно удовлетворяет следующей рекурентной формуле:<p>
<tex dpi="145">P(n,k) = \left \{ \begin{array}{ll} P(n,k - 1) + P(n - k, k), & k \leqslant n \\ P(n, n), & k > n \end{array} \right.</tex>
</p>
При этом считается , что
<tex>P(0,0) = 1</tex>
<tex> P(i, 0) = 0, \forall i>0</tex>
Рассмотрим самый быстрый алгоритм нахождения количества разбиений числа <tex>n</tex> на слагаемые, который работает за <tex> O(n \sqrt{n}) </tex>.
Итак, обозначим количество таких разбиений за <tex> p(n) </tex>.
Рассмотрим следующее бесконечное произведение:
<tex> (1 + x + x^2 + \dots)(1 + x^2 + x^4 + \dots) \dots (1 + x^k + x^{2k} + \dots) \dots </tex>
После раскрытия скобок каждый член произведения получается в результате умножения мономов (одночленов), взятых по одному из каждой скобки. Если в первой скобке взять <tex dpi = "145">x^{m_1}</tex>, во второй — <tex dpi = "145">x^{2m_2}</tex> и т.д., то их произведение будет равно <tex dpi = "145">x^{m_1 + 2m_2 + 3m_3 + \dots}.</tex> Значит, после раскрытия скобок получится сумма мономов вида <tex dpi = "145">x^{m_1 + 2m_2 + 3m_3 + \dots}</tex>.
Можно увидеть, что <tex> x^n </tex> встретится в полученной бесконечной сумме столько раз, сколькими способами можно представить <tex>n</tex> как сумму <tex>m_1 + 2m_2 + 3m_3 + ...</tex> Каждому такому представлению отвечает разбиение числа <tex>n</tex> на <tex>m_1</tex> единиц, <tex>m_2</tex> двоек и т.д. Таким образом, очевидно, получаются все разбиения, так как из первой скобки мы можем взять любое <tex>x^{m_i}</tex>, где <tex>m_i \in [0 \dots \infty),</tex> то есть произвольное количество единиц в нашем разбиении. Аналогично, мы можем взять произвольное количество двоек и т.д. Но при раскрытии скобок мы находим произведения всех возможных комбинации множителей из разных скобок. Поэтому коэффициент при <tex>x^n</tex> равен числу разбиений <tex>p(n)</tex>. Для доказательства следующего факта обратимся к [[Производящая функция|методу производящих функций]].
Посмотрим теперь на выражения в скобках. Каждое из них — бесконечная геометрическая прогрессия. По Полагая <tex>0 \leqslant x < 1</tex>, по формуле ее суммирования:
: <tex dpi = "145">1 + x + x^2 + x^3 + ... = \frac{1}{1 - x}</tex>,
: <tex dpi = "145">1 + x^2 + x^4 + x^6 + ... = \frac{1}{1-x^2}</tex>
: <tex dpi = "145">p(0) + p(1) x + p(2) x^2 + p(3) x^3 + \dots = \frac{1}{(1-x)(1-x^2)(1-x^3)\dots}</tex> <tex>(1)</tex>
Пентагональная теорема Эйлера
| statement =
<texdpi="145"> \prod\limits_{k=1}^\infty \left(1-x^k \right) = \sum\limits_{q=-\infty}^\infty (-1)^q x^{\frac{3q^2+q}{2}}</tex>
| proof =
<tex>A = (1-x)(1-x^2)(1-x^3)...(1-x^n)...</tex>
Раскроем в этом произведении первые <tex>22 </tex> скобки. Мы получим выражение
<tex>A = [1-x-x^2+x^5+x^7-x^{12}-x^{15}+x^{22}+...]\times(1-x^{23})(1-x^{24})...(1-x^n)...</tex>,
где в квадратной скобке точками обозначены слагаемые, содержащие <tex>x</tex> в более высокой степени, чем <tex>22</tex>. Не будем выписывать эти члены, так как после умножения квадратной скобки на <tex>1-x^{23}</tex>, <tex>1-x^{24}</tex> и т.д. они изменятся. Выписанные же члены больше меняться не будут. Поэтому, если раскрыть все скобки, то получится бесконечный ряд, первые члены которого имеют вид
<tex>A=1-x-x^2+x^5+x^7-x^{12}-x^{15}+x^{22}+...</tex>
<tex>A = (1-x)(1-x^2)(1-x^3)...(1-x^n)...</tex>
в ряд, то в этом ряду отличны от нуля лишь слагаемые, вида <texdpi="145">(-1)^q x^{\frac{3q^2+q}{2}}</tex>, где <tex>q</tex> - — натуральное число.
При раскрытии скобок в исходном произведении слагаемое <tex>\pm x^N</tex> встретится столько раз, сколькими способами можно разбить число <tex>N</tex> на различные слагаемые. При этом, если число слагаемых четно, то появляется <tex>x^N</tex>, а если это число нечетно, то появляется <tex>-x^N</tex>. Например, разбиению <tex>12=5+4+2+1</tex> - соответствует слагаемое <tex>(-x^5)(-x^4)(-x^2)(-x^1)=x^{12},</tex> а разбиению <tex>12=5+4+3</tex> - — слагаемое <tex>(-x^5)(-x^4)(-x^3)=-x^{12}</tex>. Таким образом, коэффициент при <tex>x^N</tex> в разложении <tex>A</tex> в ряд равен разности между количеством разбиений <tex>N</tex> на четное число различных слагаемых и количеством разбиений <tex>N</tex> на нечетное число различных слагаемых. Тогда теорему можно переформулировать следующим образом:
{{Теорема
Переформулировка пентагональной теоремы
| statement =
Если число <tex>N</tex> не может быть представлено в виде <texdpi="145">N =\frac{3q^2+q}{2}</tex>, то оно имеет одинаковое количество разбиений на четное и на нечетное число различных слагаемых. А для чисел вида <texdpi="145">N =\frac{3q^2+q}{2}</tex> разность между этими количествами равна <tex>(-1)^q</tex>.
}}
Иными словами, если <tex>q</tex> четно, то на одно больше разбиений на четное число слагаемых, а если <tex>q</tex> нечетно, то на одно больше разбиений на нечетно нечетное число слагаемых. Докажем эту теорему с помощью диаграмм Юнга<ref name="young" />. Покажем несколько способов превращения диаграммы с четным числом строк диаграмму из стольких же точек с нечетным числом строк и обратно.
Так как мы рассматриваем лишь разбиения на различные слагаемые, то диаграммы таких разбиений состоят из нескольких трапеций, поставленных друг на друга. Обозначим число точек в верхней строке диаграммы через <tex>m</tex>, а число строк нижней трапеции через <tex>n</tex> (на рис. 1 слева изображена диаграмма, для которой <tex>m=2</tex>, <tex>n=3</tex>).
[[Файл:Diagramma_2_3.jpg|thumb|right|500px|Рис. 1. <tex>m = 2, n = 3</tex>]]
'''Преобразование 1.''' Предположим, что диаграмма содержит не менее двух трапеций, причем <tex>m \leqslant n</tex>. В этом случае отбросим первую строку из <tex>m</tex> точек, но удлиним последние <tex>m</tex> строк нижней трапеции на одну точку(рис. 1). После этого общее число точек не изменится, все строки окажутся различной длины, но четность числа строк изменится.
Точно такое же преобразование можно сделать, если диаграмма состоит из одной трапеции, но <tex>n \geqslant m+1</tex>. Стираем верхнюю строку и к нижним <tex>n-1</tex> строчкам приписываем <tex>m</tex> точек.
Осталось выяснить, какие же диаграммы не допускают ни одного из описанных преобразований. Ясно, что эти диаграммы состоят из одной трапеции, причем для них либо <tex>m=n</tex>, либо <tex>m=n+1</tex>. В первом случае диаграмма содержит
<texdpi="145">n+(n+1)+(n+2)+\dots+(2n-1)=\frac{3n^2-n}{2}</tex>
точек, а во втором - на <tex>n</tex> точек больше, т.е. <texdpi="145">\frac{3n^2+n}{2}</tex>.
Приведенные рассуждения показывают, что если <tex>N</tex> не является числом вида <texdpi="145">\frac{3n^2 \pm n}{2}</tex>, то оно имеет поровну разбиений на четное и нечетное число различных слагаемых.
Очевидно также, что равенства <texdpi="145">N=\frac{3n^2+n}{2}</tex> и <texdpi="145">N=\frac{3l^2-l}{2}</tex> одновременно выполняться не могут, поэтому если <texdpi ="145">N=\frac{3n^2\pm n}{2}</tex>, то без пары останется ровно одна диаграмма, не допускающая преобразования и имеющая <tex>n</tex> строк (слагаемых <tex>N</tex>). Если <tex>n</tex> - четное число, то разбиений на четное число слагаемых окажется на <tex>1 </tex> больше, чем на нечетное число слагаемых. Если же <tex>n</tex> - нечетное число, то на <tex>1 </tex> больше будет разбиений на нечетное число слагаемых. Теорема доказана.
}}
Умножим обе части равенства <tex>(1) </tex> на <tex>\prod\limits_{k=1}^\infty \left(1-x^k \right) </tex> и воспользуемся пентагональной теоремой:: <tex> ( p(0) + p(1) x + p(2) x^2 + \dots)(1 - x - x^2 + x^5 + x^7 - x^{12} - x^{15} + \dots) = 1 </tex> <tex>(2)</tex>
Начнем раскрывать скобки, для наглядности мономы с одинаковыми степенями <tex>x </tex> пишем друг под другом:
: <tex> p(0) + p(1)x + p(2) x^2 + p(3) x^3 + p(4) x^4 + p(5) x^5 + p(6) x^6 + \dots </tex>
:: <tex> - p(0)x - p(1) x^2 - p(2) x^3 - p(3) x^4 - p(4) x^5 - p(5) x^6 - \dots </tex>
::::::::::::: <tex>+ p(0) x^5 + p(1) x^6 + \dots </tex>
Так как <tex>p(0) = 1</tex>, то оно сокращается с единицей справа. Так что, чтобы выражение <tex>(2) </tex> было удовлетворено при любом <tex>x</tex>, все коэффициенты должны быть равны <tex>0</tex>. Поэтому:
: <tex> p(1) = p(0) </tex>
: <tex> p(2) = p(1) + p(0) </tex>
: <tex> p(4) = p(3) + p(2) </tex>
: <tex> p(5) = p(4) + p(3) - p(0) </tex>
Получаем формулу Эйлера, позволяющую последовательно находить числа <tex>p(n)</tex>:
: <texdpi="145"> p(n) = p(n-1) + p(n-2)+ ... + (-1)^{q+1} (p\left(n - \frac{3q^2 - q}{2}\right) + p\left(n - \frac{3q^2 + q}{2}\right)) </tex>.
Асимптотика <tex> O(n \sqrt{n}) </tex> получается следующим образом. Так как <texdpi="145"> n - \frac{3q^2 + q}{2} \geqslant 0 </tex> , то получаем <tex> q </tex> порядка <tex> \sqrt{n} </tex>, а так как находим <tex>n</tex>-е число, то получаем <tex> O(n \sqrt{n}) </tex>. == Примечания ==<references><ref name="amount">[http://oeis.org/A000041 Последовательность 000041 на OEIS]</ref><ref name="pentagon">[http://ru.wikipedia.org/wiki/%D0%A4%D0%B8%D0%B3%D1%83%D1%80%D0%BD%D1%8B%D0%B5_%D1%87%D0%B8%D1%81%D0%BB%D0%B0#.D0.9F.D1.8F.D1.82.D0.B8.D1.83.D0.B3.D0.BE.D0.BB.D1.8C.D0.BD.D1.8B.D0.B5_.D1.87.D0.B8.D1.81.D0.BB.D0.B0 Википедия {{---}} Пятиугольные числа]</ref><ref name="young">[http://ru.wikipedia.org/wiki/%D0%A0%D0%B0%D0%B7%D0%B1%D0%B8%D0%B5%D0%BD%D0%B8%D0%B5_%D1%87%D0%B8%D1%81%D0%BB%D0%B0#.D0.94.D0.B8.D0.B0.D0.B3.D1.80.D0.B0.D0.BC.D0.BC.D1.8B_.D0.AE.D0.BD.D0.B3.D0.B0 Википедия {{---}} Диаграммы Юнга]</ref></references>
== Источники информации ==
* [http://en.wikipedia.org/wiki/Pentagonal_number_theorem Wikipedia {{---}} Pentagonal number theorem]
* [http://kvant.mirror1.mccme.ru/1988/11/razbienie_chisel.htm Вайнштейн Ф., Разбиение чисел. Журнал "Квант" № 11, 1988 год]
* Н.Я. Виленкин "Комбинаторика", 2007. стр. 138-141.
[[Категория: Дискретная математика и алгоритмы]]
[[Категория: Комбинаторика]]
