Вычислительная геометрия

Теоретический материал: основной и более обширный и подробный.

Полезно помнить про документацию на модуль math (ужатая версия на русском).

Для тех, кто допускает, что его жизнь будет так или иначе связана с программированием, рекомендуется создавать классы для базовых примитивов (точка, вектор, прямая, луч, окружность и т.п.) и определять соответствующие операции с ними. Например, разность двух точек может давать вектор, вектора можно скалярно и векторно перемножать. Сами вектора можно умножать на числа, складывать и вычитать. Прямую можно собрать по двум точкам, по точке и вектору, которому прямая должна быть параллельна или перпендикулярна. Функция "Пересечение" может выдавать список точек (возможно, пустой) И т.п. Здесь хороший простор для продумывания удобной архитектуры решения.

Русско-английский тематический словарик

точка point
вектор vector
прямая line
луч ray
отрезок segment
угол angle
окружность circle
треугольник triangle
прямоугольник rectangle
квадрат square
многоугольник polygon
окружность circle
медиана median
биссектриса bissector
высота altitude
пересечение intersection
длина length
периметр perimeter
площадь area
касательная tangent
скалярное произведение dot product
векторное произведение cross product
вектор нормали normal vector

Упражнения на векторные и скалярные произведения

01: Расстояние между двумя точками

Даны координаты двух точек. Найдите расстояние между ними.

Ввод Вывод
0 0
1 1
1.4142135623730951

02: Полярный угол точки

Даны два числа – координаты точки, не совпадающей с началом координат. Выведите ее полярный угол (величину от 0 до $2\pi$).

Ввод Вывод
2 3
0.98279

03: Угол между векторами

Даны четыре числа: координаты двух невырожденных векторов.

Выведите величину неориентированного угла между ними.

Ввод Вывод
2 1
3 5
0.56673

04: Площадь треугольника

Даны шесть чисел: координаты трех вершин треугольника.

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

Ввод Вывод
1 0
2 4
5 2
7.0

05: Классификация векторов

Даны четыре числа: координаты двух ненулевых векторов. Если эти вектора коллинеарны, выведите 1. Если эти вектора перпендикулярны, выведите 2. Иначе выведите 0.

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

Ввод Вывод
1 1
2 2
1
0 1
1 0
2
1 2
2 1
0

06: Три точки

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

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

Ввод Вывод
0 0
1 1
2 2
YES

07: Принадлежность точки лучу

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

Программа должна вывести YES, если точка принадлежит лучу, или NO в противном случае.

Ввод Вывод
1 6
3 7
5 8
NO

08: Принадлежность точки отрезку

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

Программа должна вывести YES, если точка принадлежит отрезку, или NO в противном случае.

Ввод Вывод
3 3
1 2
5 4
YES

09: Расстояние от точки до луча

Дано шесть чисел: координаты точки, координаты начала и конца вектора.

Программа должна вывести единственное число: расстояние от точки до луча, заданного вектором.

Ввод Вывод
2 1
1 1
0 2
1.0

10: Расстояние от точки до отрезка

Дано шесть чисел: координаты точки и координаты двух концов отрезка.

Программа должна вывести единственное число: расстояние от данной точки до данного отрезка.

Ввод Вывод
0 4
2 3
2 5
2.0

11: Принадлежит ли точка углу

Дан угол AOB (O - вершина угла, A и B - точки на сторонах) и точка P. Определите, принадлежит ли точка P углу AOB (включая его стороны: лучи OA и OB).

Программа получает на вход координаты точек A, O, B, P. Все координаты - целые, не превосходят 100 по модулю. Точки A, O, B не лежат на одной прямой.

Программа должна вывести слово YES или NO.

Ввод Вывод
0 1
0 0
1 0
1 1
YES
1 0
0 0
0 1
-1 -1
NO

12: Пересекаются ли два луча

Даны два луча: AB и CD (A и C - вершины лучей, B и D лежат на лучах). Проверьте, пересекаются ли они.

Программа получает на вход координаты точек A, B, C, D. Все координаты - целые, не превосходят 100 по модулю.

Программа должна вывести слово YES или NO.

Ввод Вывод
0 1
1 2
1 -1
1 0
YES
0 0
1 0
0 1
1 2
NO

13: Пересекаются ли два отрезка

Даны два отрезка AB и CD. Проверьте, пересекаются ли они.

Программа получает на вход координаты точек A, B, C, D. Все координаты - целые, не превосходят 100 по модулю.

Программа должна вывести слово YES или NO.

Ввод Вывод
5 1
2 6
1 1
7 8
YES

Упражнения на уравнение прямой

14: Уравнение прямой

Прямая задана двумя точками. Выведите коэффициенты A, B, C нормального уравнения прямой.

Ввод Вывод
1 2
3 1
-1 -2 5

15: Перпендикулярная прямая

Дано уравнение прямой и координаты точки. Выведите коэффициенты уравнения прямой, перпендикулярной данной прямой и проходящей через данную точку.

Ввод Вывод
0 1 -1
0 0
1 0 0

16: Принадлежность точки прямой

Даны координаты точки и уравнение прямой. Определите, принадлежит ли точка прямой, выведите YES или NO.

Ввод Вывод
3 7
-2 1 -1
YES

17: Взаимное расположение двух точек

Даны две точки и уравнение прямой, точки не лежат на прямой. Выведите YES, если точки лежат по одну сторону от прямой и NO в противном случае.

Ввод Вывод
0 0
2 4
2 -1 -1
YES

18: Классификация прямых

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

Программа должна вывести 1, если эти прямые совпадают, 2 – если параллельны, 3 – если перпендикулярны и 0 во всех остальных случаях.

Ввод Вывод
1 2 3
-1 -2 -3
1

19: Расстояние от точки до прямой

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

Программа должна вывести одно число: расстояние от данной точки до данной прямой.

Ввод Вывод
1 1
1 1 -1
0.70711

20: Параллельная прямая

Даны четыре числа: коэффициенты нормального уравнения прямой и величина d.

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

Ввод Вывод
0 -1 1
1
0 -1 2

21: Основание перпендикуляра

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

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

Ввод Вывод
1 1
1 1 -1
0.5 0.5

22: Пересечение двух прямых - 1

Дано шесть чисел: коэффициенты нормальных уравнений двух непараллельных прямых.

Программа должна вывести два числа: координаты точки пересечения данных прямых.

Ввод Вывод
1 1 -1
1 -1 0
0.5 0.5

23: Пересечение двух прямых - 2

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

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

Если прямые не пересекаются, выведите одно число 0. Если прямые совпадают, выведите 2. Если прямые пересекаются ровно в одной точке, то выведите сначала число 1, а затем два вещественных числа - координаты точки пересечения.

Ввод Вывод
1 1
2 2
1 10
2 11
0
0 0
1 1
1 0
-1 2
1 0.5 0.5

Упражнения на треугольники

Треугольник задается координатами трех точек - вершин треугольника.

24: Биссектриса угла

Даны координаты трех точек O, A, B. Постройте уравнение прямой, являющейся биссектрисой угла ∠AOB.

Ввод Вывод
1 1
1 0
0 1
-1.0 1.0 0.0

25: Точка пересечения медиан

Дан треугольник. Найдите точку пересечения его медиан.

Ввод Вывод
0 0
1 0
0 1
0.333333 0.333333

26: Точка пересечения биссектрис

Дан треугольник. Найдите точку пересечения его биссектрис.

Ввод Вывод
0 0
1 0
0 1
0.292893 0.292893

27: Точка пересечения высот

Дан треугольник. Найдите точку пересечения его высот.

Ввод Вывод
0 0
1 0
0 1
0.0 0.0

28: Точка пересечения срединных перпендикуляров

Дан треугольник. Найдите точку пересечения его срединных перпендикуляров.

Ввод Вывод
0 0
1 0
0 1
0.5 0.5

29: Вписанная окружность

Выведите радиус окружности, вписанной в данный треугольник.

Ввод Вывод
0 0
1 0
0 1
0.2928932188134525

30: Описанная окружность

Выведите радиус окружности, описанной вокруг данного треугольник.

Ввод Вывод
0 0
1 0
0 1
0.7071067811865476

31: Принадлежность точки треугольнику

Дан треугольник и точка. Выведите слово YES, если точка принадлежит треугольнику (она может лежать на границе треугольника) и слово NO в противном случае.

Ввод Вывод
0 0
1 0
0 1
1 1
NO

32: Минимальная окружность

Дан треугольник. Найдите минимальную окружность, содержащую внутри себя треугольник.

Выведите три числа: координаты центра и радиус данной окружности.

Ввод Вывод
0 0
1 0
0 1
0.5 0.5 0.7071067811865476

33: Расстояние от прямой до окружности

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

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

Ввод Вывод
0 0 1
1 0 -2
1.0

34: Пересечение прямой и окружности

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

В первой строке одно число K, равное количеству точек пересечения прямой с окружностью. Далее в K строках координаты самих точек.

Ввод Вывод
1 1 1
1 -1 0
2
1.70711 1.70711
0.29289 0.29289

35: Угол обзора

Даны пять чисел: координаты центра окружности, ее радиус (в первой строке), координаты точки, лежащей вне окружности (во второй строке).

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

Ввод Вывод
1 1 1
2 2
1.5707963267948966

36: Поворот точки

Даны три числа: координаты точки и угол (в радианах, задан в виде действительного числа).

Выведите координаты точки, полученной поворотом данной точки вокруг начала координат на данный угол в положительном направлении.

Ввод Вывод
1 0
1.5707963267948966
0.0 1.0

37: Поворот прямой

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

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

Ввод Вывод
1 1 -1
1.5707963267948966
1.0 -1.0 1.0

38: Точки касания

Даны пять чисел: координаты центра окружности, ее радиус (в первой строке), координаты точки (во второй строке).

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

Ввод Вывод
1 1 1
2 2
2
1.0 2.0
2.0 1.0

39: Длина дуги

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

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

Ввод Вывод
1 1 1
2.0 1.0
1.0 2.0
1.5707963267948966

40: Пересечение окружностей

Даны шесть чисел – координаты центров и радиусы двух окружностей окружности.

В случае если количество общих точек окружностей конечно, в первой строке вывести одно число K, равное этому количеству, далее в K строках координаты самих точек. Если указанных точек бесконечно много, вывести единственное число "3".

Ввод Вывод
2 3 1
3 2 1
2
3.0 3.0
2.0 2.0

Задачи на многоугольники

41: Площадь многоугольника

Задан многоугольник: сначала задано количество вершин многоугольника N, не превосходящее 100.000, затем N вершин многоугольника в порядке обхода. Выведите единственное число: площадь многоугольника.

Ввод Вывод
3   
1 0
0 1
1 1
0.5

42: Выпуклость многоугольника

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

Ввод Вывод
3   
1 0
0 1
1 1
YES

43: Принадлежность точки многоугольнику

Первая строка входных данных содержит число N, 3≤N≤105. Далее идет N точек – координаты вершин многоугольника. Затем идут координаты проверяемой точки. Программа должна вывести YES, если точка лежит внутри многоугольника или на его границе и NO в противном случае.

Координаты всех точек могут быть дробными числами!

Ввод Вывод
3  
4
0 0
100 0
1 1
0 100
2 2
NO

44: Принадлежность точки выпуклому многоугольнику

Первая строка входных данных содержит число N, 3≤N≤105. Далее идет N точек – координаты вершин выпуклого многоугольника. Затем записано число K , 3≤N≤105 – количество проверяемых точек. Затем в K строчках записаны координаты K проверяемых точек.

Для каждой проверяемой точки выведите YES, если точка лежит внутри многоугольника или на его границе и NO в противном случае.

Координаты всех точек могут быть дробными числами!

Ввод Вывод
3  
0 0
2 0
0 2
2
1 1
2 2
YES
NO

Разные задачи

45: Ловушка для Слонопотама

Пятачок и Винни-Пух каждое утро ходят пить чай в гости к Кролику. Естественно, самым коротким путем.

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

Напишите программу, которая посчитает длину самого короткого безопасного пути от домика Винни-Пуха до домика Кролика.

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

В первой строке входных данных записаны координаты домика Винни-Пуха $X_B$, $Y_В$, во второй строке — координаты домика Кролика $X_K$, $Y_K$, в третьей строке — координаты центра и радиус ловушки $X_L$, $Y_L$, $R_L$. Все координаты — целые числа из диапазона от –32000 до 32000. Радиус ловушки — натуральное число, не превышающее 32000.

Домики Винни-Пуха и Кролика не могут находиться внутри ловушки, но могут находиться на ее границе.

Выведите в выходной файл одно число — длину самого короткого безопасного пути от домика Винни-Пуха до домика Кролика с тремя знаками после точки.

Ввод Вывод
0 0
0 1
10 10 1
1.000
5 0
0 5
0 0 5
7.854
-5 0
5 0
0 0 3
11.861

46: Воздушный шарик

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

Поскольку Винни Пух очень любит покушать, то в данной задаче (да и не только в задаче) примем его за сферу радиуса $P$. Центр медведя находится на высоте $H_p$ над уровнем земли. Строго над медведем , находится еще одна сфера, радиуса $S$ — воздушный шарик; центр шарика находится на высоте $H_s$ над уровнем земли. Центры обеих сфер находятся на одной вертикальной прямой. По понятным причинам гарантируется, что сферы не пересекаются, однако могут касаться.

Считая, что ружье стреляет строго по прямой, вычислите минимальное расстояние $L$, на которое Пятачок должен отойти от места взлета, чтобы успешно поразить шарик. Шарик считается пораженным, если траектория пули хотя бы касается его поверхности; при этом если траектория пули касается медведя, то он считается невредимым.

Вводятся положительные целые числа $P$, $H_p$, $S$ и $H_s$, не превосходящие 10000.

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

Ввод Вывод
1 9 10 21
0.0

47: Космический кегельбан

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

Каждая кегля представляет собой высокий цилиндр с основанием в виде круга радиусом $r$ метров. Все кегли одинаковые. Кегли расставлены по следующим правилам. Кегли образуют $n$ рядов, в первом ряду стоит одна кегля, во втором — две, и так далее. В последнем $n$-м ряду стоит $n$ кеглей. Введем на плоскости систему координат таким образом, чтобы единица измерения была равна одному километру. Центр единственной кегли в первом ряду находится в точке $(0, 0)$. Центры кеглей во втором ряду находятся в точках $(–1, 1)$ и $(1, 1)$. Таким образом, центры кеглей в $i$-м ряду находятся в точках с координатами $(–(i – 1), i – 1)$, $(–(i – 3), i – 1)$, ..., $(i – 1, i – 1)$.

Игра происходит следующим образом. Используется шар с радиусом $q$ метров. Игрок выбирает начальное положение центра шара $(x_c, y_c)$ и вектор направления движения шара $(v_x, v_y)$. После этого шар помещается в начальную точку и двигается не останавливаясь в направлении вектора $(v_x, v_y)$. Считается, что шар сбил кеглю, если в процессе движения шара имеет место ситуация, когда у шара и кегли есть общая точка. Сбитые кегли не меняют направления движения шара и не сбивают соседние кегли при падении.

На рисунке приведен пример расположения кеглей для $r = 500$, $n = 4$ и шара для $q = 1000$, $x_c = –2$, $y_c = –2$, $v_x = 1$, $v_y = 1$.

Требуется написать программу, которая по заданным радиусу кегли $r$, количеству рядов кеглей $n$, радиусу шара $q$, его начальному положению $(x_c, y_c)$ и вектору направления движения $(v_x, v_y)$ определяет количество кеглей, сбитых шаром.

Первая строка входных данных содержит два целых числа: $r$ и $n$, разделенных ровно одним пробелом ($1 \le r \le 700$, $1 \le n \le 200000$).

Вторая строка входных данных содержит целое число $q$ ($1 \le q \le 10^9$).

Третья строка входных данных содержит два целых числа $x_c$ и $y_c$, разделенных ровно одним пробелом ($–10^6 \le x_c \le 10^6$, $–10^6 \le y_c$, $1000y_c \lt –(r + q)$).

Четвертая строка входных данных содержит два целых числа $v_x$ и $v_y$, разделенных ровно одним пробелом ($–10^6 \le v_x \le 10^6$, $0 \lt v_y \le 10^6$).

Программа должна вывести одно целое число — количество сбитых кеглей.

Ввод Вывод
500 4
1000
-2 -2
1 1
7

48: Кольцевая автодорога

К 2110 году город Флэтбург, являясь одним из крупнейших городов мира, не имеет обходной автомагистрали, что является существенным препятствием для его развития как крупнейшего транспортного центра мирового значения. В связи с этим ещё в 2065 году при разработке Генерального плана развития Флэтбурга была определена необходимость строительства кольцевой автомобильной дороги.

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

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

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

Программа получает на вход четыре строки. Каждая из них содержит по два целых числа: $x_i$ и $y_i$ — координаты места расположения достопримечательности. Никакие две достопримечательности не находятся в одной точке.

Все числа не превосходят 100 по абсолютной величине.

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

Во второй строке требуется вывести координаты центра дороги минимальной длины и её радиус. Если существует несколько разных способов построения дороги минимальной длины, необходимо вывести координаты центра и радиус любой из них. Координаты центра и радиус должны быть выведены с точностью не хуже $10^5$ и не должны превышать $10^9$ . Гарантируется, что существует хотя бы один план с такими параметрами.

Ввод Вывод
0 0
0 1
1 0
2 2
7
1.5 0.5 1.14412281
0 0
0 1
1 0
1 1
Infinity
0.5 0.5 0.0