CC BY
33
О ВЫЧИСЛИТЕЛЬНОЙ ЭФФЕКТИВНОСТИ АЛГОРИТМА СПЕКТРАЛЬНОГО РАСЩЕПЛЕНИЯ ПРОВЕРКИ ИЗОМОРФИЗМА ГРАФОВ
А.В. Пролубников, Р.Т. Файзулли Омский государственный университет
Мы рассматриваем эвристический алгоритм проверки изоморфизма графов. Алгоритм основан на последовательном расщеплении кратных собственных значений матриц, сопоставленных графам: матрицы, с которыми работает алгоритм, представляют собой модификации матриц смежности графов. Алгоритм строится на последовательном возмущении матриц и решении связанных с ними систем линейных алгебраических уравнений. Матрицы видоизме-нёны до положительно определенных, что позволяет решать системы линейных уравнений, связанные с ними. Доказывается вычислительная эффективность алгоритма в одной из принципиально тяжелых для решения задачи проверки изоморфизма графов ситуаций.
Алгоритм спектрального расщепления проверки изоморфизма графов
Задача проверки изоморфизма графов принадлежит к задачам, относительно которых нет ясности: являются ли они полиномиально разрешимыми или нет [1], тогда как известно, что она является полиномиально разрешимой для некоторых классов графов. В частности для планарных, регулярных графов, графов с ограниченной степенью вершин, графов с ограниченной кратность собственных значений из спектра матрицы смежности и некоторых других[4], [5], [6] построены эффективные алгоритмы решения этой задачи.
В задаче проверки изоморфизма графов даны два неориентированных графа ОА = (УА, ЕА) и Ов = (Ув, Ев), где УА , Ув - множества вершин графов, ЕА , Ев - множества ребер. Везде далее графы Оа и Ов - связные. Предполагается, что
У = У
A B
\ea\ \ч
Задача изоморфизма графов формулируется следующим образом: существует ли биективное отображение р: УА ^ Ув , такое, что (/', у) е ЕА, то (р(/), р(у)) е Ед ?
Алгоритм решения задачи определения изоморфизма графов работает с видоизмененными матрицами смежности графов. Пусть А0 - матрица
смежности графа0А, то есть А0 = (а..), где
и
ГШ у) е еа , у [ 0, иначе.
В0 = (Ь,) - матрица смежности графа Ов. По матрице А0 строим матрицу ПА :
(d + d-.
0 ^
d + d
0 ... d + d
n
Da - диагональная матрица, где:
di = 2 aij ,
j *i
то есть di - степень вершины i графа GA
d = тах di.
1<г'<п
Аналогично строится матрица по матрице
смежности графа Ов. Рассматриваемые далее матрицы
А = А0 + Да0, в = в0 + Бво (1)
- матрицы, с которыми будет работать алгоритм, являются симметрическими положительно определенными матрицами.
Если графы О а = (Уа , Еа ) и О в = (Ув, Ев) изоморфны, то соответствующие им матрицы описанного типа могут быть получены одна из другой перестановкой строк с одновременной перестановкой столбцов с теми же номерами. Таким образом, задача проверки изоморфизма двух графов может быть сформулирована как частный случай задачи Фробениуса: может ли быть получена матрица в из матрицы А последовательной перестановкой строк с одновременной перестановкой столбцов?
Для произвольной матрицы перестановка строк с номерами i и у эквивалентна ее умножению на матрицу перестановки Р, . Перестановка столбцов матрицы эквивалентна ее умножению справа на ту же матрицу. При умножении вектора на матрицу Р, как справа,
так и слева происходит перестановка компонент вектора с номерами i и у .
Рассмотрим две системы уравнений следующего
вида:
Ах = е,, ву = ек, (2)
где векторы ei = (0,... ,0,1,0,... ,0) - базисные векторы
в пространстве Яп, матрицы А ив - описанного выше вида. Обе системы уравнений имеют решение, и решение единственно, так как А и в - матрицы с диагональным преобладанием, и, следовательно, их
0
0
0
0
определители не равны нулю. Пусть далее хр -
решение системы линейных уравнений Ах = в р, ук - решение системы линейных уравнений Ву = вк.
Отметим, что решив системы уравнений (2), мы получим обратные матрицы для матриц А и В . Так, для I -ой компоненты вектора хр- верно:
Хр = Ау/1А |, где Ау - алгебраическое дополнение элемента а ^ матрицы А . То есть векторы решений Хр, 1 < р < п , являются столбцами обратной к А матрицы.
Если В = РркАРрк , то для решений систем
уравнений (2) должно выполняться равенство: хРрк = у . Действительно:
(Ах = вр) ~(РркАх = Ррквр) ~ ~(РркАхРрк = РркврРрк) ~(РркАРркх = вр)~ ~ (РркАРркхРрк = врРрк) ~(ВхРрк = вк). То есть хРрк = у .
Если матрица В получена из матрицы А многократной одновременной перестановкой строк и столбцов, то:
В = Рр к ... Рр к АРр к ... Рр к , и, соответственно,
ргкг р1к1 р1к1 р1к1
хР рк ...Ррк = у . р1к1 р1к1
Таким образом, если в системе уравнений (2) при фиксированном р индекс к будет пробегать значения от 1 до п, то векторы хр_и_ук - соответствующие друг другу решения полученных систем (2) будут совпадать с точностью до перестановки компонент только в том случае, если строке матрицы А соответствует строка к
матрицы В . То есть элементы строки к матрицы В есть переставленные элементы строки / матрицы А . То же верно и для столбцов матриц.
При нахождении в матрице В строки и столбца, соответствующих строке и столбцу с номером I матрицы А , на каждой итерации алгоритма будем последовательно производить возмущения ее диагональных элементов.
Алгоритм работает с симметрическими матрицами, которые могут быть приведены к диагональной форме с помощью ортогональных преобразований. То есть:
А = иААПА, В = и ВВи В,
где А , В - диагональные матрицы, с собственными значениями на диагонали, и А , и В - матрицы ортогональных преобразований, и 'А, и 'В -транспонированные к ним матрицы. Диагональные элементы А , В - собственные числа матриц А и В .
Спектры матриц смежности изоморфных графов совпадают[2]. Спектры матриц вида (1) изоморфных графов будут совпадать тоже, так как описанная выше процедура замены нулевых диагональных элементов матрицы смежности приводит лишь к сдвигу спектра матрицы. Причем, если матрицы соответствуют изоморфным графам, то спектры матриц будут совпадать и после сдвига.
Если спектр каждой матрицы является простым, или кратность собственных значений мала, то задача проверки изоморфизма графов разрешима однозначно
путем сопоставления строк матриц и А
иВ
А , В [3]. Матрицы иА, следующий вид:
и В могут иметь, например,
' V» ^2 V13 ^4 V. ^
V 22 V 23 V 24 V 25
N А и V31 V32 ^3 V34 V35
V41 V 42 V 43 V 44 V 45
V V51 ^2 V53 V54 V55 ,
' ^31 V32 V33 V34 V35 ^
V 21 V 22 V23 V 24 ^5
^2 V13 V14 V15 ,
^51 V52 ^3 V54 V55
ч V 41 V 42 V43 V 44 V45 )
иВ =
где {V>П=1 - собственные векторы матрицы А . То
есть строки матриц, составленных из собственных векторов, те же с точностью до перестановки, задающей изоморфизм.
Основные трудности возникают при рассмотрении графов, спектры которых содержат кратные собственные значения. Возмущая матрицы в ходе итераций алгоритма, мы будем получать возмущение спектра матриц, при котором происходит расщепление кратных собственных значений, что позволит установить однозначное соответствие между строками и столбцами матриц.
В результате, после того как кратные собственные значения матриц А и В в ходе работы алгоритма будут расщеплены, будет возможно однозначное определение перестановки, задающей биекцию ф, устанавливающую изоморфизм графов
Оа = (УА, Кл > и Ов =<Ув , Ев > .
Приведенный ниже алгоритм, однако, работает не с собственными векторами, а с векторами, являющимися столбцами обратной матрицы к матрице А , и в качестве характеристик, на основе которых строится соответствие между вершинами графов, рассматриваются компоненты и нормы этих векторов. Норма:
х2к
Ы = хк
V к=1
В этом случае невозможность установления однозначного соответствия возникает при наличии среди этих векторов группы векторов, которые могут
быть получены друг из друга с помощью перестановки их компонент, и, следовательно, имеют одинаковые нормы. При этом компоненты этих векторов, образующие диагональ обратной матрицы, совпадают. Однако, при возмущении матрицы, влекущем расщепление ее собственных значений, происходит как возмущение ее собственных векторов, так и возмущение вектор-столбцов обратной к А матрицы. В результате описанные выше группы векторов с одинаковыми нормами расщепляются.
Пусть х/ и хк - вектор-столбцы обратной матрицы такие, что xi = Рхк для некоторой матрицы перестановки Р , и их компоненты, лежащие на диагонали обратной матрицы, равны: х/ = х¡к. В
ходе работы алгоритма на каждой итерации будем работать не с исходными, но уже с возмущенными в ходе предыдущих итераций алгоритма матрицами. Получив соответствие строки у матрицы А у строке ¡ матрицы в у на у -ой итерации, а также столбцами с теми же номерами, рассматриваем далее возмущенные матрицы А1+1 и в1+1:
Ау+1 = Ау +Б]Е}, в;+1 = в +Б]ЕК.
Возмущения производим с помощью матриц
Еп с элементами еу :
ч и
1, i = ] = ч,
0, иначе.
Численные эксперименты показывают, что расщепление спектра матриц, необходимое для определения соответствующих друг другу строк и столбцов матриц, происходит значительно раньше заключительной итерации. К примеру, для графов, представляющих собой решетку на торе с числом вершин п от 9 до 400, уже на итерациях с номером 4п не требовалось дальнейших возмущений диагональных элементов матриц А и в для установления однозначного соответствия между вершинами графов.
В результате, если матрица в может быть получена из матрицы А последовательными одновременными перестановками строк и столбцов, то в ходе работы алгоритма, пока у пробегает значения от 1 до п , будет получена перестановка строк и столбцов матрицы А - перестановка Р . Перестановка Р является, вообще говоря, одной из возможных перестановок. То есть, как подстановка Р выглядит следующим образом:
7} ву+1 = ву
Р =
1 2
, ¡1 к 2
кп
где ¡ у - номер строки матрицы в , полученный на у -й итерации. Перестановка Р задает такую перенумерацию вершин графаОА , то есть отобра-
жение р: УА ^ Ув, при которой матрицы А и в , представляющие графы, совпали бы, что возможно только в случае изоморфности графов.
Алгоритм спектрального расщепления проверки изоморфизма графов
Шаг 0. А0: = А, у: = 1. Р(/): = 0, 1 < i < п.
Шаг 1. Если у < п, то Шаг 1.1., иначе - работу алгоритма завершить.
Шаг 1.1.
(d+^ )2 + п (II у+2+2 I-^+^).
Шаг 1.2. Ау: = Ау+1 +еЕ .
Шаг
2.
Решение
системы
уравне-
ний А1 х = е :. х . - полученное решение. }
Шаг 3. к: = 1. Если к < п , то - Шаг 3.1, иначе перейти на Шаг 4.
Шаг 3.1. вк: = ву+1 +е]Ек .
Шаг 3.2. Решение системы уравнений
вку = ек . ук - полученное решение.
Шаг 3.3. к: = к +1. Перейти на Шаг 3.
Шаг 4. Сравнение норм векторов ху и ук ,
к : л/ < у : Р(0 ф к .
Если •к : ||ху | ф ||ук ||, то графы ОА и Ов неизоморфны. Работу алгоритма завершить. Если Эк : ||ху || = ||ук || и ху = уЦ, то Р(у): = к (вершине графа О А ставим в соответствие вершину к графа Ов), вк: = ву+1 +еуЕк .
Шаг 5. ]: = у +1.
Шаг 6. Решение систем уравнений А1 х = е1: V/' Р(1) ф i. х1 - полученные решения.
Шаг 7. Если 'р, I: ||х^| ф ||хр ||, то перейти на Шаг
8, иначе перейти на Шаг 1.
Шаг 8. Если у < п , то решение систем уравнений для всех i таких, что Р(0 ф р .
Шаг 9. Сравнение норм полученных решений:
Если Э р, I: ||хг || = ||ур 11, л/ Р(1) ф i, Л/ Р(/) ф р , то Р(1): = р .
В действительности в связи с погрешностями вычислений проверяются близости норм, а не их равенства, и хотя алгоритм спектрального расщепления проверки изоморфизма графов хорошо зарекомендовал себя в ходе вычислительных экспериментов, следует отметить, что для произвольных графов нельзя гарантировать получения решения задачи изоморфиз-
п
ма графов с помощью представленных схем алгоритма. Так, представляется допустимой ситуация, когда проверки векторов на равенство (близость) норм на шаге 4 алгоритма и проверки равенства
(близости) компонент хр и ук недостаточно для ответа на вопрос, может ли вектор ук быть получен из вектора хр последовательной перестановкой его компонент.
При установлении соответствия между вершинами графов, алгоритм спектрального расщепления:
1) проверяет равенство, а в связи с погрешностями вычислений, в действительности проверяется близость норм решений систем линейных уравнений Ах = вр и Ву = вр, 1 < р < п ;
2) в случае равенства норм векторов х р и
ук проверяется равенство (близость) компонент
рк этих векторов хр и ук .
Таким образом, для ответа на вопрос об эффективности алгоритма спектрального расщепления для проверки изоморфизма графов ОА иОВ необходимо оценить насколько будут отличаться после последовательных возмущений матрицы А нормы векторов х ^ и хр, соответствующих разным вершинам графа Оа .
Эффективность алгоритма
Везде далее будем использовать следующие обозначения:
А ^, 1 < I < п - вектор, компоненты которого являются элементами / -й строки матрицы А, и -обратная матрица к матрице А , и и 1 < / < п - вектор, компоненты которого являются элементами / -го столбца матрицы и (и) = х1), а1 - угол между вектором А^ и ортом в^,
в^ - угол между вектором и^ и ортом в^, - угол между вектором А^ и ортом и^ .
Будем также считать, что максимальная степень вершин графа (ё) не менее 3. Ситуация, когда максимальная степень вершин графа менее 3, не представляет интереса с точки зрения проверки изоморфизма графов.
При возмущении диагонального элемента а рр матрицы А возрастает р -я компонента вектора Ар. То есть при возмущении диагональных элементов матрицы происходит приближение векторов {Ар}пр=1 к соответствующим ортам вр
декартовой системы координат в Яп.
Действительно, пусть а р - угол между вектором Ар и ортом вр, тогда:
cos ар =
(Ap > ep > = арр
ир
арр = dj + d+sj ■
A К a ]
ким образом,
+ I a2]k =J(dp + d + Sp )2 + dp. Та-
k*j
cos a =
dp + d + Sp
,J(dp + d + sp)] + d p 1
1 + -
d
(dp + d + Sp)
d = max dp, и, следовательно, при возмущении диа-
i<j<«
тональных элементов матрицы cosap-, j = 1 < j< n, будет приближаться к единице, а сами углы а р, соответственно, к нулю.
Столбцы обратной матрицы к A матрицы U -
векторы {U р }П=1 образуют систему векторов, каждый вектор U из которых ортогонален всем векторам Ak, k Фр, и (Ap, Up) = 1, так как UA = I.
Векторы {U p }П=1 близки по направлению к векторам {Ap}n=1. Действительно, так как Up L Ak,k Ф p , то Up определяет гиперплоскость в пространстве Rn, в которой лежат векторы Ak, k Ф i. В процессе возмущения матрицы векторы {Ak}n=1 приближаются к ортам {ek}n=\, а значит и гиперплоскость, натянутая на векторы Ak, k Ф p, приближается к гиперплоскости Rn\Lp, гдеLp - одномерное линейное пространство, натянутое на орт ep. Таким образом, и вектор Up будет приближаться к реперу e р, и,
следовательно, угол между векторами Up и Ap будет уменьшаться
Лемма 1. Для произвольной матрицы вида (1), представляющей граф, для любого р
3d
cosy ■ >
р Л0\\л,.
Доказательство
Угол вр равен углу между гиперплоскостью Яп\Ьр и гиперплоскостьюЯп\Ь(ир), где Ь(ир) - одномерное линейное пространство, натянутое на вектор ГиперплоскостьЯп\Р(и р) натянута на векто-
U
Uj L\R"\L(Uj)),
ры Ак, к * j, так как e, i{r"\Lj ). Следовательно Ук, к * j, б к > в j, а
dk + d
значит
cos в , > max
к* j
V(dk + d )2
+ dк
(Uj, e,) Ujj
cos в i = |. .... .. = ttAt , значит
Uj
U,
\\u,
■ > max
к * j
dk + d
1
3
M + d )2
+ dк
cos у, =■
(U,, A,)
i+9 ^
3. 10 jj'
Известен следующий результат Фань-Цзы
[7].
Пусть А - матрица с диагональным преобладанием, и si = \ац1 - ^ ак . Если ^ > 0, то
к ф/
существует обратная к матрице А матрица А', и
К1 <-, 1 <i<п.
В нашем случае, матрица А , ставящаяся в соответствие графу ОА , - матрица с диагональным преобладанием, и
si = d + di - di = d ; и у\ = и у , следовательно,
u,, , 11 d
cosy, >
3d
10l A/Uii
10 A,
Лемма 1 доказана.
Пусть далее векторы и углы со штрихом -векторы и углы, в которые переходят соответствующие им нештрихованные векторы и углы в ходе итерации алгоритма.
Лемма 2
| cos - cos y' | ^
Доказательство
cos Yi - cos =
ии ||a;|||u;||
1 Г 1
Ai
ikiiiu; II nUNI и; II
(u,'|| -II
Ai Ui
По правилу треугольника из Ui - U' + U' = Ui следует:
\\Ui - Ui 11 + |\U' 1|> Ц |, или
Ui -U; 11> HI - lUi'l |.
То есть
| cos Yi - cos/i | ^
1
U'. - U ,, ii
. (4)
(Ui, A j) = 0, тогда как
(U., A j) = (U., A, + в, e,) = (U., A,,) + (U', s ,e,) = 0, сл едовательно:
(\;, ai ) = -в\.
И поэтому
(Ui - и;, Aj) = SjU;.
Из (4) и (5) следует:
(5)
s . U'..
j ■j
cos/■ - cos/ '■ ^ -.
Ai pJ pi AJ|cosn|
где ^ - угол между вектором U' - Ui и A,.
У к * j : (U; - Ui, Ак) = 0, а так как векторы {Ак }пк=1 образуют базис вRn, то, учитывая (5), получаем, что при возмущении вектора А, происходит изменение Ui в одном из двух возможных направлений - или в направлении вектора А,, и тогда cos п = 1, или в направлении противоположном вектору А,, и тогда cos п = -1, и значит,
cos y■ - cos/'А ^ ]
s U'..
J и
hllh ||ui hi
s U'..
J и
|ui III A. Ill i 1
cos/j
s . j
I
То есть cos/■ - cos/' ^ n-jr. Лемма 2 доказана.
Лемма 3
1
1
1
1
1
>
и
3
е
1
1
cos Yj =
1 + J JJ
Л'
cosYj
cos y j
1
sPAp--
И
cosYj
Доказательство 1
cos y j = |—jj—
, cosy,- = и—
1 л'
AW J
cos Yj - cos Y J = iaiiiu, 1 Л,|и (Uj|| - U j||)(6) Векторы U - и U J имеют следующие ком-
поненты:
Up =
Лемма 3 доказана.
Вернемся теперь к исходной задаче. Для того, чтобы определить эффективность алгоритма, нам необходимо оценить снизу разность норм векторов U' и U р после увеличения р -го диагонального элемента матрицы на некоторое s р, при условии равенства норм векторов Ui и Up до возмущения, и, соответственно, cos yi = cos ур.
f Aj1 A AJn , UJ = f AJ1 AJn M-b|= ||AJ||cos Yj -|| AJII cos y)
l И " И V U A)" '" И V ||AJ|| cosYj||A J|| cos yJ
где А'рк , |А', соответственно, - алгебраическое дополнение к элементу а кр и определитель возмущенной матрицы. Но, так как на р -ой итерации возмущается только элемент а рр матрицы А , то У к А'рк = А рк . Значит
U J =
fА А ^
AJ1 Ajn
.И'""' и
Следовательно,
А 2 У ■ -
УИ2 1
А 2
У Ajk
У И12
1
У и
k=1
I___1
И" Й1Л
У Л
k=1
И-| и
. \Л'\\И .
Л J cos y J -
cosy¿
(7)
До
возмущения матрицы: ||л, || = ||Лг-1|. так как возмущение не затрагивает j-ю
+ s pi.
строку матрицы. Пусть spi такое, что ||Aj|| = Тогда (7) эквивалентно
(IAi|| + sp) cos y'j -IЛг|| cos y'
л;|||л; ii
Рассмотрим числитель выражения выше. Пусть S = Ai || + spi) cos y) —IAi || cos y' = |A ||(cos y' -- cos y') + sji cos y'} = ||Ai ||(cos y'} - cos yp) + Щ || X X (cos yi - cos y') + Si cos y'p ,
Возмущение определителя: |Л' | = |Л| + s-A-И таким образом:
sPAP
У A2J
^ И' И
Возвращаясь к (6), получаем:
2 SjAJj
cos yJ - cosy, = ñ—¡п.—1—iñ—ir У Apk —¡¡—г
|AJII|UJ|A;|U;Nу Ia 'И
И1=и ц
УЛ
k=1
УЛ
k=1
следовательно,
cos y) - cos Yj = -
А, значит,
sjAjj
ИМ У A2
sjAjj
cos y ,■
| k =1
S =||Лг||(cos y'j - cos Yp),
52 =|И ||(cos Yi - cos y) ),
53 =sji cos y) , (S = Sj + S2 + S3).
Лемма 4
N > |s 2I.
Доказательство
По лемме 3
S1 =| И
по лемме 2
SpApp
1 +„ .. cosy •
A' ■'
S 2 <
IA
Следовательно
= s
p •
1
и
>
>
и
s
( \
S1 + S2 >IKII1 -eJ =IKIIC0S^j +
+ s
ИИ
A'
cosy , -s,.
cos y, >
3d
10 A,
по лемме 1.
Таким образом,
3d 3d S1 + S 2 >^= + s
<1 A-
Vrn 1 л/Ю Al |A|
- s.
Оценим, каким должно быть е у для того, чтобы выполнялось необходимое нам неравенство 1
е >-.
п
Пусть
е ц = + - = — + 8
оценки снизу модуля разности |U' 11- Ц || учитывать только S3 . То есть
S = S1 + S2 + S3 > S3 = s. cos Yj. По лемме 1
cosy- >
Если
sH > -, то
n
3d
10 a,
S3 >
3d
3d
V1o| Ai|n VI0^V4
n 2 + n
Тогда
U' 11- U >
(||4|| + ^ )cos у) -||4||coy
A
= |S|
3d
\\A1
|ain||a; ii vi0w4n2+n\\a'
где £ > 0. Откуда получаем, что
S1 =
■^A(n)2 + 2A(n)^(d + d, )2 + d, + (d + dj)
1
- (d + d,), где A(n) = — + S, а S - мало.
si =
J(d + d, )2 + 2A(n)yj(d + d, )2 + d, + A~(ñ)2
- (d + d ,) < (d + d , \ 1 + 2A(n) —— +
л/19 A(n)2
-(d + d, ).
Таким образом,
1
16 16
1-1-3
■ 1*1 о 1 • 1 • 5 i
s, <-S'(n)--S (n) +-S3(n) +
1 2 2 - 4 2 - 4 - 6
1-1-3 - 5 S 4 ( )
■S (n),
2 - 4 - 6 - 8
где S' (n) = 2A(n)^ +A(n)2
16
16
поэтому
3d
Итак, s .■ <__
лД0
e e 3d 3d ЦП An S
S1 + S2 >^= + s1—j= -^-гц-TT-si > S0,
■Л0 1 лЯО ||a Ц |a|
где So > 0 .
То есть
S1 + S2 >So, >|S2I + So
> 21. Лемма 4 доказана.
В итоге, доказав леммы выше, мы можем отбросить слагаемые и £2, и для построения
л/10п(4п2 + n к/ 4n
то есть
| i-|u;||
V10n(4n2 + n)n2 + n
Таким образом, нами доказано следующее утверждение.
Утверждение
такое что:
Если |uJ = ЦЦ и s,
а/К +sí)2 + dí 4a2 + dí > n ■
то после -й итерации алгоритма
N I- N
9
V10n(4n2 + n)
Таким образом, если погрешность вычислений будет, к примеру, на порядок меньше разницы точных
значений норм векторов 11и' || и |и' 11, то на каждой итерации алгоритма возможно выделение одного вектора из группы векторов с одинаковыми нормами.
Пусть решение систем линейных уравнений осуществляется с помощью метода Зейделя. Приближения к точному решению на т -й итерации метода Зейделя могут быть оценены как
и; - и\ (m)|| <¿ums°
U[■ - ¿7'/m)|| <MmS0,
] j
здесь 8° ,8у - погрешности начальных приближений и', и), и'(т), и,]{т) - приближения и', и} на т -й итерации метода Зейделя, а
1
S
3
>
>
9
>
n
2
9
>
n
>
и
U = max
dk - dk
1<k<n d + dk - dk
k-1
где
dk =I
i=1
и, следовательно, u < —, и
U' - U'(m)
] m
up -Up(m)
<-L s0. ]m
б0 = шах{б° ,бр0} - погрешность начального
приближения. То есть необходимо выполнение неравенства
9
— S0 < ] m
л/iQn ] (4n ] + n)]
Так, это неравенство будет выполнено уже
при
«1 л бо , 2542 т = 5^ п + 1оя2 б + 1оЯ2——.
Поэтому можно утверждать о разрешимости проблемы расщепления групп векторов с одинаковыми нормами и вычислительной устойчивости алгоритма в данной ситуации.
Полученная оценка на число итераций, необходимых для эффективной работы алгоритма, получила подтверждение в ходе вычислительного эксперимента для дешифрования шифра двойной перестановки с помощью алгоритма спектрального расщепления [8]. При решении задачи проверки изоморфизма не взвешенных неориентированных графов алгоритм работает с целочисленными, за исключением диагональных элементов, матрицами. При решении задачи дешифрования шифра
двойной перестановки алгоритм работает с нецелочисленными матрицами, что ведет к возрастанию погрешностей вычислений. Несмотря на это, оценка числа итераций, необходимого для приближения решений, достаточного для установления однозначного соответствия, подтвердилась.
Литература
1. Гэри М., Джонсон Д. Вычислительные машины и труднорешаемые задачи // М.: Мир, 198]. ]. Цветкович Д. и др. Спектры графов. Теория и применение // Киев: Наукова думка, 1984.
3. Кикина А.Ю., Файзуллин Р.Т. Алгоритм проверки изоморфности графов // деп. ВИНИТИ 21.06.95 1789-В95.
4. Hopcroft, Wong. A linear time algorithm for isomorphism of planar graphs// Proceedings of the Sixth Annual ACM Symposium on Theory of Computing, p. 172-184, 1974.
5. Luks. Isomorphism of graphs of bounded valence can be tested in polynomial time // Proceedings of the 21-st IEEE FOCS Symp., 4], 49, 1980.
6. Hoffmann. Group-Theoretic Algorithms and Graph Isomorphism// Lecture Notes in Computer Science (Chapter V). P.127-138, 198].
7. Беллман Р. Введение в теорию матриц // Москва: Наука, 1976.
8. Пролубников А.В., Файзуллин Р.Т. Эвристический алгоритм дешифрования шифра двойной перестановки // Математические структуры и моделирование: Сб. научн. тр. Под ред. А.К. Гуца. Омск: Омск. гос. ун-т, ]00]. Вып. 9.
1
