CC BY
63
Научный журнал КубГАУ, №115(01), 2016 года
1
УДК 519.157.1
01.00.00 Физико-математические науки
ИНТЕРВАЛЬНЫЕ ЗАДАЧИ ОБ ОСТОВНЫХ ДЕРЕВЬЯХ С ТОПОЛОГИЧЕСКИМ КРИТЕРИЕМ
Шапошникова Ольга Ивановна к.ф.-м.н., доцент
Северо-Кавказская государственная гуманитарнотехнологическая академия, Черкесск, Россия [email protected]
Темирова Лилия Гумаровна к.ф.-м.н., доцент SPIN-код: 5002-3657
Северо-Кавказская государственная гуманитарнотехнологическая академия, Черкесск, Россия blg1961 @rambler. ru
В статье представлена задача об остовных деревьях с топологическими критериями и интервальными весами. Введены отношения предпочтения и несравнимости для нахождения полного множества альтернатив в случае интервальных весов. Базу математических расчетов составляет интервальная математика
Ключевые слова: ГРАФ, ОСТОВНОЕ ДЕРЕВО, ИНТЕРВАЛ, ПАРЕТОВСКИЙ ОПТИМУМ, ПОЛНОЕ МНОЖЕСТВО АЛЬТЕРНАТИВ, ЦЕЛЕВАЯ ФУНКЦИЯ
UDC 519.157.1 Physics and Math
INTERVAL SPANNING TREE PROBLEM ON A TOPOLOGICAL CRITERION
Shaposhnikova Olga Ivanovna Cand.Phys.-Math.Sci., associate professor North-Caucasian State Humanities and Technology Academy,Cherkessk, Russia [email protected]
Temirova Lilija Gumarovna Cand.Phys.-Math.Sci., associate professor RSCI SPIN-code: 5002-3657
North-Caucasian State Humanities and Technology Academy,Cherkessk, Russia blg1961 @rambler.ru
The article presents the problem of spanning trees with topological criteria and interval scales. We have introduced relationship preferences and incomparability to find the complete set of alternatives in the case of interval scales. The base for mathematical calculations is interval mathematics
Keywords: GRAPH, SPANNING TREE, INTERVAL, PARETO OPTIMUM, FULL SET OF ALTERNATIVES, OBJECTIVE FUNCTION
В процессе практического использования каких-либо
математических моделей необходимо учитывать неоспоримый факт того, что все или большинство числовых параметров этих моделей определены приближенно. Как правило, достоверным является тот факт, что значение рассматриваемого параметра принадлежит определенному интервалу или отрезку числовой оси [1].
Фундаментальная проблема определения наилучшего решения интервальной задачи возникает как для многокритериальных постановок, так и в случае, когда рассматриваемая задача является оптимизационной, когда качество решений оценивается целевой функцией с единственным критерием.
http://ej.kubagro.ru/2016/01/pdf/23.pdf
Научный журнал КубГАУ, №115(01), 2016 года
2
Вопрос нахождения наилучшего решения интервальной задачи сводится к следующим этапам. Вначале определяется состав критериев целевой функции, после чего находится полное множество альтернатив (ПМА) X0[2]. На заключительном этапе с помощью процедур теории выбора и принятия решений из ПМА X0 выбирается искомый «компромиссный» оптимум х0.
Рассмотрим общую постановку многокритериальной задачи об остовных деревьях [2]. Дан n -вершинный N -взвешенный граф G = (V, E), в котором каждому ребру eе E приписаны веса wv(e)> 0, v = 1,2,...,N. Допустимым решением является остовное дерево х = (V, Ех) графа G; X = X (G )={х} - множество допустимых решений (МДР). На МДР X определена векторная целевая функция (ВЦФ):
f (х )=(f1 (х \ f2 (х х..^ fn (х)), (1)
состоящая в общем случае из критериев вида MINSUM
Fv (х)= 2 Wv (e)® min, v = 1,2,..., N,, N < N (2)
eeE V /
и критериев вида MINMAX
Fv (х ) = max wv (e)® min, v = N +1,..., N
ee Ех
(3)
При N > 2 ВЦФ (1)-(3) определяет на МДР X паретовское множество (ПМ) X с X [2], состоящее из паретовских оптимумов (ПО). Элемент ~ е X называется ПО, если X не содержит такого элемента х*, для которого выполняются неравенства F (х* )< F (~), v = 1, N, среди которых хотя бы одно является строгим. Искомым решением данной векторной задачи является так называемое полное множество альтернатив (ПМА) X0. Подмножество X0 с X называется ПМА, если его мощность X0
минимальна и при этом выполняется равенство F (X0) = F (X), где F(X*) = {F(х): хе X*} "X* с X .
http://ej.kubagro.ru/2016/01/pdf/23.pdf
Научный журнал КубГАУ, №115(01), 2016 года
3
Существует, однако, немало реальных ситуаций, в которых многообразие критериев не исчерпывается весовыми критериями вида (2) или (3). Сравнивать получаемые решения необходимо с учетом также топологических критериев. Среди них наиболее известными являются: радиус р(x), диаметр d(x) остовного дерева; число висячих вершин, т.е.
мощность |V1(x)|, где V1(x) = {v e V;degx v = 1} - множество вершин степени 1; степень дерева W(x) = max degx v и др.
veV
В настоящей статье рассматривается такой класс 2-критериальных задач об остовных деревьях с топологическим критерием, ВЦФ которых содержит первый критерий F1 (x) ® min вида MINSUM (2) или MINMAX
(3), а второй критерий является топологическим. Таким образом получаем двукритериальную задачу с ВЦФ
F (x )=(F1 (x), F2 (x)) (4)
у которой F1 (x) представляет собой минимизируемый метрический критерий
F (x)e \ 2 К (e), max wv (e)|
I eeEt
(5)
а второй критерий является топологическим
F2 (x) e {d(x) ® max, W(x) ® min, | V1 (x)| ® min}. (6)
Известно, что проблема нахождения ПМА для векторной задачи об остовных деревьях является труднорешаемой, т.е. имеет
экспоненциальную вычислительную сложность в случае, если ВЦФ (1) включает в себя N1 > 2 критериев вида MINSUM (2). Эта проблема является полиномиально разрешимой при N1 £ 1 и N = 2 [2]. Проблема нахождения ПМА для каждой из 2-критериальных задач с ВЦФ (4)-(6) является NP-трудной проблемой [3].
http://ej.kubagro.ru/2016/01/pdf/23.pdf
Научный журнал КубГАУ, №115(01), 2016 года
4
Рассмотрим интервальную задачу об остовных деревьях с топологическим критерием.
Дан n -вершинный граф G = (V, E), в котором каждому ребру e е E приписан интервальный вес w(e) = [w1 (e), w2 (e)]. Допустимым решением является остовное дерево x = (V, Ex), Ex c E. Через X = {x} обозначим множество всех остовных деревьев, образующих множество допустимых решений (МДР). На X определена интервальная целевая функция (ИЦФ)
w(x )= 2 w(e)® min, (7)
где под знаком суммы (7) осуществляется сложение интервалов w(e) = [w1(e), w2(e)], согласно аксиомам интервального исчисления [1,4], в частности, если w(e') = [ w1(e/), w2(e')], w(e') = [w1(eA’), w2(e')], то
w(e') + w(e") = [ w1 (e', e"), w2 (e', e^’)], где w1 (e', e^’) = w1 (e') + w1 (e*) и
w2(e, e№) = w2(e) + w2(e№).
В качестве топологического критерия возьмем степень остовного дерева xе X
F2 (x) = W(x) ® min .
(8)
Зафиксируем значения критерия (8) в виде равенства
W(x) = к,2 < к < n -1, (9)
и решаем оптимизационную 1-критериальную задачу с целевой функцией (7). В качестве конкретного значения возьмем к = 3 в условии (9). Согласно определению операции сложения интервалов [4] получим значение ИЦФ w( x) = [w1 (x), w2 (x)], где значения ее критериев
w, (x) = 2 w, (e) ® min, i = 1,2.
eeEx
(10)
Требуется найти такой элемент x0 е X, на котором значение ИЦФ (7) достигает требуемого экстремума - минимума.
http://ej.kubagro.ru/2016/01/pdf/23.pdf
Научный журнал КубГАУ, №115(01), 2016 года
5
В случае задания интервальных весов решение задачи, вообще говоря, не единственно, и для определения множества решений необходимо ввести, следуя терминологии [1, 4] отношения предпочтения и несравнимости.
Определение 1. Пусть x1, x2 е X. Тогда решение x1
предпочтительнее решения x2 (x1 p x2), если wt (x1) < wt (x2), i = 1,2, при этом хотя бы одно неравенство строгое.
В случае, когда w (x1) с w (x2) либо w (x2) с w (x1), решения
x1, x2 считаем несравнимыми, x1 » x2.
Таким образом, сравнимы лишь интервалы, сдвинутые один относительно другого вдоль числовой оси, причем сдвинутый вправо (влево) интервал, является большим (меньшим) [5].
Введенные на МДР X отношения предпочтения и несравнимости порождают Паретовское множество X с X [2]. Паретовское множество состоит из паретовских оптимумов.
Определение 2. Решение ~ е X называется паретооптимальным для задачи (7), если не существует x е X, такого, что x p ~ .
Следовательно, паретовское множество X состоит из несравнимых решений.
В качестве решения поставленной задачи (7)-(9) можно рассматривать как паретовское множество X, так и используемое в многокритериальной оптимизации полное множество альтернатив X0 [2]. ПМА X0 определяется как подмножество X0 с X минимальной мощности, содержащее по одному представителю на каждое значение ИЦФ (7).
Определение 3. Интервальная задача с ИЦФ (7)-(9) называется полной (обладает свойством полноты), если для ее МДР X существуют такие веса w(e), e е E, что выполняются равенстваX = X = X0.
http://ej.kubagro.ru/2016/01/pdf/23.pdf
Научный журнал КубГАУ, №115(01), 2016 года
6
Как правило, вычисление мощности МДР представляет меньшую трудность по сравнению с вычислением мощности ПМА. Если рассматриваемая задача обладает свойством полноты, то с учетом равенств X = X = X0 снижается сложность нахождения максимальной мощности ПМА.
Согласно [6] всякая интервальная задача на графах с ИЦФ (7) является полной в случае, если ее множество типовых графов является однородным, т.е. когда мощности множеств вершин и ребер одинаковы для всех типовых графов. Под множеством типовых графов понимаем множество допустимых решений. В нашей задаче допустимым решением на графе G = (V, E) является остовное дерево x = (V, Ex), Ex c E, и для
каждого выполняется условие \V\ = n, |Ex| = n -1. Следовательно, МДР X = {x} рассматриваемой задачи об остовных деревьях с топологическим критерием однородно.
Мощность МДР X = {x} задачи об остовных деревьях с топологическим критерием растет экспоненциально с ростом n в случае, когда G - полный граф. При этом максимум мощности МДР X = {x}
определяется формулой |X| = nn-2. С учетом этих замечаний, интервальная
задача об остовных деревьях с топологическим критерием является труднорешаемой.
Рассмотрим конкретную индивидуальную интервальную задачу об остовных деревьях степени 3 на 6-вершинном графе G = (V, E), представленном на рисунке 1.
http://ej.kubagro.ru/2016/01/pdf/23.pdf
Научный журнал КубГАУ, №115(01), 2016 года
7
Рисунок 1. 6-вершинный граф
Каждому ребру e е E графа G = (V, E) приписаны интервальные веса, таким образом, как они представлены в таблице 1.
Таблица 1. Интервальные веса ребер графа G = (V, E)
ei, i = 1,9 e1 e2 e3 e4 e 5 e6 e7 e8 e9
w(e) [4,10] [0,5] [1,8] [6,6] [3,3] [3,4] [2,7] [3,6] [3,6]
МДР X рассматриваемой индивидуальной задачи состоит из шести остовных деревьев xk = (V, E4), к = 1,6, степени 3, где множества ребер
Exj = {e1, e2, e3, e7, e8}, Ex2 = {e2, e1, e6, e9, e7 }, Ex3 = {e2, e3 , e4, e8 , e9 }, Ex4 = {e3 , e4, e5, e7, e9}, Ex5 = {e5 , e6, e1, e9, e8 }, Ex6 = {e4, e5 , e6 , e7 , e8}.
Значения ИЦФ каждого решения xk =(v, E4), к = 1,6 согласно аксиомам интервального исчисления представлены в таблице 2.
Таблица 2. Значения интервальной целевой функции для МДР X
xk к = 1,6 xt x2 x3 x4 x5 x6
w( x) [10,36] [12,32] [13,31] [15,30] [16,29] [17,26]
http://ej.kubagro.ru/2016/01/pdf/23.pdf
Научный журнал КубГАУ, №115(01), 2016 года
8
Согласно введенным выше отношениям предпочтения и несравнимости все решения в данной индивидуальной задаче несравнимы, т.е. X = X = X0, и рассматриваемая задача обладает свойством полноты.
Одним из подходов к определению решения оптимизационной задачи на графах с интервальными параметрами является понятие реализации весов ребер графа [7, 8], что означает выбор вещественных значений весов ребер, находящихся внутри заданных интервалов.
В статье представлены три формулы для реализации весов ребер графа G = (V, E).
Первая реализация, в случае , когда веса ребер графа могут отвечать
Гw1(e), если e е Ex
следующим требованиям: w(e) = < X .
[w2 (e), если e € Ex
Применив эту реализацию в нашем примере, мы получим МДР X со следующими значениями (см. табл. 3), среди которых, определяем, согласно условию экстремума (7), оптимальное решение, равное x1 , со значением ИЦФ w( x1) = 10.
Таблица 3. Численные значения множества допустимых решений при
первой реализации
Xkk = 1,6 X1 X2 X3 X4 X5 X6
w( x) 10 12 13 15 16 17
Вторая реализация весов ребер графа могут быть веса
w(e) = 1w1 (e) + (1 - 1)w2 (e) при 1 е (0,1). Пусть 1 = 0,5. При этой
реализации рассматриваемая задача имеет значения ИЦФ, представленные таблицей 4.
http://ej.kubagro.ru/2016/01/pdf/23.pdf
Научный журнал КубГАУ, №115(01), 2016 года
9
Таблица 4. Численные значения множества допустимых решений
при второй реализации
xkk = 1,6 xl x2 x3 x4 x5 x6
w(x) 23 22 22 22,5 22,5 21,5
Среди которых, оптимальным решением является x6, со значением ИЦФ w( Хб) = 21,5.
Третий случай реализации весов ребер графа могут быть
\w2(e), если eе Ex
представлены весами (см. табл.5): w(e) = < Х .
[w1 (e), если e <£ Ex
Таблица 5. Численные значения множества допустимых решений при третьей реализации
xk k = 1,6 x1 x2 x3 x4 x5 x6
w(x) 36 32 31 30 29 26
Оптимальным решением при этой реализации является x6, со значением ИЦФ w( x6) = 26.
Таким образом, проведенные расчеты показывают, что решение xk будет оптимальным, если оно является оптимальным по большинству числа реализаций.
В статье рассмотрена интервальная задача об остовных деревьях с топологическим критерием. Исследована вычислительная сложность задачи, что привело к ее труднорешаемости. Введены отношения предпочтения и несравнимости для нахождения полного множества альтернатив в случае интервальных весов. Предложен ряд конкретных реализаций весов графа для нахождения оптимального решения.
http://ej.kubagro.ru/2016/01/pdf/23.pdf
Научный журнал КубГАУ, №115(01), 2016 года
10
Литература
1. Вощинин А.П., Сотиров Г.Р. Оптимизация в условиях неопределенности. Изд-во МЭИ, Москва, изд-во Техника, София, 1989.
2. Емеличев В.А., Перепелица В.А. Сложность дискретных
многокритериальных задач // Дискретная математика.- 1994.-Т.6, вып.1.- С.3 - 33.
3. Шапошникова О.И. Алгоритмы с оценками для векторной задачи об остовных деревьях. Деп. рукопись в ВИНИТИ, рег. № 3618-1398 от 09.12.98.
4. Алефельд Г., Херцбергер Ю. Введение в интервальные вычисления. - М.: Мир, 1987.- 360с.
5. Левин В.И. Сравнение интервальных величин и оптимизация
неопределенных систем // Информационные технологии. 1998. №7. С.22-32.
6. Перепелица В.А., Тебуева Ф.Б. Задачи дискретной оптимизации с интервальными параметрами // Журнал вычислительной математики и математической физики. 2010. Т.50. №5. С.836-847.
7. Yaman H, Karasan O.E, Pinar M.C. Minimum SpanningTree Problem with Interval Data. - Technical Report 9909, Department of Industrial Engineering, Bilkent University, Ankara, Turkey, 1999.
8. Козина Г.Л. Правило выбора решений оптимизационных задач на графах с интервальными параметрами. -Труды XIII Байкальской международной школы-семинара «Методы оптимизации и их приложения», Иркутск: ИСЭМ СО РАН, 2005, с. 51-54.
References
1. Voshhinin A.P., Sotirov G.R. Optimizacija v uslovijah neopredelennosti. Izd-vo MJeI, Moskva, izd-vo Tehnika, Sofija, 1989.
2. Emelichev V.A., Perepelica V.A. Slozhnost' diskretnyh mnogokriterial'nyh zadach // Diskretnaja matematika.- 1994.-T.6, vyp.1.- S.3 - 33.
3. Shaposhnikova O.I. Algoritmy s ocenkami dlja vektornoj zadachi ob ostovnyh derev'jah. Dep. rukopis' v VINITI, reg. № 3618-1398 ot 09.12.98.
4. Alefel'd G., Hercberger Ju. Vvedenie v interval'nye vychislenija. - M.: Mir, 1987.360s.
5. Levin V.I. Sravnenie interval'nyh velichin i optimizacija neopredelennyh sistem // Informacionnye tehnologii. 1998. №7. S.22-32.
6. Perepelica V.A., Tebueva F.B. Zadachi diskretnoj optimizacii s interval'nymi parametrami // Zhurnal vychislitel'noj matematiki i matematicheskoj fiziki. 2010. T.50. №5. S.836-847.
7. Yaman H, Karasan O.E, Pinar M.C. Minimum SpanningTree Problem with Interval Data. - Technical Report 9909, Department of Industrial Engineering, Bilkent University, Ankara, Turkey, 1999.
8. Kozina G.L. Pravilo vybora reshenij optimizacionnyh zadach na grafah s interval'nymi parametrami. -Trudy XIII Bajkal'skoj mezhdunarodnoj shkoly-seminara «Metody optimizacii i ih prilozhenija», Irkutsk: ISJeM SO RAN, 2005, s.
http://ej.kubagro.ru/2016/01/pdf/23.pdf
