Краткие общие сведения

Общепринятым методом решения транспортной задачи является так называемый "метод потенциалов". Приведем структуру соответствующего алгоритма.

1. Поиск начального решения, т.е. набора { X_ij }, удовлет-воряющего ограничениям (5.3).

2. Если{X_ij} -оптимальное решение,токонец.

3. Модификация решения, не увеличивающая его стоимости, и пе­реход на 2.

Задача (5.2), (5.3) является частным случаем задачи линейного программирования, а вышеприведенный алгоритм укладывается в общую схему решения этой задачи.

Выделим специфические черты, отличающие метод потенциалов.

1. В силу (5.1) условие (5.3) содержит ровно m+n-1незави­симых ограничений,и, следовательно, любое базисное решение транс­портной задачи имеет не болееm+n-1ненулевых компонент.

2. Транспортная задача допускает несколько простых способов поиска начального базисного решения. Один из возможных алгоритмов реализован в лабораторной работе:

а) выбрать: в табл.5.1 незаполненную и непомеченную клетку с минимальной стоимостью (пусть, для определенности, это клетка с индексами (k, l));

б) заполнить выбранную клетку числом

X_kl = min { α_k, B_l },

т.е. максимально загрузить соответствующий путь;

в) модифицировать величины спроса и предложения:

1) α_k=α_k-X_kl

B_l = B_l – X_kl.

2) еслиα_k= 0,

топометить все незаполненные и непомеченные клеткиk-й строки,

иначепометить все незаполненные и непомеченные клеткиl-го столбца;

д) еслиостались незаполненные или непомеченные клетки,

идтина а).

Заполненные клетки будем называть базисными (некоторые из них могут содержать нули).

- 21 -

Заполненные клетки будем называть базисными (некоторые из них могут содержать нули).

3. Проверку оптимальности полученного решения удобно осущест­влять, переходя к двойственным переменным. Не входя в подробности обоснования, приведем соответствующий алгоритм:

а) найти "потенциалы" и, такие, что

U_i + V_j = C_ij (5.5)

для всех пар (i, j), выделяющих базисные клетки. Легко видеть, что одну из переменныхU_iилиV_j(по нашему выбору) можно поло-жить равной чему угодно, после чего система (5.5) решается одно­значно;

б) для всех небазисных клеток вычислить так называемые "ха­рактеристические разности"

X_ij = U_i + V_j - C_ij

в) есливсеX_ijнеотрицательны,

торешение оптимально,

иначесм. следующий пункт.

При решении задачи вручную проверку оптимальности удобно произво­дить с использованием табл.5.2.

Таблица 5.2

C11	...	C1n	Ui
X11		X1n
...	...	...	...
Cm1	...	Cmn	Um
Xm1		Xmn
V1	...	Vn

4. В случае, если полученное решение неоптимально, можно от-носительно просто улучшить его:

а) выбрать какую-либо клетку, для которой X_ij < 0;

б) найти последовательность клеток (цикл)

(i₁=i, j₁=j), (i₂,j1), (i₂,j₂), ..., (i_p,j_p-1), (i_p,j_p), (i_p+1=i₁,j_p), в которой все клетки, кроме клеток(i₁,j₁),являются базисными;

в) перераспределить грузопоток по найденному циклу:

пропустить дополнительно Xединиц груза по каждому из пу- тей(i_k, j_k),

уменьшить на Xколичество груза, перевозимого по путям

(i_k+1, j_k).

- 22 -

Здесь

(5.6)

г) включить в базис клетку (i, j)и исключить из базиса одну из клеток, поток в которой после модификации, проведенной в п. в), оказывается равным 0.

Вообще говоря, п. в) вышеприведенного алгоритма следует сводить к задаче поиска цикла в ориентированном графе (связанном с модифици­руемым решением), однако при небольших размерах задач, решаемых в ходе лабораторной работы, цикл можно найти путем простого направ- ленного перебора в глубину. Следует запомнить, что величина Х, вычисляемая в (5.6), может оказаться равной 0. При использовании такого цикла меняется не грузопоток, а лишь набор базисных пере­менных.