Производящие функции нескольких переменных

Обычные Производящие функции хранили одну последовательность чисел. Для задач, которые требуют больше информации применяются производящие функции нескольких переменных (англ. multivariable generating function). Очень часто применяются функции от двух переменных(далее они и будут рассматриваться), которые в общем случае принимают вид

[math]G(x, y) = \sum\limits_{n, k \geqslant 0} g_{n, k} x^n y^k[/math]

Такие последовательности удобно задавать в виде треугольника. Рассмотрим несколько примеров :

Треугольник Паскаля

Рис.[math]1[/math]

Треугольник Паскаля (рис.[math]1[/math]). Элементы треугольника перечисляют пути, идущие из его вершины в соответствующую клетку. Пути имеют вид ломаных, составленных из векторов единичной длины двух видов: идущих вправо-вниз и идущих влево-вниз.

Теорема:

Числа, стоящие в треугольнике - биномиальные коэффициенты [math] c_{n,k} = \begin{pmatrix} n \\ k \end{pmatrix}[/math].

Доказательство:

[math]\triangleright[/math]

Числа, стоящие в треугольнике - биномиальные коэффициенты [math] c_{n,k} = \begin{pmatrix} n \\ k \end{pmatrix}[/math]. Кратко докажем (индукция по n) . Предположим, что числа в [math]n[/math]-ой строчке треугольника совпадают с коэффициентами разложения многочлена [math](1 + s)^{n}[/math]. число различных путей, ведущих в точку [math](n + 1, k)[/math], равно сумме числа путей, ведущих в точку [math](n, k - 1)[/math], и числа путей, ведущих в точку [math](n, k)[/math] , [math]c_{n+1,k} = c_{n,k-1} + c_{n,k}[/math]. Поэтому число [math]c_{n+1, k}[/math] совпадает с коэффициентом при [math]s^k[/math] в многочлене [math](1+s) \cdot (1+s)^n = (1+s)^{n+1}[/math].

[math]\triangleleft[/math]

Производящая функция может быть сопоставлена треугольнику Паскаля несколькими способами. Например, можно рассмотреть производящую функцию

[math]\sum\limits_{n,k = 0}^{\infty} c_{n,k} x^k y^n = \sum\limits_{n,k = 0}^{\infty} \begin{pmatrix} n \\ k \end{pmatrix} x^k y^n = \sum\limits_{n = 0}^{\infty}\Big(\sum\limits_{k = 0}^{n} \begin{pmatrix} n \\ k \end{pmatrix} x^k\Big) y^n = \sum\limits_{n = 0}^{\infty} (1 + x)^n y^n = \dfrac{1}{1 - y - xy}[/math]

Рис.[math]2[/math]

Второй способ соответсвует нумерации элементов треугольника числом отрезков каждого типа на путях, ведущих в соответствующую точку (рис.[math]2[/math]) [math]C_{n,m} = c_{n+m, n} = \begin{pmatrix} n + m \\ m \end{pmatrix}[/math]. Тогда производящая функция будет иметь вид

[math]\sum\limits_{n,m = 0}^{\infty} C_{n, m} x^n y^m = \sum\limits_{n,m = 0}^{\infty} \begin{pmatrix} n + m \\ m \end{pmatrix} x^n y^m = \sum\limits_{k = 0}^{\infty} \sum\limits_{n + m = k} \begin{pmatrix} n + m \\ n \end{pmatrix} x^n y^m = \sum\limits_{k = 0}^{\infty} (x + y)^k = \dfrac{1}{1 -x - y}[/math]

Также существует еще один способ: сопоставить треугольнику Паскаля экспоненциальную производящую функцию. Экспоненциальная производящая функция отличается от обычной тем, что в качестве коэффициентов степенного ряда берутся не элементы последовательности [math]a_n[/math], а числа [math]\dfrac{a_n}{n!}[/math]

Экспоненциальные производящие функции

Зафиксируем произвольную последовательность [math]\{ \alpha_n \}[/math]. Каждой последовательности [math]\{ \alpha_n \}[/math] мы можем сопоставить производящую функцию

[math]\{ \alpha_n \} \mapsto \sum\limits_{n = 0}^{\infty} a_n \alpha_n s^n[/math]

определяемую последовательностью [math]\{ \alpha_n \}[/math]. Если в последовательности [math]\{ \alpha_n \}[/math] отсутствуют нулевые элементы, то такое сопоставление взаимно однозначно. До сих пор мы пользовались обычными производящими функциями - отвечающими последовательности [math]\{ \alpha_n \} \equiv 1[/math]. В зависимости от преследуемых целей целей могут принести и другие последовательности .

Определение:

Производящие экспоненциальные функции(англ. exponential generating function) - функции, соответсвующие последовательности [math]\{ \alpha_n \} = \dfrac{1}{n!}[/math].

Определение:

Экспоненциальные производящие функции для целочисленных последовательностей называют функциями Гурвица(англ. Hurwitz function).

Чем отличаются экспоненциальные производящие функции от обычных? Посмотрим на поведение экспоненциальных производящих функций при выполнении операции над ними. Сумма ведет себя обычным образом:

[math]\sum\limits_{n = 0}^{\infty} \dfrac{a_n}{n!}s^n + \sum\limits_{n = 0}^{\infty} \dfrac{b+n}{n!}s^n = \sum\limits_{n = 0}^{\infty} \dfrac{(a_n + b_n)}{n!} s^n[/math]

а с произведением по-другому :

[math]\bigg(\dfrac{a_0}{0!} + \dfrac{a_1}{1!}s + \dfrac{a_2}{2!}s^2 + \ldots\bigg)\bigg(\dfrac{b_0}{0!} + \dfrac{b_1}{1!}s + \dfrac{b_2}{2!}s^2 + \ldots\bigg) = \dfrac{a_0}{0!} \dfrac{b_0}{0!} + \bigg(\dfrac{a_0}{0!} \dfrac{b_1}{1!} + \dfrac{a_1}{1!} \dfrac{b_0}{0!}\bigg)s + \bigg(\dfrac{a_0}{0!} \dfrac{b_2}{2!} + \dfrac{a_1}{1!} \dfrac{b_1}{1!} + \dfrac{a_2}{2!} \dfrac{b_0}{0!}\bigg)s^2 + \ldots[/math]

Коэффициенты [math]\dfrac{c_n}{n!}[/math] произведения вычисляются по формуле

[math]c_n = \begin{pmatrix} n \\ 0 \end{pmatrix} a_0 b_n + \begin{pmatrix} n \\ 1 \end{pmatrix} a_1 b_{n - 1} + \ldots + \begin{pmatrix} n \\ n \end{pmatrix} a_0 b_0[/math]

Еще одно существенное отличие экспоненциальных производящих функций от обычных наблюдается при взятии производных (и интегрировании). Дифференцирование или интегрирование экспоненциальной производящей функции приводит к сдвигу последовательности ее коэффициентов без изменения их велечины :

[math]\bigg( \dfrac{a_0}{0!} + \dfrac{a_1}{1!}s + \dfrac{a_2}{2!}s^2 + \ldots \bigg)^{'} = \dfrac{a_1}{0!} + \dfrac{a_2}{1!}s + \dfrac{a_3}{2!}s^2 + \ldots [/math]

[math]\int \bigg(\dfrac{a_0}{0!} + \dfrac{a_1}{1!}s + \dfrac{a_2}{2!}s^2 + \ldots \bigg) = \dfrac{a_0}{1!}s + \dfrac{a_1}{2!}s^2 + \dfrac{a_2}{3!}s^3 + \dfrac{a_3}{4!}s^4 + \ldots [/math]

Обычная производящая функция [math]A(s) = a_0 + a_1s + a_2s^2 + \ldots[/math] выражается через экспоненциальную [math]\mathcal{A} (t) = \dfrac{a_0}{0!} + \dfrac{a_1}{1!}t + \dfrac{a_2}{2!}t^2 + \ldots [/math] по формуле

[math]A(s) = \int_{0}^{\infty} e^{-t} \mathcal{A}(st) dt[/math]

Действительно,

[math]k! = \int_{0}^{\infty} e^{-t}t^kdt[/math]

Теперь можно выписать экспоненциальную производящую функцию для треугольника Паскаля :

[math]\sum\limits_{n, m = 0}^{\infty} \dfrac{1}{(n + m)!} \begin{pmatrix} n + m \\ m \end{pmatrix} x^n y^m = \sum\limits_{n = 0}^{\infty} \dfrac{(x + y)^n}{n!} = e^{x + y}[/math]

Треугольник Дика

Рис.[math]3[/math] Представление треугольника Дика

Треугольник Дика перечисляет пути в положительном квадранте плоскости, выходящие из начала координат и составленные из векторов [math](1, 1)[/math] и [math](1, -1)[/math] (рис.[math]3[/math]). Мы можем увидеть, что элементы [math]d_{ij}[/math] треугольника Дика отличны от нуля в том и только том случае, если [math]i \geqslant j [/math] и [math]i + j[/math] четно. Обозначим через [math]D(x, y)[/math] производящую функцию Дика двух пременных

[math]D(x, y) = \sum\limits_{i,j =0}^{\infty} d_{ij} x^i y^j[/math]

Правило построения треугольника Дика подсказывает нам уравнение на эту производящую функцию

[math]xyD(x, y) + (D(x,y) - D(x,0))\dfrac{x}{y} = D(x, y) - 1[/math]

Действительно, коэффициент при любом мономе [math]x^iy^j[/math] , отличном от единичного, представляет собой сумму коэффициентов при мономах [math]x^{i-1}y^{i-1}[/math] и [math]x^{i-1}y^{j+1}[/math] . Функция [math]D(x, 0) = \dfrac{1 - \sqrt{1 - 4x^2}}{2x^2}[/math] нам известна, и ряд [math]D(x, y)[/math] находится решением линейного уравнения,

[math]D(x, y) = \dfrac{1 - \sqrt{1 - 4x^2} - 2xy}{2x(xy^2 + x - y)}[/math]

См. также

• Производящая функция

Источники информации

• Ландо С.К. Лекции о производящих функциях — МЦНМО, 2007 - С. 57