close

Вход

Забыли?

вход по аккаунту

?

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

код для вставкиСкачать
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ И СИСТЕМЫ АВТОМАТИЗАЦИИ ПРОЕКТИРОВАНИЯ
Заключение. Таким образом, построив
множество вариантов конфигурации КД4, на
практике можно моделировать инженерные
сети такими кратчайшими деревьями, и построить сеть, удовлетворяющей наперед заданным требованиям. Очевидно, при различных значениях q получит другие топологии.
Для построения кратчайших связывающих
заданное множество точек линий разработана
программа POLSET на языке Турбопаскаль.
Библиографический список
Definition of optimum topology of the
shortest tree for four points of the plane with
polar metrics
K.A. Kuspekov, V.J. Volkov
In article the technique of construction of an
optimum configuration of the shortest tree for four
points of a plane with the polar metrics is
considered. The weight is enclosed to each point
– the factor considering indicators of an
engineering network.
2
1. Куспеков К.А. Минимальное деревья на Е с
полярной метрикой. Материалы 6-Международной
науно-практической конференции, посвященной
125-летию
Национального
технического
университета
«Харьковский
политехнической
институт» и 10-летию Украинской ассоциации по
прикладной геометрии. 21-24 апреля 2009 г. –
Харьков, С.93-97.
2. Есмухан Ж.М., Куспеков К.А. Проблемы
Штейнера и её прикладной алгоритм. Научный
журнал «ПОЙСК» №1, 2006, С.227-231.
Волков Владимир Яковлевич – д-р техн. наук,
профессор, зав. кафедрой «Начертательной геометрии, инженерной и машинной графики» Сибирской государственной автомобильно-дорожной академии. Основное направление научных исследований геометрическое моделирование многокомпонентных
многофакторных процессов. Имеет более 200 опубликованных работ. E-mail: volkov_vy39@mail.ru
Куспеков Кайырбек Амиргазыулы – канд. техн.
наук, доцент, зав. кафедрой «Начертательной геометрии и графики» Казахского национального технического университета. Основное направление
научных исследований - геометрическое моделирование инженерных объектов. Имеет более 66 опубликованных работ.
E-mail: kuspekov_k@mail.ru
УДК 621.87; 681.5
РЕЗУЛЬТАТЫ СРАВНИТЕЛЬНОГО АНАЛИЗА АЛГОРИТМОВ
ПЛАНИРОВАНИЯ ТРАЕКТОРИИ ДВИЖЕНИЯ ОБЪЕКТА С УЧЕТОМ
ЕГО УГЛОВЫХ КООРДИНАТ В ТРЕХМЕРНОМ ПРОСТРАНСТВЕ
С ПРЕПЯТСТВИЯМИ
В.С. Щербаков, М.С. Корытов
Аннотация. Приводятся некоторые результаты сравнительного анализа алгоритмов планирования оптимальной траектории перемещения объекта произвольной
формы с учетом его угловой ориентации в трехмерном пространстве с произвольными препятствиями, заданными в дискретном виде.
Ключевые слова: планирование оптимальной траектории, поиск пути, трехмерное пространство, препятствия, угловые координаты
Введение
Проблема оптимизации траектории перемещения объекта в трехмерном пространстве с
препятствиями является актуальной. В качестве
такого объекта может рассматриваться груз, перемещаемый грузоподъемными машинами, роботами и манипуляторами.
68
Авторами были разработаны алгоритмы оптимизации траектории перемещения объекта с
учетом его угловой ориентации на основе наиболее эффективных современных вычислительных
подходов: роевого интеллекта, генетического,
вероятностной дорожной карты, декомпозиции
линейных и угловых координат, направленного
волнового [1,2,3,4,5,6].
Вестник СибАДИ, выпуск 1 (19), 2011
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ И СИСТЕМЫ АВТОМАТИЗАЦИИ ПРОЕКТИРОВАНИЯ
Задача сравнения алгоритмов имеет практическое применение: необходимо определить,
насколько быстро алгоритмы реализуются при
помощи средств вычислительной техники, какой
объем памяти для них требуется, и насколько
они точны.
Сравнительный анализ алгоритмических
и
программных
реализаций
методик
планирования траектории
Для сравнительной оценки использовалось
3 частных критерия оценки эффективности ал-
горитмов и методик расчета на их основе: Tр –
математическое ожидание значения времени
расчетов, с; Me – информационная длина исследуемого алгоритма (требуемый объем памяти, занимаемой всеми массивами и переменными алгоритма при его реализации),
байт; δ Lусл – математическое ожидание значения относительной погрешности к условному глобальному оптимуму целевой функции,
% (таблица 1).
Таблица 1 – Принятые критерии оценки эффективности алгоритмов и методик
Название критерия
1. Математическое ожидание
значения времени расчетов Tр
2. Информационная длина исследуемого алгоритма Me (требуемый объем памяти, занимаемой всеми массивами и переменными алгоритма при его реализации)
3. Математическое ожидание
значения относительной погрешности к условному глобальному
оптимуму δLусл
Аналитическое
выражение
Tр =
1
⋅
ne
ne
∑ (T )
р i
i =1
me
pe
j =1
k =1
Me = ∑(mm) j + ∑(mp )k
δLусл =
1
⋅
ne
Указанные три частных критерия использовались для сравнительного анализа предложенных методик путем сравнения численных значений критериев при программной реализации соответствующих алгоритмов, т.е.
при применении метода эталонных тестов и
проведении вычислительных экспериментов
[7,8].
Было проведено несколько серий вычислительных экспериментов со случайным расположением препятствий в пределах области
рассматриваемых перемещений объекта. Положение объекта в пространстве описывалось
3 линейными координатами и 2 углами поворота. Данное сочетание координат, рассмотренное в качестве примера, описывает довольно распространенный частный случай
положения объекта в форме тела вращения
(например, цилиндра).
Наиболее часто при работе грузоподъемных машин на строительных площадках, в цехах и складских помещениях встречаются
препятствия в виде стен, перекрытий, стеллажей, контейнеров, столбов и свай, которые
Вестник СибАДИ, выпуск 1 (19), 2011
ne
∑ (δL )
усл i
i =1
Параметры
Tр – время расчетов для отдельного
эксперимента, с; ne – число независимых экспериментов
mm – объем памяти, занимаемый отдельным массивом данных, байт; me
– количество массивов данных алгоритма; mp – объем памяти, занимаемый отдельной переменной, байт; pe
– количество переменных алгоритма
δLусл – относительная погрешность к
условному глобальному оптимуму
для отдельного эксперимента, %
могут быть описаны в форме параллелепипедов различных размеров [9].
Методом эталонных тестов производилось
сравнение разработанных методик поиска оптимальной траектории объекта на основе: 1 –
направленного волнового алгоритма; 2 – алгоритма роевого интеллекта; 3 – генетического подхода; 4 – алгоритма декомпозиции линейных и угловых координат; 5 – алгоритма
вероятностной дорожной карты; 6 – распараллеленного алгоритма роевого интеллекта;
7 – распараллеленного алгоритма декомпозиции линейных и угловых координат. Для этого
были проведены несколько серий вычислительных экспериментов, моделирующих процесс поиска оптимальной по значению целевой функции траектории перемещения объекта в среде с полидистантными поверхностями, построенными вокруг реальных поверхностей препятствий [10].
В качестве целевой функции использовалось среднее взвешенное длин линейных и
угловых перемещений, выражение которого
для отдельной траектории перемещения имеет вид:
69
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ И СИСТЕМЫ АВТОМАТИЗАЦИИ ПРОЕКТИРОВАНИЯ
L=
imax
∑  (x
− x i −1 ) + (y i − y i −1 ) + (zi − zi −1 ) + cγω ⋅
2
i
2
2
(γ i − γ i −1 )2 + (ωi − ωi −1 )2 ,

i =2
где imax – число точек отдельной дискретной
траектории (imax=20); cγω – весовой коэффициент угловых координат (cγω=1); x, y, z, γ, ω –
соответственно три линейные и две угловые
координаты, характеризующие положение
объекта.
а)
Исходные поверхности препятствий формировались случайным образом из сочетания
нескольких параллелепипедов, каждый из которых имел случайные размеры. Для всех
экспериментов рассматривалось рабочее
пространство в виде куба с размерами
10×10×10 условных линейных единиц (УЛЕ)
(рисунок 1).
б)
в)
Рис. 1. Рабочее пространство размером 10×10×10 УЛЕ с детерминированным (а)
и стохастическим (б, в) расположением препятствий
Линейные координаты начальной и конечной точек положения условного центра груза
(начала локальной системы координат груза
OgХgYgZg) принимались постоянными для всех
экспериментов. Начальная точка имела линейные координаты (УЛЕ): [xн0; yн0; zн0]=[0; 2;
5], конечная – [xк0; yк0; zк0]=[10; 2; 5]. Начальные и конечные значения угловых координат
также принимали постоянные нулевые значения (рад): [γн0; ωн0]=[0; 0], [γк0; ωк0]=[0; 0].
Количество разбиений каждой траектории
на кусочно-линейные участки вдоль оси Х0
неподвижной системы координат принималось равным: imax=20. Соответственно, максимальные значения индексов остальных координат груза при поиске траектории принимались равными: jmax=20; kmax=20; lmax=5; mmax=17
(индексы i, j, k, l, m соответствовали координатам x, y, z, γ, ω).
Сочетание линейных и угловых координат
груза позволяет представить область перемещений груза в виде гиперкуба с размерами
равномерной
сетки
индексов
(РСИ):
20×20×20×5×17. То есть, полный граф состояний объекта при поиске траектории состоял из
680000 вершин.
В то же время для задания исходной поверхности препятствий, построения полидистантных поверхностей вокруг препятствий и
для проведения дискретной локальной оптимизации найденных траекторий, использовалось более подробное описание рабочей области в трехмерном пространстве, в виде куба
70
с размерами (РСИ): 100×100×100. То есть,
использовался иерархический подход, при
котором задание исходной поверхности и локальная оптимизация траектории проводились
с более мелким шагом описания препятствий
(0,1 УЛЕ), а собственно поиск траектории с
использованием разработанных методик – с
более крупным шагом описания препятствий
(0,5 УЛЕ). Это позволило более точно провести локальную оптимизацию.
Применение иерархического подхода вызвано большой размерностью задачи, сложность которой экспоненциально возрастает
при уменьшении шага дискретизации координат. Переход от исходного к укрупненному
описанию линейных координат выполнялся
методом редукции. Использовался коэффициент масштабирования kм линейных координат, значение которого принималось равным:
kм=5.
Для методик, использующих элементы вероятностного выбора (на основе алгоритмов
№ 2, 3 и 5), предварительно были проведены
несколько серий экспериментов, по 1200 независимых экспериментов для каждого сочетания значений варьируемых параметров, на
рабочей области с одинаковым, детерминированным расположением препятствий, задаваемых высотами поверхности препятствий
Yпр(i,k) и описываемых следующими условиями на РСИ 100×100 (см. рисунок 1,а):
пр(i,k)=2,5
УЛЕ – при (29≤i≤31) и (20≤k≤80);
Вестник СибАДИ, выпуск 1 (19), 2011
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ И СИСТЕМЫ АВТОМАТИЗАЦИИ ПРОЕКТИРОВАНИЯ
Yпр(i,k)=8,0 УЛЕ – при (49≤i≤51) и (50≤k≤100);
Yпр(i,k)=8,0 УЛЕ – при (79≤i≤81) и (1≤k≤50);
Yпр(i,k)=0,0 УЛЕ – в остальных случаях.
Груз в форме цилиндра с габаритными
размерами: габаритный диаметр 0,5 УЛЕ, и
высота 2,0 УЛЕ, был представлен в виде набора точек (cг=12) на поверхности объемного теv
ла с координатами в УЛЕ: { Rig }={[0,25;0;1;1];
[0,25;0;0;1];
[0,25;0;–1;1];
[–0,25;0;1;1];
[–
0,25;0;0;1];
[–0,25;0;–1;1];
[0;0,25;1;1];
[0;0,25;0;1]; [0;0,25;–1;1]; [0;–0,25;1;1]; [0;–
0,25;0;1]; [0;–0,25;–1;1]}.
Для каждой траектории объекта, найденной
в отдельном независимом эксперименте (рисунок 2),определялись собственные значения
целевой функции (L до локальной оптимизации, L* после локальной оптимизации), информационной длины исследуемого алгоритма Me
и времени поиска Tр.
Проведение вычислительных экспериментов на тестовом примере с детерминированным расположением препятствий позволило
обосновать рациональные значения основных
внутренних параметров соответствующих алгоритмов с элементами вероятностного выбора
для решения поставленной задачи и для последующего сравнения с детерминированными
алгоритмами. Ввиду ограниченного объема
статьи результаты вычислительных экспериментов на тестовом примере с детерминированным расположением препятствий не приводятся.
Была проведена серия из 10000 компьютерных экспериментов, моделирующих процесс поиска оптимальной траектории перемещения объекта в среде с полидистантными
поверхностями, построенными вокруг реальных поверхностей препятствий [Yпр], сформированных случайным образом из сочетания
нескольких параллелепипедов, каждый из которых имел случайные размеры.
Число параллелепипедов n в каждом эксперименте генерировалось по равномерному
закону распределения случайной величины в
интервале от 1 до 15. Размеры каждого параллелепипеда формировались в пределах
(x×y×z) от 0×0×0 УЛЕ до 4×5×4 УЛЕ также по
равномерному закону распределения. Допускалось перекрытие объемов и поверхностей
параллелепипедов при их наложении.
Для каждого эксперимента определялись путем непосредственных вычислительных измерений, реализованных программно, значения Tр и
Вестник СибАДИ, выпуск 1 (19), 2011
а)
б)
Рис. 2. Примеры найденной алгоритмом на основе
генетического подхода (№ 3)
траектории: а) до локальной оптимизации (L); б)
после локальной оптимизации (L*)
Me, и рассчитывалось по результатам вычислительных измерений значение δLусл:
*
δLусл= (L /Lусл)·100%,где Lусл – условный глобальный оптимум решения задачи, определяемый как
минимальное значение из множества значений
целевой функции, полученных при использовании различных методик, но для одних и тех же
численных значений исходных данных задачи
[Yпр]:
*
Lусл=min{(L )i}, i ∈ [1; 5],где i – номер применяемой методики на основе: 1 – направленного
волнового алгоритма; 2 – алгоритма роевого
интеллекта; 3 – генетического подхода; 4 – алгоритма декомпозиции линейных и угловых
координат; 5 – алгоритма вероятностной дорожной карты.
Некоторые результаты сравнительного
анализа методик и алгоритмов по принятым
критериям оценки эффективности приведены
на рисунке 3. Исследования проводились на
программных реализациях методик и алгорит-
71
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ И СИСТЕМЫ АВТОМАТИЗАЦИИ ПРОЕКТИРОВАНИЯ
мов в среде MS Visual C++ на ПК средней производительности (AMD Athlon 64 X2 Dual Core
Processor 5600+ 2.90 GHz, RAM 4,00 ГБ).
Поскольку использовалось 3 самостоятельных критерия оценки алгоритмов и методик, ре-
Tр ,с
ниям каждого критерия. Векторное представление разработанных алгоритмов и методик
по принятым критериям эффективности как
компонентам векторов, а также проекции векторов на три плоскости ([ Tр Me], [ Tр δ Lусл ], [Me
24,897
23,196
25
зультаты анализа были интерпретированы графически в удобном для восприятия и наглядном виде, как совокупность векторов в трехмерном пространстве координат [ Tр Me δ Lусл ].
20
15
δ Lусл ]) приведены на рисунке 4.
10
5
0
2,611
1
3,31
2
0,782
3
Tр ,с
5
δ Lусл ,
2
5
4
3
20
1
0
2,611
3,31
1
2
0,782
3
Tр ,с
4
5
24,897
23,196
25
2,93 1,782
20
15
10
5
0
2,611
3,31
1
2
0,782
3
2,93 1,782
4
5
Рис. 3. Результаты сравнительного анализа методик и алгоритмов по принятым критериям оценки
эффективности: 1 – направленный волновой алгоритм; 2 – алгоритм роевого интеллекта; 3 – алгоритм на основе генетического подхода; 4 – алгоритм декомпозиции линейных и угловых координат;
5 – алгоритм вероятностной дорожной карты;
–
математическое ожидание значения времени расчетов распараллеленных алгоритмов 2 и 4, верхний пределДля этого было проведено масштаби-
3
1
2
1
10
5
25
3
15
5
4
4
24,897
23,196
25
4
δ Lусл ,
2
2,93 1,782
Tр
δ Lусл
1
5
Me,
Мб
3
Me
Tр
Рис. 4. Векторное представление
критериев эффективности
разработанных алгоритмов и методик
При сравнении алгоритмов по принятым
критериям оценки эффективности, в трехмерном пространстве положений векторного
критерия эффективности могут быть выделены отдельные области, соответствующие оптимальному значению какого-либо одного, двух
или всех трех компонент вектора, определяемые условиями вида:
Tр ≤( Tр )max;
Me≤(Me)max;
δ Lусл ≤( δ Lусл )max,где
( Tр )max; (Me)max; ( δ Lусл )max – максимальные предельно допустимые (заданные) значения критериев.
Соответственно можно выделить область
условно быстрых алгоритмов, область
условно точных алгоритмов и область условно компактных (т.е. занимающих малый объем
в памяти ПК) алгоритмов (рисунок 5, а, б, в).
рование отображаемых результатов экспериментов по максимальным полученным значе-
72
Вестник СибАДИ, выпуск 1 (19), 2011
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ И СИСТЕМЫ АВТОМАТИЗАЦИИ ПРОЕКТИРОВАНИЯ
Условно быстрые алгоритмы
δ Lусл , %
δ Lусл , %
2
Условно точные алгоритмы
5
5
4
4
3
1
2
3
Tр ,с
Me, Мб
Tр ,с
1
Me, Мб
б)
а)
δ Lусл , %
δ Lусл , %
5
2
Условно
компактные
алгоритмы
2
5
3
3
1
Tр ,с
1
Me, Мб
Me, Мб
4
Tр ,с
в)
4
г)
Рис. 5. Области условно быстрых (а), условно точных (б), условно компактных (в) алгоритмов и их пересе– последовательные алгоритмы, попадающие в соответствующую
чения (г):
область оптимальности
В качестве примера, при ( Tр )max=5 с;
(Me)max=2 Мб; ( δ Lусл )max=5% к условно быстрым
алгоритмам могут быть отнесены 1, 3 и 5 алгоритмы, к условно точным 1 и 4 алгоритмы, к
условно компактным 2, 3 и 5 алгоритмы. В случае распараллеливания алгоритмов 2 и 4, они
также попадают в область условно быстрых.
Заключение
Ни один из разработанных алгоритмов не
удовлетворяет условию оптимальности по
всем трем критериям одновременно, т.е. не
входит в область пересечения условно быстрых, условно точных и условно компактных алгоритмов (см. рисунок 5, г), однако отдельные
алгоритмы оптимальны по двум критериям. В
частности, алгоритмы 3 и 5 одновременно
быстрые и компактные, причем алгоритм 3 в
большей степени. Алгоритмы 1 и распараллеленный 7 одновременно быстрые и точные,
причем алгоритм 1 в большей степени.
Результаты сравнения алгоритмов подтверждают предположение о том, что существует взаимосвязь между оптимальностью
алгоритма, которая в данном случае выража-
Вестник СибАДИ, выпуск 1 (19), 2011
ется его точностью δ Lусл , и сложностью алгоритма (временной или пространственной), которая в данном случае выражается через критерии Tр и Me. Т.е. невозможно упростить алгоритм (временную либо пространственную
сложность, которые также взаимосвязаны
между собой), не жертвуя при этом его оптимальностью, что хорошо иллюстрируется алгоритмом на основе генетического подхода (№
3).
Учитывая, что при современном уровне
развития компьютерной техники требование
компактности не является критичным в диапазоне значений всех рассматриваемых алгоритмов (не более 6 Мб), следует отметить
направленный волновой алгоритм (№ 1) как
наиболее точный и в то же время достаточно
быстрый. В пользу его перспективности для
решения поставленной задачи говорит и тот
факт, что это детерминированный алгоритм в
отличие от всех остальных, кроме алгоритма
№ 7. Этот алгоритм находит единственно возможное и постоянное решение задачи при од-
73
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ И СИСТЕМЫ АВТОМАТИЗАЦИИ ПРОЕКТИРОВАНИЯ
них и тех же численных значениях исходных
данных.
Библиографический список
1. Щербаков, В.С. Использование алгоритмов
поиска пути перемещения груза автокраном на графах / В.С. Щербаков, М.С. Корытов // Вестник Воронежского государственного технического университета, 2009. – Т.5. – № 5. – С. 37-41.
2. Щербаков, В.С. Поиск оптимальной траектории груза, перемещаемого автокраном, в среде с
произвольными препятствиями, с учетом координат
угловой ориентации в трехмерном пространстве /
В.С. Щербаков, М.С. Корытов // Вестник Брянского
государственного технического университета, 2009. –
№ 4 (24). – С. 48-51.
3. Щербаков, В.С. Методика поиска субоптимальной траектории движения объекта в трехмерной
среде с произвольными препятствиями с учетом координат угловой ориентации / В.С. Щербаков, М.С.
Корытов // Вестник СибАДИ, 2009. – Вып. 4 (14). – С.
5-10.
4. Корытов, М.С. Декомпозиция обобщенных координат при решении задач оптимизации траектории
перемещения груза // Вестник МАДИ (ГТУ), 2010. –
Вып. 3(22). – С. 32-35.
5. Щербаков, В.С. Использование генетических
алгоритмов для поиска оптимальной траектории перемещения груза / В.С. Щербаков, М.С. Корытов //
Вестник КГТУ им. А.Н. Туполева, 2010. – № 3. – С.
155-158.
6. Щербаков, В.С. Об одной модификации алгоритма муравьиных колоний для планирования траектории перемещения груза в пространстве с препятствиями с учетом угловой ориентации / В.С. Щербаков, М.С. Корытов // Известия высших учебных заведений. Поволжский регион. Технические науки. –
2010. – № 3. – С. 143-151.
7. Кормен, Томас X. Алгоритмы: построение и
анализ: пер. с англ. / Томас X. Кормен, Чарльз И. Лейзерсон, Рональд Л. Ривест, Клиффорд Штайн. – М.:
Изд. дом «Вильямс», 2005. – 1296 с.
8. Рассел, С. Искусственный интеллект: современный подход: пер. с англ. / Стюарт Рассел, Питер
Норвиг. – М.: Изд. дом «Вильямс», 2006. – 1408 с.
74
9. Щербаков, В.С. Влияние стохастических параметров пространства с препятствиями на длину
траектории груза, перемещаемого грузоподъемным
краном / В.С. Щербаков, М.С. Корытов // Вестник СибАДИ, 2009. – Вып. 3 (13). – С. 13-17.
10. Корытов, М.С. Использование полидистантных поверхностей в задаче поиска пути перемещения
груза в среде с препятствиями // Материалы 64-й
научно-технической конференции ГОУ «СибАДИ». –
Омск: СибАДИ, 2010. – Кн. 1. – С. 302-306.
A comparative analysis of scheduling algorithms trajectory of an object given its angular coordinates in three-dimensional space
with obstacles
V.S. Shcherbakov, M.S. Korytov
Some results of comparative analysis of algorithms for planning optimal trajectories move an
object of arbitrary shape, given its angular orientation in three-dimensional space with arbitrary
obstacles, given in a discrete form.
Щербаков Виталий Сергеевич – доктор техн.
наук, профессор, декан факультета «Нефтегазовая и строительная техника» Сибирской государственной автомобильно-дорожной академии.
Основное направление научных исследований –
совершенствование систем управления строительных и дорожных машин, общее количество
публикаций – более 200, адрес электронной почты – sherbakov_vs@sibadi.org.
Корытов Михаил Сергеевич – канд. техн.
наук, доцент, докторант Сибирской государственной автомобильно-дорожной академии. Основное направление научных исследований – автоматизация планирования рабочих процессов
мобильных грузоподъемных машин, общее количество публикаций – более 80, адрес электронной
почты – kms142@mail.ru.
Вестник СибАДИ, выпуск 1 (19), 2011
1/--страниц
Пожаловаться на содержимое документа