计算几何凸包的算法 Andrew和Melkman算法

最新推荐文章于 2025-07-01 11:49:29 发布

Akahieveman

最新推荐文章于 2025-07-01 11:49:29 发布

阅读量5k

点赞数 3

CC 4.0 BY-SA版权

分类专栏：计算几何文章标签：凸包算法计算几何

本文链接：https://blog.youkuaiyun.com/Akahieveman/article/details/52416476

先介绍下二者的时间复杂度：

Andrew算法是葛立恒扫描法的变种，但是更快，时间O(nlogn)。

Melkman算法是采用双端队列，时间O(n)。

第一种是经典算法，第二种则是在解决时间要求高的问题上的一个也是目前我所知道最快的。

Andrew代码如下：

int Andrew(Point *p,int n,Point *q)
{
  sort(p,p+n);
  int m=0;
  for(int i=0;i<n;i++)
  {
    while(m>1&&Cross(q[m-1]-q[m-2],p[i]-q[m-2])<=0) m--;
    q[m++]=p[i];
  }
  int k=m;
  for(int i=n-2;i>=0;i--)
  {
    while(m>k&&Cross(q[m-1]-q[m-2],p[i]-q[m-2])<=0) m--;
    q[m++]=p[i];
  }
  if(n>1) m--;
  return m;
}

Melkman代码如下：

double Cross(Vector A,Vector B) { return A.x*B.y - A.y*B.x;}
double side(Point a,Point b,Point p)
{
  Vector A=Vector(b.x-a.x,b.y-a.y);
  Vector B=Vector(p.x-a.x,p.y-a.y);
  return Cross(A,B);
}
void Melkman(int n,int &head,int &tail)
{
  int i;
  head=tail=n;
  ch[tail++]=p[0];
  for(i=0;i<n;i++)
  {
    ch[tail]=p[i];
    if(side(ch[head],ch[head+1],p[i+1]