О способах оценки числа планов транспортной задачи
О способах оценки числа планов транспортной задачи
Аннотация
Код статьи
S042473880012408-7-1
Тип публикации
Статья
Статус публикации
Опубликовано
Авторы
Погодин И. Е. 
Аффилиация: ВМПИ ВУНЦ ВМФ ВМА
Адрес: Российская Федерация, Санкт-Петербург
Выпуск
Страницы
116-120
Аннотация

В связи со значительным расширением круга проблем, образующих класс так называемых транспортных задач, включая выход на нелинейные, целесообразно иметь оценки числа всех допустимых планов, среди которых выбираются оптимальные. Рассматриваются оценки количества возможных планов решения замкнутых транспортных задач для матриц различной размерности и структуры. В качестве базы исследования проанализировано около пятидесяти примеров различных транспортных задач. Обнаружено, что при перераспределении между собой, к примеру, только мощностей производителей и постоянстве их суммы число возможных планов задачи монотонно уменьшается с увеличением их относительного среднеквадратического отклонения. Из анализа частных эмпирических зависимостей для матриц различной структуры получены аналитические обобщения для простейших ситуаций. Получено, что для задач: (а) с размерами (2×M) и (N×2) возможен прямо аналитический подсчет числа планов; (б) с матрицами произвольных размеров и структур требуется компьютерный алгоритм из (M) вложенных циклов; (в) с равными мощностями и равными емкостями возможен вероятностный способ расчета оценок, результаты которого коррелируют с точными на уровне 0,8. Приведены конкретные алгоритмы оценки числа допустимых планов. Задача может представлять интерес при оценке эффективности различных методов оптимизации.

Ключевые слова
планы транспортной задачи, поставщик, заказчик, поставка, емкость, ограничение перевозки.
Классификатор
Получено
01.12.2020
Дата публикации
16.12.2020
Всего подписок
14
Всего просмотров
1066
Оценка читателей
0.0 (0 голосов)
Цитировать   Скачать pdf
1 ВВЕДЕНИЕ
2 Транспортная задача , задача о назначениях и другие распределительные задачи обычно являются частными задачами линейного программирования , у истоков создания которого еще в 1930-е годы стоял академик Л.В. Канторович (Канторович, 1939), который, в частности, начал делать метрологические оценки количества вычислений, необходимых для решения транспортных задач классическими методами. Однако даже для таких относительно простых задач характерно большое число переменных решения . Для транспортных задач , имеющих практическое значение, применение общих методов может стать неэффективным. Кроме того, ряд современных модификаций транспортных задач (например, (Ассаул, 2019)) вообще выводят их из класса задач линейного программирования. Платой за универсальность методов и задач часто является снижение их эффективности, заключающееся в медленной сходимости, высоком объеме вычислений и т.п. (Худокормов, 1994).
3 Рассматриваемая задача оценки числа всех возможных планов замкнутой транспортной задачи может представлять интерес для оценки и сравнения степени эффективности различных, порой весьма трудоемких и приближенных, методов оптимизации при решении широкого класса транспортных и сводимых к ним задач с различными критериями оптимальности (Ассаул, 2019), например по сравнению со способом прямого перебора всех планов без какой-либо подготовки. Здесь становится уместным знать и сравнить трудоемкость прямого перебора всех планов, т.е. их число, подобно тому, как в статистических задачах обычно проводится сравнение с исключительно случайными условиями или с нормальным законом распределения. Кроме того, на практике встречаются чрезвычайные ситуации поиска исчезнувших из поля зрения в нескольких районах (одушевленных или неодушевленных) существ или предметов по местам их возможной концентрации также в различных районах, и когда требуется оценить число подлежащих рассмотрению (опросу) связей.
4 МЕТОДЫ ИССЛЕДОВАНИЯ И РАСЧЕТА
5 Если рассматривать только базисные планы транспортной задачи, то в качестве грубой верхней (существенно завышенной) оценки их числа можно принять, например, число сочетаний из (N×M) по (N+M–1) CMNM+N-1 (здесь N – число производителей, M – число заказчиков), которое зависит только от размеров задачи.
6 Как показало практическое рассмотрение более 50 замкнутых транспортных задач с размерами от ((N = 2)(M = 2)) до ((N = 3)(M = 4)) (табл. 1), полное число (Y) всех их возможных планов варьируется на 2 порядка (от 2 до 186), существенно зависит как от размеров таблицы N×M, так и от значений мощностей и емкостей. Однако попытка построить по массиву рассмотренных 54 матриц общую регрессионную зависимость вида YNαMβSγDδ, где S = N×A = M×B — полный поток (сумма значений мощностей или емкостей), D — сумма дисперсий мощностей и емкостей, оказалась несостоятельной (например, оценки α и β различаются вдвое и т.п.).
7 Таблица 1. Примеры фактического числа допустимых планов (Y) для задач с матрицами размером N×M и произвольным распределением мощностей ( Ai) и емкостей ( Bj )
8
{Ai}/{Bj} Y {Ai}/{Bj} Y {Ai}/{Bj} Y {Ai}/{Bj} Y {Ai}/{Bj} Y {Ai}/{Bj} Y
N×M=2×2 N×M=2×2 N×M=2×2 N×M=2×2 N×M=2×2 N×M=2×2
8,2/8,2 3 6,4/6,4 5 7,3/7,3 4 9,1/9,1 2 9,1/8,2 2 9,1/7,3 2
9,1/6,4 2 8,2/7,3 3 8,2/6,4 4 5,5/5,5 6 9,1/5,5 2 8,2/5,5 3
7,3/5,5 4 6,4/5,5 5 6,4/8,2 3 3,3/1,5 2 3,3/2,4 3 3,3/3,3 4
6,6/1,11 2 6,6/2,10 3 6,6/3,9 4 6,6/4,8 5 6,6/5,7 6 6,6/6,6 7
N×M=2×3 N×M=2×3 N×M=2×3 N×M=2×3 N×M=2×3 N×M=2×3
3,3/1,2,3 5 3,3/1,1,4 3 3,3/2,2,2 8 1,5/2,2,2 3 2,4/2,2,2 6 6,6/1,1,10 4
6,6/1,2,9 6 6,6/1,3,8 8 6,6/1,4,7 10 6,6/1,5,6 12 6,6/2,3,7 12 6,6/2,4,6 15
6,6/2,5,5 15 6,6/2,8,2 9 6,6/4,4,4 21 4,8/4,4,4 15 3,9/4,4,4 10 2,10/4,4,4 6
N×M=2×3 N×M=2×3 N×M=2×3 N×M=2×3 N×M=2×3 N×M=2×3
1,11/4,4,4 3 1,7/1,3,4 3 2,6/1,3,4 5 3,5/1,3,4 7 4,4/1,3,4 8 2,3/1,2,2 5
3,3/2,2,2 13 N×M=3×3 N×M=3×3 N×M=3×3 N×M=3×3 N×M=3×4
2×4 1,2,3/4,1,1 8 6,3,3/3,4,5 80 2,3,4/1,2,6 18 3,5,7/1,6,8 39 3,4,7/1,2,5,6 186
10,18/5,6,8,9 262
9 Однако для каждого подкласса матриц одинакового размера и полного потока искомое число планов Y обнаруживает уверенную статистическую тенденцию к уменьшению по мере увеличения различий (неравномерности, выражаемой, например, суммой их дисперсий D) среди групп значений мощностей и емкостей (рис. 1).
10 Рис. 1. Зависимость числа планов от суммы дисперсий мощностей и емкостей (для задач: а) 2×2 при S = 10 (r = –0,86); б) 2×3 при S = 12 (r = –0,74))
11 Это обстоятельство позволяет предположить, что в случае получения необходимых оценок числа планов для равномерных задач с одинаковыми значениями в группах мощностей (A) и емкостей (B) переход к конкретной задаче уже без этого условия можно получить по соответствующим регрессионным номограммам (см. рис. 1), предварительно скорректировав (интерполяцией) задачу от равномерной к реальным размерам матрицы соответствующего вида, либо использовать эту величину просто в качестве верхней оценки количества планов конкретной задачи. Во всех этих расчетах удобнее производителей и заказчиков располагать в порядке возрастания их мощностей и емкостей.
12 Для задач хотя бы с одним малым размером матриц число планов рассчитывается точно и достаточно просто.
13 1. При (2×2): Y = U = min(A1B1) + [1]; прибавлять единицу в квадратных скобках следует только если поставка Х11 может принимать значение «0», а это, в свою очередь, возможно, если А1 = D12 и В1 = D21.
14 2. При (3×3) (N = M = 3) число планов U определяется всевозможными сочетаниями значений четырех поставок (i, j, k, l) в левом верхнем углу таблицы:
15 max(0, A1D13, B1D31) ≤ i D11,max(0, A1i, A1D13) ≤ j D12,
16 max(0, B1i, B1D31) ≤ k D21,max(0, B2jD32, A2kD23) ≤ l D22
17 и может быть найдено сразу (т.е. без применения более сложной процедуры) в виде
18 Y=U=ijk(D22-max(0,B2-j-D32,A2-k-D23)).
19 3. При (N×2) (или (2×M)) в равномерных задачах также может быть получена точная оценка U = Y (пример для M = 4), которая будет получена как частный случай
20 Y=i=0Dj=max(0,(A-2D-i))min(D,(A-i))k=max(0,(A-D-j-i))min(D,(A-j-i))1.
21 4. Для матриц с произвольными размерами N×M и без условия равномерности распределений мощностей {Ai} и емкостей {Bj} можно обобщить п.п. (2) и (3) и точно рассчитать число планов. Однако для этого потребуется создать достаточно громоздкий (индивидуальный для каждой структуры N×M) алгоритм в виде вложенных однотипных N1×M1 циклов с накоплением единиц в первоначально обнуленном сумматоре S=S+1 в самом внутреннем цикле. При этом для индекса Ii, j, связанного с перевозкой в ячейке с адресом (i, j), следует организовать цикл по этому индексу, начиная с max(0,(Ai-s=j+1MDi,s-s=1j-1Ii,s),(Bj-s=i+1NDs,j-s=1i-1Is,j)) ) до min(Di,j,(Ai-s=1jIi,s),(Bj-s=1iIs,j))) , где Di,jmin(Ai,Bj) . Построенный таким образом алгоритм пригоден для различных значений мощностей {Ai} и емкостей {Bj}замкнутой транспортной задачи.
22 5. Для матриц с совершенно произвольными размерами N×M и равномерными распределениями мощностей {A} и емкостей {B} начнем с распределения заявок одного из участников какой-нибудь группы (например, поставщика) по соответствующим участникам из другой группы (заказчикам), что эквивалентно известной задаче о числе способов совершенно свободного (без ограничений) размещения S шаров по M ящикам (Генкин, 1994; Сергеев, 2011), в которой условно вводится (M–1) разделителей между M ящиками, расположенными в линию. После этого находят номера (M–1) позиций разделителей среди (S+M–1) мест, вариантом которых будет CS+M-1M-1 (число сочетаний из (S+M –1) по (M –1)).
23 Заметим, что если в каком-то одном произвольном из M ящиков оказалось U>DminA,B шаров, то для размещения остальных (SU) шаров вариантов будет CA-U+M-2M-2 , поскольку U шаров уже не требуется размещать и разделителей между остальными ящиками (кроме выбранного) потребуется (M –2). Такой недопустимо переполненной ячейкой может оказаться любая из M в рассматриваемой строке поставщика, а значит, можно предположить, что из величины CS+M-1M-1 следует вычесть M таких величин (CA-U+M-2M-2) . Однако, к примеру, переполнение может появиться более чем в одной ячейке строки, а следовательно, среди вычитаемых могут появиться повторяющиеся планы и вычитать предстоит несколько иное (большее) число планов. То же искажение дублированием некоторых планов может влиять на величину CS+M-1M-1 : Z=CA+M-1M-1-Mq=D+1ACA-q+M-2M-2 . Поскольку посчитать это измененное из-за дублирований уменьшенное число вычитаемых планов представляется затруднительным, это обстоятельство является недостатком предлагаемой вычислительной модели. Оценку для одной такой линейки можно получить простым перебором с ограничениями, как это сделано выше в (п. 3).
24 Наличие другой такой линейки (поставщика), если она не последняя (см. выше п.3), требует повторения этой процедуры, проделанной для первой рассмотренной линейки.
25 Далее рассмотрим создание планов для всей таблицы из построенных строчных планов. На примере одного столбца видно, что из n возможных планов по столбцам условию равенства суммы по столбцам величине B будут удовлетворять только m планов, т.е. можно говорить о вероятности p=m/n того, что построенные планы по столбцам проходят через фильтр условия заказчика. Снова используем модель распределения q шаров (q = 0, ..., ND) по N ящикам, а также (ND – j) шаров (j = (D+1), ..., ND) по N1 ящикам и получим
26 n=q=0NDCq+N-1N-1, m=q=DNDCq+N-1N-1-Nj=D+1NDC(ND)-j+N-2N-2.
27 Поскольку такие вероятности отбора планов по столбцам с суммой, равной B, считаем независимыми во всех M столбцах, то для оценки числа K планов всей задачи получаем K=Z(m/n)M.
28 Таблица 2. Примеры фактического числа допустимых планов (Y) и рассчитанного вероятностным способом (K) для задач с матрицами размером N×M и равномерным распределением мощностей ( Ai=A) и емкостей ( Bj=B )
29
Обозначения и параметры транспортной задачи N=2 N=2 N=2 N=2 N=2 N=2 N=2 N=3 N=3 N=3 N=3
M=2 M=2 M=3 M=4 M=2 M=3 M=4 M=3 M=3 M=3 M=4
S=4 S=6 S=6 S=8 S=8 S=12 S=16 S=6 S=9 S=12 S=12
A=2 A=2 A=3 A=4 A=4 A=6 A=8 A=2 A=3 A=4 A=4
B=2 B=3 B=2 B=2 B=4 B=4 B=4 B=2 B=3 B=4 B=3
D=2 D=2 D=2 D=2 D=4 D=4 D=4 D=2 D=3 D=4 D=3
Реальное число допустимых планов Y 3 7 7 19 5 19 85 8 18 80 415
Полное число расчетных (фиктивных) планов K 4, 8 16, 9 7, 4/ 9, 6 9, 8/ 23, 4 16, 9 64, 6/ 365 3302/ 306 57, 4 298 1253 2490/ 5938
30 РЕЗУЛЬТАТЫ
31 Такие модельные расчеты были проведены для 16 примеров равномерных задач (см. табл. 2) и, по-видимому, по указанной выше причине (присутствие неисключенных дублирующих планов) обнаружены, как правило, завышение оценок K над точным числом планов Y, а также несовпадение оценок для задач с транспонированными матрицами. Тем не менее можно построить регрессионную зависимость вида Yapr=1,88K0,62 , прогнозы Yapr числа планов по которой ожидаются коррелирующими с точным числом планов на уровне r = 0, 84 (рис. 2).
32 Рис. 2. Связь точного числа планов Y с оценкой, прогнозируемой по формуле Yapr=1,98K0,58 в логарифмическом (слева) и линейном (справа) масштабах
33 Анализ оценок числа возможных планов решения замкнутых транспортных задач для матриц различной размерности и структуры обнаружил, что при перераспределении между собой только мощностей производителей и постоянстве их суммы число возможных планов задачи монотонно уменьшается с увеличением их относительного среднеквадратического отклонения. Для матриц различной структуры получены аналитические обобщения в простейших ситуациях, а также более универсальный компьютерный алгоритм. Предложен вероятностный путь оценки числа планов равномерных задач. Сформирована обучающая база для дальнейших исследований.

Библиография

1. Ассаул В.Н., Погодин И.Е. (2019). О транспортной задаче с экологическим критерием // Экономика и математические методы. Т. 55. Вып. 1. C. 87–93.

2. Генкин С.А., Итенберг И.В., Фомин Д.В. (1994). Ленинградские математические кружки. Киров: АСА.

3. Канторович Л.В. (1939). Математические методы организации и планирования производст-ва. Ленинград: Изд. ЛГУ.

4. Сергеев А.Н., Петросян Г.А., Тулякова Е.В. (2011). Теория вероятностей и математическая статистика. СПб.: АФ. СПб.

5. Худокормов А.Г. (1994). История экономических учений. Ч. II. Учебник. Под ред. А.Г. Ху-докормова. М.: Изд-во МГУ.

Комментарии

Сообщения не найдены

Написать отзыв
Перевести