题目描述:
覆盖平面上N个点的最小凸多边形
形象来讲就是平面上有n个柱子,把一根封闭的弹性绳套上这些柱子,绳子绷紧以后形成的多边形就是我们要求的凸包.
算法实现:
上一篇文章介绍了Jarvis算法
复杂度为 O(n*h)
但如果Jarvis的时间效率不够高,我们可以选择另一种更高效的算法–Graham算法.
Gramham算法其实是基于对Jarvis算法的一个优化,让我们可以更高效的找出凸包。
在Jarvis里,我们每次找寻下个点都需要 n 的去扫描
如果我们可以提前预处理出这个序列,那么就可以O(n)的找寻平面凸包。
具体流程如下:
首先找出最靠左的点(相同最靠下),然后将每个点按照到第一个点的逆时针极角排序。
然后搞一个栈,每次判断栈顶的两个元素和新加入的元素形成的角,若向右转则需弹出栈顶元素(因为是凸多边形),到最后留在栈里的元素即为凸包元素。
代码实现:
#include <cstdio>
#include <iostream>
#include <cmath>
#include <algorithm>
const int maxm=10005;
const double eps=1e-10;
struct Point{
double x,y;
Point operator -(Point &s)
{return (Point){x-s.x,y-s.y};}
};
Point d[maxm];
double operator * (Point p1,Point p2)
{
return (p1.x*p2.y-p2.x*p1.y);
}
int stk[maxm],top,n;
double distance(Point p1,Point p2)
{
return (p1.x-p2.x)*(p1.x-p2.x)+(p1.y-p2.y)*(p1.y-p2.y);
}
bool judgeleft(Point p0,Point p1,Point p2)
{
double s=(p1-p0)*(p2-p0);
return s<0||(s==0&&distance(p1,p0)>=distance(p2,p0));
}
bool comp(Point p1,Point p2)
{
double s=(p1-d[1])*(p2-d[1]);
return s<0||(s==0&&distance(p1,d[1])>=distance(p2,d[1]));
}
int main()
{
scanf("%d",&n);
int Fis=1;
for(int i=1;i<=n;i++)
{
scanf("%lf%lf",&d[i].x,&d[i].y);
if(d[i].x<d[Fis].x||(d[i].x==d[Fis].x&&d[i].y<d[Fis].y)) Fis=i;
}
std::swap(d[Fis],d[1]);
std::sort(&d[2],&d[n+1],comp);
d[n+1]=d[1];
stk[1]=1,stk[2]=2,top=2;
for(int i=3;i<=n+1;i++)
{
while(top>1&&judgeleft(d[stk[top-1]],d[i],d[stk[top]])) top--;
stk[++top]=i;
}
top--;
printf("%d\n",top);
for(int i=1;i<=top;i++)
printf("%.2lf %.2lf\n",d[stk[i]].x,d[stk[i]].y);
//for(int i=1;i<=n;i++) printf("%.2lf %.2lf\n",d[i].x,d[i].y);
return 0;
}