Ориентированная площадь как найти

Ориентированные площади и объёмы

Ориентированная площадь параллелограмма

Ориентированной площадью параллелограмма, построенного на неколлинеарных векторах и , называется его площадь , взятая со знаком плюс, если ориентация пары векторов и правая , и со знаком минус, если ориентация — левая

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

1) ;

2) ;

3) .

4) Если векторы и в правом ортонормированием базисе имеют координаты и соответственно, то внешнее произведение этих векторов находится по формуле

Если — координатные столбцы векторов в стандартном базисе на плоскости, то их внешнее произведение находится по формуле

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

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

при (рис.1.49,а)

при (рис.1.49,b)

так как пара в первом случае правая (рис.1.49,а), а во втором случае — левая (рис.1.49,б).

Пусть теперь векторы и не коллинеарны (рис.1.49,в). Построим проекции и на прямые, содержащие базисные векторы: . Из концов векторов и опустим перпендикуляры и соответственно на прямую, содержащую вектор . Учитывая, что векторы и противоположно направлены, а также подобие прямоугольных треугольников с гипотенузами и , находим абсциссу вектора :

так как пара — правая, а пара — левая. Аналогично находится ордината (векторы и одинаково направлены)

Таким образом, вектор имеет следующее разложение по базису на плоскости:

Рассмотрим применение формулы (1.19) для решения системы двух линейных уравнений с двумя неизвестными:

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

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

что совпадает с правилом Крамера.

---

Ориентированный объем параллелепипеда

Ориентированным объемом параллелепипеда, построенного на некомпланарных векторах , называется его объем , взятый со знаком плюс, если ориентация тройки векторов правая и со знаком минус, если ориентация — левая .

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

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

так как определитель транспонированной матрицы равен определителю исходной матрицы.

При помощи ориентированных объемов задача разложения вектора по базису в пространстве решается так же как и на плоскости с использованием ориентированных площадей. В результате получаем

где

Формула (1.21) соответствует правилу Крамера решения системы трех линейных уравнений с тремя неизвестными.

---

Пример 1.23. Заданы координатные столбцы

векторов в стандартном базисе. Разложить вектор по векторам .

Решение. По формуле (1.20) находим смешанные произведения

Коэффициенты разложения определяем по формуле (1.21):

Следовательно, .

Математический форум (помощь с решением задач, обсуждение вопросов по математике).

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

Определение. Множество концов
векторов вида,
гдеипри
условии, что начала векторовисовпадают с концом вектора,
называетсяпараллелограммом, натянутым
на вектораи.

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

Определение. Ориентированная
площадь параллелограмма, натянутого
на вектораиесть
число (как положительное так и
отрицательное), причём модуль (абсолютная
величина этого числа) есть обычная
площадь параллелограмма. Знак этого
числа определяется сравнением ориентации
пары векторовис
ориентацией пары базисных векторовиплоскости векторови(для
простоты можно считать, чтоиобразуют декартов базис, то есть эти
вектора имеют единичную длину и взаимно
ортогональны (перпендикулярны)). Если
ориентация парыисовпадает с ориентацией парыи,
то ориентированная площадь положительна,
если не совпадают, то ориентированная
площадь отрицательна.

Аналогичноопределяется ориентированный
объём параллелепипеда, натянутого на
тройку векторов,и:
это есть число, модуль которого – обычный
объём этого параллелепипеда, а знак
числа определяется по совпадению
ориентаций тройки векторов,ис ориентацией тройки базисных векторов
пространства,и.

Определение. Грассманово произведениедвух векторов есть операция, обладающая
следующими свойствами.

1. Эта операция линейна относительно
   своих сомножителей:
   .

2. Эта
   операция антикоммутативна относительно
   своих сомножителей:

.
Следствие антикоммутативности:.

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

Далее,
,
поскольку ориентация парыпротивоположна ориентации пары.

Рассмотрим грассманово произведение
двух векторов, разложенных по векторам
идекартова базиса:и.

Видно, что:

.

В данном случае использованы свойства
линейности и антикоммутативности
грассманова произведения и его следствие:

Отсюда видно, что
,
где— определитель второго порядка.
Рассматриваякак именованную ориентированную площадь,
а величинукак некоторую единицу измерения этой
площади, видим, что определитель как
отвлечённое число есть отношение
именованных площадей:.

Задача. Проверить формулу вычисления
определителя второго порядка
непосредственно на декартовой плоскости.

Решение вытекает из рисунка.

Именно:
.

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

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

Время на прочтение
8 мин

Количество просмотров 125K

Здравствуйте, уважаемые хабравчане! Это моя вторая статья, и мне хотелось бы поговорить о вычислительной геометрии.

Немного истории

Я являюсь студентом уже 4 курса математического факультета, и до того как я начал заниматься программированием, я считал себя математиком на 100 процентов.

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

Мне очень нравится подход моего преподавателя: «разберись с этой темой, а потом расскажи нам, да так чтоб мы все поняли».

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

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

Вступление

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

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

Немного теории о векторах

Отрезок, для которого указано, какой из его концов считается началом, а какой — концом, называется вектором. Любая точка пространства также может рассматриваться как вектор. Такой вектор называется нулевым. Начало и конец нулевого вектора совпадают, и он не имеет какого-либо определенного направления.

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

Скалярное произведение векторов

Скалярное произведение векторов — это число, равное произведению длин этих векторов на косинус угла между ними.
(a, b) = |a||b|cos∠(a, b)

Если векторы заданы своими координатами a(x₁, y₁), b(x₂, y₂) то скалярное произведение (a, b) = x₁x₂ + y₁y₂.

Косое произведение векторов

Псевдоскалярным или косым произведением векторов на плоскости называется число
[a, b] = |a||b|sinθ
где — угол вращения (против часовой стрелки) от a к b. Если хотя бы один из векторов a и b нулевой, то полагают [a, b] = 0.
Если векторы заданы своими координатами a(x₁, y₁), b(x₂, y₂) то косое произведение [a, b] = x₁y₂ — x₂y₁.
Геометрически косое произведение векторов представляет собой ориентированную площадь параллелограмма, натянутого на эти вектора.

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

А теперь займемся практикой

Начнем с треугольников

Задача №1

Задача очень простая, а именно: по введенным трем числам a, b, c определить существует ли треугольник с такими сторонами.

Решение
Понятно, что здесь нужно только проверить неравенство треугольника: a + b > c, a + c > b, b + c > a. Интересно, при изучении неравенства треугольника только ли у меня возник вопрос: не могут ли отрицательные числа тоже удовлетворять этим трем неравенствам? Оказывается, нет! Если мы сложим каждое неравенство, то получим a > 0, b > 0, c > 0. Поэтому неравенство треугольника является необходимым и достаточным условием существования треугольника.

Задача №2

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

Решение
С первого взгляда решение кажется очевидным: вычислить стороны треугольника и свести задачу к предыдущей. Однако поскольку расстояние между двумя точками A(x₁, y₁), B(x₂, y₂) вычисляется по формуле √(x₁-x₂)²+(y₁-y₂)² то при извлечении корня возможна потеря точности, что плохо скажется на проверке неравенства треугольника. Оказывается, что если треугольник задан координатами своих вершин, то вычислять длины его сторон и проверять неравенство треугольника не требуется. В этом случае треугольника не существует тогда и только тогда, когда данные три точки лежат на одной прямой. А это легко проверяется через косое произведение векторов. Если оно равно нулю, то векторы коллинеарные, то есть все три точки лежат на одной прямой.

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

Задача №3

Треугольник задан своими сторонами. Определить тип треугольника: тупоугольный, прямоугольный или остроугольный.

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

Из курса геометрии известно, что напротив большей стороны лежит больший угол (он нам и нужен). Поэтому если мы выясним чему равен больший угол, то поймем тип треугольника:

1. Угол больше 90° – треугольник тупоугольный
2. Угол меньше 90°– треугольник остроугольный
3. Угол равен 90°– треугольник прямоугольный

Воспользуемся теоремой косинусов:

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

• Если cosα > 0, то a² < b² + c² – треугольник остроугольный
• Если cosα = 0, то a² = b² + c² – треугольник прямоугольный
• Если cosα < 0, то a² > b² + c² – треугольник тупоугольный

где a – большая сторона.

Задача №4

Задача аналогична предыдущей задаче, только треугольник задан не своими сторонами, а координатами вершин.

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

Задача №5

По данным сторонам треугольника найти его площадь.

Решение
Очевидно решение, заключается в применение формулы Герона.

Кстати, никого не интересовало доказательство этой формулы?

Задача №6

Вычислить площадь треугольника заданного координатами своих вершин.

Решение
Не будем говорить о решении, которое сводится к предыдущей задачи, а попробуем воспользоваться геометрическим смыслом косового произведения. Геометрически косое произведение двух векторов определяет ориентированную площадь параллелограмма натянутого на эти вектора. Поскольку диагональ параллелограмма разбивает его на два равновеликих треугольника, то можем найти площадь нашего треугольника, как половину площади параллелограмма.
Для векторов a(x₁, y₁), b(x₂, y₂)

S = (x₁y₂ — x₂y₁) / 2 — ориентированная площадь треугольника

Задача №7

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

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

Метод площадей

Если сумма площадей треугольников AKB, AKC, BKC (не ориентированных, а «обычных») больше площади треугольника ABC точка лежит вне треугольника. Если же сумма первых трех площадей равна четвертой, то нужно проверить, не равна ли нулю одна из трех площадей. Если равна, то точка лежит на границе треугольника, иначе – внутри.
Вычислять площади треугольников, естественно, надо через косое произведение векторов. Этот метод не очень хороший. Поскольку здесь используются сравнение чисел с плавающей точкой, а это в свою очередь может привести к принятию неверного решения при сравнении. Второй метод опять таки опирается на вектора, он намного эффективнее во всех отношениях.

Проверка полуплоскостей

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

В первом примере сторона AB разводит вершину C и точку K по разным полуплоскостям, поэтому точка лежит снаружи.

Задача №8

Вычисление площади многоугольника заданного координатами своих вершин.

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

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

Метод трапеций

Для того чтобы посчитать площадь многоугольника нужно разбить его на трапеции, так как это показано на рисунке, а затем сложить ориентированные площади полученных трапеций это будет ориентированной площадью исходного многоугольника.
S = S_{A1 A2 B2 B1} + S_{A2 A3 B3 B2} + S_{A3 A4 B5 B3} + S_{A4 A5 B6 B5} + S_{A5 A6 B4 B6} + S_{A6 A1 B1 B4}
Площади трапеций считаем по известной формуле: полусумма оснований на высоту
S_{A1 A2 B2 B1} = 0.5 * (A₁B₁ + A₂B₂) *(B₂ — B₁)

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

Метод треугольников

Аналогично предыдущему методу можно разбивать многоугольник не на трапеции, а на треугольники, как показано на рисунке. В результате, сложив ориентированные площади этих треугольников, мы получим опять-таки ориентированную площадь многоугольника.
S = S_OA1A2 + S_OA2A3 + S_OA3A4 + S_OA4A5 + S_OA5A6 + S_OA6A1

Как вы видите задача вычисления площади многоугольника достаточна проста. Не знаю, почему, но мне больше нравится решать эту задачу методом разбиения на трапеции (наверно потому, что на всех олимпиадах я ее так решал). Тем более, что при втором решении площади треугольников надо вычислять через косое произведение. О формуле Герона надо забыть!!!

Задача №9

Многоугольник задан координатами своих вершин в порядке его обхода. Необходимо проверить является ли многоугольник выпуклым.

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

Задача опять сводится к вычислению косового произведения векторов, а именно у выпуклого многоугольника знаки косых произведений [A_i A_i+1, A_i+1 A_i+2] либо положительны, либо отрицательны. Поэтому если мы знаем направление обхода, то знак косых произведений для выпуклого многоугольника одинаков: он неотрицателен при обходе против часовой стрелки и неположителен при обходе по часовой стрелки.

Задача №10

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

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

Для любого многоугольника с целочисленными координатами вершин справедлива формула Пика: S = n + m/2 — 1, где S – площадь многоугольника, n – количество целых точек лежащих строго внутри многоугольника, m – количество целых точек лежащих на границе многоугольника. Поскольку площадь многоугольника мы знаем как вычислять, то S известно. Так же мы можем вычислить количество целых точек лежащих на границе многоугольника, поэтому в формуле Пика остается лишь одна искомая неизвестная которую мы можем найти.
Рассмотрим пример:

S = 16 + 4 + 4,5 + 6 + 1 + 2 = 33,5
m = 15
n = 33,5 – 7,5 +1 = 27 — точек лежит строго внутри многоугольника
Вот так вот решается эта задачка!

Вот и все! Надеюсь, Вам понравилась статья, и я напишу ее вторую часть.

Знаковая площадь треугольника и предикат «По часовой стрелке»

Определение

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

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

Вычисление

Воспользуемся понятием косого (псевдоскалярного) произведения векторов. Оно как раз равно удвоенной знаковой площади треугольника:

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

(Модуль косого произведения двух векторов равен модулю векторного произведения их.)

Косое произведение вычисляется как величина определителя, составленного из координат точек:

Раскрывая определитель, можно получить такую формулу:

Можно сгруппировать третье слагаемое с первыми двумя, избавившись от одного умножения:

Последнюю формулу удобно записывать и запоминать в матричном виде, как следующий определитель:

Реализация

Функция, вычисляющая удвоенную знаковую площадь треугольника:

Функция, возвращающая обычную площадь треугольника:

Функция, проверяющая, образует ли указанная тройка точек поворот по часовой стрелке:

Функция, проверяющая, образует ли указанная тройка точек поворот против часовой стрелки:

Ориентированная площадь треугольника формула

Ориентированной площадью треугольника ABC называется величина (ABC), равная его площади, взятой со знаком плюс, если обход треугольника в порядке A–B–C–A совершается против часовой стрелки и со знаком минус, если по часовой стрелке (рис. 1). Таким образом, строго говоря, ориентированная площадь (ABC) определена не для треугольника как такового, а для «ориентированного треугольника», т.е. треугольника с заданным порядком вершин, причем (ABC) = (CAB) = (BCA) = –(ACB) = –(BAC) = –(CBA).

С помощью ориентированной площади во многих ситуациях можно избавиться от необходимости рассматривать разные расположения точек. Например, если точка D лежит внутри треугольника ABC (рис. 2), то справедливо равенство:

(которое, в частности, используется при выводе формулы площади треугольника через радиус вписанной окружности); если D лежит вне треугольника, но внутри угла ABC, то S_ABC = S_DAB + S_DBC – S_DCA, для других расположений получатся другие наборы знаков в левой части. Равенство для ориентированных площадей
(ABC) = (DAB) + (DBC) + (DCA) (1)
справедливо при любом расположении точек A, B, C, D на плоскости. Также имеет место векторное равенство
,
из которого следует, что барицентрические координаты точки D относительно треугольника ABC пропорциональны ориентированным площадям треугольников DBC, DCA, DAB.

Если координаты точек A и B равны (x_A;y_A) и (x_B; y_B), то ориентированная площадь треугольника OAB, где O – начало координат, равна
(OAB) =( x_A y_B – y_A x_B)
(при этом треугольник с вершинами (0; 0), (1; 0), (0; 1) считается положительно ориентированным). Эта формула позволяет написать выражение для площади произвольного треугольника ABC через координаты его вершин: для этого надо записать в координатах правую часть формулы (1), взяв в ней D = O). Приведем обобщение этой формулы для произвольного n-угольника A₁A₂. A_n:

(A₁A₂. A_n) = ((x₁y₂–y₁x₂) + (x₂y₃–y₂x₃) + . + (x_{n – 1}y_n – y_{n – 1}x_n) + (x_ny₁ – y_nx₁)), (2)

где (x_i; y_i) – координаты точки A_i. Можно показать, что модуль этой величины равен (обычной) площади n-угольника, а знак определяется направлением обхода вершин в порядке A₁ – A₂ – … – A_n – A₁, как и в случае треугольника. При этом правая часть (2) не зависит от того, с какой вершины начинать обход – существенна только последовательность вершин. Если же поменять порядок на обратный – A_nA _{n – 1}. A₁, то и знак правой части равенства (2) изменится на обратный (а модуль не изменится).

«Направление против часовой стрелки» не является математически строгим понятием, поэтому и данное в начале определение нестрогое. Формализовать наглядное представление о «направлении обхода», или «ориентации», многоугольника довольно сложно. Наиболее простой путь – непосредственно использовать формулу (2) как определение ориентированной площади; тогда многоугольник можно считать положительно или отрицательно ориентированным в зависимости от знака его ориентированной площади. Нетрудно проверить, что при этом будут выполняться интуитивно понятные утверждения об ориентации, например: два треугольника в общей стороной – ABC и ABD будут одинаково ориентированы тогда и только тогда, когда вершины C и D лежат по одну сторону от AB, при осевой симметрии ориентация многоугольника меняется на противоположную и др.

Вычислительная геометрия, или как я стал заниматься олимпиадным программированием.Часть 1

Здравствуйте, уважаемые хабравчане! Это моя вторая статья, и мне хотелось бы поговорить о вычислительной геометрии.

Немного истории

Я являюсь студентом уже 4 курса математического факультета, и до того как я начал заниматься программированием, я считал себя математиком на 100 процентов.

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

Мне очень нравится подход моего преподавателя: «разберись с этой темой, а потом расскажи нам, да так чтоб мы все поняли».

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

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

Вступление

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

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

Немного теории о векторах

Отрезок, для которого указано, какой из его концов считается началом, а какой — концом, называется вектором. Любая точка пространства также может рассматриваться как вектор. Такой вектор называется нулевым. Начало и конец нулевого вектора совпадают, и он не имеет какого-либо определенного направления.

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

Скалярное произведение векторов

Скалярное произведение векторов — это число, равное произведению длин этих векторов на косинус угла между ними.
(a, b) = |a||b|cos∠(a, b)

Если векторы заданы своими координатами a(x₁, y₁), b(x₂, y₂) то скалярное произведение (a, b) = x₁x₂ + y₁y₂.

Косое произведение векторов

Псевдоскалярным или косым произведением векторов на плоскости называется число
[a, b] = |a||b|sinθ
где — угол вращения (против часовой стрелки) от a к b. Если хотя бы один из векторов a и b нулевой, то полагают [a, b] = 0.
Если векторы заданы своими координатами a(x₁, y₁), b(x₂, y₂) то косое произведение [a, b] = x₁y₂ — x₂y₁.
Геометрически косое произведение векторов представляет собой ориентированную площадь параллелограмма, натянутого на эти вектора.

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

А теперь займемся практикой

Начнем с треугольников

Задача №1

Задача очень простая, а именно: по введенным трем числам a, b, c определить существует ли треугольник с такими сторонами.

Решение
Понятно, что здесь нужно только проверить неравенство треугольника: a + b > c, a + c > b, b + c > a. Интересно, при изучении неравенства треугольника только ли у меня возник вопрос: не могут ли отрицательные числа тоже удовлетворять этим трем неравенствам? Оказывается, нет! Если мы сложим каждое неравенство, то получим a > 0, b > 0, c > 0. Поэтому неравенство треугольника является необходимым и достаточным условием существования треугольника.

Задача №2

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

Решение
С первого взгляда решение кажется очевидным: вычислить стороны треугольника и свести задачу к предыдущей. Однако поскольку расстояние между двумя точками A(x₁, y₁), B(x₂, y₂) вычисляется по формуле √(x₁-x₂) 2 +(y₁-y₂) 2 то при извлечении корня возможна потеря точности, что плохо скажется на проверке неравенства треугольника. Оказывается, что если треугольник задан координатами своих вершин, то вычислять длины его сторон и проверять неравенство треугольника не требуется. В этом случае треугольника не существует тогда и только тогда, когда данные три точки лежат на одной прямой. А это легко проверяется через косое произведение векторов. Если оно равно нулю, то векторы коллинеарные, то есть все три точки лежат на одной прямой.

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

Задача №3

Треугольник задан своими сторонами. Определить тип треугольника: тупоугольный, прямоугольный или остроугольный.

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

Из курса геометрии известно, что напротив большей стороны лежит больший угол (он нам и нужен). Поэтому если мы выясним чему равен больший угол, то поймем тип треугольника:

1. Угол больше 90° – треугольник тупоугольный
2. Угол меньше 90°– треугольник остроугольный
3. Угол равен 90°– треугольник прямоугольный

Воспользуемся теоремой косинусов:

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

• Если cosα > 0, то a 2 2 + c 2 – треугольник остроугольный
• Если cosα = 0, то a 2 = b 2 + c 2 – треугольник прямоугольный
• Если cosα 2 > b 2 + c 2 – треугольник тупоугольный

где a – большая сторона.

Задача №4

Задача аналогична предыдущей задаче, только треугольник задан не своими сторонами, а координатами вершин.

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

Задача №5

По данным сторонам треугольника найти его площадь.

Решение
Очевидно решение, заключается в применение формулы Герона.

Кстати, никого не интересовало доказательство этой формулы?

Задача №6

Вычислить площадь треугольника заданного координатами своих вершин.

Решение
Не будем говорить о решении, которое сводится к предыдущей задачи, а попробуем воспользоваться геометрическим смыслом косового произведения. Геометрически косое произведение двух векторов определяет ориентированную площадь параллелограмма натянутого на эти вектора. Поскольку диагональ параллелограмма разбивает его на два равновеликих треугольника, то можем найти площадь нашего треугольника, как половину площади параллелограмма.
Для векторов a(x₁, y₁), b(x₂, y₂)

S = (x₁y₂ — x₂y₁) / 2 — ориентированная площадь треугольника

Задача №7

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

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

Метод площадей

Если сумма площадей треугольников AKB, AKC, BKC (не ориентированных, а «обычных») больше площади треугольника ABC точка лежит вне треугольника. Если же сумма первых трех площадей равна четвертой, то нужно проверить, не равна ли нулю одна из трех площадей. Если равна, то точка лежит на границе треугольника, иначе – внутри.
Вычислять площади треугольников, естественно, надо через косое произведение векторов. Этот метод не очень хороший. Поскольку здесь используются сравнение чисел с плавающей точкой, а это в свою очередь может привести к принятию неверного решения при сравнении. Второй метод опять таки опирается на вектора, он намного эффективнее во всех отношениях.

Проверка полуплоскостей

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

В первом примере сторона AB разводит вершину C и точку K по разным полуплоскостям, поэтому точка лежит снаружи.

Задача №8

Вычисление площади многоугольника заданного координатами своих вершин.

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

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

Метод трапеций

Для того чтобы посчитать площадь многоугольника нужно разбить его на трапеции, так как это показано на рисунке, а затем сложить ориентированные площади полученных трапеций это будет ориентированной площадью исходного многоугольника.
S = S_{A1 A2 B2 B1} + S_{A2 A3 B3 B2} + S_{A3 A4 B5 B3} + S_{A4 A5 B6 B5} + S_{A5 A6 B4 B6} + S_{A6 A1 B1 B4}
Площади трапеций считаем по известной формуле: полусумма оснований на высоту
S_{A1 A2 B2 B1} = 0.5 * (A₁B₁ + A₂B₂) *(B₂ — B₁)

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

Метод треугольников

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

Задача №9

Многоугольник задан координатами своих вершин в порядке его обхода. Необходимо проверить является ли многоугольник выпуклым.

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

Задача опять сводится к вычислению косового произведения векторов, а именно у выпуклого многоугольника знаки косых произведений [A_i A_i+1, A_i+1 A_i+2] либо положительны, либо отрицательны. Поэтому если мы знаем направление обхода, то знак косых произведений для выпуклого многоугольника одинаков: он неотрицателен при обходе против часовой стрелки и неположителен при обходе по часовой стрелки.

Задача №10

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

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

Для любого многоугольника с целочисленными координатами вершин справедлива формула Пика: S = n + m/2 — 1, где S – площадь многоугольника, n – количество целых точек лежащих строго внутри многоугольника, m – количество целых точек лежащих на границе многоугольника. Поскольку площадь многоугольника мы знаем как вычислять, то S известно. Так же мы можем вычислить количество целых точек лежащих на границе многоугольника, поэтому в формуле Пика остается лишь одна искомая неизвестная которую мы можем найти.
Рассмотрим пример:

S = 16 + 4 + 4,5 + 6 + 1 + 2 = 33,5
m = 15
n = 33,5 – 7,5 +1 = 27 — точек лежит строго внутри многоугольника
Вот так вот решается эта задачка!

Вот и все! Надеюсь, Вам понравилась статья, и я напишу ее вторую часть.

{«root»:»text»,»url»:»lopshits_vychislenie_ploshhadey_orientirovannyh_figur_1956″,»surl-package»:»/text/%PACKAGE%/?query=%QUERY%»,»surl-page»:»/text/%PACKAGE%/p%PAGE%/?query=%QUERY%»,»query»:»»»»,»section»:»library»,»mode-gfx»:true,»mode-html»:true,»mode-prefer»:»gfx»,»layout-prefer»:»1×1″,»zoom»:{«1×1»:{«level»:100,»_w»:false,»_h»:true},»2×1″:{«level»:100,»_w»:true,»_h»:false},»html»:{«level»:100}},»textsize-prefer»:»2″,»textfont-prefer»:»a»,»tree-type»:»ajax»,»tree-state»:»visible»,»printbox-state»:»hidden»,»print-allowed»:»1″,»searchbox-state»:»hidden»,»searchbox-type»:»inline»,»goto-pageno»:null,»goto-page»:-1,»defw»:»800″,»defh»:»1307″,»minh»:1307,»maxh»:1307,»fixeven»:null,»package»:»left»,»sitemode»:»live»,»user»:{«uuid»:»»}}