Алгоритм Форда-Фалкерсона

Размещено на http://www.allbest.ru/

2

Размещено на http://www.allbest.ru/

Курсовой проект

на тему

Алгоритм Форда-Фалкерсона

Реферат

алгоритм форд фалкерсон граф

Пояснительная записка к курсовой работе: с., рис., табл., разделов, приложения.

АЛГОРИТМ, ГРАФ, РЕБРО, ВЕРШИНА, МАКСИМАЛЬНЫЙ ПОТОК, ПАРАЛЛЕЛЬНЫЙ, ПСЕВДОКОД, БЛОКСХЕМА, РЕАЛИЗАЦИЯ

Объект исследования: алгоритм Форда-Фалкерсона

Цель работы: закрепление знаний и умений, полученных в процессе теоретического обучения.

Метод работы: закрепление теоретических знаний с помощью разработки программы, реализующий данный алгоритм.

Введение

Объектом исследования курсовой работы стала реализация алгоритма Форда-Фалкерсона.

Целями работы являлись:

1) ознакомление с алгоритмом Форда-Фалкерсона, его историей;

2) реализация алгоритма, для построения нахождения максимального потока сети;

3) анализ трудоёмкости алгоритма;

4) тестирование алгоритма.

В теории графов транспортная сеть -- ориентированный граф, в котором каждое ребро имеет неотрицательную пропускную способность и поток. Выделяются две вершины: источник и сток такие, что любая другая вершина сети лежит на пути из источника в сток. Транспортная сеть может быть использована для моделирования, например, дорожного трафика. Целочисленная транспортная сеть -- транспортная сеть, все пропускные способности ребер которой -- целые числа.

Поток -- функция со следующими свойствами для любых вершин:

1) Ограничение пропускной способности;

2) Поток не может превысить пропускную способность;

3) Поток из вершины в другую должен быть противоположным в другом направлении;

4) Сохранение потока. для всех , кроме источника и стока.

Первый отчет «Максимальный поток в сети» датируется 19 ноября 1954 года. Авторы отчета - Форд и Фалкерсон, доказали теорему о максимальном потоке и минимальном разрезе для неориентированных графов: значение максимального потока в сети равно минимальной пропускной способности разреза. (Разрезом в сети называется разбиение множества ее вершин на два непересекающихся класса, таких что источник и сток лежат в разных классах. Пропускной способностью разреза называется сумма пропускных способностей ребер, концы которых лежат в разных классах.) В 1955 году в работе Робахера было упомянуто, что впервые эту теорему доказал Фалкерсон.

Граф можно задавать несколькими путями, однако самый просто способ задания - задание графа матрицей смежности. Тем не менее, этот способ подходит исключительно для таких сетей, которые имеют начальную пропускную способностью равную длине ребра.

Алгоритм широко применяется в сфере перевозок, для того чтобы найти оптимальные маршруты перевозки грузов. Также широкое применение он находит в компьютерных технологиях, коммуникационных сетях, электрических сетях.

Суть распараллеливания алгоритма заключается в уменьшении времени работы с большими объемами входных данных. Стоит заметить, что сложность и трудоемкость алгоритма напрямую зависит от способа реализации. Пример описанных в проекте, работает при помощи алгоритма обхода графа в ширину. Этот алгоритм нужен для того, чтобы искать пути в заданном графе из источника в сток.

1. Графы

1.1 Основные характеристики графов

Граф G - это математический объект, состоящий из множества вершин X = {x₁,x₂,...,x_n} и множества ребер A={a₁, a₂,..., a_m}. Таким образом, граф полностью определяется совокупностью множеств X, A и обозначается G(X, A).

Если ребрам графа приданы направления от одной вершины к другой, то такой граф называется ориентированным. Ребра ориентированного графа называются дугами. Соответствующие вершины ориентированного графа называют началом и концом. Если направления ребер не указываются, то граф называется неориентированным (или просто графом).Граф, имеющий как ориентированные, так и неориентированные ребра, называется смешанным.Неориентированное ребро графа эквивалентно двум противоположно направленным дугам, соединяющим те же самые вершины.

Множество ребер графа может быть пустым. Множество вершин графа не может быть пустым. Граф G (X, A)- полный, если для любой пары вершин x_iи x_j существует ребро (x_i, x_j).

1.2 Пути и маршруты

Маршрутом или цепью в неориентированном графе G называется такая последовательность (конечная или бесконечная) ребер a₁, a₂,...,a_n,…, что каждые соседние два ребра a_i и a_i+1имеют общую инцидентную вершину. Одно и то же ребро может встречаться в маршруте несколько раз. В конечном маршруте (a₁,a₂,...,a_n) имеется первое ребро a₁и последнее ребро a_n. Вершина x₁, инцидентная ребру a₁, но не инцидентная ребру a₂, называется началом маршрута, а вершина x_n, инцидентная ребру a_n, но не инцидентная ребру a_n-1, называется концом маршрута.

Длиной (или мощностью) маршрута называется число ребер, входящих в маршрут, причем каждое ребро считается столько раз, сколько оно входит в данный маршрут.

Замкнутый маршрут называется циклом. Маршрут (цикл), в котором все ребра различны, называется простым маршрутом или простой цепью (циклом). Маршрут (цикл), в котором все вершины, (кроме первой и последней), различны, называется элементарным маршрутом или элементарной цепью (циклом).

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

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

Путь (контур), в котором все дуги различны, называется простым. Путь (контур), в котором все вершины, кроме первой и последней, различны, называется элементарным.

Иногда дугам графа сопоставляются числа называемые весом, или длиной дуги. Тогда граф называется графом со взвешенными дугами. Иногда веса приписываются вершинам графа, и тогда получается граф с взвешенными вершинами. Если в графе веса приписаны и дугам, и вершинам, то он называется просто взвешенным. Вес пути в графе определяется как сумма весов ребер этого пути. Кратчайшим путем между двумя вершинами называется путь наименьшего веса, соединяющий эти вершины.

2. Алгоритм Форда-Фалкерсона

2.1 Краткое описание алгоритма

Идея алгоритма заключается в следующем. Выбираем такой путь от источника с к стоку, чтобы для каждого ребра остаточная пропускная способность была строго больше нуля. При этом ребра на данном пути могут проходиться как в прямом, так и в обратном направлении. Выбираем минимальное значение среди остаточных пропускных способностей ребер данного пути. Увеличиваем поток на каждом из ребер данного пути на выбранное минимальное значение. Далее ищем следующий аналогичный путь.

Работа алгоритма продолжается до тех пор, пока удается находить данные пути.

Сразу отметим, что данный алгоритм относится к классу недетерминированных, т.е. каждый следующий шаг алгоритма определен неоднозначно. И время работы (количество шагов) алгоритма зависит от того, как будут выбираться шаги.

Неформальное описание алгоритма:

ѕ Обнуляем все потоки. Остаточная сеть изначально совпадает с исходной сетью;

ѕ В остаточной сети находим любой путь из источника в сток. Если такого пути нет, останавливаемся;

ѕ Пускаем через найденный путь (он называется увеличивающим путём или увеличивающей цепью) максимально возможный поток;

ѕ На найденном пути в остаточной сети ищем ребро с минимальной пропускной способностью ;

ѕ Для каждого ребра на найденном пути увеличиваем поток на , а в противоположном ему -- уменьшаем на ;

ѕ Модифицируем остаточную сеть. Для всех рёбер на найденном пути, а также для противоположных им рёбер, вычисляем новую пропускную способность. Если она стала ненулевой, добавляем ребро к остаточной сети, а если обнулилась, стираем его;

ѕ Возвращаемся на шаг 2.

Важно, что алгоритм не конкретизирует, какой именно путь мы ищем на шаге 2 или как мы это делаем. По этой причине алгоритм гарантированно сходится только для целых пропускных способностей, но даже для них при больших значениях пропускных способностей он может работать очень долго. Если пропускные способности вещественны, алгоритм может работать бесконечно долго, не сходясь к оптимальному решению. На рисунке 1 представлена общая схема работы программы.

Более формально, мы имеем ориентированный граф G = (V; E), со множеством вершин V и множеством ребер E. Во множестве V выделены две вершины: s - источник и t - сток. При этом, у вершины s есть только ходящие ребра, у вершины t есть только входящие ребра.

Если ребро имеет максимальную пропускную способность x и потокy, то остаточной пропускной способностью ребра называется число xy. Смысл остаточной пропускной способности заключается в том, на сколько еще можно увеличить поток на данном ребре.

Следующим этапом выполнения программы является нахождение маршрутов в графе. Поиск маршрутов реализуется с помощью алгоритма обхода графа в ширину. На рисунке 2 представлена блок-схема алгоритма обхода графа в ширину.

Алгоритм Форда-Фалекрсона работает пока можно найти хотя бы один путь из источника в сток, итеративно увеличивая поток на минимальную величину. На рисунке 3 представлена блок-схема алгоритма Форда-Фалкерсона.



Рисунок 1 - Общая схема

Рисунок 2 - Блок-схема алгоритма поиска в ширину



Рисунок 3 - Алгоритм Форда-Фалкерсона

2.2 Анализ сложности алгоритма

На каждом шаге алгоритм добавляет поток увеличивающего пути к уже имеющемуся потоку. Если пропускные способности всех рёбер -- целые числа, легко доказать по индукции, что и потоки через все рёбра всегда будут целыми. Следовательно, на каждом шаге алгоритм увеличивает поток, по крайней мере, на единицу, следовательно, он сойдётся не более чем за O(f) шагов, где f -- максимальный поток в графе. Можно выполнить каждый шаг за время O(E), где E -- число рёбер в графе, тогда общее время работы алгоритма ограничено O(Ef).

Если величина пропускной способности хотя бы одного из рёбер -- иррациональное число, то алгоритм может работать бесконечно, даже не обязательно сходясь к правильному решению. На рисунке 4 представлен псевдокод алгоритма Форда-Фалкерсона.

2.3 Псевдокод

Ford_Fulkerson(G, s, t)

1 for каждого ребра (u, v) є E[G]

2 do f[u,v] «-0

3 f[v,u] «- 0

4 while существует путь p из s в t в остаточной сети G_f

5 do сf(р) «-- min {c_f(u, v) : (u, v) принадлежит р}

6 for каждого ребра (u, v) in p

7 do f[u,v] «- f[u] + c_f(p)

8 f[v,u] «- -f[u,v]

Рисунок 4 - Псевдокод алгоритма Форда-Фалкерсона

3. Распараллеленный вариант алгоритма

3.1 Библиотека omp.h

Разработчики OpenMP хотели предоставить простой способ создания потоков в приложениях, не требуя от программиста знаний о создании, синхронизации и уничтожении потоков, а также необходимости определения, сколько потоков следует создать. Для этого разработчики OpenMP создали независимый от платформы набор прагм, директив, вызовов функций и переменных среды, которые явным образом указывают компилятору, как и где именно следует вставить потоки в приложение. Большинство циклов можно распараллелить, вставив всего одну прагму непосредственно перед циклом. Более того, оставив исполнение рутинных функций компилятору и OpenMP, вы можете больше времени уделить определению того, какие циклы следует распараллелить, и как наилучшим образом изменить структуру алгоритмов для достижения максимальной производительности. Максимальная производительность OpenMP реализуется при использовании этого инструмента для распараллеливания "горячих точек", - наиболее трудоемких циклов в приложении.

Написание параллельного варианта сводится к написанию прагм. Прагмы - это директивы препроцессора, игнорируемые компиляторами, которые не поддерживают данную технологию. То есть программа, написанная с применением данной технологии, могут быть скомпилированы как последовательные, при том среда разработки не выдаст никакой ошибки. Список директив, использованных в данной работе:

1) #pragma omp parallel - открывает новый регион для распараллеливания внутренних инструкций. Программа делится на потоки, число которых можно указать с помощью функции omp_num_threads(int n), где n - число потоков.

2) #pragma omp for - применяется для распараллеливания итерационных циклических операций. То есть, несколько итераций выполняется параллельно на нескольких процессорах.

Так как потоки выполняют свои операции параллельно, то заканчивать свою работу будут в разное время, поэтому данный вариант распараллеливания будет бесполезен при использовании на задачах класса P = NP.

3.2 Краткое описание распараллеливания алгоритма

Процесс распараллеливания алгоритма Форда-Фалкерсона сводится к распараллеливанию обхода графа в ширину. Если же при реализации алгоритма, для поиска путей в графе использовать, например, алгоритм обхода в глубину - распараллеливание не представляется возможным.

При распараллеливании используется библиотека omp.h, использующая инструменты OpenMP технологии. Алгоритм делится на несколько потоков, и на каждом из потоков выполняется своя часть программы. Параллелизм реализуется за счет написания прагм в соответствующих фрагментах кода. На рисунке 5 представлен фрагмент кода, с использованием нижеописанных прагм. Полный листинг представлен в приложении Б.

Основной целью распараллеливания был анализ ускорения работы алгоритма. На практике, ускорение заметно только на очень больших наборах данных, но совсем незначительное. Ожидаемое ускорение, вычисленное по формуле Амдаля, при 25% операций выполненных параллельно на 2 процессорах составляет 300% от последовательного выполнения, однако на практике такого ускорения достигнуто не было, в виду целого ряда факторов, влияющих на быстроту выполнения операций.



Рисунок 5 - Фрагмент распараллеленного кода программы

Заключение

алгоритм форд фалкерсон граф

В процессе создания данного проекта была рассмотрена теория графов в общем, изучены понятия путей и маршрутов. Разработана программа, реализующая алгоритм Форда-Фалкерсона в MicrosoftVisualStudio 2010.Алгоритм может найти свое применение в программах для транспортных и коммуникационных сетей, таких как коммутация информационного пакета в Internet, где вершины - роутеры, а ребра это связи между роутерами. Таким образом, чем короче будет путь, тем быстрее пакет достигнет пункта назначения и меньше будут информационные потери.

Также были углублены знания, полученные в процессе выполнения лабораторных работ по предмету «Программирование».

Список использованных источников

1. ГОСТ 7.32-2001. Отчёт о научно-исследовательской работе. Структура и правила оформления [Текст]. - Взамен ГОСТ 7.32-91 ;введ. 2001-07-01. - Минск :Межгос. совет по стандартизации, метрологии и сертификации ; М. : Изд-во стандартов, 2001. - 16 с. - (Система стандартов по информации, библиотечному и издательскому делу).

2. ГОСТ 7.1-2003. Библиографическая запись. Библиографическое описание. Общие требования и правила составления [Текст]. - Взамен ГОСТ 7.1-84, ГОСТ 7.16-79, ГОСТ 7.18-79, ГОСТ 7.34-81, ГОСТ 7.40-82 ;введ. 2004-07-01. - М. : Изд-во стандартов, 2004. - 116 с. - (Система стандартов по информации, библиотечному и издательскому делу).

3. Канева, Ольга Николаевна. Дискретная математика [Текст] : учеб.пособие / О.Н. Канева, 2009. - 87 с.

4. Левин, А. Самоучитель работы на компьютере в вопросах и ответах[Текст]: учеб. пособие / А. Левин -М.: «Аквариум», 2007г. - 644с.

5. Алгоритм Форда-Фалкерсона [Электронный ресурс]: максимальный поток транспортной сети. - http://habrahabr.ru/- Загл. с экрана.

6. Алгоритм Распараллеливание поиска в ширну [Электронный ресурс]: оптимизация алгоритмов. -http://neerc.ifmo.ru/mediawiki/index.php?title.redirect=no/- Загл. с экрана.

Приложение

Приложение А

(обязательное)

Текст основного файла

#include"ff_alg.h"

int main()

{

clock_t time = clock();

read_input_file();

printf("%d\n",max_flow(0,n-1));

time = clock() - time;

cout<<time<<endl;

system("pause");

return 0;

}

Приложение Б

(обязательное)

Текст включаемого файла

#include<stdio.h>

#include<iostream>

#include<time.h>

using namespace std;

#ifdef _DEBUG

#undef _DEBUG

#include <omp.h>

#define _DEBUG

#else

#include <omp.h>

#endif

#define WHITE 0

#define GRAY 1

#define BLACK 2

#define MAX_NODES 100

#define INFINITY 10000

int min (int x, int y);

void enqueue (int x);

int dequeue () ;

int bfs (int start, int target) ;

int max_flow (int source, int sink);

int n;

int e;

int capacity[MAX_NODES][MAX_NODES];

int flow[MAX_NODES][MAX_NODES];

int color[MAX_NODES];

int pred[MAX_NODES];

int head,tail;

int q[MAX_NODES+2];

int min (int x, int y) {return x<y ? x : y;}

void enqueue (int x)

{

q[tail] = x;

tail++;

color[x] = GRAY;

}

int dequeue ()

{

int x = q[head];

head++;

color[x] = BLACK;

return x;

}

int bfs (int start, int target)

{

int u,v;

for (u=0; u<n; u++)

{

color[u] = WHITE;

}

head = tail = 0;

enqueue(start);

pred[start] = -1;

omp_set_num_threads(32);

#pragma omp parallel

while (head!=tail)

{

u = dequeue();

//#pragma omp for

for (v=0; v<n; v++)

{

if (color[v]==WHITE && capacity[u][v]-flow[u][v]>0)

{

enqueue(v);

pred[v] = u;

}

}

}

return color[target]==BLACK;

}

int max_flow (int source, int sink) {

int u;

int max_flow = 0;

for (int i = 0; i < n; i++)

{

for (int j = 0; j < n; j++)

{

flow[i][j] = 0;

}

}

omp_set_num_threads(32);

//#pragma omp parallel

while (bfs(source,sink))

{

int increment = INFINITY;

//#pragma omp for

for (u=n-1; pred[u]>=0; u=pred[u])

{

increment = min(increment,capacity[pred[u]][u]-flow[pred[u]][u]);

}

#pragma omp for

for (u=n-1; pred[u]>=0; u=pred[u])

{

flow[pred[u]][u] += increment;

flow[u][pred[u]] -= increment;

}

max_flow += increment;

}

return max_flow;

}

void read_input_file()

{

int a,b,c;

FILE* input = fopen("C:\\data.txt","r");

fscanf_s(input,"%d %d",&n,&e);

for (int i = 0; i < n; i++)

{

for (int j = 0; j < n; j++)

{

capacity[i][j] = 0;

}

}

for (int i = 0; i < e; i++)

{

fscanf_s(input,"%d %d %d",&a,&b,&c);

capacity[a][b] = c;

}

fclose(input);

}

Размещено на Allbest.ru