Оптимизация маршрута с использованием теории графов в пакетах прикладных программ

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

Е.Ю. Вардомацкая, В.Л. Шарстнев, УДК 004.9 : 658

Я.А. Алексеева

РЕФЕРАТ

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

В статье приведен анализ методов решения задачи о кратчайшем пути с использованием теории графов в пакетах прикладных программ. Для возможности применения системы компьютерной алгебры данная проблема была представлена в виде математической модели на графе. Подробно рассмотрены несколько вариантов решения искомой задачи с помощью СКА Maple. Так, решение данной задачи реализовано методом Дейкстры, методом имитации отжига и методом муравьиного алгоритма. СКА Mathematica также обладает расширенной поддержкой графов, необходимой для решения задачи о кратчайшем пути. \___

В современных условиях развития рыночных отношений всё более актуальным становится всестороннее обеспечение конкурентоспособности как экономики Республики Беларусь в целом, так и её отдельных отраслей. В частности, легкой промышленности - одного из важнейших секторов экономики, который производит большое количество товаров народного потребления. В последние годы в сфере товарного обращения ряда стран произошли существенные преобразования. В хозяйственной практике стали использоваться новые методы и технологии доставки товаров, базирующиеся на концепции логистики. Важную роль в создании объективных возможностей для развития логистики сыграл технический прогресс в средствах связи и информатики. Он позволил на более высоком уровне проводить отслеживание всех основных и вспомогательных процессов товародвижения.

ABSTRACT

GRAPH THEORY, MATHEMATICAL MODELS, FLOW NETWORKS, ROUTE OPTIMIZATION THE SHORTEST PATH, TRANSPORTATION EXPENSES, COMPUTER ALGEBRA SYSTEM

The analysis of methods for solving the using graph theory in application packages is represented in this article. For possibility of use of computer algebra system this problem was presented in the form of а mathematical model on the graph. Some versions of the solution of the problem by means of various CAS are in detail considered. Indisputable leaders in in this area are CAS Maple and CAS Mathemat-ica. The solution of this problem by means of CAS Maple is realized by Dijkstra's algorithm, Simulated annealing and Ant colony optimization algorithms. CAS Mathematica has expanded support graphs necessary for solution of the shortest path problem. _)

Так, постепенно на первый план выдвигается поиск возможностей сокращения производственных затрат, в том числе транспортных расходов.

Цель исследования - разработка средств компьютерного моделирования процесса управления цепями поставок.

Объект исследования - логистические системы предприятий легкой промышленности.

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

Методы исследования - абстрагирование, математическое моделирование, анализ.

Инструментарий исследования - система компьютерной математики (далее - СКМ) Maple.

Научная новизна исследования заключается

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

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

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

Таким образом, транспортные сети удобно представить в виде графа.

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

Транспортная сеть может быть представлена только связным графом.

Моделирование транспортной сети начинают

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

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

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

Для снижения размерности задачи и ускорения расчетов для транспортных сетей больших городов или районов применяют микро- и макрорайонирование. При микрорайонировании транспортной сети в качестве вершин используют не пересечения улиц (дорог) и конкретные пункты отправления и назначения, а центры микрорайонов (районов получения или назначения грузов). Макрорайонирование транспортной сети заключается в разбиении ее на отдельные подсети, расчеты по которым выполняются раздельно, а затем объединяются для получения общего результата. При изменениях дорожной обстановки производится пересчет отдельной подсети, в которой произошли изменения, но не по всей транспортной сети [5, с. 171-172].

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

роко используется для нахождения путей между физическими объектами на картографических сервисах, в GPS-навигаторах, при определении наименьшего расстояния в сети дорог.

Выделяют несколько вариантов решения задачи о кратчайшем пути на графах с помощью различных пакетов прикладных программ. Бесспорным лидером в данной области является СКА Maple. Рассмотрим решение исходной проблемы на примере использования данного программного пакета, имеющего специализированную библиотеку networks для работы с графами.

Выделяют несколько способов решения данной задачи - это, прежде всего, метод Дейкстры, метод отыскивания всех гамильтоновых циклов, алгоритм ближайшего соседа (Nva), муравьиный алгоритм, метод имитации отжига и другие. Остановимся на некоторых из них более подробно.

Наиболее точным и при этом простым методом решения транспортной задачи в СКА Maple на графах является метод Дейкстры - пошаговый алгоритм определения кратчайшего расстояния от вершины А до В.

Реализацию метода Дейкстры в СКА Maple проиллюстрируем решением следующей практической задачи. Предприятие легкой промышленности города Витебска наладило поставки своей продукции в магазины регионов Респуб-

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

В программе, представленной на рисунке 2, реализован данный метод на примере орграфа с неотрицательными весами, смоделированного на основе путей сообщения между пятью городами Республики Беларусь (рисунок 1).

Ответ заносится в переменную MinPath. Алгоритм заканчивает свою работу, когда flag принимает значение true, то есть конечная вершина (target) приобретает постоянную метку. После завершения работы программы список постоянных меток можно посмотреть, раскрыв переменную V: evalm(V) [3, с. 111].

Также кратчайший путь в орграфе можно найти, используя стандартные процедуры СКА Maple. Рассмотрим заданный в алгоритме Дейкстеры граф и вычислим кратчайший путь. В программе на рисунке 3 независимо используются два оператора, shortpathtree и allpairs.

В алгоритме ближайшего соседа (рисунок 4) выбор дальнейшего пути производится, исходя из минимального расстояния до очередной вершины. Обход графа начинается с произвольной вершины, от выбора которой часто зависит результат. Движение по графу происходит по исходящим (Out) дугам минимального веса с по-

г

Рисунок 1 - Ориентированный взвешенный граф

У

> restart. wiife [networks) .

> : п !— 5 :

> adebertex(itt = 1 щ G):

> added&{[[\42\ [г й [1,31 [3, 111X21 [5, Щ [4.31 [5,41 [4, 2\ [% А\ [*:5}. [5. 1Ц weights = [ 110,111, 35,34, 210, 140, 280, 130, 215.5, 215, 230, 200]. С?) :

> У:= Vectori 1 ..й) :

> for; to п do P[i] ~ infinity, od

> s ~ [. target ~ 5 : к '= s : V\k) ¡= 0 :U~ | 0$3-n ] :

> flag s— false :

> for ? while not/7<ig do

U[i ] ¡= к :

d «— outdegree[k, G) : z G) :

for j to d do CW1 := eweighi

«( (n^V]]« o) or №Ш] > + ПК)) then V[z[]]] := eveigkl{op(edgeH[k,z[j]j, (?)), G) + V\k\: fl; od, Next •— r4 for/ from 2 to fi do

iSnatmember(j, U) and V[}] < V\N$Xt] then Next к == Next,

flag is(k - target)', od

> evab(t)i

| 169 110 Z5 325 340 |

"> Mr,Path ;= V[4\t

Mir,Path := 325

Рисунок 2 - Программная реализация метода Дейкстры

К.

г

> restart: with(networks):

> iiew(G) : ft 5 :

> V == [ I 2, X 4, 51: n nops(F):

> c/ddvirrfezi i' G] :

оЫиВДШ.П [2, 1}, [1,3]. [ЭИ1 ($21 [5,3], [4,3], [5,41 [4.2], [2,4], [2, Я [5, 1]], weights =*[110, 111, Ё5, 84, 210, 140, 230, 130, 215.5, 215, 230, 200], G) :

> T s= shortpathtrss {G, 1) :

> vweight(T),

W ~tabb

> MmPatk W{4];

MinPath := 325

> alipairs[G)[ 1,4];

325

Рисунок 3 - Программа нахождения кратчайшего пути с использованием операторов shortpathtree и allpairs

У

г

_> restart: witk{netu!orks) :

> nsw[G} : п := 5 :

> addvertex{i%i - 1£?);

1, 2. 3.5

> addedge ([[1,21, [2, 1],[U 3], [3, 1 ]. [3, 2], [5, 3], [4, 3]. (5, 4], [4.2\ [X А\[2, S\ [5,1]],

wet&tts= [ПО, 111. 05, 84, 210, 140, 280, ПО, 215 5, 215, 2Ж 200], £?) : L> Gl ¡= dupiicate(G) ;

> aO — 5 : # Начальная вершина

> ai := ¿fO : rw i= 0 ; J_> i == faO]:

> foi' k to n -1 do

for v in incidental, GI, Out) do

if eweightiv, Gi) = inm{e>p[eweight(\op(incident{.a!, GI, Out)) J, Gl))) (lien и := v, iw «— i* + eweight(v, GI); break; Ё, od, a2 — ends(uy Gi )[2 |, delete {a I, Gl ] : ai = a2 s [op(s)7 a2\:

od:

> jm fw -f- eweigkt(op (edges (\ a2, aO |, O)), G)", # Сумма

sw := 541.5 + table

> s Контур

5.4, 2, 1,31

V.

Рисунок 4 - Программная реализация метода ближайшего соседа

J

следующим удалением пройденной вершины. Поэтому перед началом процесса необходимо запастись дубликатом графа, из которого затем будем извлекать информацию о длине пройденного пути. Оператор break досрочного выхода из цикла введен для случая двух или более одинаковых дуг минимального веса, выходящих из текущей вершины [3, с. 145]. Если дать возможность алгоритму случайно выбирать решение, оптимальное на каждом шаге, то можно пропустить ход, не лучший локально, но дающий в результате более эффективное решение, что, несомненно, является существенным недостатком данного алгоритма. Однако несмотря на предельную простоту алгоритма, для небольших графов он дает достаточно близкие к точному решению ответы.

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

присваиваем СЦ,]]:=РЦ,}\, так как простое присваивание С:=Р, предусмотренное для работы с матрицами, работает ненадежно и часто является причиной ошибок [3, с. 143\.

Для решения задачи о кратчайшем пути необходимо вычислить количество n1:=nops(H) гамильтоновых циклов и найти минимальный из них. Для этого разработана процедура вычисления веса дуги и-у в цикле к и использована стандартная функция отыскания минимума последовательности (рисунок 6) [3, с. 146\.

Время вычисления при этом будет в несколько раз больше, и с увеличением порядка графа разность в быстродействии программ будет рас-

Несомненное достоинство алгоритма - его точность и простота.

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

г

> rsstart: witk(networks) : rmw{G) :

> V = [ а, Ь, с, d, g |. » nops{ У) :

> addvertex

[> addedgei ([а, с}, {a,d}, fb.g), {b.c). {c, d}, {d, g}. [a, gl [a. ЬЩ G).

> Б = Matrix(n):

> -4 adjacency(G):

> foritondi)

forj to n do if A[i, j] = J then Б|i. j] == Vjj]; fi, od. od

'> С >- A :

> forkton-ldo

for i to n do

forj ton do P[i,ih= evp:jnd[sijbs(V[i] = Q, V[j] = Q, addr B; i, ro ]tS*C] m, j], m = 1 .,11)}); od; od;

for i to 11 do

forj to n do Ciüjl í= P|i,j3; od: od; od:

> {>:

> foritondn ifCfviH Otilen S S unión {¡^afld(C[i Í]«fe*V[i]))p II; od

> F{}:

> for qin S do if whattype(q) =' +' tJipn F = Funión {op(q)} else F = Funión  q) : fi: od;

> S ™ map{x —» convertíx, Iisl\, F)

> for qin S do iftiopsjqi Ф л then S S niinus (q); fi; od;

> к RQps[S) :

:> я - Í j -

> forj to к do

for i to H do LT^iIIi] = V\1 ] tJienm i, fi: od,

for i ta n da z [i ] j II (í + m - 2 mod n) + 1 ]; od; H # unión {eortvert(Z, ftsí)}:

od

> H,

[[a.b, c, g. d.f\ [o. b, g, c. d, fi [et,f,d, c. b.gi [a,ftd,c,g,b\)

Рисунок 5 - Программа отыскивания всех гамильтоновых циклов

Рисунок 6 - Продолжение программы отыскивания гамильтоновых циклов

J

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

Рассмотрим программу минимизации транс-

портных расходов для п городов, координаты х[г] и у[г] которых заданы (рисунки 7, 8). В начале программы вводятся константы задачи. Улучшение сходимости во многом зависит от их значений. Число Ьтгп требуется для сравнения при выборе минимального маршрута, масштабная константа Q порядка длины маршрута выбирается пропорциональной порядку графа.

> tor кО to 100 do U Основной цикл. for ant to n do # Цикл no муравьям s {$ i } ; # Список непосвященных вершин j — ant : U Начальная вершина для маравъя ant for jl to n — 1 do s s minus {j}; # Tabu list увеличился sp [ ] : Ш i— 0 :

for / in s do # Каждой вершине - свой вес sp — [op{sp)7\í W[j, i* Wt[л if^q] : kl + 1 : №î[jti] := í; od:

jO := nn\P{sp) j; # Выбор направления

j jO : # Начало дуги - конец предыдущей od: # Цикл jl по всем вершинам для муравья ant L add{ /[/],/=! ,.n - 1 ) -h W [ op (s), ant]; vfO] := ant, v\ n ] := ant7 # Начало и конец дуги vi -— seqг( [ v[m]7 v[m + 1 ] ]_, m = 0 L ) : U Дуги for / to n do U Пометка дуг феромоном

DWt[op{vl{i])\ := DWt[ op{vl[i]) ] + Q! L\ od:

ifL < Lmin t hen Lmin :== L :fi: # Выбор min od: # л/?/

И7? := Wt + DWtwp : # Добавление новых следов Wt •= Wt* (1 -p) : U Испарение феромона od.tikO

Рисунок 8 - Основной фрагмент программы муравьиного алгоритма

Ч.

функция Р вероятности перехода имеет в качестве аргумента список из чисел х[г]. Вычисля-

ется сумма sm= Yxftl Отрезок прямой от 0 до 1 разбивается на n участков [Beg[i], End[i]] с

длинами x[i]/sm, где Beg[1] =0, End [n] = 1. Затем случайное число 0 < s < l указывает выигрышный номер - номер вершины для дальнейшего движения [3, с. 147].

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

В методе отжига (рисунок 9) очередной порядок следования по маршруту между городами выбирается случайно, небольшим изменением предыдущего решения, предположительно оптимального. Самый простой вариант изменения - перестановка двух случайно выбранных городов в маршруте следования. Если полученный маршрут лучше всех существовавших ранее, то этот маршрут берется за очередной, если маршрут хуже, то самый простой вариант - не брать его. Следует отметить, что решение зависит от нескольких параметров. Меняя, например, число циклов, можно получать различные ответы [3, с. 89].

Таким образом, рассмотрев возможные решения задачи о кратчайшем пути на графах, нельзя однозначно выделить один из методов решения, каждый из них имеет как очевидные достоинства, так и недостатки, такие как, необходимость знания Maple-языка, являющегося функционально полным процедурным языком программирования четвертого поколения (4GL), ручной ввод исходных данных и другие. Однако согласно проведенному анализу, более предпочтительным методом решения исходной проблемы в СКМ Maple является метод Дейкстры.

Экономический эффект от внедрения средств

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

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

Рассмотрим на конкретном примере возможный экономический эффект от внедрения разработанной методики оптимизации транспортных расходов по оценочным данным на примере предприятия легкой промышленности города Витебска, занимающегося производством трикотажных изделий, в котором автоматизируется поиск оптимального маршрута транспортировки произведенной продукции в фирменные магазины. Средняя стоимость пакета СКМ Мар1е составляет 10,25 млн. руб. Предполагаем, что внедрением, технической поддержкой и возможными доработками разработанной методики в данных СКМ будет заниматься дополнительно обученный штатный работник. Стоимость дополнительного обучения примем равной 1 млн. руб.

В итоге затраты на внедрение нового программного обеспечения составят:

Рисунок 9 - Фрагмент программы "Метод имитации отжига"

-У

рассчитан по формуле

Зм*Р,е = 10,25 + 1 = 11,25 млн. руб.

Для оценки экономии от разработанной методики в качестве критерия можно выбрать:

• значение коэффициента использования пробега;

• оптимизация маршрута (его протяженность, пропускная способность, скорость, загруженность);

• затраты топлива на 1 т перевезенного груза.

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

Рассчитаем экономию за счет сокращения транспортных расходов на топливо. По оценочным данным, средняя величина сокращения протяженности пройденного за счет оптимизации маршрута на одну поездку составляет около 60 км, норма потребления грузовым автомобилем топлива на 100 км в загруженном состоянии может достигать 15 л. Проанализировав текущую ситуацию на рынке автомобильного топлива, можно сделать заключение, что средняя цена топлива составляет 11 000 руб. Экономия, связанная с сокращением транспортных расходов на топливо на одну грузоперевозку, составляет:

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

Тогда экономический эффект от внедрения нового программного обеспечения может быть

где Э - годовая суммарная экономия, млн.руб.; З - годовые суммарные затраты, связанные с внедрением нового программного обеспечения, млн. руб.; к - количество лет реализации проекта плюс один год, лет.

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

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

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

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

Р = 99000 х

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

предложенной методики на предприятиях промышленности и бытового обслуживания населения всех форм собственности Республики Беларусь.

СПИСОК ИСПОЛЬЗОВАННЫХ REFERENCES

ИСТОЧНИКОВ

1. Аладьев, В.З., Бойко, В.К., Ровба, Е.А. (2011), Программирование в пакетах Maple и Mathematica: Сравнительный аспект, Гродно, 2011, 516 с.

2. Березина, Л.Ю. (1979), Графы и их применение, Москва, Просвещение, 1979, 143 с.

3. Кирсанов, М.Н. (2007), Графы в Maple. Задачи, алгоритмы, программы,Москва, 2007, 168 с.

4. Кристофидес, Н. (1978), Теория графов. Алгоритмический подход, Москва, Мир, 1978, 432 с.

1. Alad'ev, V.Z., Bojko, V.K., Rovba, E.A. (2011), Programmirovanie v paketah Maple i Mathematica: Sravnitel'nyj aspekt [Programming in the Maple and Mathematica packages: A comparative aspect], GGU, Grodno, 2011, 516 p.

2. Berezina, L. Ju. (1979), Grafy iihprimenenie [Graphs and their application], Prosveshhenie, 1979, 143 p.

3. Kirsanov, M.N. (2007), Grafy vMaple. Zadachi, algoritmy, programmy [Graphs in Maple. Tasks, programs, algorithms], Moscow, 2007, 168 p.

5. Хлевной, И.И. (2006), Грузовые перевозки, СПб., 2006, 290 с.

6. Расчет экономического эффекта от внедрения системы автоматизации [Электронный ресурс] / Компания «Antegraconsulting». - Режим доступа: http://www.antegra.ru. - Дата доступа: 14.05.2014.

7. Sharstniou, U. L., Vardamatskaja, A.U. (2007) Computer information technology: software packages for modeling and analysis of problems in economics: a tutorial [Компьютерные информационные технологии: пакеты прикладных программ для моделирования и анализа задач экономики, Vitebsk EE «VSTU», 2007. 138str

4. Kristofides, N. (1978) Teorija grafov. Algoritmicheskij podhod [Graph theory. An algorithmic approach], Moscow, Mir, 1978, 432 p.

5. Hlevnoj, I.I. (2006), Gruzovye perevozki [Freight transport], SPb., 2006, 290 p.

6. The calculation of the economic effects of the introduction of automation systems [electronic resource] [Raschet jekonomicheskogo jeffekta ot vnedrenija sistemy avtomatizacii] [Elektronnyjresurs] / Kompanija «Antegraconsulting». - Rezhim dostupa: http://www. antegra.ru. - Data dostupa: 14.05.2014.

7. Sharstniou, U. L., Vardomatskaja, A.U. (2007) Computer information technology: software packages for modeling and analysis of problems in economics: a tutorial [Kompjuternye informacionnye tehnologii: pakety prikladnyh programm dla modelirovanija i analiza zadach economiki], Vitebsk EE «VSTU», 2007. 138 str.

Статья поступила в редакцию 18. 01. 2016 г.