Задача состоит в том, чтобы перечислить все точки, принадлежащие границе выпуклой оболочки заданного множества точек, в порядке ее обхода, например, против часовой стрелки (в некоторых задачах требуется перечислить только угловые точки). Для эффективного решения этой задачи существует несколько различных алгоритмов (см. [1]-[4]). Приведем наиболее простую реализацию одного из них — алгоритма Джарвиса.
{ следующий абзац проиллюстрировать рисунком из номера 16/2000, стр. 11 }
Перечисление точек искомой границы выпуклого многоугольника начнем с правой нижней точки, которая заведомо принадлежит границе выпуклой оболочки. Обозначим ее координаты (x 0, y 0). Следующей, при указанном порядке обхода, будет точка (x 1, y 1), для которой угол между осью OX и вектором (x 0, y 0)‑(x 1, y 1) минимален. Если таких точек несколько, то угловой в многоугольнике станет точка, для которой длина вектора (x 0, y 0)‑(x 1, y 1) максимальна, а следующей точкой, принадлежащей выпуклой оболочке — та, длина вектора у которой минимальна (таким образом наш выбор будет зависеть от конкретной постановки задачи). Для нахождения следующей точки значения углов между векторами вычислять необязательно. Мы опять можем воспользоваться понятием знака векторного произведения. Пусть, далее, мы уже нашли i ‑ю вершину выпуклой оболочки (xi, yi). Тогда, (i + 1)-я точка есть такая точка, еще не вошедшая в выпуклую оболочку, для которой угол между вектором (xi ‑1, yi ‑1)‑(xi, yi) и вектором (xi, yi)‑(xi +1, yi +1) минимален, при минимальной длине вектора (xi, yi)‑(xi +1, yi +1) среди всех векторов с таким углом. Заметим, что для всех остальных точек (x, y), вектор (xi, yi)‑(x, y) будет лежать вне угла, образованного указанными векторами, левее него. Тогда векторное произведение (xi +1 – xi)(y – yi) – (yi +1 – yi)(x – xi) ³ 0, для любой точки (x, y), еще не вошедшей в границу выпуклой оболочки. Следовательно, мы можем сначала считать следующей, (i + 1)‑ой, любую, еще не вошедшую в выпуклую оболочку, точку, а затем, вычисляем указанное выражение для остальных “свободных” точек (х, y). Если для одной из них (xi +1 – xi)(y – yi) – (yi +1 – yi)(x – xi) < 0, считаем следующей ее и продолжаем проверку остальных точек (аналогично алгоритму поиска минимального элемента в массиве). Если же значение выражения равно 0, то сравниваем квадраты длин векторов, а именно (xi +1 – xi)2 + (yi +1 – yi)2 и (x – xi)2 + (y – yi)2.
Таким образом, при решении данной задачи в случае изначально целочисленных координат мы полностью можем избежать применения вещественной арифметики, а, следовательно, избавиться от потери точности вычислений. В противном случае, в решение могут быть включены “лишние” точки, близкие к границе выпуклой оболочки или не учтены некоторые из “нужных” точек. Сложность данного алгоритма составит O (kN), где k — количество точек в выпуклой оболочке, в худшем случае равное N. Существуют алгоритмы решения этой задачи, основанные на предварительной сортировке точек исходного множества, например, по значению угла в полярной системе координат с центром в одной из точек выпуклой оболочки, с вычислительной сложностью O (N log N) (алгоритм Грехема). То есть наиболее трудоемкой задачей оказывается именно сортировка исходных точек.
Приведем программу решения данной задачи алгоритмом Джарвиса:
var a, b: array [1..100] of record
x,y: integer;
f: boolean
end;
min, m, j, k, n: integer;
Begin
readln(n);
for i:=1 to n do
Begin
read(a[i].x, a[i].y);
a[i].f:= false
end;
{ищем нижнюю правую точку}
m:=1;
for i:=2 to n do
if a[i].y < a[m].y then m:=i else
if (a[i].y = a[m].y) and
(a[i].x > a[m].x) then m:=i;
b[1]:=a[m];
a[m].f:= true;
k:=1;
Repeat
min:=1;
{ищем первую еще свободную точку}
while a[min].f do inc(min);
{ищем очередную вершину выпуклой оболочки}
for j:=1 to n do
if (not a[j].f) and
((a[min].x-b[k].x)*(a[j].y-b[k].y)-
(a[j].x-b[k].x)*(a[min].y-b[k].y)<0)
then min:=j else
if (not a[j].f) and
((a[min].x-b[k].x)*(a[j].y-b[k].y)-
(a[j].x-b[k].x)*(a[min].y-b[k].y)=0) and
(sqr(a[min].x-b[k].x)+sqr(a[min].y-b[k].y) >
sqr(a[j].x-b[k].x)+sqr(a[j].y-b[k].y))
then min:=j;
k:=k+1;
a[min].f:= true;
b[k]:=a[min] {записана очередная вершина}
until min=m; {пока ломаная не замкнется}
for j:=1 to k do {печать результата}
writeln(b[j].x,' ',b[j].y);
end.
Приведем примеры задач, при решении которых используется построение выпуклой оболочки.
Задача 1. Стена. (В англоязычном варианте задача предлагалась на студенческих командных соревнованиях по программированию в Санкт-Петербурге в 2001 г.)
Однажды жадный король приказал своему архитектору построить стену вокруг дворца. Король был настолько жадный, что не стал слушать предложения архитектора о построении стены совершенной формы. Вместо этого король приказал потратить на строительство стены определенной высоты и произвольной формы (не обязательно в виде ломаной!!!) как можно меньше кирпичей, но потребовал, чтобы стена отстояла от дворца не меньше, чем на L футов. В случае невыполнения условия или перерасхода средств архитектор может лишиться головы.
Ваша задача помочь бедному архитектору. Ваша задача написать программу, которая определит минимально возможную длину стены, которой можно окружить дворец и при этом выполнить все требования короля.
Первая строка входного файла содержит 2 числа N (3 £ N £ 1000) — количество углов в здании дворца и L (1 £ L £ 1000) минимальное расстояние на которое стена может приближаться ко дворцу. Следующие N строк описывают координаты на поверхности земли углов дворца, в порядке их обхода по часовой стрелке. Каждая строка содержит два целых числа — Xi и Yi, разделенных пробелом (-10000 £ Xi, Yi £ 10000), которые описываю координаты углов дворца в футах. Все углы дворца различны, а стороны не имеют других общих точек, кроме угловых.
Запишите в выходной файл одно число, определяющее минимальную длину дворца в футах, удовлетворяющую условию задачи. Ответ должен быть записан в виде целого числа, так как с вещественными числами король не знаком, однако округлять его следует так, чтобы целое число футов отличалось от настоящего ответа менее, чем на 8 дюймов (в 1 футе 12 дюймов).
| Пример входного файла | Выходной файл |
| 9 100 200 400 300 400 300 300 400 300 400 400 500 400 500 200 350 200 200 200 |
Решение. Ответом на данную задачу будет длина выпуклой оболочки, увеличенная на длину окружности с радиусом L и округленная до ближайшего целого.
Задача 2. На плоскости заданы N точек. Построить замкнутую ломаную без самопересечений и самокасаний, вершинами которой должны стать все заданные точки. (См., например, [5], аналогичная задача предлагалась на кировской областной олимпиаде 2002г.).
Решение. Следующий рисунок проиллюстрирует идею одного из возможных способов решения данной задачи:
