凸包(graham算法)--转载

凸包是指覆盖平面坐标系内若干点的面积最小的凸多边形。求凸包的第一步是确定:凸包的定点都在给定的点中。通过几何方法反证很容易得到这一结论。所以,只要从所有点中挑选若干正确的点,按顺序(顺时针或逆时针)排列,就相当与求得了凸包。一种简便的方法是"包裹法"(Gift-Wrapping)。将y坐标最低的点作为凸包的第一个顶点H1(易证,所有点中x或y坐标取到最大或最小值的点一定是凸包顶点之一)找到满足条件的点,该点与水平线的叉积为正且夹角最小,作为凸包的第二个点H2。再求与线段H1H2叉积为正且夹角最小的点,。。。。依此类推。
数据从文件读取 文件名data.txt 文件格式:

//计算几何中的凸包问题程序(graham算法)
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#define MAXN 10000

/**//* 顶点的类型定义 */
typedef struct {
    double x;
    double y;
    double arCos;
}Point1;
int n;     // 顶点的个数
Point1 points[MAXN];     // 顶点数组
int stack[MAXN];  // 栈

/**//*主函数*/
int main()
{
    void Init();
    void Make();
    Init();    //程序数据的读入
    Make();    //程序算法过程
    return 0;
}

/**//*数据读入函数*/
void Init()
{
    FILE *in;   // 采用读文件的方式,读入数据
    int i;
    in = fopen("data.txt", "r");
    fscanf(in, "%d", &n);
    for(i = 0;i<n;++i)
        fscanf(in, "%lf%lf", &points[i].x, &points[i].y);
    fclose(in);   // 关闭文件
}

/**//*算法实现函数*/
void Make()
{
    int Multi(Point1, Point1, Point1);   //计算两个向量的积
    double Angle(int);      //计算其余顶点与第一顶点的角度,为排序做准备
    void QSort(int, int);     //对顶点进行快速排序
    void Swap(int, int);     
    int i, j, t;
    double min = 32767.0;
    for(i = 0; i < n; ++i){     //找第一个顶点,做为算法的起始顶点
        if(points[i].y < min) {
            j = i;
            min = points[i].y;
        }
    }
    Swap(0, j);
    for(i = 1; i < n; ++i){     //计算除第一顶点外的其余顶点到第一点的线段与x轴的夹角
        points[i].arCos = Angle(i);
    }
    QSort(1, n-1);       //根据所得到的角度进行快速排序.
    for(i = 0; i <= 2; ++i) stack[i] = i;   //将前3个顶点压栈
    t = 2;
    while(i < n) {
        /**//*如果新的点,与最近入栈中的2点构成了一个"凹"角, 则将栈顶元素出栈. 直到把栈检查完*/
        while(Multi(points[stack[t-1]], points[stack[t]], points[i]) && t >= 1) 
            t--;    
        t++;        // 将新点压栈
        stack[t] = i;
        i++;
    }
    /**//*打印结果*/
    for(i = 0; i <=t; ++i)
        printf("<%.2lf, %.2lf> ",points[stack[i]].x, points[stack[i]].y);
}
int Multi(Point1 px, Point1 py, Point1 pz)
{
    double k;
    k = (py.x-px.x)*(pz.y-py.y) - (pz.x-py.x)*(py.y-px.y);  // 计算两个向量的向量积,
    // 判断3个点所成的角是不是一个"凹"角.
    if(k < 0) return 1;
    return 0;
}

/**//*角度计算函数*/
double Angle(int i)
{
    double j, k, m, n;
    j = fabs(points[i].x - points[0].x);
    k = fabs(points[i].y - points[0].y);
    m = sqrt(j*j+k*k);     //得到顶点i 到第一顶点的线段长度.
    n = acos(j/m);      //得到该线段与x轴的角度
    //强悍
    return n;
}

void QSort(int top, int bot)
{
    //快排
    int Loc(int, int);
    int pos;
   
    if(top < bot) {
        pos = Loc(top, bot);
        QSort(top, pos-1);
        QSort(pos+1, bot);
    }
}

int Loc(int top, int bot)
{
    void Swap(int, int);
    double x = points[top].arCos;
    int j, k;
    j = top+1;
    k = bot;
    while(1) {
        while(j < bot && points[j].arCos < x) j++;
        while(k > top && points[k].arCos > x) k--;
       
        if(j >= k) break;
       
        Swap(j, k);
    }
    Swap(top, k);
    return k;
}

void Swap(int px, int py)
{
    Point1 k;
    k = points[px];
    points[px] = points[py]; //注意
    points[py] = k;
}

 

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
好的,以下是Graham-Scan算法的C++实现。 ```cpp #include <bits/stdc++.h> using namespace std; struct Point { int x, y; }; // 叉积 int cross(const Point &a, const Point &b, const Point &c) { return (b.x - a.x) * (c.y - b.y) - (b.y - a.y) * (c.x - b.x); } // 求凸包,返回点集 vector<Point> convexHull(vector<Point> &points) { int n = points.size(); // 先按照x坐标从小到大排序,x坐标相同按y坐标从小到大排序 sort(points.begin(), points.end(), [](const Point &a, const Point &b) { return a.x == b.x ? a.y < b.y : a.x < b.x; }); // 求下凸壳 vector<Point> lower; for (int i = 0; i < n; ++i) { while (lower.size() >= 2 && cross(lower[lower.size() - 2], lower.back(), points[i]) <= 0) lower.pop_back(); lower.push_back(points[i]); } // 求上凸壳 vector<Point> upper; for (int i = n - 1; i >= 0; --i) { while (upper.size() >= 2 && cross(upper[upper.size() - 2], upper.back(), points[i]) <= 0) upper.pop_back(); upper.push_back(points[i]); } // 合并下凸壳和上凸壳 vector<Point> ans(lower); ans.insert(ans.end(), upper.begin() + 1, upper.end() - 1); return ans; } int main() { int n; cin >> n; vector<Point> points(n); for (int i = 0; i < n; ++i) cin >> points[i].x >> points[i].y; vector<Point> hull = convexHull(points); cout << "Convex Hull:" << endl; for (const Point &p : hull) cout << p.x << " " << p.y << endl; return 0; } ``` 在上述代码中: - `Point` 结构体表示一个点,包含 `x` 和 `y` 坐标; - `cross` 函数用于计算向量 $\overrightarrow{AB}$ 和 $\overrightarrow{AC}$ 的叉积,即 $(\overrightarrow{AB} \times \overrightarrow{AC})$,结果为正表示 $\overrightarrow{AB}$ 在 $\overrightarrow{AC}$ 的逆时针方向,结果为负表示 $\overrightarrow{AB}$ 在 $\overrightarrow{AC}$ 的顺时针方向,结果为 $0$ 表示 $\overrightarrow{AB}$ 与 $\overrightarrow{AC}$ 共线; - `convexHull` 函数用于求解凸包,输入为点集 `points`,输出为凸包点集; - 在函数内部,先按照 x 坐标从小到大排序,x 坐标相同按 y 坐标从小到大排序; - 接着求下凸壳,利用单调栈维护; - 再求上凸壳,同样利用单调栈维护; - 最后将下凸壳和上凸壳合并,得到最终的凸包点集。 希望能对你有所帮助。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值