Принадлежность точки N-угольнику. - C++ - Ответ 1831977

#include <iostream>
#include <vector>
 
class Point
{
    double x, y;
 
public:
    Point(): x(0.), y(0.) {   }
    Point(double X, double Y): x(X), y(Y) {   }
    Point(const Point& p): x(p.x), y(p.y) {   }
 
    double X () const {return x;}
    double Y () const {return y;}
 
    bool operator == (const Point& p) const { return p.x == x && p.y == y; }
    void operator = (const Point& p) { x = p.x; y = p.y; }
};
 
class Polyline
{
    std::vector<Point> vertex;
    void swap(Point segment[2])
    {
        if(segment[0].X() > segment[1].X())
        {
            Point p = segment[0];
            segment[0] = segment[1];
            segment[1] = p;
        }
    }
public:
    Polyline() {   }
    void addVertex(const Point& p) { vertex.push_back(p); }
 
    bool isVertex(const Point& p)
    {
        for(std::vector<Point>::const_iterator it = vertex.begin(), end = vertex.end(); it != end; ++it)
            if(*it == p)
                return true;
        return false;
    }
 
 
    bool onPlant(const Point& p)
    {
        /*Идея такая, если в полиноме найдутся два ребра, одно из которых выше точки, а другое ниже - точка наша!*/
        Point segment[2];
        segment[0] = *vertex.begin();
        segment[1] = *(vertex.end() - 1);
        swap(segment);
        bool under = false,
             on = false;
        if(segment[0].X() > p.X() || segment[1].X() < p.X() && segment[1].X() != segment[0].X() )
        {
            under = p.Y() + p.Y() < ( segment[0].X() - p.X() ) * ( segment[1].Y() - segment[0].Y() ) / ( segment[1].X() - segment[0].X() );
            on    = p.Y() + p.Y() > ( segment[0].X() - p.X() ) * ( segment[1].Y() - segment[0].Y() ) / ( segment[1].X() - segment[0].X() );
        }
 
        for(std::vector<Point>::const_iterator it = vertex.begin(), end = vertex.end() - 1; it != end; ++it)
        {
            segment[0] = *it;
            segment[1] = *(it + 1);
            swap(segment);
            if(segment[0].X() > p.X() || segment[1].X() < p.X() && segment[1].X() != segment[0].X() )
            {
                if(!under)
                    under = p.Y() < p.Y() + ( p.X() - segment[0].X() ) * ( segment[1].Y() - segment[0].Y() ) / ( segment[1].X() - segment[0].X() );
                if(!on)
                    on    = p.Y() > p.Y() + ( p.X() - segment[0].X() ) * ( segment[1].Y() - segment[0].Y() ) / ( segment[1].X() - segment[0].X() );
            }
        }
        return under && on;
    }
};
std::istream& operator >> (std::istream& cin, Point p)
{
    double m[2];
    cin >> m[0] >> m[1];
    p = Point(m[0],m[1]);
    return cin;
}
 
int main()
{
    setlocale( LC_ALL, "Russian" );
    Polyline polyline;
    Point p;
    std::cout << "Введите количество вершин: ";
    unsigned size;
    std::cin >> size;
    for(unsigned i = 0; i != size; ++i)
    {
        std::cout << i + 1 << ": ";
        std::cin >> p;
        polyline.addVertex(p);
    }
    std::cout << "Введите координаты проверяемой точки: ";
    std::cin >> p;
    std::cout << polyline.isVertex(p) ? "Это вершина многоугольника.\n" : "";
    std::cout << polyline.onPlant(p) ? "Точка лежит внутри многоугольника" : "Точка лежит вне многоугольника.";
}