凸包——Andrew(安德鲁)算法
2014年07月12日 17:37:20
阅读数:1439
凸包(ConvexHull)就是把给定点包围在内部的,面积最小的凸多边形。
原理:首先把所有点按照从小到大排序(如果x相同,按照y从小到大排
序),删除重复点后得到序列p1; p2; :::,然后把p1和p2放到凸包
中。从p3开始,当新点在凸包“前进”方向的左边时继续,否则
依次删除最近加入凸包的点,直到新点在左边。
从左到右和从右到左各扫描一次
时间复杂度为O(n)
1.计算凸包,输入点数组p,点个数为n,输出点数组ch。函数返回凸包顶点数
2.如果不希望在凸包的边上有输入点。把<=改成<
bool cmp(Point x,Point y)
{
return x.x<y.x||(x.x==y.x&&x.y<y.y);//x从小到大排序,如果x相同则y从小到大排序
}
int Cross(Point x,Point y,Point z)
{
int x1=x.x-y.x;
int y1=x.y-y.y;
int x2=z.x-y.x;
int y2=z.y-y.y;
if((x1*y2-x2*y1)<=0) return 0;//在右边(关于 2. )
return 1;
}
int ConvexHull(Ponit *p,int n,Point *ch)
{
sort(p,p+n,cmp);
int m=0;
for(i=0;i<n;i++)//从左到右扫描
{
while( m>1 && !Cross(ch[m-1],ch[m-2],p[i])) m--;
ch[m++]=p[i]; //ch和p一样是定义的结构体
}
int k=m;
for(i=n-2;i>=0;i--)//从右到左扫描
{
while( m>k && !Cross(ch[m-1],ch[m-2],p[i])) m--;
ch[m++]=p[i];
}
if(n>1) m--;//凸包有m个顶点
return m;
}