Типовой расчет графов

Данная работа является типовым расчетом N2 по курсу "Дискретная математика" по теме "Графы", предлагаемая студентам МГТУ им. Баумана. (Вариант N 17).

Сразу хочу сказать для своих коллег: Граждане! Имейте терпение и совесть, поймите, что я это делаю для Вас с целью помочь разобраться в этой теме, а не просто свалить очередной предмет. Мне известно, как непросто сейчас с литературой, и с информацией вообще. Поиски неизвестно какой книги занимают много времени, поэтому в конце я привел небольшой список литературы, составленный мной из различных источников в дополнение к списку, написанному ранее в работе по графам (о постановке лаб. работ по алгоритму Прима и Дейкстра), которая, я надеюсь, есть в сети.

Содержание работы:

Типовой расчет состоит из 11-ти задач:

1, 2 и 3 задачи относятся к способам задания графов и опредению их характеристик, таких как диаметр, радиус и т.д.

4 и 5 задачи соответственно на алгоритм Прима и Дейкстра. Здесь я снова отсылаю Вас к более ранней работе (см. выше).

6-я задача о поиске максим ального потока в сети (метод Форда-Фалкерсона).

7-я задача - Эйлерова цепь (задача о почтальоне).

8-я задача - Гамильтонова цепь.

9-я задача - метод ветвей и границ применительно к задаче о коммивояжере.

10-я задача - задача о назначениях; венгерский алгоритм.

11-я задача - тоже методом ветвей и границ.

Gор(V,X)

Рис. 1

Задача1 Для неориентированного графа G, ассоциированного с графом Gор выписать (перенумеровав вершины) :

а) множество вершин V и множество ребер X, G(V,X);

б) списки смежности;

в) матрицу инцидентности;

г) матрицу весов.

д) Для графа Gор выписать матрицу смежности.

Нумерация вершин - см. Рис 1

а) V={0,1,2,3,4,5,6,7,8,9}

X={{0,1},{0,2},{0,3},{1,2},{1,4},{1,5},{1,6},{1,7},{2,3},{2,5},{3,8},{3,9},{4,5},{4,6},{5,3},{5,6},{5,8},{6,9},{7,8},{7,9},{8,9}}

В дальнейшем ребра будут обозначаться номерами в указанном порядке начиная с нуля.

б) Г0={1,2,3};

Г1={0,2,4,5,6,7};

Г2={0,1,3,5};

Г3={0,2,5,8,9};

Г4={1,5,6};

Г5={1,2,3,4,6,8};

Г6={1,4,5,9};

Г7={1,8,9};

Г8={1,3,5,7,9};

Г9={3,6,7,8};

в) Нумерация вершин и ребер соответственно п. а)

0

1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

17

18

19

20

0

1

1

1

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

1

1

0

0

1

1

1

1

1

0

0

0

0

0

0

0

0

0

0

0

0

0

2

0

1

0

1

0

0

0

0

1

1

0

0

0

0

0

0

0

0

0

0

0

3

0

0

1

0

0

0

0

0

1

0

1

1

0

0

1

0

0

0

0

0

0

4

0

0

0

0

1

0

0

0

0

0

0

0

1

1

0

0

0

0

0

0

0

5

0

0

0

0

0

1

0

0

0

1

0

0

1

0

1

1

1

0

0

0

0

6

0

0

0

0

0

0

1

0

0

0

0

0

0

1

0

1

0

1

0

0

0

7

0

0

0

0

0

0

0

1

0

0

0

0

0

0

0

0

0

0

1

1

0

8

0

0

0

0

0

0

0

0

0

0

1

0

0

0

0

0

1

0

1

0

1

9

0

0

0

0

0

0

0

0

0

0

0

1

0

0

0

0

0

1

0

1

1

г) Показана верхняя половина матрицы, т.к. матрица весов неориентированного графа симметрична относительно главной диагонали.

0

1

2

3

4

5

6

7

8

9

0

¥

8

3

5

¥

¥

¥

¥

¥

¥

1

¥

1

¥

2

2

4

5

¥

¥

2

¥

2

¥

5

¥

¥

¥

¥

3

¥

¥

1

¥

¥

1

6

4

¥

4

2

¥

¥

¥

5

¥

2

¥

1

¥

6

¥

¥

¥

2

7

¥

1

1

8

¥

6

9

¥

д) Матрица смежности для графа Gор.

0

1

2

3

4

5

6

7

8

9

0

¥

1

1

1

¥

¥

¥

¥

¥

¥

1

-1

¥

1

¥

1

1

1

1

¥

¥

2

-1

-1

¥

1

¥

1

¥

¥

¥

¥

3

-1

¥

-1

¥

¥

-1

¥

¥

1

1

4

¥

-1

¥

¥

¥

1

1

¥

¥

¥

5

¥

-1

-1

1

-1

¥

1

¥

1

¥

6

¥

-1

¥

¥

-1

-1

¥

¥

¥

1

7

¥

-1

¥

¥

¥

¥

¥

¥

1

1

8

¥

¥

¥

-1

¥

-1

¥

-1

¥

1

9

¥

¥

¥

-1

¥

¥

-1

-1

-1

¥

Задача 2 Найти диаметр D(G), радиус R(G), количество центров Z(G) для графа G ; указать вершины, являющиеся центрами графа G.

D(G)=2

R(G)=2

Z(G)=10

Все вершины графа G(V,X) являются центрами.

Задача 3 Перенумеровать вершины графа G, используя алгоритмы:

а) "поиска в глубину";

б) "поиска в ширину".

Исходная вершина - a .

а)

б)

Задача 4 Используя алгоритм Прима найти остов минимального веса графа G. выписать код укладки на плоскости найденного дерева, приняв за корневую вершину a .

Вес найденного дерева - 14.

Код укладки дерева: 000011000001111111.

Задача 5 Используя алгоритм Дейкстра найти дерво кратчайших путей из вершины a графа G.

Вес найденного пути - 8.

Задача 6 Используя алгоритм Форда - Фалкерсона, найти максимальный поток во взвешенной двуполюсной ориентированной сети {Gор , a , w}. Указать разрез минимального веса.

Последовательность насыщения сети (насыщенные ребра отмечены кружечками):

1-й шаг

2-й шаг

3-й шаг

4-й шаг

5-й шаг

6-й шаг

7-й шаг

Окончательно имеем:

Как видно из рисунка, ребра {6,9},{7,9},{3,9}, питающие вершину w , насыщенны, а оставшееся ребро {8,9}, питающееся от вершины 8, не может получить большее значение весовой функции, так как насыщенны все ребра, питающие вершину 8. Другими словами - если отбросить все насыщенные ребра, то вершина w недостижима, что является признаком максимального потока в сети.

Максимальный поток в сети равен 12.

Минимальный разрез сети по числу ребер: {{0,1},{0,2},{0,3}}. Его пропускная способность равна 16

Минимальный разрез сети по пропускной способности: {{6,9}, {7,9}, {3,9}, {3,8}, {5,8}, {7,8}}. Его пропускная способность равна 12.

Задача 7 (Задача о почтальоне) Выписать степенную последовательность вершин графа G.

а) Указать в графе G Эйлерову цепь. Если таковой цепи не существует, то в графе G добавить наименьшее число ребер таким образом, чтобы в новом графе можно было указать Эйлерову цепь.

б) Указать в графе G Эйлеров цикл. Если такого цикла не существует, то в графе G добавить наименьшее число ребер таким образом, чтобы в новом графе можно было указать Эйлеров цикл.

Степенная последовательность вершин графа G:

(3,6,4,5,3,6,4,3,4,4)

а) Для существования Эйлеровой цепи допустимо только две вершины с нечетными степенями, поэтому необходимо добавить одно ребро, скажем между вершинами 4 и 7.

Полученная Эйлерова цепь: 0,3,2,0,1,2,5,1,4,5,6,1,7,4,6,9,7,8,9,3,8,5,3.

Схема Эйлеровой цепи (добавленное ребро показано пунктиром):

б) Аналогично пункту а) добавляем ребро {3,0}, замыкая Эйлерову цепь (при этом выполняя условие существования Эйлерова цикла - четность степеней всех вершин). Ребро {3,0} кратное, что не противоречит заданию, но при необходимости можно ввести ребра {0,7} и {4,3} вместо ранее введенных.

Полученный Эйлеров цикл: 0,3,2,0,1,2,5,1,4,5,6,1,7,4,6,9,7,8,9,3,8,5,3,0.

Схема Эйлерова цикла (добавленные ребра показаны пунктиром):

Задача 8

а) Указать в графе Gор Гамильтонов путь. Если такой путь не существует, то в графе Gор изменить ориентацию наименьшего числа ребер таким образом, чтобы в новом графе Гамильтонов путь можно было указать.

б) Указать в графе Gор Гамильтонов цикл. Если такой цикл не существует, то в графе Gор изменить ориентацию наименьшего числа ребер таким образом, чтобы в новом графе Гамильтонов цикл можно было указать.

а) Гамильтонов путь (ребра с измененной ориентацией показаны пунктиром):

б) Гамильтонов цикл (ребра с измененной ориентацией показаны пунктиром):

Задача 9 (Задача о коммивояжере) Дан полный ориентированный симметрический граф с вершинами x1, x2,..xn.Вес дуги xixj задан элементами Vij матрицы весов. Используя алгоритм метода ветвей и границ, найти Гамильтонов контур минимального (максимального) веса. Задачу на максимальное значение Гамильтонова контура свести к задаче на минимальное значение, рассмотрев матрицу с элементами ,где . Выполнить рисунок.

Исходная таблица.

x1

x2

x3

x4

x5

x6

x1

¥

3

7

2

¥

11

x2

8

¥

06

¥

4

3

x3

6

05

¥

7

¥

2

x4

6

¥

13

¥

5

¥

x5

3

3

3

4

¥

5

x6

8

6

¥

2

2

¥

Таблица Е 14

x1

x2

x3

x4

x5

x6

x1

¥

1

5

01

¥

7

2

x2

8

¥

01

¥

4

1

x3

6

00

¥

7

¥

00

x4

1

¥

8

¥

01

¥

5

x5

01

00

00

1

¥

00

3

x6

6

4

¥

00

00

¥

2

2

Дробим по переходу x2-x3:

Таблица 23 å =14+0=14

x1

x2

x4

x5

x6

x1

¥

1

01

¥

7

x3

6

¥

7

¥

06

x4

1

¥

¥

01

¥

x5

01

01

1

¥

00

x6

6

4

00

00

¥

Таблица 23 å =14+1=15

x1

x2

x3

x4

x5

x6

x1

¥

1

5

01

¥

7

x2

7

¥

¥

¥

3

03

1

x3

6

00

¥

7

¥

00

x4

1

¥

8

¥

01

¥

x5

01

00

05

1

¥

00

x6

6

4

¥

00

00

¥

Продолжаем по 23. Дробим по переходу x3-x6:

Таблица 23E36 å =14+0=14

x1

x2

x4

x5

x1

¥

1

01

¥

x4

1

¥

¥

01

x5

01

01

1

¥

x6

6

¥

00

00

Таблица 2336 å =14+6=20

x1

x2

x4

x5

x6

x1

¥

1

01

¥

7

x3

01

¥

1

¥

¥

6

x4

1

¥

¥

01

¥

x5

00

01

1

¥

07

x6

6

4

00

00

¥

Продолжаем по 2336. Дробим по переходу x4-x5:

Таблица 23E3645 å =14+0=14

x1

x2

x4

x1

¥

1

01

x5

01

01

1

x6

6

¥

00

Таблица 233645 å =14+1=15

x1

x2

x4

x5

x1

¥

1

01

¥

x4

00

¥

¥

¥

1

x5

01

01

1

¥

x6

6

¥

00

00

Продолжаем по 233645. Дробим по переходу x5-x1:

Таблица 23364551 å =14+1=15

x2

x4

x1

1

¥

1

x6

¥

00

Таблица 23364551 å =14+6=20

x1

x2

x4

x1

¥

1

01

x5

¥

01

¥

x6

0

¥

00

6

Окончательно имеем Гамильтонов контур: 2,3,6,4,5,1,2.

Прадерево разбиений:

Задача 10 (Задача о назначениях) Дан полный двудольный граф Knn с вершинами первой доли x1, x2,..xn.и вершинами другой доли y1, y2,..yn.Вес ребра {xi,yj} задается элементами vij матрицы весов. Используя венгерский алгоритм, найти совершенное паросочетание минимального (максимального веса). Выполнить рисунок.

Матрица весов двудольного графа K55 :

y1

y2

y3

y4

y5

x1

2

0

0

0

0

x2

0

7

9

8

6

x3

0

1

3

2

2

x4

0

8

7

6

4

x5

0

7

6

8

3

Первый этап - получение нулей не нужен, т. к. нули уже есть во всех строк и столбцах.

Второй этап - нахождение полного паросочетания.

y1

y2

y3

y4

y5

x1

2

0

0

0

0

x2

0

7

9

8

6

x3

0

1

3

2

2

x4

0

8

7

6

4

x5

0

7

6

8

3

Третий этап - нахождение максимального паросочетания.

y1

y2

y3

y4

y5

x1

2

0

0

0

0

X

x2

0

7

9

8

6

X

x3

0

1

3

2

2

x4

0

8

7

6

4

x5

0

7

6

8

3

X

X

Четвертый этап - нахождение минимальной опоры.

y1

y2

y3

y4

y5

x1

2

0

0

0

0

x2

0

7

9

8

6

5

x3

0

1

3

2

2

1

x4

0

8

7

6

4

2

x5

0

7

6

8

3

3

4

Пятый этап - возможная перестановка некоторых нулей.

y1

y2

y3

y4

y5

x1

3

0

0

0

0

x2

0

6

8

7

5

5

x3

0

0

2

1

1

1

x4

0

7

6

5

3

2

x5

0

6

5

7

2

3

4

Решение с ненулевым значением. Переход ко второму этапу.

Полное паросочетание:

y1

y2

y3

y4

y5

x1

3

0

0

0

0

x2

0

6

8

7

5

x3

0

0

2

1

1

x4

0

7

6

5

3

x5

0

6

5

7

2

Максимальное паросочетание:

y1

y2

y3

y4

y5

x1

3

0

0

0

0

X

x2

0

6

8

7

5

X

x3

0

0

2

1

1

x4

0

7

6

5

3

x5

0

6

5

7

2

X

X

Минимальная опора:

y1

y2

y3

y4

y5

x1

3

0

0

0

0

6

x2

0

6

8

7

5

7

x3

0

0

2

1

1

1

x4

0

7

6

5

3

2

x5

0

6

5

7

2

3

4

5

Перестановка нулей:

y1

y2

y3

y4

y5

x1

3

0

0

0

0

6

x2

0

6

8

7

5

7

x3

0

0

2

1

1

1

x4

0

7

6

5

3

2

x5

0

6

5

7

2

3

4

5

Полное паросочетание:

y1

y2

y3

y4

y5

x1

3

0

0

0

0

6

x2

0

6

8

7

5

7

x3

0

0

2

1

1

1

x4

0

7

6

5

3

2

x5

0

6

5

7

2

3

4

5

Максимальное паросочетание:

y1

y2

y3

y4

y5

x1

3

0

0

0

0

X

x2

0

6

8

7

5

x3

0

0

2

1

1

X

x4

0

7

6

5

3

X

x5

0

6

5

7

2

X

X

X

Минимальная опора:

y1

y2

y3

y4

y5

x1

3

0

0

0

0

x2

0

6

8

7

5

1

x3

0

0

2

1

1

x4

0

7

6

5

3

x5

0

6

5

7

2

2

3

Перестановка нулей:

y1

y2

y3

y4

y5

x1

5

0

0

0

0

x2

0

4

6

5

3

1

x3

2

0

2

1

1

x4

2

7

6

5

3

x5

0

4

3

5

0

2

3

Полное паросочетание:

y1

y2

y3

y4

y5

x1

5

0

0

0

0

x2

0

4

6

5

3

x3

2

0

2

1

1

x4

2

7

6

5

3

x5

0

4

3

5

0

Максимальное паросочетание:

y1

y2

y3

y4

y5

x1

5

0

0

0

0

X

x2

0

4

6

5

3

X

x3

2

0

2

1

1

X

x4

2

7

6

5

3

x5

0

4

3

5

0

X

X

X

X

X

Минимальная опора:

y1

y2

y3

y4

y5

x1

5

0

0

0

0

x2

0

4

6

5

3

x3

2

0

2

1

1

x4

2

7

6

5

3

1

x5

0

4

3

5

0

Перестановка нулей:

y1

y2

y3

y4

y5

x1

5

0

0

0

0

x2

0

4

6

5

3

x3

2

0

2

1

1

x4

0

5

4

3

1

1

x5

0

4

3

5

0

Полное паросочетание:

y1

y2

y3

y4

y5

x1

5

0

0

0

0

x2

0

4

6

5

3

x3

2

0

2

1

1

x4

0

5

4

3

1

1

x5

0

4

3

5

0

Максимальное паросочетание:

y1

y2

y3

y4

y5

x1

5

0

0

0

0

X

x2

0

4

6

5

3

X

x3

2

0

2

1

1

X

x4

0

5

4

3

1

x5

0

4

3

5

0

X

X

X

X

X

Минимальная опора:

y1

y2

y3

y4

y5

x1

5

0

0

0

0

x2

0

4

6

5

3

3

x3

2

0

2

1

1

x4

0

5

4

3

1

1

x5

0

4

3

5

0

2

Перестановка нулей:

y1

y2

y3

y4

y5

x1

6

0

0

0

0

x2

0

3

5

4

2

3

x3

3

0

2

1

1

x4

0

4

3

2

0

1

x5

1

4

3

5

0

2

Полное паросочетание:

y1

y2

y3

y4

y5

x1

6

0

0

0

0

x2

0

3

5

4

2

3

x3

3

0

2

1

1

x4

0

4

3

2

0

1

x5

1

4

3

5

0

2

Максимальное паросочетание:

y1

y2

y3

y4

y5

x1

6

0

0

0

0

X

x2

0

3

5

4

2

X

x3

3

0

2

1

1

X

x4

0

4

3

2

0

x5

1

4

3

5

0

X

X

X

X

X

Минимальная опора:

y1

y2

y3

y4

y5

x1

6

0

0

0

0

x2

0

3

5

4

2

4

x3

3

0

2

1

1

x4

0

4

3

2

0

1

x5

1

4

3

5

0

5

2

3

В результате имеем:

y1

y2

y3

y4

y5

x1

6

0

0

0

0

x2

0

1

3

2

2

4

x3

3

0

2

1

1

x4

0

2

1

0

0

1

x5

1

4

3

5

0

5

2

3

Исходный граф

Полученный граф:

Вес найденного совершенного паросочетания = 12.

Задача 11 Решить задачу 10, используя алгоритм ветвей и границ (отождествив вершины xi и yj).

Таблица Е (исходная). Строки - xi , столбцы - yj. å =0

1

2

3

4

5

1

2

01

03

02

02

2

06

7

9

8

6

3

01

1

3

2

2

4

04

8

7

6

4

5

03

7

6

8

3

Дробим по переходу x2 - y1:

Таблица Е21 å =0+8=8

2

3

4

5

1

00

02

01

00

3

01

2

1

1

1

4

4

3

2

02

4

5

4

3

5

03

3

Таблица 21 å =0+6=6

1

2

3

4

5

1

2

01

03

02

00

2

¥

1

3

2

01

6

3

01

1

3

2

2

4

04

8

7

6

4

5

03

7

6

8

3

Продолжаем по 21:

Дробим по переходу x4 - y1:

Таблица 21Е41 å =6+4=10

2

3

4

5

1

00

02

01

00

2

1

3

2

01

3

01

2

1

1

1

5

4

3

5

03

3

Таблица 2141 å =6+4=10

1

2

3

4

5

1

2

01

03

02

00

2

¥

1

3

2

01

3

01

1

3

2

2

4

¥

4

3

2

02

4

5

03

7

6

8

3

Продолжаем по Е21:

Дробим по переходу x5 - y5:

Таблица Е21Е55 å =8+2=10

2

3

4

1

00

01

00

3

01

2

1

4

2

1

01

2

Таблица Е2155 å =8+3=11

2

3

4

5

1

00

02

01

00

3

01

2

1

1

4

4

3

2

02

5

1

01

2

¥

3

Продолжаем по Е21Е55:

Дробим по переходу x3 - y2:

Таблица Е21Е55Е32 å =10+0=10

3

4

1

01

00

4

1

01

Далее решение очевидно: x1 - y3 и x4 - y4. Это не увеличит оценку.

В итоге имеем совершенное паросочетание с минимальным весом:

Прадерево разбиений:

Литература

1. Грешилов А.А. Как принять наилучшее решение в реальных условиях:-М.:Радио и связь, 1991.-320с.:ил.

2. Беллман Р. Динамическое программирование: Пер. с англ./Под ред. Н.Н. Воробьева.-М.: ИЛ, 1960.-400 с.

3. Беллман Р., Дрейфус С. Прикладные задачи динамического программирования: Пер с англ./Под ред. А.А. Первозванского.-М.: Наука, 1965.-458 с.

4. Вентцель Е.С. Исследование операций.-М.: Сов. радио, 1972.-551 с.

5. Вильямс Н.Н. Параметрическое программирование в экономике (методы оптимальных решений):-М.:Статистика, 1976.-96с.

6. Гольштейн Е.Г., Юдин Д.Б. Новые направления в линейном программировании:-М.: Сов радио, 1966.- 524 с.

7. Зангвилл У.И. Нелинейное программирование: Пер. с англ./Под ред. Е.Г. Гольштейна.-М.: Сов радио, 1973.- 312 с.

8. Зуховицкий С.И., Авдеева Л.И. Линейное и выпуклое программирование (справочное руководство).-М.: Наука, 1964.-348 с.

9. Исследование операций. Методологические основы и математические методы: Пер. с англ./ Под ред. И.М. Макарова, И.М. Бескровного.-М.: Мир, 1981.- Т.1.-712 с.

10. Исследование операций. Модели и применение: Пер. с англ./ Под ред. И.М. Макарова, И.М. Бескровного.-М.: Мир, 1981.- Т.1.-712 с.

11. Лазарев В.Г., Лазарев Ю.В. Динамическое управление потоками информации в сетях связи.-М.: Радио и связь, 1983.- 216 с.

12. Мартин Дж. Системный анализ передачи данных.: Пер с англ./ Под ред. В.С. Лапина.-М.: Мир, 1975.- М.2.- 431 с.

13. Монаков В.М., Беляева Э.С., Краснер Н.Я. Методы оптимизации. Пособие для учителя.-М.: Просвещение, 1978.- 175с.

14. Муртаф Б. Современное линейное программирование: Теория и практика. Пер. с англ./Под ред. И.А. Станевичуса.- М.: Мир, 1984.- 224 с.

15. Рокафеллор Р. Выпуклый анализ: Пер. с англ./Под ред. А.Д. Иоффе, В.М. Тихомирова.-М.: Мир, 1973.- 469 с.

16. Сухарев А.Г., Тимохов А.В., Федоров В.В. Курс методов оптимизации.- М.:- Наука, Физматгиз, 1986.- 326 с.

17. Ху Т. Целочисленное программирование и потоки в сетях: Пер. с англ./Под ред. А.А. Фридмана.- М.: Мир, 1974.-419 с.

18. Фиакко А., Мак-Кормик Г. Нелинейное программирование. Методы последовательной безусловной минимизации: Пер. с англ./Под ред. Е.Г. Гольштейна. -М.:- Мир, 1972.- 240 с.

19. Филлипс Д., Гарсиа-Диас А. Методы анализа сетей: Пер. с англ./ Под ред. Б.Г. Сушкова.- М.: Мир, 1984.- 496 с.

20. Юдин Д.Б., Гольштейн Е.Г. Линейное программирование. Теория и конечные методы,- М.:- Физматгиз, 1963.- 775 с.


Вместе с этим смотрят:


"Инкарнация" кватернионов


*-Алгебры и их применение


10 способов решения квадратных уравнений


РЖнварiантнi пiдпростори. Власнi вектори i власнi значення лiнiйного оператора


РЖнженерна графiка