CC BY
48
удк 517.987
ЗАДАЧИ ПЕРЕЧИСЛЕНИЯ ГРАФОВ С ПОМЕЧЕННЫМИ ВЕРШИНАМИ ENUMERATION PROBLEMS OF GRAPHS WITH LABELED VERTICES
Ю.П. Вирченко, Л.П. Остапенко Yu.P. Virchenko, L.P. Ostapenko
2Белгородский национальный исследовательский университет, Россия, 308015, г.Белгород, ул. Победы, 85 Belgorod National Research University, 85 Pobedy St, Belgorod, 308015, Russia
E-mail: virch @bsu.edu.ru;
Аннотация. Изучаются комбинаторные задачи, связанные с перечислением графов с помеченными вершинами. Они состоят в вычислении таких функций, как Nп и Мп , П G П , где Nп - число всех связных
графов, а м - число всех связных графов без вершин сочленения в зависимости от числа П вершин. На основе специального алгебраического метода дается решение этих задач в терминах производящих функций чисел Nn и
Resume. Combinatorial problems which concerns the enumeration of graphs with labeled vertices are studied. They are made in the calculation of these functions, as Nn well as M n , He , where Nn is the number of
all connected graphs, and Mn is the number of all connected graphs with no articulation vertices depending on the number of vertices. On the basis of a special algebraic method provides a solution to these problems in terms of generating functions of numbers N and M .
Ключевые слова: разложения конечного множества, производящая функция, алгебра симметричных функций, связный граф, граф с помеченными вершинами, вершина сочленения.
Key words: partition of finite set, generation function, algebra of symmetric functions, connected graph, graph with labeled vertices, articulation vertex.
Введение
В настоящей работе изучается комбинаторная задача, связанная с перечислением графов с помеченными вершинами, в зависимости от числа вершин.Здесь и далее, мы используем этот термин или, эквивалентно, термин «помеченные графы», следуя терминологии монографии [1]. В этой монографии, а также в [2],под задачами перечисления графов понимаются такие, решение
которых состоит в описании определенного класса Зп графов c п вершинами, элементыкоторого
объединяются по признаку наличия у них того или иного свойства. В настоящей работе мы занимаемся не самим перечислением графов в указанном смысле, а решением более узкой задачи, которая состоит в определении мощности | §п | некоторых из таких классов графов с фиксированным числом вершин п без конкретного описания элементов, их составляющих. Основной целью нашей работы является определение значений | Зп | как функции от п е □ в том случае, когда класс %
состоит из графов без вершин сочленения (см. [1], [3]). Важность решения такой задачи связана с ее прикладной значимостью, так как связные графы с помеченными вершинами интенсивно используется в равновесной статистической механике, в частности, в статистической механике таких
простейших, с физической точки зрения, систем, какими являются разреженные одноатомные газы (см., например, [4], [5]), где графы такого типа называются графиками Майера, каждому из которых сопоставляется определенное аналитическое выражение. Далее, все графы, о которых идет речь в настоящей работе, являются графами с помеченными вершинами.
Итак, согласно выше сказанному, в настоящей работе мы занимаемся решением чисто
комбинаторной задачи о вычислении арифметической функций | Зп |, п е □ . Как правило, когда
такая функция не «элементарна», то к решению задачи привлекается подход, основанный на применении производящих функций (см., например, [6, 7]). В рамках этого подхода, вместо непосредственного вычисления значений арифметической функции | ^ |, ставится задача о вычислении соответствующей производящей функции ¥(г), аналитической в некоторой окрестности нуля
комплексной плоскости г е □ и подходящим образом определенной для возможности вычисления на основе ее коэффициентов разложения около нуля значений искомой арифметической функции. При таком подходе, целью решения задачи становится установление явного вида функции ¥ (г), и для поиска этого решения становится возможным применение методов классического математического анализа. При этом если даже возникает затруднение при установлении явного вида функции ¥(г) , но имеется возможность формулировки однозначно разрешимой задачи для ее определения и возникает возможность анализа этого решения в окрестности нуля, то можно говорить о «неявном» решении соответствующей задачи определения арифметической функции. В этом смысле, такую переформулировку исходной комбинаторной задачи в терминах производящей функции ¥(г) можно рассматривать как ее решение. Именно такие решения мы изучаем в настоящей работе.
Во втором разделе приводятся необходимые сведения из теории графов и ставится задача, которой посвящена работа. В третьем разделе мы вводим специальные производящие функции, зависящие от дополнительного параметра £, в терминах которых решается поставленная задача. В четвертом разделе мы даем введение в алгебраическую технику, посредством которой далее дается решение основной задачи о вычислении арифметических функций | §п |, где §п - некоторый
класс связных помеченных графов с п вершинами. Далее, в пятом разделе мы приводим примеры применения описанной алгебраической техники для решения комбинаторных задач, связанных с решением основной задачи. В шестом подготовительном разделе выводится рекуррентная формула, с помощью которой в седьмом разделе мы решаем ту основную задачу, которой посвящена работа, а именно, мы получаем функциональное уравнение для производящей функции О(г;£) числа всех помеченных графов с п вершинами без вершин сочленения.
2. Помеченные графы
Пусть имеется множество V элементов, которые мы будем называть далее вершинами и обозначать малыми латинскими буквами. Это предполагает, что вершины, как элементы множества V различимы. Обозначим посредством V(2) множество всех пар (л, у} с V . Помеченнъм графом над множеством V мы называем упорядоченную пару & = (V, ¥), в которой ¥ является
подмножеством из V(2), ¥ с К(2). Множество ¥ называется множеством смежности графа, а его элементы называются ребрами графа (V, ¥) . С точки зрения топологического подхода к понятию графа, определенная нами конструкция соответствует понятию неориентированного простого графа с ребрами одного типа (не содержащего петель и кратных ребер), вершины которого отмечены какими-либо метками (номерами). Наличие меток у вершин графа (V, ¥) приводит к
тому, что какая-то из перестановок номеров вершин из V может приводить к другому графу, отличному от & . Указанное обстоятельство связано с тем, что биекция множества V на себя приводит к соответствующим ей одновременным перестановкам строк и столбцов матрицы смежности графа б . Если такие перестановки не изменяют матрицы смежности, то биекция V не изменяет помеченного графа, в противном случае, она приводит к графу, отличному от исходного. Таким образом, помеченные графы различны тогда и только тогда, когда они имеют различные матрицы смежности, в отличие от графов в обычном смысле, у которых вершины не помечены. Для последних их топологическое различие определяется тем, что их матрицы смежности не могут быть получены друг из друга при какой-либо одновременной перестановке строк и столбцов.Далее, мы будем рассматривать только случай, когда множество V конечно.
Данное выше определение помеченного графа предполагает, что два графа &1 = (V, ¥х) и
&2 = (V, ¥2) совпадают в том и только в том случае, когда V = V и ¥ = ¥2 . Отсюда, в частности, следует, что два графа над одним и тем же множеством V различны, если они отличаются своими множествами смежности ¥ и ¥2 .
Множество смежности ¥ с V(2) графа б определяет бинарное симметричное отношение смежности на V . Каждое такое бинарное отношение порождает бинарное отношение связности на V, которое обладает свойствами симметричности и транзитивности. Оно конструируется на основе понятия пути на графе б = (V, ¥). Путем на графе б = (V, ¥) называется каждая последовательность (х,х1,х2,...,хп_х,у) вершин из V такая, что (х.,х.+,}е¥, j = 0,1,...,п — 1, х0 = х, Хп = у. При этом пара (х, у} с V является связной на графе б , если существует путь, состоящий
из вершин, которые входят в состав V. Получаемое таким образом подмножество пар из V(2) определяет бинарное отношение. Если же это подмножество связных пар (х, у) с х Ф у дополнить парами (х, х), х е V, то отношение связности становится рефлексивным. Тогда такое отношение превращается в отношение эквивалентности.
Так как каждое отношение эквивалентности на V разбивает это множество на непересекающиеся подмножества эквивалентных между собой вершин, то отношение связности разбивает множество V на непересекающиеся подмножества (компоненты) связных между собой вершин.
Граф (V, ¥) называется связным, если у него, при указанном разбиении, имеется только одно подмножество эквивалентности, то есть он состоит из одной связной компоненты. Если в
указанном разбиении имеется более одной компоненты Ух, У2, ..., Ут так, что имеет место дизъ-
т
юнктивное разложение У} = У, то такое разбиение индуцирует разбиение на компоненты
7=1
¥2, ..., ¥т , т > 1 множества смежности ¥ так, что
т
хр = ихр/' ^при
7=1
и ^ у(2), j = 1 + т . В этом смысле можно говорить, что граф б состоит из несвязанных между собой графов б . = <у, ¥.), j = 1 ^ т, каждый из которых называется связной компонентой исходного графа.
Путь (X,x1,х2,...,хп_1,у) на графе с п > 3 такой, что X = у, называется циклом. Если на графе (У, ¥) отсутствуют циклы, то такой граф называется древесным.
Граф б ' = (У',¥ '), у которого У' ^ У, ¥ ' = ¥п>У'(2) называется подграфом графа б . При этом если граф б - связный, то граф б' не обязательно является связным.
Пусть & = (У, ¥х) и б 2 = (У2, ¥2 ) - связные подграфы связного графа б = (У, ¥), причем
У0У2 = {х}, У = УиУ2; ¥ = ¥хп¥2=0, ¥ = ¥хи¥2. (1)
Тогда X называется вершиной сочленения графа б , а сам этот граф называется графом скленным по этой вершине. При этом граф б , по определению, является склейкой графов бх и б2.
Если вершина X является вершиной сочленения в графе б , и этот граф представляется в виде склейки графов б . = У, ¥. ), где j = 1 ^ р,
б = У б ] (2)
j=1
таких, что
Уj = {х}, ¥у ^¥к =0, j Ф к; j, к = 1 р,
то число р называется степенью вершины сочленения X.
Если вершина X не является вершиной сочленения то, по определению, будем считать, что ее степень равна 1.
Если граф не обладает вершинами сочленения, то его, в дальнейшем, будем называть блоком. По определению, любой одновершинный граф {X} не имеет вершины сочленения.
Если граф б содержит такой подграф б ' = (У', ¥ ' ), что в У' существует подмножество вершин А ^ У', являющихся вершинами сочленения в б , и набор непустых связных попарно несвязных друг с другом под графов б(£) = (У (£), ¥(х)), 7 е А такой, что возникает биекция между элементами набора и вершинами из А , и при этом граф б представим в виде склейки
б = б^ [уб (г)],
то подграф б' будем называть блоком в графе б .
Очевидно, что если подграф б' = (V', ¥') графа б = (V, ¥) является блоком, то вершины
множества V" = V V' являются вершинами сочленения в графе б и степень каждой из вершин этого множества не превосходит числа | V" | — | V" |.
Вершина сочленения в связном графе б характеризуется следующим свойством.
Теорема 1. Для того чтобы вершина х связного графа б = (V, ¥) была вершиной сочленения, необходимо и достаточно чтобы существовала такая пара вершин у1 е V, у2 еУ,
у. Ф х, j = 1,2, для которой любой путь у(у1, у2) из у1 в у2 обязательно содержит вершину х.
Ш Необходимость. Если имеет место разложение (1) графа б на два подграфа бх , ¥) и б 2 = (V, ¥2) , то для любых двух вершин у е V и у2 еV2, отличных от х, во всяком пути (ух,...,у2) содержится вершина х. В самом деле, рассмотрим произвольный путь у = (у1,...,у2).
Найдем в этом пути вершину г , которая является первой из не принадлежащих V вершин, и допустим, что вершина V, предшествующая г в пути у, не совпадает с х. Тогда, с одной стороны, связь (V, г} е¥, а, с другой стороны, эта связь не может принадлежать ни одному из множеств ¥ ■, j = 1,2 .В самом деле, по построению, V е V2 и, следовательно, (V, г} е ¥2. Кроме того, г е V и поэтому (V, г} её ¥. В результате, получаем противоречие с равенством ¥ = ¥х и ¥ . Достаточность. Пусть у , j = 1,2 - указанные в условии теоремы вершины графа б . Нужно дока-
у2 = (V , V) и(х} и ¥ = ((", V} е¥:(м, V} с V}, j = 1,2. Тогда, по построению V uУ2 = V , У1^У2 = (х}, ¥х п¥2 =0. Покажем, что ¥ и ¥2 = ¥. Допустим противное, что существует связь (^, w2 } её ¥ ■, j = 1,2 и , } е ¥. Это возможно только в том случае, когда wj, содержатся в разных множествах и не совпадают с х .В силу определения V , существует путь у(у, w1)
зать существование разложения (1). Определим
и
на графе б , полностью расположенный в V . Если пара , w2} содержится в ¥ , то существует путь у(у, ) v( w1, w2), который не проходит через вершину х, что противоречит тому, что её V . Определив графы б . = (V, ^^), j = 1,2, получим, что бх V б2 и склейка происходит по вершине х. +
Докажем, что для любого связного графа имеет место разложение (2).
Теорема 2. Пусть X - вершина сочленения связного графа б = (У, ¥). Тогда существует число $ > 2 и однозначным образом определяемый набор связных графов б. = (У, ¥■),
j = 1 ^ $, в каждом из которых вершина X уже не является вершиной сочленения, и при этом имеют место следующие соотношения
У = У иУ2 и... иУ$, У ПУ;. = {X}, I Ф i,j = $;
¥=¥и¥2и .., ¥,.п¥,.= {X},iфj;i,j = $
и Построим указанное в условиях теоремы разложение. Так как X - вершина сочленения графа б , то имеется два графа б(Х) =(У/(1), ¥(1)), j = 1,2 таких, что У(1) и У2(1) = У, У(1) п У2(1) = {X} и
¥<ц = 0, п¥(1) = ¥.
Если в графах б(1), j = 1,2 вершина X не является вершиной сочленения, то построение закончено и, согласно определению, $ = 2. Если же, по крайней мере, в одном из этих графов вершина X является вершиной сочленения, то для каждого такого графа c номером у , равным 1
или 2, существует пара графов б=(У,(2/, ¥(2у), j = 1,2 таких, что У(2) и у(22 = У^1, У™ п У^ = {X} и ¥(2) п¥(2)2 = 0, ¥(2) и¥(2'2 = ¥(1). Если же у графа б(1) вершина X с номером у не является вершиной сочленения, то графу б сопоставим пару графов б ^ = би б¿2 = ({X}, 0) . Таким образом, для каждого у = 1,2 имеется пара графов б, у2 = 1,2.
Далее, если у вновь построенных графов ¥ . ■ ; у'1, у'2 = 1,2 вершина X не является вершиной сочленения, то построение заканчивается и, по определению, $ равно числу тех графов б (2^ , у которых множество связности ■ не пусто. Если же среди всех вновь построенных графов найдутся такие, у которых вершина X является вершиной сочленения, то построение продолжается. Каждому графу б(2;, обладающему вершиной сочленения, сопоставим пару графов
б у =У(Хл ) • у" 1.2 таких, что V. и ^ = V. , V. п Г»; = {X} и
¥(3)лД п¥Ъ,2 = 0 , ¥(3)ЛД и¥(3^2 = ¥(2^ . Для тех же графов б(2)2, у которых вершина X не
является вершиной сочленения, положим б(3^ 1 = б(2^ . и б(3^ 2 = ({X}, 0) . Если у всех построенных графов вершина X не является вершиной сочленения, процедура построения новых графов останавливается. Если же это не так, то переходим к новому шагу построения.
Описанный алгоритм построения все более мелких графов, у которых вершина X является вершиной сочленения, должен завершиться на каком - то шаге т , так как множество вершин исходного графа б конечно. При завершении описанного процесса построения набора подграфов
исходного графа б, выберем из всех полученных графов б(т) , = (У(m) , ¥(т) . ); 7, = 1,2,
к = 1 * т те, у которых ¥(т) , Ф 0 . Пусть число таких графов равно 5 . В результате, получен
набор неодновершинных графов, каждый из которых содержит вершину х , которая у каждого из них не является вершиной сочленения. Обозначим эти графы, занумеровав их произвольным образом, посредством б(1) = (V(1), ¥(1)), ..., б(5) = (V(5), ¥(5)) . Для этого набора графов, по построению, выполняется
V(1) uУ(2) и... иУ(5) = V, ¥(1) и¥(2) и... и¥(5) =¥; (4)
V^ п V(к) = (х}, ¥°') п¥(к) =0 при j Ф к; j,к = 1 *5. (5)
Покажем, что набор графов б(1), ..., б(5) выбирается однозначным образом. Допустим противное, что имеется два разных набора графов, удовлетворяющих (4) и (5): уже указанный набор
б(Л = (V^¥(])) , j = 1,...,^ и набор графов б^' =(У(J)',¥(Л'), j = 1,...,5 ', удовлетворяющий аналогичным условиям, а именно
V(1) ' uУ(2) ' и... иУ(5 ' ) ' = V, ¥ ( 1 ) ' и ¥( 2 ) ' и. . . и ¥ (5 ' ) ' = ¥ ; 4 '
VО') ' п V(к) ' = (х}, ¥(^) ' п ¥(к) ' =0 при j Ф к;!, к = 1 * 5'. 5 '
Не ограничивая общности будем считать, что 5 < 5 '. В противном случае, поменяем штрихованный и нештрихованный наборы графов ролями в нижеследующих рассуждениях.
Возьмем любую вершину у из V(1), отличную от х, и найдем номер того графа среди б(]) , j = 1 * 5 ', в котором она содержится. Не ограничивая общности, можно считать, что таковым является граф б(1) . Тогда V(1) п V(1) з (х, у1}. Допустим, что V(1) Ф V(1) . Следовательно,
(V(1), V(1) ') и (V(1) ', V(1)) Ф0
. Положим для определенности, что не пуста первая компонента объединения и г е V(1) , V(1) . Тогда, так как объединение всех множеств V(^), j = 1 * 5 совпадает с объединением всех множеств V(') , j = 1 * 5 ', то найдетсятакой номер к Ф 1, для которого г е V(^ . Так как г1 и у принадлежат одному и тому же множеству V(1), то существует путь у(у1, г ), расположенный в V и не проходящий через вершину х. Но это противоречит тому, что у и г содержатся в различных множествах V(1) и V(к1) , склеенных в вершине х. Это противоречие указывает на то, что У(1) = У(1) . Удалив эти блоки V(1) и V(1) , каждый из своего набора, применим всю описанную для блока V(1) процедуру к блоку V(2), то есть найдем вершину у2 Ф х,
у2 е V; найдем по этой вершине соответствующий граф в наборе штрихованных графов и обозначим его посредством б(2) . В результате аналогичных рассуждений покажем, что У(2) = У(2) . Продолжим описанную процедуру вплоть до исчерпания всего набора графов б(11, j = 1 * 5 . Та-
ким образом, в результате проведенных рассуждений мы установили, что каждый из графов б(]) , j = 1 * 5 совпадает с соответствующим графом б(1) . Так как при этом исчерпается весь запас вершин из множества V, то 5 = 5 ' и оба набора графов совпадают.
Теорема 3. Пусть х - вершина графа б = (V, ¥), не являющаяся вершиной сочлене-
ния.Тогда в графе б найдется единственный блок б = (V', ¥ ' ), в котором она содержится.
□ Пусть V - множество вершин сочленения графа б их е V - вершина, о которой идет речь в
формулировке теоремы. Если V = 0, то б - блок и в этом случае теорема доказана. Если V Ф 0, то выделим тот блок, о котором идет речь в формулировке теоремы, посредством алгоритма дробления
графа б . Опишем первый шаг этого алгоритма. Возьмем первую из вершин у со степенью 5 , принадлежащих V . Тогда представим граф б в виде склейки б = б(1) V б(1) двух графов
б(1) = (У1, ¥) и б (1) =^2, ¥2) так, что V пУ2 = (у1}, V = V иУ2. Выберем из этих двух графов тот, в котором содержится вершина х . Так как х не является вершиной сочленения, то только один граф из представленных двух графов обладает таким свойством. Положим, что этим свойством обладает граф б(1) .При этом множество V является собственным подмножеством V . При этом вершина сочленения у в графе б(1) обладает меньшей степенью чем 5 . Мы описали первый шаг алгоритма дробления. Если граф б(1) является блоком, то есть не имеет вершин сочленения, то теорема доказана. Если же он не является блоком, то перейдем к следующему шагу дробления.
Применим описанный алгоритм к графу б(1). При этом мы будем использовать вершину сочленения у , если она осталасьтаковой в графе б(1), либо возьмем новую вершину сочленения
у2 со степенью 52, которая содержится в графе б(1) В результате применения алгоритма дробления, получим новый граф б (2) , в котором содержится вершина х и в котороммножество его вершин сочленения V2 является собственным подмножестсвом множества V и степень вершины сочленения у2 в нем меньше 52, причем граф б(1) представим в виде склейки б(1) = б(2) V б(2) со связным графом б(2) .Таким образом, исходный граф б представим в виде склейки б = б(1) V б(2) V б(1).
Продолжим применение алгоритма дробления до тех пор пока на каком-то шаге вновь полученный граф б(^ не имеетвершин сочленения, однако при этом он содержит вершину х. Кроме того, так как имеет место б = б) V (б ^ 1) V... V б(1)). +
3. Алгебра А коэффициентов производящих функций
Настоящая работа посвящена решению комбинаторной задачи о вычислении функции Мп , значением которой при каждом п е □ является число всех помеченных графов без вершин сочленения с числом вершин п е □ , а также связанной с ней задачей вычисления чисел Nп, п е □ ,
равных при каждом йеО числу всех связных помеченных графов с числом вершин п .
При решении комбинаторной задачи очень редко удается найти явный вид функции нату-рального(ых) аргумента(ов), представляющие такие решения. Чаще представляется возможность найти рекуррентную(ые) формулу(ы), определяющую(ие) значения этой функции последовательно.
Именно такая ситуация реализуется при вычислении функции , п е Ш .
В более сложных ситуация, когда невозможно найти явные рекуррентные формулы, имеется возможность найти решение задачи в виде явного выражения для подходящим образом определенной производящей функции (см., например, [8-10]), которая(ые) является(ются) аналитическими от комплексного(ых)переменного(ых) 2 в окрестности нуля комплексной плоскости их изменения и при этом решение исходной комбинаторной задачи выражается через коэффициенты их степенных рядов. Наконец, в еще более сложной ситуации, явные формулы для производящих функций не удается найти на основе элементарных аналитических функций, а возможно только лишь найти уравнение, которому они подчиняются и являются при подходящих условиях его единственным решением. Тогда найденное уравнение можно рассматривать в качестве искомого решения комбинаторной задачи. Именно такая ситуация реализуется для комбинаторной функции Мп , П £ □ .
Далее, мы будем рассматривать производящие функции вида
œ n
z
¥д (2) = £-К!
п=о п!
чисел ! §п ! элементов в классах §п .Для вычисления таких производящих функций, а также для
установления уравнений связывающих производящие функции, которые соответствуют связанным друг с другом комбинаторным функциям, мы применим специальный алгебраический метод. Этот метод основан на применениибесконечномерной, коммутативной алгебре с единицей для последовательностей коэффициентов функций ^ (2) •
Обозначим посредством /п ={1,2,...,и}- стандартное п -элементное множество, йеО . Пусть Ь линейное пространство бесконечных последовательностей 3 = (ат;/йбО + ), в которых каждая компонента при любом теП + является симметричной измеримой функцией ат на [ОД]'" со значениями в □ , ат(Хт) = ат(х 1,х2,...,хт), Хт = {ху,х2,...,хт) е [0,1]'". Для значения т = 0 такие функции я0 £ [0,1], полагаются константами. Линейное пространство Ь превращается в коммутативную алгебру А при определении на нем следующей бинарной операции * , ко-
торую будем называть умножением. Каждой паре В = {ат; т е □ + ) и Ь = {Ьт; т е □ + ) последова-тельностейиз А сопоставим последовательность а* Ь - результат умножения последовательностей а и Ь , компоненты которой определяются формулой
(а*Ь)„(X(1п)) = ^ а (X(А))Ъ (X(1п, А)), (6)
А х 1п. А
Ас/,
в которой использованы следующие обозначения:
X(А) = (х1, х2,. ., х ) для каждого А = (/1, /2,..., ¡5} с 1п , где | А | = 5 - число элементов в А так, что Хп = Х(1п). Очевидно, что функции (а*ь)п(х1,х2,...,х;)) симметричны при любом йёО и
результатом применения операции * к двум последовательностям а и Ь снова является последовательность симметричных функций. Точно также, очевидным образом, операция * , действительно, является умножением, так как она обладает всеми свойствами, которые предъявляются к операции умножения в общеалгебраическом смысле (см, например, [11]), так как она согласована с операциями линейного пространства. Введенная операция умножения ассоциативна и коммутативна.
Таким образом, линейное пространство А , снабженное операцией * является бесконечномерной коммутативной алгеброй. В этой алгебре имеется единица, которая представляет последовательность е = (1,0,0,...) .
Рассмотрим в алгебре А линейное подпространство А , которое состоит из таких последовательностей а = (ат (Хт); т е □ + ), у которых а0 = 0. Это линейное пространство представляет собой максимальный идеал алгебры А , то есть оно замкнуто относительно умножения и не содержится ни в каком другом множестве, обладающем таким свойством.
Напомним, что разложением стандартного множества !п называется неупорядоченный набор А непустых множеств А., j = 1 * 5, которые составляют дизъюнктивное разложение множества /„:
слАи Ф0, /фАг; /', А: = 1-^5,
1=1
где 5 - число элементов разложения (мощность). Класс всех разложений множества 1п обозначим посредством S п.
Точно также обоначим класс всех разложений А , каждое из которых содержит m подмножеств множества !п , обозначим посредствомS (m). Для каждого ае А0 рассмотрим степень
а"' с любым иеО . При этом, по определению, положим а'! = е . Тогда компоненты (а™)и последовательности аm равны нулю при п > m . Если же п < m., то эти компоненты даются формулой
(аШ)я(Х(1п)) = ш\ V П«^(А)).
А
А&у (ш) АеА
Определим, теперь, на алгебре А экспоненциальную функцию:
ад ^
ехр* а = V— аш.
ш=0 Ш\
Принимая во внимание замечание относительно степеней элементов из А, находим, что для а Е А0 имеет место формула
[ехр* а]п(X(1п)) = [V — аШ](Х(1Я)) = V П«и (X(А)), (7)
ш=0 Ш \ Ае5 п АеА
где суммирование производится по всем разложениям множества !п с любым допустимым числом ^ = 1,..., п компонент.
Определим, теперь, на алгебре А линейный функционал 1(•) , зависящий от параметра
г е□ . Для каждого элемента Э = {ат,т еО }еА значение /.(г;а) этого функционала определяется формулой
ад Я
1(2; а) = £ —. | «п (Хп)dXn, (8)
<? п \
п=0 [0,1]п
где dXn = dx]dx2...оХи - мера Лебега на [0,1]п . Область определения функционала .(•;•) содержит все последовательности а, в которых все функции равномерно ограничены и при этом его значения конечны при всех г е □ .
Теорема 4.Функционал I мультипликативен, то есть для каждой пары таких последовательностей а = (ат; т е □ ) и Ь = (Ьт ;теП) из А, у которых | /.(г; а) |< оо, | /.(г; Ь) |< оо, имеет место равенство
1(2; а * Ь) = 1(2; а).(г;Ь). (9)
Ш Доказательство проводится приводимым ниже прямым вычислением. Согласно определению (8),
ад п
Цг;а*Ь) = £ —. | (а*ь)п))dXn =
<? п \
п=0 [0,1]п
ад п
V п / V ^' ^^=
п=0 п \ [0 1]п ю^п
ад п
V "7 V / «а(XА^,/X^ ' ^^ .
п=0 п\ Ас 1п [0,1]
После замены переменных интегрирования X(А) , |А |= ш на X(1т) и X(/и , А) на X(!„ , ) в каждом из слагаемых и используя правило перестройки суммирований
е ..., Е 1 =
Лс 1п т=0 Лс1п :\Л\=т Лс1п :\Л\=т
ГпЛ
V т у
получаем
ад п
¿(-; а*Ь) ^ Е I «т X )Ь-т (^ (^п , 4 )Жп =
п=0 п ! Л с 1п [0
" ГП п (П1 г
= Е - Е I | ^^ Х Ж-т (X(1п , 1т )№ =
п=0 п • т=0 Vту [0 1]п
ад т ад п-ти
Е — «» (хи Ж(/и) { Е 7-X(1П , 4 )^Х(/п-т^ =
т=0 т ! [0,1]т [0,1]п-т п=т (п т) !
ад _т
= (Е — I (Хт )ЛХ{/т )) • ( \ Е — (Х(/п ))^Х(/п)) = ¿(-; а)/(—; Ь).
т=0 т ![0,1Г [0,1Г п=п!
Перестановки порядка суммирований обосновывается абсолютной сходимостью рядов для значений ¿(-; а), /(-; Ь) функционала в некотором круге ненулевого радиуса в комплекскной
плоскости - . +
Следствие.Еслг/ последовательность Э = (ат; т е □ ) е А0 г/ состоит из равномерно
ограниченньх в совокупности функций, то
I (г; ехр* а) = ехр(/. (г; а)), (10)
где функция комплексного переменного - в правой части является аналитической в окрестности нуля.
и Формула (10) является следствием линейности и мультипликативности функционала /(-; •) +
Продемонстрируем эффективность использования алгебры А и с вязанных с ней понятий на примере решения комбинаторной задачи о числе Сп всех разложений множества /п с фиксированным значением п е □ . При малых значениях п = 1,2,3 соответствующее число Сп легко устанавливается перебором всех разложений:
1. При п = 1, / = {1} имеется только одно разложение с 5 = 1, Л1 = / = {1} . С = 1.
2. При п = 2, /2 = {1,2} имеется два разложения с 5 = 1, Л1 = /2 и с 5 = 2, Д = {1}, Л2 = {2} .
С2 = 2.
3. При п = 3, /3 = {1,2,3} имеется пять разложений с 5 = 1, Л1 = /3 и с 5 = 2 : Лх = {1}, Л2 = {2,3} >
Л, = {1,2}, Л2 = {3}; Л, = {1,3}, Л2 = {2}; с 5 = 3, Л = {1}, Л2 = {2}, Лз = {3}. С3 = 5. Для вычисления числа С для любого значения п е заметим, что
с = Е
Аеу „
Поэтому, полагая в (7) ап (X (/и )) = 1, получим, что
С„=(ехр»а)„(Х(/„)), а = <0,1,1,....), ие0+ (и)
для любой точки Хи е [0,1]п , где последовательность а является элементом идеала А • Эта формула является основой для вычисления производящей функции О (г) чисел Сп
ОД = £ -ся, (12)
П~0> п!
в которой полается, по определению, С0 = 0 • На основании Следствия Теоремы 4, находим, что справедлива
Теорема 5. Производящая функция (12) является аналитической функцией, определяемой в достаточно малой окрестности нуля на вещественной оси г, радиус которой определяется наименьшим решением уравнения
О(г) = ехр(г -1). (13)
□ Так как последовательность Э = (ап(Х(/;)));яе0 + ), где ап(Х(1п)) = 1, п е I , находится в А0, то
ад „п
г
0( г) = Ц г; а) = £ - = ег-1 п!
и, согласно формулам (10) и (11),
I.(г; ехр* а) = ехр (/.(г; а)) = ехр (ег -1). +
Произведем расчет первых пяти значений Сп на основе формулы (13), где Си = Оп)(0). Заметим, что О") (г) = О(г)Рп (ег), где Рп (х) - полином п -й степени, где Р (х) = х • При этом полиномы Рп, п е вычисляются рекуррентно по формуле
Рп+1( х) = х( Р( х) + Рп (х)), что проверяется непосредственным дифференцированием
0(п+1) = Я(г)Ря (ег) + Щг)Р, (ег )ег = ОД^ (ег) + Рп. (ег ))ег. Так как 0(0) = 1 и С = 0(п)(0) = 0(0)^(1), Р+1(1) = Ря(1) + Р(1), то
Р( х) = х С = 1;
Р (х) = х + х2 С = 2;
Р (х) = х + 3х2 + х3 С = 5;
Р (х) = х + 7 х2 + 6 х3 + х4, С = 15;
Р (х) = х + 15х2 + 25х3 +10 х4 + х5, С = 52;
то есть эти значения совпадают со значениями С , С2, С , вычисленными выше перебором имеющихся вариантов.
4. Применение алгебраического метода
В этом пункте мы продемонстрируем применение алгебраической техники, описанной в предыдущем пункте, на примере вычисления производящей функции для числа Р всех связных древесных графов с п помеченными вершинами(см.[9]). в простейших случаях при п = 1,2,3,4 числа р даются следующими значениями.
1. При п = 1, V = {1} имеется только один граф ({1}, 0) , состоящий из единственной вершины 1. Он, по определению, является связным и древесным, то есть Р = 1.
2. При п = 2, V = {1,2} имеется только один связный граф ({1,2}, ¥ = {1,2}0), который, из-за наличия только двух вершин, не содержит циклов. Следовательно, он является древесным, то есть Р = 1.
3. При п = 3, V = {1,2,3} имеется три связных графа без циклов, соответственно, со следующими множествами смежности: ¥ = {{1,2},{2,3}} , ¥ = {{1,2},{1,3}}, ¥ = {{1,3},{2,3}} . Все они различны как графы с помеченными вершинами. Следовательно, Р = 3 .
4. При п = 4, V = {1,2,3,4} имеется 16 различных древесных графов. Соответствующие им множества смежности которых мы распределим на две группы топологически совпадающих друг с другом графов в каждой из них:
4 графа вида ¥ = {{/, '}; ' = 4, ' Ф /}, где I = 1 ^ 4;
12 графов вида ¥ = {{/, к},{¡},{к, ¡}}, где {/, 7} с {1,2,3,4}, {к, 1} с {1,2,3,4}, {/,'} о{к, ¡} = 0;
Следовательно, Р = 16 .
Теорема 6. Для функции Ьп, п е имеет место следующая рекуррентная формула
4+1=2 иеП. (14)
АеУ п ЛеА
□ Для каждого древесного графа б = (V, ¥) над множеством вершин V = /и+1 выполним следующее построение. Сопоставим графу б совокупность древесных графов б , ' = 1,..., ^, получаю-
щихся удалением из него вершины (п +1) вместе со всеми ребрами, инцидентными этой вершине, то есть сопоставляется, вообще говоря, несвязный граф (V , {п +1}, ¥ , {х е V, {п +1}: {х, п +1} е ¥}) .Связные компонентах этого графа представляют совокупность графов бJ, j = 1,...,5 . Каждая из этих связных компонент является древесным графом, так как если бы в какой-то из них имелся цикл, то он должен присутствовать в исходном графе. Каждая из связных компонент образована некоторым подмножеством вершин А. ^ /и и
множеством ребер ¥у.
Отметим вершины х. е А, j = 1 ^ 5 такие, что ребро (ху, п +1} содержится в ¥ , которое удаляется вместе с удалением вершины (п +1) .
Для всей совокупности множеств {А.; j = 1 ^ 5} и связанных с ними ребер {¥^; j = 1 ^ 5} выполняется
3=1 7=1
Прик ^ I; к, I = 5.
Таким образом, каждомуграфу б , на основе описанной конструкции, сопоставляется, сначала, разложение А еЗ „., а затем, на основе этого разложения, - множество отмеченных вершин {х. е А. : А. е А, j = 1 ^ 5}. и набор древесных графов б . = (А., Ф.), j = 1 ^ 5. Они характеризуют граф б, то есть, при фиксации указанных математических объектов, граф б однозначно восстанавливается, что достигается присоединением к вершинам /п вершины (п +1) и ребер
{х], п +1}, / = 1 л - к совокупности ребер Ф .
У=1
В связи с доказанной характеризацией каждого из деревьев на множестве вершин V = /и+1, распределим все деревья над V по классам эквивалентности^ (А ), занумерованным разложениями А еЗ „ так, что каждому древесному графу из фиксированного класса З (А ) соответствует одно и то же разложение А еЗ п . Тогда
4+1 =£ З (А) |. (15)
АеЗ п
Найдем число элементов в каждом классе З (А), А еЗ „. Для этого заметим, что, для фиксированного А = {А^; j = 1 ^ 5} , каждая из вершин х - и каждый древесный граф З ■ над множеством вершин А., j = 1 ^ 5 могут быть выбраны независимо от всех остальных вершин х е Ак и всех графов З к над этими множествами, к е / , {^ . Следовательно, в силупринципа
умножения (см., например, [7]), |\ (A) |= ^ s(A), где s(A) - число способов выбора отмеченной
AeA
вершины x в A и, одновременно, выбора древесного графа над множеством вершин A . Так как каждый из этих объектов может быть выбран независимо от другого, то, в силу принципа умножения, s(A) равно произведению числа вершин в A на число древесных графов, которые можно
построить на множестве A , s(A) =| A | L^ . Тогда из (15) следует, что \ (A) |= ^ |A . +
AeA
Полученная рекуррентная формула (14), в принципе, позволяет последовательно вычислять все числа Ln, и g □ на основе начального значения Л, = 1 , выполняя многократные суммирования в этой формуле. Но мы, здесь, с методической целью, вычислим производящую функцию
ад n
H (z) = Ê z~Ln+l, (16)
n=o n!
в которой положено L0 = 0, для 'этого бесконечного набора чисел, точно также как и для набора
чисел С-п, n G □ .
Теорема 7.Производящая функция H(z) является аналитической функцией в достаточно малой окрестности нуля на вещественной оси z и определяется какнаименьшее решение уравнения
H ( z) = exp( zH ( z)). (17)
□ Определим последовательность функций Э = {о(Хп )\п <Е + ), где о(Хп ) = nLn, n G □ и
a(0) = 0, ае A, то есть все функции этой последовательности являются постоянными. Согласно формуле (14),
Ln+i = \ (exp* a)n(Хя)dXn. [0,1]n
Тогда определение (16) производящей функции H (z) может быть представлено в виде
H ( z) = L( z; exp* а). Воспользовавшись формулой (10), запишем это равенство в виде
H ( z) = exp L( z; а). Наконец, преобразуем форму, стоящую в экспоненте,
ад n ад n
L(z;а) = Xz { an(Xn)dXn =^-zLn = zH(z). +
n= ![0J1]n n=o(n -1)!
Пример: Правильность полученного уравнения для производящей функции легко проверяется вычислением на его основе значений Ln при n = 1,2,3,4.
Положим H(z) = exp R(z) , R(z) = zH(z) . Тогда при n G □ имеем
R(n) = zH(n) + nH(n-1). (18)
Определим функции 5И (2) равенством Н(п)(^) = еЯ(2)£и (2) . Из этого определения следует рекуррентное соотношение
пеп, (19)
причем 5 (2) = 2Н'(2) + Н(2) . Равенство (18) при п = 2,3,... запишем в виде
Я(п)(2) = 2Бп (2) + п^и_1(2). (20)
Вычислим теперь Ъп при п = 1,2,3,4, исходя из неявно заданной уравнением (17) производящей функции, Ц = Н(п-1) (0) . Из (18) следует, что Я"(0) = 1, а из (20) - Я(п) (0) = ^0). Тогда ц = н (0)=1, 4 = н ' (0)=я ' (0)=1.
Для вычисления значений Ц , Ц нужно найти 52 (0), 53 (0) на основе рекурренции (19).
Для этого нужно вычислить значения производных (0), 5" (0). Они получаются из рекуррентного соотношения
(2) = ад + ^(2т2) + ^?(2№), ИеП (21)
и отдельно
5 ( 2) = гН 2 ( 2) + 2Н 2 ( 2). Отсюда следует при п = 1, что 5" (0) = 52(0) + 52 (0)Я '(0) + 5 (0)Я " (0) . Так как 5(0) = 2Н "(0) = 2 и из (20) следует Я"(0) = 2^ (0) = 2Н(0) = 2, то 52 (0) = 5 (0) + 5"(0) = 3 и, следовательно, Ц = 3 .
Далее, из вычисленных значений также следует, что 5" (0) = 5"(0) + 4, и так как
5"(2) = 2Нт(2) + 3Н " (2) , то 5"'(0) = 3Н" (0) = 352 (0) = 9. Поэтому 5" (0) = 13 и, следовательно, из
(20) следует Ц = 53 (0) = 52 (0) + 5" (0) = 16. Полученные значения совпадают со значениями, вычисленными в начале раздела.
Уравнение (14) определяет аналитическую функцию в круге с радиусом е- .В самом деле, особенность неявно заданной уравнением (14) функции Н (2) находится на вещественной оси. Это
следует из (14), так как Н " ( 2) = Н 2( 2)(1 - 2Н (2)) и особенности функции лежат там, где 2Н(2) = 1. Тогда, ввиду вещественности 2Н(2) в этих точках, из (14) следует, что Н (2) и, следовательно, 2 также вещественны. Далее, на вещественной оси, ввиду выпуклости функции вещественной переменной вхН относительно Н, ее график имеет два пересечения с графиком линейной функцией от Н, если х > 0 и достаточно мало. График имеет одно пересечение при х < 0 . Следовательно, уравнение (14) относительно Н имеет два решения при вещественных положительных и достаточно близких к нулю х . В этом случае мы выбираем, в качестве значения неявно заданной функции Н (х), наименьший корень уравнения. Именно он, по непрерывности, перехо-
дит в единственный корень уравнения (14) при X < 0 . Исчезновение же корней уравнения (14) при X > 0 происходит в точке бифуркации X, в которой исчезает пересечение линейной по Н функции с функцией вхН . Для этого необходимо и достаточно, чтобы в этой точке, кроме равенства Н = ех*Н, имело место равенство производных по Н, х»ех,Н = 1. Из полученной системы двух
уравнений находим хН = 1, то есть бифуркация реализуется в точке, где Н = в при х„ = в"1 и поэтому радиус сходимости степенного по X ряда для функции Н(х) равен е-1.
5. Уравнение для производящей функции ¥(г; е) .
Рассмотрим теперь комбинаторную задачу об определении числа Nn всех связных графов с помеченными п вершинами, в смысле определения этого понятия, данного в разд.~2 (см.~[10]). Легко понять, что при п = 1,2,3 функция N принимает значения N = 1, N = 1, N = 4. Непосредственное вычисление ее значений, посредством перечисления возможных графов,при п > 3, например, N = 38, N = 728 уже становится довольно рутинным.В этом разделе мы вычислим функцию N . Решение такой задачи будем понимать как получение рекуррентной формулы, связывающей число со всемичислами N. , j = 1 ^ п при п > 1. Эта формула позволяет последовательно определять значений N . Ее вывод основан на вычислениипроизводящей функции.
Пример:В качестве примера найдем несколько первых значений числа Nn непосредственным перебором имеющихся возможностей. Такой перебор, однако, начиная с п = 4 перебор всех возможных связных деревьев становится довольно рутинным.
1. При п = 1, V = {1} имеется только один граф ({1}, 0) , состоящий из единственной вершины 1. Он же, по определению, является связным.
2. При п = 2, V = {1,2} имеется два графа и один из них связный (V = {1,2}, ¥ = { {1,2} }) .
Следовательно, N = 1.
3. При п = 3, V = {1,2,3} имеется четыре связных графа, соответственно, со следующими множествами смежности: ¥ = {{1,2},{2,3}} , ¥ = {{1,2},{1,3}}, ¥ = {{1,3},{2,3}} , ¥ = {{1,3},{2,3},{1,2}} . Все они различны в смысле данного выше определения. Следовательно, N3 = 4.
4. При п = 4, V = {1,2,3,4}.
Среди связных графов с четырьмя вершинами имеется 16 различных древесных графов (см. [2]) и связные графы, содержащие, по крайней мере, один цикл. Эти недревесные графы разобьем на группы согласно различному топологическому типу, не учитывая наличия нумерации вершин. Таких групп четыре. Каждую из них опишем множеством смежности, из которого все остальные множества смежности рассматриваемой группы графов получаются подходящей пере-
становкой вершин. Затем для каждой группы укажем число графов, имеющихся в ней, которое легко устанавливается на основе непосредственного перебора.
Первую группу представляет множество смежности ¥ = {{1,2},{1,3},{2,3},{3,4}}. В этой группе имеется 12 графов. Вторую группу представляет множество смежности ¥ = {{1,2},{2,3},{3,4},{1,4}} и в этой группе имеется 3 графа. Далее, третью группу представляет множество смежности ¥ = {{1,2},{2,3},{3,4},{1,4},{2,3}}, в которой имеется 6 графов. Наконец, последняя группа содержит один граф с множеством смежности ¥ = {{1,2},{2,3},{3,4},{1,4},{1,3},{2,4}}. На основе такого перечисления заключаем, что имеется 22 связных недревесных графа с помеченными четырьмя вершинами. Следовательно, вместе с древесными графами, получаем, что N = 38.
Перейдем к общему решению задачи на основе разработанной в разд. 4 техники. Начнем с того, что введем в рассмотрение последовательность целых функций Кп (£) = (1 + £)и(и 1)2,
п е □ от комплексного переменного £ е .Для них в области комплексного £, где выполняется 11 + £ |< 1, абсолютно сходится степенной ряд
ад „и 2
Е (2,£) = £ ~Кп (£). (22)
и=о И!
Заметим, теперь, что любой граф (не обязательно связный) с помеченными вершинами из V = /и , однозначным образом, описывается посредством матрицы смежности , к = 1 ^ и, равной 1, если в этом графе имеется ребро, соединяющее вершины с номерами у и к и равной о, в противном случае. Введем переменные £( . е и определим в алгебре А элементы Ь(о, ;) = фп (Хл; еи); п е □ + >, Ь0 (еи) = 1 с компонентами
Ъ,ЛХп-£1])= П (1+ *,.,.), /!£□ (23)
(2)
такие, что их компоненты, функции Ъп (Хп; £ .) не зависят от точки Хи е [0,1]и .Тогда при £ ■ = £ значения этих функций совпадает с Кп £) = Ъп (Хп; £). Следовательно,
Е( г,е) = 1( г;Ь(£)). (24)
В силу своего определения, значения каждой из введенных функций совпадает при £ = 1 с числом Кп = 2"(" 1) 2 всех графов с помеченными п вершинами, п е .Последнее связано с тем, что Кп равно мощности множества всех слов длины п(п — 1) / 2 в алфавите из двух букв.
Производя перемножение всех скобок в определении (23), получим (см., например, [4, 5]),
что
К (Хя; ) = Е П «и даА); £,), п е □ , (25)
м
Ае п МеА
ап (Хп е) = У П , (2б)
¥п {1,1 }е¥„
где суммирование производится по всем связным графам (/и, ¥и ) над /п с отношением связности . Отсюда следует, что
Ъ(£у) = (ехр, а)(ег 1), (27)
если считать, что аЦ;) = <а;г (Хп; ); и е □ + > е А0.
Полагая в формуле (24) = е, мы определим функцию Ып (е) = ап (; е), а именно,
N (е) = (¥п), (28)
¥п
где I(¥п ) - число ребер в графе (/п, ¥п ). В частности, при е = 1 имеем Ып (1) =| {¥п: (/п, ¥п )} | -
число всех связных графов с помеченными вершинами над множеством вершин / .
Заметим, что функции Кп (е) и N (е) являются целыми функциями параметра е .Введем, теперь, производящую функцию
¥ (г;е) = £ % (е), Ып (е) = (^И) г=о- (29)
п=0 п! ¿г
Полагая в (10) а = а(е^ ), находим формулу, связывающую производящие функции последовательностей (КП(е);пе□ ) и (Ы„(е);пе□ ),
Е( г,е) = [/.( гМе, ))]е=е= [¿(г; ехр* а( е ^ = ехр([ц г; а( ^^ Л^) =
ш гп
= ехр У —Ып (е) = ехр ¥ (г, е). (30)
п=1 п! п
Так как при 11 + £ |< 1 ряд (22) представляет целую функцию от г е □ , то, в тех же условиях, степенной ряд, представляющий функцию
¥ (г;е) = 1п Е( г;е),
сходится в достаточно малой окрестности нуля, так как любая целая функция имеет на каждом компакте не более чем конечное множество нулей. Следовательно, при указанном ограничении на е, справедлива формула
1п ш п
К (е) = — 1п У —Кт (е)
V т=0 п! У г=0
(31)
При вычислении каждой п -й производной, согласно этой формуле, функция N (е) представляется полиномом от переменных Кт (е), п < п .В самом деле, положив Е(г; е) = 1 + Е(г; е) так, что
^ - ^
— Е(г; а) = — (а), л еП , (32)
^ 1=0
из (31) следует
N (а) = ^ Е Е (г;а)|,
ёг 1=1 1
Почленное дифференцирование согласно правилу
л +...+л=п
приводит при и (г) =... = и 1 (г) = Е (7; а), принимая во внимание (31), находим
" (-1)11 „ 1 К (а)
ЗД=«Е^ Е П—Г (33)
1=1 1 <л,.]=1 V
л +...+Л = "
Кх(е) = Кх(е) при .V е I и К0(е) = 0. Здесь мы использовали, что разложение Л ( г; а) начинается с первой степени г и, следовательно, ряд по т в правой части обрывается при т = п . Сформулированное утверждение следует из полученной формулы.
Так как каждый из коэффициентов К$ (а), .V е является полиномом относительно а , то N (а), определяемая (31), также является полиномом от а .Следовательно, формула (31), полученная нами при ееО , удовлетворяющих неравенству 11 + а |< 1, может быть однозначно аналитически продолжена на всю комплексную плоскость.В частности, эта формула справедлива при а = 1.
Так как — (1) = , то, вместе с — (а) = 0, при а = 1 из (33) приходим к заключению о справедливостиследующего утверждения.
Теорема 8. Число Ып всех связанных помеченных графов над множеством вершин 1п определяется формулой
" ( IV-1 1
N =«! ^ Е П . (34)
1=1 ' .....]=1
л +...+л =«
Расчет чисел для любого йёО на основе формулы (33) позволяет избежать рутинного перебора всех возможных связных графов с фиксированным числом помеченных вершин. Применение этой формулы легко приводит к значениям N при п = 1,2,3,4,5 , приведенных в начале этого раздела.
Заметим, что если бы мы использовали для доказательства теоремы производящую функцию, определяемую непосредственно на основе чисел , п е □ , то так как связь между этой
функцией ипроизводящей функцией для чисел п е П , которные получаются из (30) при
£ = 1, теряет смысл, ряд (22) при этом значении £ расходитсяво всех точках z Ф 0 комплексной плоскости. Поэтому возникает необходимость в обосновании получаемых решений, так как в процессе его нахождения приходится производить операции с расходящимися рядами,которые эквивалентны перестановкам бесконечного множества членов этих рядов.
6. Рекуррентная формула связи чисел Nn и М .
Приступим, теперь, к решению основной задачи вычисления числа М всех связанных графов с помеченными n вершинами без вершин сочленения. По определению, положим, что Мх = 1, М2 = 1, то есть уединенная вершина и единственный связный двухвершинный граф рассматриваются как графы без вершин сочленения. Очевидно, что М3 = 1 и непосредственным перебором имеющихся возможностей несложно устанавливается, что М4 = 10 . Однако уже вычисление следующего значения М5 = 238 уже требует определенных усилий. В этом разделе мы получим формулу, связывающую последовательности чисел (Nn;/?£□) и {Мп;ие0 } в виде рекуррентной формулы. Эта формула связывает число Nn+1, n > 1 c числами N. и М., j = 1 ^ n.
Она аналогична, по своему смыслу, формуле, которая связывает последовательности (Кп;ие0 ) и (Nп;иеD ), однако, по сравнению с ней, нуждается в более сложном обосновании.
Используя ее, можно вычислять числа Мп, n > 1 последовательно друг за другом, начиная со значения М = 1.
Теорема 9. Числа Nn, n = 2,3,... и Мп = 1,2,... связаны следующим рекуррентным соотношением
N/J+1 =1 П(I Мв+1 I ПNn(35)
AeS n AeA B.BcA (V(z);zeB)eS B+b zeB
начиная с n = 1.
U Доказательство состоит из следующих пунктов.
1. Обозначим посредством F(A ^{х}) - класс всех связных графов над множеством вершин A ^ {х}, в котором имеется одна отмеченная вершина х. Положим A = /и , | /и | = n и вершина х является вершиной сочленения степени s e{1,...,п}в связном графе G =(/и ^{х},¥) из
этого класса. При этом если s = 1, то вершина х, фактически, не является вершиной сочленения. Однако, чтобы не рассматривать этот случай отдельно, мы включаем ее в рассмотрение как вершину сочленения со степенью 1.
Рассмотрим графы, у которых s = 1 положим A = {Al = /п}, A eS „ и обозначим ¥ = ¥ и G = G . Если s Ф 1, то, согласно утверждению Теоремы 2, для выделенной вершины х имеется
однозначным образом определяемый набор связных графов G . = < A. u {x}, ), j = 1 ^ s, в каждом из которых вершина x уже не является вершиной сочленения, и при этом имеютместо следующие соотношения:
/„ = A u A u ••• u A, A ° A =0; ^^u^u... , ¥го¥у=0, i Ф j; i, j = 1 s.
Тогда набор {A,..., A} образует разложение A ES „множества ln так, что каждый из графов Gj, j = 1 ^ s построен над множеством вершин A u {x} и в каждом из них вершина x уже не является вершиной сочленения. При этом граф G представим в виде склейки графов Gj, j = 1 ^ s :
s
G = VG;. (36)
j=1
Введем для каждого разложения A = {A,...,A}eS „, состоящего из непустых множеств A , j = 1,...,s класс F[A ] = F(A,.., A ) связных графов G , которые определяются на основе формулы (36) из связных графов G ., j = 1 ^ s, каждый из которых является графом с помеченными вершинами над соответствующим ему множеством A u {x} , j = 1,..., s и в каждом из них отмеченная вершина x не является вершиной сочленения. Заметим, что для двух несовпадающих разложений A1 и A2 из S n классы F[AJ и F[A2] не пересекаются.
Тогда из приведенных выше рассуждений следует, что класс F(ln u{x}) связных графов с помеченными вершинами представим в виде дизъюнктивного разложения
и
F(i„^{x}) = (J и Л4»-Л)= U ЛА]. (37)
•5=1 Асу,1/1 Асу,,
Ввиду дизъюнктивности разложения (37), имеет место
n
N„+i =| F(ln u{x})| = Z Z I F(A,...,As)l = Z I F[A]| . (38)
s=1 AeS ns) AeS n
2.Обозначим посредством F (A u{x}) = F[A], A = {A} множество всех связанных графов с помеченными вершинами над множеством A u{x} с одной выделенной вершиной {x} , которая не является в каждом из этих графов вершиной сочленения и N4+1 =l F (A u{x}) |.
Заметим, что, ввиду справедливости формулы (36), для каждого набора {Д,...,As}, образующего разложение A ES п множества ln с непустыми составляющими, имеет место
s
F(л,..., As) = 0F' (AJ u {x}), то есть
j=1
F [A ] = (F' (A u {x}). (39)
AeA
Тогда, применяя комбинаторный принцип умножения, получаем
и [А ]|=П N А+1,
АеА
N„+1 = Е П NА|+1. (40)
Ае п АеА
3. Рассмотрим теперь графы каждого из классов ^ (А о{х}), А Ф0 , А е А . Пусть в = < А о {х}, ¥) - произвольный граф из этого класса. Согласно Теореме 3, каждому такому графу взаимно-однозначным образом, сопоставляется подграф вВ = (В о{х}, ¥В) графа , вВ Ф 0
которыйявляется в нем блоком. Тогда класс ^ (А о{х}) представляется в виде дизъюнктивного объединения
где ^ (А, В о{х}) - класс связных графов с помеченными вершинами, являющийся подклассом
класса ^ (А о{х}) таких графов, у которых вершина х содержится в блоке над множеством В о{х}. Тогда
и' (А о{х})|= Е '(А, В о{х})|. (41)
В:0фВс А
4. Вершины множества В не являются вершинами сочленения в графе в в, но они могут быть вершинами сочленения в графе в . Сопоставим каждой вершине г е В блока в в число р(г) , равное степени сочленения этой вершины в графе в . Если г не является вершиной сочленения, то, как и ранее, полагаем р(г) = 1. Тогда каждой вершине г е В, согласно Теореме 2, однозначным образом, сопоставляется набор связных графов в . (г) = (г) о {г}, (г)), : = 1 — р(г) , р(г) - степень вершины сочленения вершины г , р(г) <| А |. Среди всех таких графов содержится блок вв = (В о {х}, ). Не ограничивая общности, можно считать, что блок вв имеет номер р (г) в указанном наборе связных графов.
Для каждой вершины г е В и для каждой пары связных графов в . (г) и вк (г) из указанного набора имеет место V. (г) п V (г) = 0, ¥ ■ (г) п ¥ (г) = 0 ] Ф к, где :, к = 1 - р(г) — 1. Кроме того, для любых двух различных вершин ^ и из В все графы в ■ (^) , : = 1 — р(г) — 1 и в; (г ), : = 1 — р(г2 ) — 1 являются несвязанными друг с другом и при этом УЛ (г) о V: = 0, (г) пТ: (ъ) = 0, : = 1 -р(-1,: = 1 -р(-1.
Тогда, согласно формуле (3), все графы в класса Т (А, В о{х}) представимы в виде склейки
Р ( z )"1
G" = GB v ( V V GJ (z)).
zeB j=1
Определим теперь для каждого z e B связные графы
Р ( z )-1
G(z) = V G (z) (42)
j=1
такие, что G( z) = (V( z) ^{z}, ¥( z)>, где
1 P(Z)~ 1
V(z) = U Vj(z), 4>(z) = (J 4>,(z),
j=1 j=1
где, согласно Теореме 2,
p( z )-1
q(z) =| V(z) |= ^ Vj (z) |, ^ q(z) =| A , B |, (43)
j=1 zeB
где каждый из этих графов имеет выделенную вершину z . Заметим, что множества V(z) могут быть пустыми. Такая ситуация реализуется при p (z) = 1.
Следовательно, каждому графу G e F (A, B о{х}), однозначным образом, сопоставлены: упорядоченный набор (V(z); z e B> (согласно порядка вершин z в множестве B) множества
А , В с числом компонент |5| такой, что ^J V(z ) = А , В,и соответствующий ему набор графов
zgB
G(z) = (V(z) ^{z}, ¥(z)>, z e B такие, что имеет место формула склейки
G' = Gв v ( V G(z)). (44)
zeB
(Заметим, что здесь используются именно упорядоченные наборы (V(z); z e В>, а не разложения множества A , В . В частности, некоторые компоненты V(z) в этих наборах могут быть пустыми и, тем самым, они могут быть повторяющимися.)
Обратно, каждому упорядоченному набору (V(z); z e В>, являющемуся разбиением единицы для множества A , В , и набору графов G(z) = (V(z) ^{z}, ¥(z)> с отмеченными вершинами z e В, однозначным образом, по формуле (44) сопоставляется граф G из класса F (A, В о{х}) . Ввиду описанного взаимно однозначного соответствия, вводя обозначение^
для класса таких разложений множества A с числом компонент s , у которых допустимы пустые компоненты, имеет место дизъюнктивное объединение
г (А, В^{х})= и Г\А, в и (4; V(z), z Е в\
где F"(A, В {х}; V(z), z е В) - класс графов, содержащийся в F"(A, В ^ {х}), у которых зафиксирован упорядоченный набор (V(z); z e В>, определяющий разбиение единицы множества
А , В, на каждой компоненте которого строится связный граф (V(г) о {г}; ¥(г)) с выделенной вершиной г , которая приклеивается к блоку В при построении графа из класса Ц"(А, В о{х}). Из дизъюнктивности объединения следует формула
|Ц' (А,В о{х})|= Е Ц" (А,В о{х}^(г), г е В) |. (45)
У;егС(г)=Л, В
В свою очередь, из формулы (44) следует, что класс Ц"(А, В о1 {х}; V(г), г е В) представляется в виде прямого произведения
У \А, В о {х}; V(г), г е В) = 9(В о {х}) ® (V(г) о {г})), (46)
геВ
где §(В о {х}) - класс всех связных графов с помеченными вершинами над множеством вершин В о {х}, не имеющих вершин сочленения.
Из формулы (45) следует
| Г (А В о{х}; V (г), г е В) |=| §( В о{х}) | ПЦ (г) о{г}) |. (47)
геВ
Учитывая формулу (45), находим
и\А,В о{х}) |=| §( В о{х})| Е ПЦ(г) о{г})|. (48)
5.Из формул (41), (48) следует
\у'(au{*})\= Е G(Bu{*})\ Е П^(*)
В:0# Be A
<F(z):zeB): zgB
или, вводя обозначение \ G(B ^ {х}) \=М]в]+
NA+ = Е м\в\+1 Е ПNv(Z)\+i-
В:0*В-~А (V(Z):zgB): ZGB
(49)
Подстановка этой формулы в (4о)приводит к формуле (35) +
Утверждение теоремы переформулируется следующим образом.
Следствие. Числа N и М при п е подчинены следующему уравнению
(
n __S
N.+1=шее-1 Е П
k +...+ks=n
Л
kj __m 1
Е>+1 Е
+1
<7i,.JmM h+-+L=ki
(50)
□ Положим, что элемент а = (0,а,а,•••) из максимального идеала А алгебры А определяется компонентами ап(Хп) = А'' , /? е I , где А1'' |+1 определяется формулой (49). Из (35) следует, что
= [ехр, а](Хп) . Тогда, согласно (7), имеем
Г Т Т—т П 1 Т—т
[ехр* а]„ (X(1п)) = ^ П (X(А)) = Е Е П ^ (X(А)) =
Ае5 п ЛеА «=1 * ! Ае (*) АеА
п 1 п! *
= Е ТУ Е лТ-! П а,, (51)
к+...+к=п
так как компоненты аи (Хп ) зависят только от п .
Преобразуем теперь выражение (49) для аи = Щ |+1. Определим для каждого т = 1 ^ п элемент С("!) = (ск (Хк ); к е □ + ) е А, компоненты которого равны
"кУ-к
„( т)
(Хк) = Е П N,)!+1, т = |В | • (52)
ск (Хк
При этом с(т) Ф 0. Тогда, на основании (35), для каждого множества А с числом элементов к ■ = | А |, которое входит в состав разложения А е5 „, получаем
к к!
ЩА|+1 = Е мв,1 сАВ|)в|(х(А, в)) = Е ... ! Мт+1 сктт,
В:0фВс А т=1!(к - т) !
так как компоненты скп) (Xк ) , к е □ зависят только от к . Таким образом, из (49) и (51) следует
п 1 * (к, 1 с(т) ^
N. ■=»1 ЕЕ т. Е ПЕ -^м к
7 V! (к:' ? ^ да! ",+1 (к, - т)!
к1 +...+к* =п
(53)
где в произведении, для каждого сомножителя , = 1 ^ т, произведена замена к ^ к. — т .Из определяющей формулы (52) следует, что
с( т)
к - ^ 11 ^ '=|У(^)|, г е В.
г(т) т 1
= Е П>. .1.
Как и выше, заменив к ^ к — т в каждом сомножителе с , = 1 ^ т и подставив полученные выражения для этих сомножителей в (53), получим формулу (50). +
На основе формулы (50) имеется возможность вычислять числа Мп, п е □ , рассматривая
ее как систему уравнений для этого числового набора. Однако такой метод определения этих чисел оказывается очень трудоемким. В связи с чем, в следующем разделе, вводится специальная производящая функция, коэффициенты которой Мп (е) непосредственно связаны с м. , и устанавли-
,,..., ¿т
вается вид функционального уравнения, которому она удовлетворяет и посредством которого она однозначно определяется.
7. Уравнение для производящей функции О (
Так как числа Мп определяются рекуррентным образом на основе чисел , и так как для последних невозможно непосредственное введение производящей функции, то таким неприятным свойством обладает и набор Мп, п е □ . По этой причине, мы введем для каждого £ е □ , I + £ |< I производящую функцию
О(г;е) = X^„¿е), Мп(е) = (^О^)2 0 (54)
п=0 П!
где коэффициенты Мп{£), п е □ определяются однозначно из следующей системы уравнений
( К/ 1 т
Т X Л
£
П 1 __£ (1 __М 1
КЛе) = 1 X П XЛми+1(е) X Л^
(55)
V' -
к^П J=l-^s■. 1 I I : т\
^+...+^=11 11+...+1т=кгт
При этом, ввиду Nп(1) = Nп, выполняется также равенство Мп{\ ) =М„ при всех /7 £ I . Так как N (е) зависят от е полиномиально, то в силу рекуррентной формулы (55), коэффициенты Мп (е) зависят от е также полиномиально.
Следующее утверждение дает нам функциональное уравнение, связывающее функцию О(z, е) с производящей функцией ¥(г, е) , и которое определяет ее однозначно при условии 0(0,е) = 1.
Теорема 10. Производящие функции ¥(г;е) и О(связаны следующим функциональным уравнениемпри изменении параметра е в области комплексной плоскости, где выполняется неравенство 11 + е |< 1:
¥ '(г; е) = ехр [О( ' (г; е); е) -1]. (56)
и Запишем ряд для производной производящей функции ¥(г; е) в виде
<» zn
¥' (г;е) = 1+ Х ^п+1(е) П=1 п!
В силу выполнимости рекуррентной формулы (55) при п е □ заменим Nп , согласно этой формуле
га п __£ ( ! ] __т ^
¥(г;е) = 1 + ]Т X 1 X ПХ X П ^ +1)
п=1 £=1 £ ! < К ,..А >, /=1т=1т! а ,..,1т>, 1=14 !
. . : Ч- J
/¿е +,/=1+тп:
кг+...+к.=п 11+...+1т=к1-ш
Переставляя соответствующим образом суммирования, произведем преобразования рядов. Перестановка суммирований обоснована тем, что в области изменения параметра е функция ¥(г; е) является целой, и поэтому все ряды, с которыми производятся преобразования, также являются целыми функциями. Справедливо следующее:
ад 1 ад __* | к, ] __т ^
¥(*;*) = 1 + Е Е П^ ЕЕ"Ч,! Е Пт+1(е)
¿1+...+к=п
, 1
(А-!,.„Л), ] = к +...+к =п
= 1+ Е-1 Е ПЕ Мт+1 (е) Е П^Ч)
*=1 (к ,...,к), ,=1 т=1 т! (I ,...,С), ¿=1 Ч !
<4,..,'„>,
, Л=\+п
ад 1 * ад
^ к,
=1+Е^ПЕ Е^+1) Е Пт^ч
к^ут^т- <4,-Л), 7=1 /
(е)
/(: ./ I :
ад
ад 1 т—г ад тт ж—I т ч
= 1+ Е ^ П|Е ^гМт+1(е) Е Е П^Л „(е)
- I Л':/• А . /// </!,..,/„>, 7=1 ( !
•' 7 1 +.'=1+»':
ад 1 5 /ад _т
г
■г'
= 1+ Е "7 П1Е ^ми+1(е) Е П уг^+1(е)
- I Л' : I I '=1 7
т ад I
л
= 1+ Е1 Е^Мт+^ПЕ+1(е) = 1+ Е1 Е^Мт+^ХПге)]
*=1 * ! ,=1 ^т=1 т! '=1 /¿=0 I' ! ) *=1 ,=1 ут=1 т!
ад 1 *
=1+ Е П 0(' (г; е)) -1] = ехр [0( г¥ (г; е))) -1]. + *=1 *!,=1
Полученное нами функциональное соотношение (56) можно рассматривать как уравнение, определяющее однозначно неявную функцию О (г;е) при дополнительном условии
0(0; е) = М =1
. Она определена в области комплексной плоскости изменения е, где имеет место
11 + е |< 1. Вычисленные на основе этой функции коэффициенты Мп (е) допускают аналитическое продолжение на всю комплексную плоскость е , ввиду их полиномиальной зависимости от этого параметра. В частности, полученные таким образом формулы дают возможность вычислить м (1) .
8. Заключение
Доказанная Теорема 8 позволяет дифференцированием последовательно находить числа Ми+1 по известным числам = 1 — (п +1) и по предварительно вычисленным числам М.,
¡1 +...+Чт =к,-
ад
¡1 +...+'т =к,
у
ад
/ = 1 ^ п. При этом нам не приходится решать алгебраическую систему уравнений, которой подчинены числа М., / = 1 ^ п +1 согласно формуле (50).
Произведем, на основе полученного уравнения (56), вычисления нескольких первых членов последовательности {МП ) с целью сравнения с приведенными ранее результатами их непо-
средственного вычисления посредством перебора. Для этого представим (56) в виде 1 + 1п ¥'(г; е) = О( ' (г; е)). Последовательно дифференцируя обе части этого равенства, находим
О '('(г; е); е)(("(г; е)) + ¥'(г; е)) = ^^,
¥ (г;е)
О " ( " ( г; е); е)(( " ( г; е)) + ¥ " ( г; е))2 + О' ( " ( г; е); е)(( " ( г; е)) + 2¥ " ( г; е)) =
_ ¥"(г;е) ¥''2(г;е) ¥'(г;е) ¥'2(г;е) '
О "" ( " ( г; е); е)(( " ( г; е)) + ¥ " ( г; е))3 +
+3О" (" (г; е); е)(( " (г; е)) + 2¥ " ( г; е))(( "( г; е)) + ¥ " ( г; е)) +
+О' ((г;е);е)(( (г;е)) + 3¥""(г;е)) = - 3 ¥ "" (^ + 2 .
Эти формулы позволяют подстановкой значений производных ¥(0; 1) = 1, ¥'(0; 1) = 1, ¥"(0; 1) = N = 4, ¥"'(0; 1) = N = 38, которые нами взяты из работы [10], найти значения для
О"(0;1) = М = 1, О"'(0;1) = М4 = 10.
Вычисление значения производной О1Г (г; е), на основе дифференцирования последнего из вышеприведенных соотношений, приводит к формуле
О!У (" (г; е); е)( " (г; е) + ¥ " (г; е))4 + +6О "" ( ' ( г; е); е)(г¥ " ( г; е) + ¥' ( г; е))2( '"( г; е) + 2¥ " ( г; е)) +
+ 3О " (' ( г; е); е)(г¥" ( г; е) + 2¥ " (г; е))2 + + 4О " (" ( г; е); е)(г¥" ( г; е) + ¥ " ( г; е))( ( г; е) + 3¥ "" ( г; е)) + + О' (г¥' (г;е);е)^¥у (г;е) + 4¥ш (г;е)) =
_¥¥ (г;е) ¥!У (г;е)¥"(г;е) ¥" 2( гге) "(гге) ¥ " 2( гге) ¥'"4( гге) ¥' (гге) ¥' 2( гге) ¥' 2( гге) ¥' 3( гге) ¥' 4( z;е),
которая позволяет вычислить число О!У (0; 1) = М5 нам основе известного уже известного числа
N = 728 связных графов с пятью помеченными вершинами. В результате, получаем М5 = 238 .
Несмотря на то, что предложенный алгоритм последовательного вычисления чисел М. предполагает рутинное, довольно громоздкое вычисление рекуррентных соотношений для произ-
водных ряда G(z;s), его достоинством является то, что он позволяет автоматизировать процесс вычисления чисел Мп, не прибегая к перечислению топологически неэквивалентных графов с
числом вершин n, а затем выявления всех тех перестановок меток вершин, которые не переводят графы каждого топологически фиксированного типа друг в друга.
Список литературы
1. Харари Ф. 2006. Теория графов. М.: КомКнига, 296. Harary F. 1996. Gbaph theory. London: Addison-Wesley Publishing Company.
2. Харари Ф., Палмер Э. 1997. Перечисление графов. М.: Мир.
Harary F., Palmer E.M. 1972. Graphical Enumeration. New-York: Academic Press.
3. Оре О. Теория графов 1980. М.: Наука,33б.
Ore Oystein, 1962. Theory of graphs American Mathematical Society. Colloquium Publications V. XXXVIII.
4. Майер Дж., Гепперт-Майер М. 1980. Статистическая механика. М.: Мир,54б.
Mayer J.E., Goeppert-Mayer M. 1977. Statistical mechanics . New York: John Wiley & Sons, Inc.
5. Рюэль Д.1971. Статистическая механика. Строгие результаты. М.: Мир, 368.
Ruelle D. 1969. Statistical mechanics, Rigorous Results. New York, Amsterdam: W.A. Benjamin, Inc.
6. Риордан Дж. 1963. Введение в комбинаторный анализ. М.: Мир.
Riordan J. 1958. An introduction to combinatorial analysis. New York: John Wiley & Sons, Inc.; London: Chapman & Hall, Limited.
7. Холл М. 1963. Комбинаторный анализ. М.: Изд. Иностр. лит., 100c.
Holl M., jr. 1958. A survey of combinatorial analysis. New York: John Wiley & Sons, Inc.
8. Вирченко Ю.П., Остапенко Л.П. 2015. Определение числа разложений конечного множества. Белгородский государственный университет Научные ведомости. Математика и Физика. 5(202); 38 : 96-100. Virchenko Yu.P., Ostapenko L.P. 2015. Number evaluation of finite set partitions. Belgorod State University Scientific Bulletin. Mathematics & Physics. o5(202); 38: 96-100.
9. Вирченко Ю.П., Остапенко Л.П. 2015. Определение числа древесных графов над конечным множеством вершин. Белгородский государственный университет Научные ведомости. Математика и Физика. 11(208); 39: 37-43. Virchenko Yu.P., Ostapenko L.P.2015. Number evaluation of tree graphs with finite vertices set. Belgorod State University Scientific Bulletin. Mathematics & Physics. 11(208); 39: 37-43.
10. Вирченко Ю.П., Остапенко Л.П.2015. Определение числа связных графов над конечным множеством вершин. Белгородский государственный университет Научные ведомости. Математика и Физика. 17(214); 40: 28-34. Virchenko Yu.P., Ostapenko L.P.2015. Number evaluation of connected graphs with finite vertices set. Belgorod State University Scientific Bulletin. Mathematics & Physics. 17(214); 40: 28-34.
11. Ван дер Ваден Б.Л.,1979. Алгебра. М. : Наука, 648.
Van Der Waerden B.L. 1971. Algebra II. Berlin : Springer-Verlag.
12. Остапенко Л.П., Вирченко Ю.П. 2016. Число связных графов без вершин сочленения. Материалы международной конференции «Воронежская зимняя математическая школа С.Г.Крейна. Воронеж: Научная книга, 310-314. Virchenko Yu.P., Ostapenko L.P. 2016. Number of connected graphs without articulation vertices. Krein's Winter Mathematical School. Voronezh State University. 310-314.
