Выпуклая оболочка

@NePascalist · Регистрация: 20.10.2018

Студворк — интернет-сервис помощи студентам

Здравствуйте!

Столкнулся с проблемой: не выходит правильно написать выпуклую оболочку

Ввод: количество точек (n), далее n координат (x, y) в декартовой системе координат
Вывод: количество точек на выпуклой оболочке и их координаты в декартовой системе координат

Вот код:

C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
#include <bits/stdc++.h>
 
using namespace std;
 
struct point
{
    long long x, y;
};
 
point operator-(const point& a, const point& b)
{
    return {a.x-b.x, a.y-b.y};
}
 
bool operator<(const point& a, const point& b)
{
    return a.x < b.x || a.x == b.x && a.y < b.y;
}
 
long long operator^(const point& a, const point& b)
{
    return a.x*b.y - a.y*b.x;
}
 
int rotation(const point& a, const point& b, const point& c)
{
    long long cross = (b-a)^(c-b);
    if (cross < 0)
        return -1;
    if (cross > 0)
        return 1;
    return 0;
}
 
void convex_hull(vector<point>& points)
{
    sort(points.begin(), points.end());
    int n = points.size();
    vector<point> top = {points[0]}, bottom = {points[0]};
    for (int i = 1; i < n; i ++)
    {
        if (rotation(points[0], points[n-1], points[i]) >= 0)
        {
            while (top.size() > 1 && rotation(top[top.size()-2], top[top.size()-1], points[i]) >= 0)
                top.pop_back();
            top.push_back(points[i]);
        }
        else
        {
            while (bottom.size() > 1 && rotation(bottom[bottom.size()-2], bottom[bottom.size()-1], points[i]) <= 0)
                bottom.pop_back();
            bottom.push_back(points[i]);
        }
    }
    reverse(bottom.begin(), bottom.end());
    bottom.pop_back();
    cout << top.size() + bottom.size() << "\n";
 
 
    for(int i = 0; i < top.size() ; i ++)
    {
        cout << top[i].x << " " << top[i].y << "\n";
    }
    for(int i = 0; i < bottom.size() ; i ++)
    {
        cout << bottom[i].x << " " << bottom[i].y << "\n";
    }
 
}
 
int main()
{
    ios_base::sync_with_stdio(0);
    cin.tie(0);
    cout.tie(0);
    int n, x, y;
    cin >> n;
    vector < point > points;
    for (int i = 0 ; i < n; i ++)
    {
        cin >> x >> y;
        points.push_back({x, y});
    }
    convex_hull(points);
}

Помогите пожалуйста, срочно

Заранее спасибо!

@NePascalist · 29.01.2019, 17:34 **[ТС]**

Тема закрыта

Если кому-то надо то вот код (правильный)

C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
struct point
{
    long long x, y;
};
 
point operator-(const point& a, const point& b)
{
    return {a.x-b.x, a.y-b.y};
}
 
bool operator<(const point& a, const point& b)
{
    return a.x < b.x || a.x == b.x && a.y < b.y;
}
 
long long operator^(const point& a, const point& b)
{
    return a.x*b.y - a.y*b.x;
}
 
int rotation(const point& a, const point& b, const point& c)
{
    long long cross = (b-a)^(c-b);
    if (cross < 0) return -1;
    if (cross > 0) return 1;
    return 0;
}
 
vector < point > convex_hull(vector<point>& points)
{
    sort(points.begin(), points.end());
    int n = points.size();
    vector<point> top = {points[0]}, bottom = {points[0]};
    for (int i = 1; i < n; i ++)
    {
        if (rotation(points[0], points[n-1], points[i]) >= 0)
        {
            while (top.size() > 1 && rotation(top[top.size()-2], top[top.size()-1], points[i]) >= 0)
                top.pop_back();
            top.push_back(points[i]);
        }
        else
        {
            while (bottom.size() > 1 && rotation(bottom[bottom.size()-2], bottom[bottom.size()-1], points[i]) <= 0)
                bottom.pop_back();
            bottom.push_back(points[i]);
        }
 
    }
 
    while (bottom.size() > 1 && rotation(bottom[bottom.size()-2], bottom[bottom.size()-1], points[n-1]) <= 0)
        bottom.pop_back();
 
    reverse(bottom.begin(), bottom.end());
 
    for (point u:bottom)
        top.push_back(u);
    top.pop_back();
    return top;
}

@Tlya · 07.02.2022, 14:14

А как насчет 3D точек?
Интернет завален решением 2D, а вот с координатами Z решение найти сложновато...

@Kuzia domovenok · 07.02.2022, 14:41

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

Новые блоги и статьи Все статьи Все блоги /
Модульная разработка через nuget packages DevAlt 07.03.2026 Сложившийся в . Net-среде способ разработки чаще всего предполагает монорепозиторий в котором находятся все исходники. При создании нового решения, мы просто добавляем нужные проекты и имеем. . .	Модульный подход на примере F# DevAlt 06.03.2026 В блоге дяди Боба наткнулся на такое определение: В этой книге («Подход, основанный на вариантах использования») Ивар утверждает, что архитектура программного обеспечения — это структуры,. . .	Управление камерой с помощью скрипта OrbitControls.js на Three.js: Вращение, зум и панорамирование 8Observer8 05.03.2026 Содержание блога Финальная демка в браузере работает на Desktop и мобильных браузерах. Итоговый код: orbit-controls-threejs-js. zip. Сканируйте QR-код на мобильном. Вращайте камеру одним пальцем,. . .	SDL3 для Web (WebAssembly): Синхронизация спрайтов SDL3 и тел Box2D 8Observer8 04.03.2026 Содержание блога Финальная демка в браузере. Итоговый код: finish-sync-physics-sprites-sdl3-c. zip На первой гифке отладочные линии отключены, а на второй включены:. . .
SDL3 для Web (WebAssembly): Идентификация объектов на Box2D v3 - использование userData и событий коллизий 8Observer8 02.03.2026 Содержание блога Финальная демка в браузере. Итоговый код: finish-collision-events-sdl3-c. zip Сканируйте QR-код на мобильном и вы увидите, что появится джойстик для управления главным героем. . . .	Реалии Hrethgir 01.03.2026 Нет, я не закончил до сих пор симулятор. Эта задача сложнее. Не получилось уйти в плавсостав, но оно и к лучшему, возможно. Точнее получалось - но сварщиком в палубную команду, а это значит, в моём. . .	Ритм жизни kumehtar 27.02.2026 Иногда приходится жить в ритме, где дел становится всё больше, а вовлечения в происходящее — всё меньше. Плотный график не даёт вниманию закрепиться ни на одном событии. Утро начинается с быстрых,. . .	SDL3 для Web (WebAssembly): Сборка библиотек: SDL3, Box2D, FreeType, SDL3_ttf, SDL3_mixer и SDL3_image из исходников с помощью CMake и Emscripten 8Observer8 27.02.2026 Недавно вышла версия 3. 4. 2 библиотеки SDL3. На странице официальной релиза доступны исходники, готовые DLL (для x86, x64, arm64), а также библиотеки для разработки под Android, MinGW и Visual Studio. . . .

@Tlya 16 / 16 / 10 Регистрация: 20.11.2015 Сообщений: 305
	07.02.2022, 14:14
	А как насчет 3D точек? Интернет завален решением 2D, а вот с координатами Z решение найти сложновато... 0

@Kuzia domovenok 4268 / 3327 / 926 Регистрация: 25.03.2012 Сообщений: 12,535 Записей в блоге: 1
	07.02.2022, 14:41
	Tlya, да тот же принцип. Перебираешь кандидатов на новые точки, оборачиваешь точками, треугольники, с которыми образуют максимальный угол с оболочкой. Всё. 0