УДК 519.1
01.00.00 Физико-математические науки
РАЗМЕЩЕНИЕ ЦЕНТРОВ НА МНОГОВЗВЕШЕННЫХ ПРЕДФРАКТАЛЬНЫХ ГРАФАХ
Кочкаров Расул Ахматович Токова Алла Аскеровна Кочкарова Асият Нерчуковна
Северо-Кавказская государственная гуманитарно-технологическая академия, Черкесск, Россия
В работе предложена многокритериальная постановка задачи размещения центров на многовзвешенном предфрактальном графе. Приведена оценка радиального критерия предфрактального графа, порожденного затравкой-звездой. Предложен полиномиальный алгоритм размещения центра предфрактального графа при сохранении смежности старых ребер. Проведена оценка вычислительной сложности алгоритма и рассмотрен пример работы алгоритма
Ключевые слова: МНОГОКРИТЕРИАЛЬНАЯ ЗАДАЧА, МНОГОВЗВЕШЕННЫЙ ГРАФ, ПРЕДФ-РАКТАЛЬНЫЙ ГРАФ, РАЗМЕЩЕНИЕ ЦЕНТРА, ВЫЧИСЛИТЕЛЬНАЯ СЛОЖНОСТЬ
Рок 10.21515/1990-4665-121-108
UDC 519.1 Physics and Math
CENTERS PLACEMENT ON MANY-WEIGHTED PREFRACTAL GRAPHS
Kochkarov Rasul Akhmatovich
Tokova Alla Askerovna
Kochkarova Asiyat Nerchukovna North-Caucasian State Humanitarian Technological Academy, Cherkessk, Russia
Multicriterial formulation for centers placement problem on many-weighted prefractal graph is proposed. Estimation of the radial criterion of prefractal graph generated by seed-star is shown. Polynomial algorithm centers placement on prefractal graph with preserving contiguity old edges is suggested. Estimation of computational complexity of the algorithm and the example of the work algorithm are considered
Keywords: MULTICRITERIAL PROBLEM, MANY-WEIGHTED GRAPH, PREFRACTAL GRAPH, CENTER PLACEMENT, COMPUTATIONAL COMPLEXITY
Постановка задачи
Рассматривается предфрактальный граф GL = (Уь, EL), порожденный множеством затравок Н = {Н}. Каждому ребру e() е Еь приписано M дей-
ствительных чисел wt(e) = wt(e(l})e (6l la,6l 1b), i = 1,M, где l = 1,L - ранг
a
ребра, a > 0, b > 0 и 0 <6 < — - коэффициент подобия.
b
Пусть x - подмножество, состоящее из p вершин множества VL предфрактального графа GL = V, EL). Через d (х, vl) обозначаем наикратчайшее из расстояний между вершинами множества x и вершиной vl, т.е. d(x, ) = min[d(у , )].
Число разделения х) для множества вершин х определяется следующим образом - х) = тах[й(х,V.)]. Множество х*, для которого
vJеV1 3
х*) = Ш1п[х)], называется р -центром предфрактального графа Оь.
На рис. 1 изображен предфрактальный граф G3 с сохранением смежности старых ребер, веса ребер взвешены единицами. Кратчайшая цепь между множеством х и вершиной V выделена пунктирной линией. Длина кратчайшей цепи или наикратчайшее из расстояний й (х, V') = 4. Число разделения р -центра 5(х*) = 2.
Всевозможные ^-центры {х} предфрактального графа GL образуют
жество допустимых решений (МДР) X = X (GL) = {х} [1, 2].
мно-
На множестве х определим векторно-целевую функцию (ВЦФ): F( х) = (Fi( х),...
, FM (х),..., F2M (х), F2M +1(х), F2M + 2 (х)), (1)
Ft(х) = st(х) ® min , i = 1,2,...,M, (2)
где st (х) - число разделения множества х;
Г
FM+1 (х) = Е P,t ® mln, i = 1,2,...,M :
t=1
Научный журнал КубГАУ, №121(07), 2016года Г t - радиус t -вершины p -центра;
F2 M +1
(X) = Н ® min , (4)
H - количество типов центров (количество центров разного ранга);
F2 M + 2
(x) = p ® min , (5)
p - количество вершин, составляющих p-центр.
Все критерии (2)-(5) ВЦФ (1) имеют конкретную содержательную интерпретацию. Веса, приписанные ребрам предфрактального графа GL, могут отражать как конкретные ограничения (время, расстояние), налагаемые на систему служб (аварийные, пожарные депо, полицейские участки, больницы), так и общие затраты, выражаемые в условных единицах.
Полученное значение числа p критерия (5) будет наименьшим числом аварийных служб, а р-центр - их оптимальным размещением, удовлетворяющим предъявляемым требованиям. Поскольку M веса одного ребра несравнимы, паретовское множество включает в себя p -центры не одного набора весов, а каждого из M . Критерий F (х) минимизирует число разделения, а критерий FM+j (х) минимизирует сумму радиусов p -центра по
i -му набору весов, i = 1,M.
Оценка радиального критерия предфрактальных графов
Рассмотрим вначале некоторые определения и обозначения. Длина кратчайшего пути, соединяющего пару вершин w,v eW, называется расстоянием между вершинами w и v и обозначается через d(w, v). Для фиксированной вершины w e W величина e(w) = max d(w, v) называется экс-
veW
центриситетом вершины w e W. Эксцентриситет, максимальный среди
всех эксцентриситетов вершин графа H = ^,О), называется диаметром графа Н и обозначается через й (Н), т.е. й (Н) = тах^У) [3].
wеW
Если пара вершин и,wеW соединяется кратчайшим путем длины й(и,w) = й(Н), то этот путь называется диаметральным. Вершина w называется периферийной, если е^) = й(Н). Радиус графа Н обозначается через р( Н) и вычисляется по формуле р(Н) = тте^).
wеW
Рассмотрим предфрактальный граф Оь, порожденный затравкой-
звездой Н , всем ребрам которого приписан вес равный единице, а коэффициент подобия равен 1. Вычислим верхние и нижние оценки радиуса предфрактального графа 01, I = 1,2,...,Ь траектории G1,G2,..., Оь.
Радиус графа G1 равен единице, так как G1 = Н есть затравка-звезда: р^) = 1. Центром графа G1 является вершина, смежная со всеми остальными вершинами. Расстояние от центра до любой другой вершины равно единице. В случае, если смежность старых ребер графа G2 сохраняется, радиус рт1п ^2) = 2 . Если старые ребра не пересекаются и инцидентны висячим вершинам подграф-затравок, радиус ртах ^2) = 4. Верхняя и нижняя оценки радиуса предфрактального графа G2 с произвольной смежностью старых ребер равна: 2 < р(02) < 4.
На рис. 2 представлены предфрактальные графы G1 и G2, порожденный 4-вершинной затравкой-звездой. Для графа G2 рассмотрены два варианта порождения: сохранение смежности и непересечение старых ребер. Центры обведены малыми пунктирными окружностями.
По аналогии с предыдущими вычислениями, верхняя и нижняя оценки радиуса графа О3 равны: 3 < р(О3) < 13 (см. рис. 3.а и 3.б).
Продолжая процесс вычисления верхних и нижних оценок радиуса, получим, что pmax(ОL) = (3Ь -1)/2 и рт,п(О1) = 1, для всех 1 = 1,2,...,Ь.
Верхние и нижние оценки радиуса графа О1, 1 = 1,2,..., Ь равны:
pmax (О1) = 1;
PmЛО2) = 3 +1;
Ри^Оь) = 32 + 3 + 1;
Лшт (О1 ) =1; ЛшП (О 2) = 2; РпаП (О3) =3;
Лmax(Оl ) = 31 -1 + 31 -2 + ... + 3 + 1 = (31 - 1)/2 ; Pmin (О1 ) = 1
Pmax(ОL ) = (3Ь - 1) / 2; РШп(Оь ) = Ь .
Тогда оценка радиуса предфрактального графа О1 выражается неравенством: 1 < р(О1) < (31 -1)/2 для всех 1 = 1,2,...,Ь. Верна следующая теорема.
Теорема 1. Для всякого предфрактального графа Gl, порождённого затравкой-звездой, радиус I < р(^) < (31 -1)/2 для всех I = 1,2,...,Ь. ◄
Вычислим далее верхнюю и нижнюю оценки диаметра предфрак-тального графа Gl, I = 1,2,...,Ь:
й max(G1) = 2;
йтах (G2) = 2(3 + 1):
йтах (Gз) = 2(32 + 3 + 1)
й тп(^) = 2;
й m1n(G 2) = 4;
йтп 3) = 6;
й тах (G I) = 2(31-1 + 31-2 +... + 3 +1) = 31 -1; ) = 2 • I;
йтах (GЬ ) = 3 Ь - 1;
йm1n(Gь ) = 2 • Ь .
Тогда оценка диаметра предфрактального графа Gl выражается неравенством 21 < й(Gl) < 31 -1 для всех I = 1,2,...,Ь. Верна следующая теорема.
Теорема 2. Для всякого предфрактального графа порождённого затравкой-звездой, диаметр 21 < <(01) < 31 -1 для всех I = 1,2,...,Ь. ◄
Примечание 1. Оценки радиуса I < р(01) < (31 -1)/2 и диаметра 21 < <(01) < 31 -1, I = 1,2,...,Ь верны для всякого предфрактального графа 01, порождённого п-вершинной затравкой-звездой. ◄
На рис. 2 и 3 цепи (выделены пунктирными линиями) между вершинами у и у" отражают значение радиуса, а между у" и V"' - значение диаметра.
Алгоритм а0 размещения центра предфрактального графа при сохранении смежности старых ребер
Рассмотрим предфрактальный граф ОЬ, порожденный множеством затравок Н = {Н}, смежность старых ребер которого сохраняется. Опишем применение алгоритма для одного фиксированного набора весов wl (е) из
М и далее обобщим для каждого г -го набора, г = 1,М [4, 5]. Алгоритм а0
начинает свою работу с подграф-затравок Ь -го ранга Ь), = 1, пЬ 1. На
последнем шаге порождения предфрактального графа GL каждая вершина графа Glбыла замещена затравкой H. Поскольку при порождении предфрактального графа действует правило сохранения смежности старых ребер, к каждой вершине GL _ склеиваются затравки одной из своих вершин.
Рассмотрим подграф-затравку z(L), так как затравка H является сильно связанной, то для всякой ее вершины можно найти путь к любой другой. Найдем кратчайшие цепи от «общей» вершины x(L) до оставшихся (n - 1) вершин v(L) подграф-затравки z(L). Среди кратчайших путей выбираем максимальный и определяем число s (x1L)) = max [ d (x1L), v(L))], которое
Л=l,n_l h
называется числом разделения вершины x(L) подграф-затравки z(L). Далее поочередно осуществляется поиск чисел разделения
s(xS(L)) = max[d(xSL), v(L))], sL = 1, nL-1, где d(xSL), xSL)) = 0 для общих вер-
L jL =1,n_1 L L L L
шин всех подграф-затравок Ь -го ранга. Поиск кратчайших путей осуществляется с помощью известного алгоритма Дейкстры, оформленного в виде процедуры.
Рассмотрим далее подграф-затравки (Ь - 1)-го ранга Ь-1),
5Ь-1 = 1, пЬ-2. Каждая из них в процессе порождения предфрактального графа ОЬ-1 была склеена к вершинам предыдущего в траектории графа ОЬ-2, так как действует правило сохранения смежности старых ребер. Тогда каждая подграф-затравка (Ь -1) -го ранга х^-1 также имеет одну общую
вершину с соответствующими затравками (Ь - 2) -го ранга Ь-2),
SL_2 = 1 П
L-3
Найдем числа разделения для каждой общей вершины ^Ь-1) подграф-
затравок (Ь -1) -го ранга х^4, 5Ь-1 = 1, пЬ 2 следующим образом:
5(х(Ь-1)) = тах[й(х(Ь-1), V(Ь-1)) + 5(х(Ь))], где й(х<Ь-1), х<Ь-1)) = 0. То есть нахо-
V £¿-1 / I- V £¿-1 Л-1 У 4 /Л 4 5Ь-1 5Ь-1 У
дим кратчайшие пути от общей вершины Ь-1) до оставшихся (п -1) вершин подграф-затравки Ь). Добавляем к длине кратчайшего пути й(х£Ь"-1),V(Ь-1)) соответствующее число разделения 5(х^)), найденное для
подграф-затравок предыдущего ранга, и среди получившихся сумм выбираем максимальную. Сумма й (х£Ь-1), х£Ь-1)) + 5 (х£Ь)) = 5 (х£Ь)), так как
d(х^L х^L = 0. Указанным способом находим числа разделения 5(xf))
°L-1 L-1
для общих вершин х^), sl = 1, nl 1 подграф-затравок z{s\) до 2-го ранга включительно, то есть для всех l = L, L -1,...,2. На последнем шаге l = 2 найдены числа разделения s(xs(2)), s2 = 1, n для общих вершин xs(2) подграф-
затравок 2-го ранга zS2) и одной подграф-затравки первого ранга zS1) = z1(1). Подграф-затравка z1(1), по сути, соответствует графу G1 траектории G1,G2,..., Gl . Тогда для каждой вершины подграф-затравки z1(1) найдено число s(xs(2)), s2 = 1, n .
Далее рассматриваем подграф-затравку z1(1) как отдельный граф и находим для каждой ее вершины xS1) число разделения следующим образом: s(x((1)) = max [d(xSX), v(1)) + s(xS2)}]. Вершина x0, для которой число раз-
s' j1=1,n-1 s j s2
деления s(x0) = min [s(xs(1))] будет минимальным, является центром пред-
s[ =1,2,..., n 1
фрактального графа GL .
Алгоритм a
Вход: Выход:
взвешенный предфрактальный граф GL вершина x0 - центр GL.
Шаг 1.
Шаг l .
Шаг l +1,
Для каждой вершины xSL) подграф-затравки L -го ранга z
(L)
Sl = 1 П
.(L)
L-1
найти
s(x; )) = max [d(х
(L) л,( L)
v.
Jl
)], где d (x
число
(L) V(L)
х
разделения ') = 0 . Поиск крат-
чайших путей между вершинами осуществляется с помощью процедуры Дейкстры. Для всех I = Ь -1, Ь - 2,...,2 выполнить:
Шаг l - l + 1. Для каждой вершины x(s) подграф-затравки l -го ранга z
Л')
найти число разделения s(х()) = max[d(x(J),vj)) + s(xj+1))]
.('+1)-
si
ji
1, П 1, где d (x
(i) x (i)
*) = 0. Поиск кратчайших путей меж-
ду вершинами осуществляется с помощью процедуры Дейкетры.
Для каждой вершины х^ подграф-затравки первого ранга
(1)
найти число разделения
) + s(xj2))], SX = 1,П, где d(xS1), xS1)) = 0.
,v(1)
s1 ' J1
Вход: Выход:
s (x^) = max[ d (x(1)
1 J1=1,п- 1
Поиск кратчайших путей между вершинами осуществляется с помощью процедуры Дейкстры.
Из всех вершин xS^, s1 = 1,2,..., п в качестве центра предф-рактального графа GL выбрать вершину x* с наименьшим числом разделения: s(x*) = min [s(xs(1))].
S[ =1,2,..., n 1
Процедура Дейкстры
взвешенный граф G = (V, E).
кратчайшие пути d(xi, v.) для всех j = 1,2,..., n.
Теорема 3. Алгоритм а0 осуществляет поиск центра предфрак-тального графа ОЬ, смежность старых ребер которого сохраняется.
Следствие 1. Алгоритма а0 выделяет центр хг 0, / = 1,М предф-рактального графа ОЬ, смежность старых ребер которого сохраняется.
Алгоритм осуществляет поиск центра для фиксированного набора
весов. Применение алгоритма по каждому набору весов / = 1,М позволяет находить центры х10.
L
s
s
1
Теорема 4. Вычислительная сложность алгоритма а0 на предф-рактальном графе Оь = (Уь, Еь), порожденном затравкой Н = (Ж, Q), равна О(4п2 • N), где \УЬ\ = N и \Ж\ = п.
Следствие 2. Вычислительная сложность алгоритма а0 поиска центров хг0, 1 = 1,М на предфрактальном графе Оь = (Уь, Еь), порожденном затравкой Н = (Ж,Q),равна О(4п2 • N • М), где \УЬ| = N и \\| = п.
Вычислительная сложность алгоритма а0 поиска центра х, 0 для одного набора весов равна О(4п2 • N). Вычислительная сложность поиска центра по каждому набору весов 1 = 1,М увеличится в М раз и составит О (4п2 • N • М).
Теорема 5. Алгоритм а0 выделяет центр х/0 на предфрактальном графе Оь, оптимальный по критериям ¥2М+1(х, 0) и Е2М+2(хг 0), и оценива-
вь — 1 вь — 1
емый по критериям а——— < ¥г(х, 0) < Ь(п — 1)-, 1 = 1,2,...,2М.
в — 1
в — 1
На рис. 4 представлен пример поиска центра предфрактального графа G3 для случая М = 1. Предфрактальный граф G3 взвешен в соответствии с правилом взвешивания ребер, где начальный отрезок [а; Ь] = [9;18], а весовой коэффициент в = 1 / 3 . Алгоритм поиска центра начинает свою
работу с подграф-затравок 3-го ранга. Внутри подграф-затравок мелким шрифтом указаны числа разделения общих вершин. Центром графа G 3 является вершина хх(1), для которой число разделения минимальное: х* = хх(1).
Радиус предфрактального графа G3 равен числу разделения вершины х*: р = 5 (х *) = 20,5. Значения критериев F1( х *) = F2( х *) = р = 5 (х *) = 20,5. Оценки критериев следующие:
F (X*) = F2(X*) < b(n - П^--1 = 18 • 2 •(1/3) -1 = 52 ; 1 ^ q-1 1/3 -1
FAх*) = F2(х*) > a^-1 = 9 • i^3-1 = 13.
1 ^ q-1 1/3 -1
Оценка критериев верна: 13 < F1(х*) = F2(х*) = 20,5 < 52 .
Исследование выполнено при финансовой поддержке РФФИ в рамках научного проекта № 16 - 07 - 00231а. The reported study was partially supported by RFBR, research project № 16 - 07 - 00231а.
Список литературы
1. Перепелица В. А. Многокритериальные модели и методы для задач оптимизации на графах. - LAP LANDERT Academic Publishing. - 333 с.
2. Кочкаров А.М. Распознавание фрактальных графов. Алгоритмический подход. - Нижний Архыз: РАН САО, 1998. - 170 с.
3. Кочкаров А. А., Малинецкий Г.Г., Кочкаров Р. А. Некоторые аспекты динамической теории графов // Журнал вычислительной математики и математической физики, 2015, том 55, № 9. - С. 1623-1629.
4. Узденов А.А., Кочкаров Р. А. Алгоритм поиска центра предфрактального графа, смежность старых ребер которого сохраняется // Вестник Северо-Кавказского государственного технического университета. - Ставрополь.: СевКавГТУ, 2011. - № 1 (26). -С. 50-53.3.
5. Кочкаров Р. А. Задачи многокритериальной оптимизации на многовзвешенных предфрактальных графах. - М.: Академинновация, 2014. - 189 с.
References
1. Perepelica V.A. Mnogokriterial'nye modeli i metody dlja zadach optimizacii na grafah. - LAP LANDERT Academic Publishing. - 333 s.
2. Kochkarov A.M. Raspoznavanie fraktal'nyh grafov. Algoritmicheskij podhod. -Nizhnij Arhyz: RAN SAO, 1998. - 170 s.
3. Kochkarov A.A., Malineckij G.G., Kochkarov R.A. Nekotorye aspekty dinamich-eskoj teorii grafov // Zhurnal vychislitel'noj matematiki i matematicheskoj fiziki, 2015, tom 55, № 9. - S. 162З-1629.
4. Uzdenov A.A., Kochkarov R.A. Algoritm poiska centra predfraktal'nogo grafa, smezhnost' staryh reber kotorogo sohranjaetsja // Vestnik Severo-Kavkazskogo gosudarstven-nogo tehnicheskogo universiteta. - Stavropol'.: SevKavGTU, 2011. - № 1 (26). - S. 50-5З.З.
5. Kochkarov R.A. Zadachi mnogokriterial'noj optimizacii na mnogovzveshennyh predfraktal'nyh grafah. - M.: Akademinnovacija, 2014. - 189 s.