Graham求凸包

最新推荐文章于 2023-02-14 17:34:43 发布

f_zyj

最新推荐文章于 2023-02-14 17:34:43 发布

阅读量1.7k

点赞数 1

分类专栏： ACM模版篇计算几何 ACM模板-f_zyj 文章标签： Graham

本文链接：https://blog.csdn.net/f_zyj/article/details/52004769

版权

ACM模版篇同时被 3 个专栏收录

156 篇文章 21 订阅

订阅专栏

ACM模板-f_zyj

147 篇文章 73 订阅

订阅专栏

计算几何

33 篇文章 0 订阅

订阅专栏

ACM模版

Graham

/*
 *  Graham 求凸包 O(N * logN)
 *  CALL: nr = graham(pnt, int n, res); res[]为凸包点集; 
 */
struct point
{
    double x, y;
};

bool mult(point sp, point ep, point op)
{
    return (sp.x - op.x) * (ep.y - op.y) >= (ep.x - op.x) * (sp.y - op.y);
}

//inline bool operator < (const point &l, const point &r)
//{
//    return l.y < r.y || (l.y == r.y && l.x < r.x);
//}

int graham(point pnt[], int n, point res[])
{
    int i, len, top = 1;
    sort(pnt, pnt + n);
    if (n == 0)
    {
        return 0;
    }
    res[0] = pnt[0];
    if (n == 1)
    {
        return 1;
    }
    res[1] = pnt[1];
    if (n == 2)
    {
        return 2;
    }
    res[2] = pnt[2];
    for (i = 2; i < n; i++)
    {
        while (top && mult(pnt[i], res[top], res[top - 1]))
        {
            top--;
        }
        res[++top] = pnt[i];
    }
    len = top;
    res[++top] = pnt[n - 2];
    for (i = n - 3; i >= 0; i--)
    {
        while (top != len && mult(pnt[i], res[top], res[top - 1]))
        {
            top--;
        }
        res[++top] = pnt[i];
    }
    return top; //  返回凸包中点的个数
}