CC BY
41
УДК 519.171.1
ИССЛЕДОВАНИЕ ПОЛИНОМИАЛЬНОСТИ МЕТОДА ВЫЧИСЛЕНИЯ ИНТЕГРАЛЬНОГО ОПИСАТЕЛЯ СТРУКТУРЫ ГРАФА
В.К. Погребной, Ан.В. Погребной
Томский политехнический университет E-mail: [email protected]
Актуальность исследования определяется большой потребностью в разработке эффективных методов инвариантного описания и анализа абстрактных структур графовых моделей.
Цель работы заключается в обосновании полиномиальности предложенного авторами метода вычисления интегрального описателя абстрактной структуры графа.
Методы исследования основываются на использовании аппарата теории графов и методов свободной и зависимой интеграции кодов структурных различий графов.
В результате исследования введено понятие устойчивых групп вершин в графе и сформулированы условия возникновения и существования таких групп в процессе интеграции кодов структурных различий при вычислении интегрального описателя структуры - Integral structure descriptor (ISD). Для устойчивых групп установлен ряд свойств, которые раскрывают правомерность применения основных правил метода ISD и его полиномиальности. На основе выделенных свойств установлено, что условия существования устойчивых групп обусловлены жесткими ограничениями, а вершины разных устойчивых групп не могут порождать новые устойчивые группы. Установлен также фактор полной обособленности устойчивых групп, что в значительной мере предопределило эффективность алгоритма вычисления интегрального описателя структуры графа. Полиномиальность метода показана для наиболее трудного случая, когда графы являются однородными и содержат устойчивые группы. Для экспериментальных исследований метода ISD на языке Java разработано программное средство GraphlSD и приведены некоторые результаты его работы.
Ключевые слова:
Интегральный описатель структуры, изоморфизм графов, абстрактная структура графа, устойчивая группа вершин, область интеграции кодов, полиномиальность алгоритма.
Введение
Предложенный в [1] метод интеграции структурных различий позволяет получать интегральный описатель структуры графа. Метод даёт возможность выявить и объединить (интегрировать в описатель) имеющиеся в графе структурные различия так, чтобы каждая вершина получила свой интегральный описатель в виде уникального числового кода. Интегральные описатели вершин графа G образуют вектор D(G), который получил название интегрального описателя структуры графа - Integral Structure Descriptor (iSD). В работе [2] метод вычисления ISD был адаптирован для решения задач определения оценок сходства структур двух графов на основе выделения в них общих частей. При этом были исследованы некоторые свойства однородных групп вершин, возникающих в процессе интеграции структурных различий.
В работах [1, 2] отмечалось также, что однородные группы при определенных условиях переходят в категорию устойчивых, наличие которых приводит к заметному усложнению работы метода определения ISD. В связи с этим актуальными становятся исследования по двум направлениям -определить влияние устойчивых групп на объемы вычислений при получении ISD и проверить, находятся ли эти объёмы в рамках полиномиальных оценок сложности [3]. Получению ответов на эти вопросы посвящена данная работа. Основное внимание уделено исследованию свойств устойчивых групп. На основе выделенных свойств установлена полиномиальность предельных оценок сложности
для наиболее неблагоприятного сценария процесса интеграции структурных различий. Экспериментальные исследования метода интеграции структурных различий проводились с помощью программы GraphlSD, написанной на языке Java.
Устойчивые однородные группы и их свойства
Рассматривается конечное множество M(S) обыкновенных связных графов, имеющих равные степенные инварианты, представленные вектором S. Множество M(S) представим совокупностью классов изоморфных графов. В одном из классов выделим граф G в качестве представителя этого класса. Соответственно, класс изоморфных графов, представленный графом G с инвариантом S, именуется как G(S).
Будем считать, что для графа G=(E,U) с помощью метода свободной интеграции кодов структурных различий (метода ISD) получен вектор D(G)={dJ интегральных описателей d вершин eteE и сформирована область интеграции We(S) для класса G(S). Пример процесса свободной интеграции представим в виде рекуррентной последовательности преобразований векторов Dk^Dk+u.
D0 ^ D ^ D2 ^ D3 ^ D64*2 ^
^ D5 ^ D46* ^ D7 ^ D8 = D(G). (1)
Из записи (1) следует, что для получения вектора D(G) потребовалось выполнить 8 шагов преобразований Dk^Dk+1. На k-м шаге метод ISD устанавливает для каждой вершины eteE код k-го уровня интеграции d‘, характеризующий особое положение данной вершины относительно других вершин
графа. В методе ISD введено два правила назначения кодов dk. Правило 1 отражает результат интеграции структурных различий, а правило 2 используется для «разрушения» устойчивых групп путём назначения вершине абстрактного (виртуального) различия.
Если в графе G на k-м шаге интеграции некоторые вершины eieE занимают одинаковое положение относительно вершин E\et, то они получат равные коды dk. Совокупность вершин, у которых значения dk в векторе Dk совпадают, образуют однородную группу Edk c dk=d. Однородные группы Edk в векторе Dk называются устойчивыми, если при переходе Dk^Dk+1 максимальное значение кода dk+1 в векторе Dk+1 не превышает максимального значения кода dk в векторе Dk. В этом случае для продолжения процесса интеграции в одну из устойчивых групп вектора Dk+1 по правилу 2 вводится виртуальное структурное различие. В записи (1) такая ситуация отмечена при переходах D°^D11, D3^D42, D5^D63. В соответствующих векторах шаги ¿=1,4,6 помечены дополнительным индексом и=1,2,3, который обозначает порядковый номер шага, где применялось правило 2.
Устойчивые группы и их свойства играют важную роль в понимании процесса интеграции кодов по методу ISD и обосновании правил формирования области интеграции WG(S). Свойства устойчивых групп позволили установить полиномиальность оценок объёмов вычислений при поиске векторов Dg(H) для графов Н из множества M(S). Поиск осуществляется на основе зависимой интеграции кодов с помощью метода ISD относительно области WG(S). Индекс G в обозначении Dg(H) показывает, что вектор ISD получен в зависимости от системы кодирования WG(S). Формат статьи не позволяет рассмотреть и обосновать всё многообразие свойств устойчивых групп, поэтому ниже рассматриваются только те из них, которые непосредственно подтверждают правомерность применения правил формирования области WG(S) и проверки наличия у графа Н из множества M(S) вектора Dg(H), который по составу кодов совпадает с вектором D(G). В этом случае граф Н принадлежит классу G(S). Совокупность свойств представлена в форме ряда утверждений. Каждое утверждение сопровождается необходимыми дополнительными пояснениями, которые основываются на теории интеграции кодов структурных различий в абстрактных структурах графов.
Утверждение 1. В графе G устойчивая группа с множеством вершин Edku^E может существовать только в виде однородного подграфа. Вершины в подграфе имеют равные значения степеней Sd из множества значений (0,1,2,...). Вершины eteEdk° в графе G, независимо от того, является он однородным или нет, также должны иметь равные степени St и при этом St>Sd. Однородный подграф на основе вершин устойчивой группы Edkv может состоять из одной или нескольких компонент связности, каждая из которых также является однородным под-
графом степени Бі. Подграфы устойчивых групп Ек"(Є) и ЕЦ"(Н), полученные в графах Є и Н на шаге к после введения виртуального различия по правилу 2 при Бй=0, являются пустыми, при Б=1 подграфы содержат четное число вершин и являются изоморфными, при Б^=2 подграфы изоморфны для числа вершин 3, 4, 5, а при 6 и более вершин подграфы могут состоять из нескольких компонент связности. Приведенные свойства легко подтверждаются на основе анализа множеств кодов й'к в инциденторах ¥(г) для вершин гієЕік" и определения устойчивой группы.
Утверждение 2. Окружением устойчивой группы ЕЦ" являются другие устойчивые группы и вершины с уникальными кодами. Однородный граф можно рассматривать в качестве устойчивой группы Ек1 на шаге к1=1, у которой ¿=1 и отсутствует окружение. Если две устойчивые группы Ек и ЕЦ° связаны между собой, то каждая вершина в{<ЕЕк с вершинами группы Е% должна иметь одинаковое число связей (ребер). Аналогичное условие должно выполняться и для вершин группы Е£, т. е. структура связей между двумя группами дол2 -жна быть однородной. Если устойчивая группа Е к связана с вершиной г, имеющей уникальный код йк*", то каждая вершина г^Ек должна быть связана с вершиной г. Из этого следует, что число вершин в группе Ек в этом случае не может превышать степень вершины г.
Утверждение 3. Все вершины устойчивой группы Е/Л полученной на шаге £„+1, принадлежат только одной из групп Ек", существовавших ранее на шаге £„, т. е. Е/"+1 сЕЦ". Это утверждение отражает принципиально важное свойство устойчивых групп и процесса интеграции кодов структурных различий. Из него следует, что на шаге кс после введения в одну из групп по правилу 2 виртуального различия число вершин ни в одной из групп на шаге к„+1 не может возрасти. При этом отдельные группы «разрушаются», остаются неизменными либо преобразуются в более мелкие. Доказательство этого утверждения строится на основе анализа и сопоставления деревьев, отражающих интеграцию кодов вершин, входящих в рассматриваемую группу.
Утверждение 4. Вершина ^єЕ/", получившая на £-м шаге интеграции уникальный код ¿к" по правилу 2, сохраняет свою уникальность до конца процесса интеграции. Это даёт возможность на очередном (£+1)-м шаге интеграции для таких вершин переносить без изменения уникальный код ¿к из вектора Бк в вектор Бк+1. Вершины, получившие уникальный код ¿к по правилу 1, также сохраняют свою уникальность, но в процессе интеграции должны участвовать на общих основаниях. Из этого следует, что если вершина гі в графе Є и вершина г1 в графе Н получили уникальные коды й*=йук, то при изоморфизме графов Є и Н вершины г{ и гі в подстановке изоморфизма образуют пару (г1,гі). Здесь доказательство строится на основе сопоставления кодов двух вершин. Если на к-м шаге вер-
шина е{ получает код йк, который отличает её от кода йквершины е, т. е. йкФйк, то это означает, что вершины е1 и е, интегрировали разные структурные различия. Эти различия в ходе последующей интеграции не теряются, они только могут пополниться новыми, т. к. код вершины интегрирует (накапливает) все структурные различия, обнаруженные на любом шаге интеграции. Поэтому у вершин с йкФйк равенство кодов й{=й- при ¿>£ невозможно.
Утверждение 5. Устойчивая группа Е/"(Н), полученная на к„-м шаге интеграции кодов в графе Н относительно области ЖО(5), и соответствующая ей устойчивая группа Е/"(О) в графе О, в случае изоморфизма графов О и Н, в подстановке изоморфизма образуют паросочетание. Данное утверждение можно рассматривать как некоторую аналогию утверждения 4, если устойчивые группы в утверждении 5 воспринимать в качестве отдельных объектов с уникальными кодами по аналогии с вершинами в утверждении 4. Важность утверждения 5 состоит в том, что оно отражает обособленность вершин устойчивой группы при определении изоморфизма графов О и Н. Другими словами, если вершины ек и е2 принадлежат устойчивой группе ЕЦО), то вершины е,г и е2, которые в подстановке изоморфизма графов О и Н с вершинами е11 и ек образуют пары, должны принадлежать устойчивой группе Е/"(Н).
Утверждение 6. При поиске в графе Н вектора БО(Н) по составу кодов й, совпадающего с кодами й вектора ДО) в условиях наличия в области ЖО(Й) шага к„, на котором в устойчивую группу Е/"(О) согласно правилу 2 вводилось виртуальное различие, необходимым и достаточным является последовательное назначение данного виртуального различия каждой вершине в соответствующей группе ЕЦ"(Н) и выполнение процесса интеграции до шага к„+1 или получение вектора БО(Н). Необходимость выполнения процессов интеграции относительно каждой вершины группы Е%"(Н) обусловлена тем, что при формировании области ЖО(5) на шаге к„ виртуальное различие назначалось произвольно для одной из вершин е1 группы ЕЦ"(О). Поэтому требуется установить, какая из вершин е, из группы Ейк"(Н) может войти в пару (е,е) в подстановке изоморфизма графов О и Н. Исходя из свойства обособленности устойчивых групп, отмеченного в утверждении 5, выполнение процессов интеграции относительно вершин группы ЕЦ"(Н) является достаточным.
Первые три утверждения в основном отражают некоторые условия возникновения и существования устойчивых групп. Здесь можно отметить, что эти условия являются весьма жёсткими. Появление устойчивых групп в структуре графа возможно лишь при наличии симметричности в отношениях между вершинами. Условия симметричности должны соблюдаться не только между вершинами внутри устойчивых групп, но и с вершинами окружения. Любые нарушения симметричности по-
рождают структурные различия и нарушают условия существования устойчивых групп.
Последние три утверждения раскрывают правомерность применения основных правил, используемых в алгоритме поиска вектора Бй(Н) при наличии устойчивых групп. В первую очередь это относится к свойству обособленности вершин с уникальными кодами и устойчивых групп. Фактор обособленности устойчивых групп дает возможность резко сократить область поиска вектора Бй(Н).
Полиномиальность алгоритма поиска вектора DG(H)
Объем вычислений, который требуется алгоритму по методу ISD, зависит от числа шагов интеграции, выполняемых при определении вектора ISD. Для графа й, представляющего класс графов й(Б), общий объем вычислений Ь(й) при формировании области интеграции Жй(Б) и определении вектора Б(й) составит
Цй)=Акт(п,ш), (2)
где т(п,ш) - объем вычислений при выполнении одного шага интеграции для графа, содержащего п вершин и ш ребер; Дк- число шагов интеграции, выполняемых при переходе от вектора Б° к вектору Бк+1=Б(й).
Значение величины т(п,ш) зависит не только от числа вершин п, но и от числа ребер ш, что отражается на заполнении матрицы смежности вершин ненулевыми элементами и определяет размерность множеств кодов инциденторов ґ(г), которая влияет на объем вычислений. В ходе интеграции с ростом числа вершин с уникальными кодами число формируемых множеств ¥(г) уменьшается и соответственно снижается значение т(п,ш). Будем считать, что для определенных интервалов значений чисел вершин и ребер установлены усредненные значения т(п,ш). Вторая величина в выражении (2) не может превышать величину п, т. е. Дк<п. Это обусловлено тем, что каждый шаг интеграции приводит к росту числа отличающихся кодов, назначаемых по правилу 1 или 2.
Для графа Н, в котором процесс интеграции осуществляется относительно области Жй(Б), объем вычислений Ь(Н) при поиске вектора Бй(Н), совпадающего по составу кодов с вектором Б(й), зависит от наличия устойчивых групп. Если в таких векторах упорядочить коды по возрастанию, то полученные векторы БДН) и Б'(й) совпадут, т. е. Бй*(Н)=Б*(й). Заметим, что если граф й неоднородный и процесс интеграции в нём не приводит к появлению устойчивых групп, то величина Ь(Н) определяется по выражению (2).
Наличие устойчивых групп в графе й усложняет поиск в графе Н вектора Бй*(Н)=Б*(й). Вместе с тем свойство обособленности устойчивых групп и отдельных вершин с уникальными кодами, которые также можно рассматривать как устойчивые группы, содержащие по одной вершине, оказалось настолько полезным, что сделало возможным задачу поиска Бй*(Н)=Б*(й) исключить из класса не-
П = 12
П2 = 6
п3“ = 2
Дко-1
д^-з
Дк2-2
Дкз-2
Рис. 1. Схема поиска вектора йс, (Н) относительно примера в записи (1)
полиномиальной сложности. Напомним, что поиск Бв*(Н)=Б*(в) соответствует решению задачи определения изоморфизма, т. к. при соблюдении этого условия графы в и Н изоморфны.
Данные вопросы рассмотрим более подробно. При этом будем придерживаться примера процесса интеграции, приведённого в записи (1). Схему поиска вектора Бв(Н) относительно области Же(й), полученную для данного примера, представим на рис. 1.
Согласно записи (1), общее число шагов Ак=8 в процессе интеграции разбивается на 4 интервала, которые в явном виде приведены на рис. 1. Начальный интервал Ак0 содержит число шагов при переходе от вектора Б° к Бк1. Если граф в и, соответственно, граф Н являются однородными, то Ак0=1. В нашем примере графы в и Н однородные, т. к. Бк1=Б1\ т. е. уже на шаге к=1 потребовалось вводить первое виртуальное различие и шаг к=1 помечать индексом ,=1. Индекс V обозначает порядковый номер вектора Б1, содержащего устойчивые группы. Последующий интервал Ак= указывает на число шагов от Б11 до Бк. Аналогично следующий интервал Ак2=2, а последний интервал Ак3=2 равен числу шагов от Б83 до вектора Б8=Бв(Н).
При подсчете общего числа шагов, которые потребуется выполнить в процессе поиска вектора Бв(Н), будем исходить из того, что на каждом интервале от Бк до Бкм виртуальное различие будет поочередно вводиться для всех вершин устойчивой группы Ек,(Н), выбранной в соответствии с Жв(й). Число вершин в этой группе обозначим величиной п“. Если на предыдущем интервале от Бк,-1 до Бк после введения виртуального различия вершинам устойчивой группы Е^-1(Н) с числом вершин п,_1 процесс интеграции относительно Жв(5) успешно завершится в векторе Бк для нескольких вершин nv из числа п“ч, п„<п“„_1, то на интервале от Бк до Бк,м процесс интеграции должен быть выполнен для каждого варианта Бк из числа п„. Таким образом, число шагов на интервале от Бк до Бкм составит величину Ак,-п,-п„ш. Например, для интервала Б42^Бв3 (рис. 1) потребуется выполнить Ак2-п2-п2“=2-3-6=36 шагов.
Исходя из этого, предельный объём вычислений Ь(Н), который может потребоваться при поиске вектора Бв(Н), складывается из объёмов вычислений по каждому интервалу и оценивается следующим полиномом:
Ь(Н) = (Ак0 + ^пП,Аку)т(п, т).
(3)
Здесь Дк0 - число шагов интеграции от вектора Б0 до вектора Бк1, содержащего устойчивые группы; Дк„ - число шагов интеграции между векторами Бк и Бк"+1, и=1,2,..., V, V- число векторов Бк°, содержащих устойчивые группы; п„ - число вершин в устойчивой группе £/"-1(Я), для которых после введения виртуального различия согласно Жд(5) и выполнения Дк„ шагов интеграции достигнуто соответствие вектора Бк" и записи в строке к„ области Же(5); пЦ - число вершин в устойчивой группе £/"(Я), для которых на шаге к„ согласно Же(й) вводятся виртуальные различия.
Полином (3) оценивает предельные объемы вычислений, когда на каждом интервале Дк„ процесс интеграции выполняется для всех вариантов, т. е. всякий раз после введения виртуального различия для каждой вершины устойчивой группы, и только после этого происходит переход от интервала Дк„ к Дк„+1. Такое представление работы алгоритма поиска вектора Бд(Я) удобно для формальной записи значения предельной оценки Ь(Я) в виде полинома.
В действительности программа GrapЫSD реализует более эффективную схему поиска, когда процесс интеграции в каждом интервале Дк„ выполняется для одного из альтернативных вариантов до тех пор, пока соблюдается соответствие кодам Же(й). При несоответствии кодов делается обратный ход к ближайшему альтернативному варианту выполнения процесса интеграции. Если вектор Бд(Я) удается получить после введения виртуального различия для очередной вершины устойчивой группы, то нет смысла продолжать поиск для последующих вершин. Программа GrapЫSD фиксирует первый найденный интегральный описатель Бд(Я) и по желанию пользователя может продолжить поиск другого описателя, в том числе и с формированием новых систем кодирования, выбирая другие устойчивые группы.
Результат работы программы GrapЫSD при поиске вектора Бд(Я) показан на рис. 2. Графы О и Я с равными степенными инвариантами представлены матрицами смежности вершин (граф О - в верхней части рис. 2, граф Я - в нижней). Для получения вектора Бд потребовалось 6 шагов интеграции. Виртуальное различие вводилось один раз на
у=1
Рис. 2. Пример поиска вектора йГз(Н) с помощью программы вгэрЫБО
4-м шаге в устойчивую группу £241(О)=(е2,е5). Для поиска вектора Бд(Я) потребовалось 9 шагов интеграции. Вектор Бд(Я) найти не удалось, т. е. граф Я неизоморфен графу О.
Очевидно, что наибольший объем вычислений предполагается в случае однородного графа и соответствует полиному (3). Эффективность алгоритма поиска Бд(Я) может быть улучшена при более полном учете свойств, приведенных в утверждениях 1-6, а также других свойств и не только устойчивых, но и однородных групп [2]. В частности, для однородных графов с целью начальной дифференциации вершин вместо степенного инварианта можно использовать другие легко вычисляемые инварианты.
Вместе с тем алгоритм поиска вектора Бд(Я) и в изложенном варианте способен успешно решать задачи разбиения множества графов на подмножества изоморфных или выбора из него таких графов, которые изоморфны заданному. При разработке метода ISD эти задачи были отмечены как частный случай более общей и существенно более сложной проблемы определения оценок сходства двух графов [2]. Поэтому изложенные выше результаты исследований по устойчивым группам, формированию области интеграции и поиску вектора Бд(Я) расширяют возможности формального описания и анализа сходства абстрактных структур графов на основе метода ISD.
Заключение
Интегральный описатель структуры (ISD) и метод его получения, названный методом ISD, стали
теоретической основой для решения ряда задач, связанных с представлением и анализом абстрактных структур графов. С помощью метода ISD, наряду с получением интегрального описателя структуры, решены задачи определения изоморфизма графов и гиперграфов, вычисления полных инвариантов графов, определения ряда оценок сходства структур. В данной работе основное внимание уделено исследованию свойств устойчивых групп вершин, которые могут выделяться в графе в процессе интеграции кодов структурных различий по методу ^. На основе этих свойств в работе показана полиномиальность алгоритма определения интегрального описателя структуры. Это означает, что все алгоритмы, основанные на методе ISD, будут иметь полиномиальные оценки сложности вычислений.
В последующих исследованиях свойств устойчивых групп и в задачах анализа абстрактных структур графов будет широко применяться программа GrapЫSD. Эта программа написана для исследовательских целей и содержит средства визуализации процесса интеграции кодов структурных различий и интерактивного участия пользователя. Программа представляет пользователю возможность наблюдать формирование системы кодирования и появление устойчивых групп, управлять процессом их разрушения, исследовать схемы поиска интегральных описателей относительно разных систем кодирования.
Работа выполнена в рамках государственного задания «Наука».
СПИСОК ЛИТЕРАТУРЫ
1. Погребной В.К. Метод интеграции структурных различий в графовых моделях и его применение для описания структур // Известия Томского политехнического университета. - 2011. -Т. 318. - №5. - С. 10-16.
2. Погребной В.К. Задача определения оценок сходства структур двух графов на основе выделения общих частей // Известия
Томского политехнического университета. - 2013. - Т. 322. -№5. - С. 194-199.
3. Зыков А. А. Основы теории графов. - М.: Вузовская книга, 2004. - 664 с.
Поступила 03.09.2013 г.
UDC 519.171.1
POLYNOMIALITY OF METHOD FOR COMPUTING GRAPH STRUCTURE INTEGRAL DESCRIPTOR
V.C. Pogrebnoy, An.V. Pogrebnoy Tomsk Polytechnic University
The relevance of the research is caused by the necessity of developing the efficient method of invariant description and analysis of abstract structures of graph models. The aim of the research is to substantiate the polinomiality of the method for computing the integral descriptor of graph abstract structure proposed by the authors. The research techniques are based on application of machinery of graph theory and methods of free and dependent integration of codes of graph structural differences. The authors have introduced the notion of stable group of vertices in graph and stated the conditions of occurrence and existence of such groups at integration of structural differences codes when computing the integral structure descriptor. A number of features which disclose the appropriateness of application of the main rules of the integral structure descriptor and its polinomiality was determined for stable groups. It was ascertained on the basis of the defined features that the conditions for stable group existing are conditioned by hard limits; the vertices of different stable groups can not generate new stable groups. The authors have also defined the factor of full isolation of stable groups that predetermined considerably the efficiency of algorithm for computing the full graph structure descriptor. Polinomiality of the technique is demonstrated for the most complex case when graphs are homogeneous and contain stable groups. The authors developed Java GraphlSD software for the experimental investigations of integral structure descriptor technique and introduced the results of its operation.
Key words:
Integral structure descriptor, graph isomorphism, abstract graph structure, stable group of vertices, code integration area, algorithm polynomiality.
REFERENCES
1. Pogrebnoy V.K. Metod integratsii strukturnykh razlichiy v gra-fovykh modelyakh i ego primenenie dlya opisaniya struktur [Integration technique for structural differences in graph models and its application to describe the structures]. Bulletin of the Tomsk Polytechnic University, 2011, vol. 318, no. 5, pp. 10-16.
2. Pogrebnoy V.K. Zadacha opredeleniya otsenok skhodstva struktur dvukh grafov na osnove vydeleniya obshchikh chastey [The task of estimating the similarity of two graphs structures on the basis of general parts definition]. Bulletin of the Tomsk Polytechnic University, 2013, vol. 322, no 5, pp. 194-199.
3. Zykov A.A. Osnovy teorii grafov [Bases of graph theory]. Moscow, Vuzovskaya kniga, 2004. 664 p.
