凸包问题的Graham扫描算法实现

凸包定义：

 

   平面点集S的凸包定义为包含S中的所有点的最小凸多边形。

凸包问题：

1. 给出点集S，求出S的凸包的顶点（极点）。
2. 以顺（逆）时针方向输出求得的凸包的顶点。                                                                                         

算法：

1. Graham扫描算法 (  O( n*logn) )
2. QuickHull算法(  O( n*logn) )

算法预备知识：

• 考虑从点p1=（x1，y1）到另一个点p2=（x2，y2）的有向线段<p1,p2>. 如果q=（x3，y3）是另外一个点。定义：q在<p1,p2>的左（右）边，若向量<p1,p2>与<p2,q>的外积大（小）于0 

算法基本思想：

• Graham扫描算法 ： 大体思路是将不是凸包顶点的点从点集中去掉。

1. 找出S中具有最小y坐标的点p（通过选取最左边的点打破平局）
2. 根据点和p的连线 与 x轴正方向所成的角度，对S中的点进行排序（由小到大），并将p放在最前面。
3. 从p点开始扫描排序后的S集合。如果这些点都在凸包上，则每三个相继的点p1,p2,p3满足以下性质：p3在向量<p1,p2>的左边.如果出现相继的三个点p1,p2,p3不满足上述性质，则p2点一定不是凸包的顶点，应立即去除。        

以下是Graham扫描算法具体实现：实现过程中使用了链表。实际上，使用双向链表会使得代码更加简洁，但同时会增加空间。  该算法实现面向过程。

 

 

/* 凸包问题：graham扫描算法实现 程序未对重复点做特殊化处理。在GCC-4.4.5 （ubuntu 10.04）测试成功 算法复杂度为O（n*log n） ，时间的主要部分是排序。 lidachao 2011 */ #include <stdlib.h> #include <stdio.h> #include <list> #include <algorithm> #include <math.h> #include <ctime> using namespace std; #define MAX 20 //给定的点的个数 //生成的点的在区域 [-RAND_max/2, RAND_max/2]X[-RAND_max/2, RAND_max/2] 内 #define RAND_max 20 typedef struct point{ float x,y; float angle; }point; /******随机生成点，此处为简单起见用整型数作为测试数据*****/ void creatdata( list<point> * plist){ srand(time(0)); point p1; int i; for ( i=0 ;i<MAX ; i++){ int r1=random()%RAND_max-RAND_max/2,r2=random()%RAND_max-RAND_max/2; p1.x=float(r1);p1.y=float(r2); plist->push_back(p1); } } /*****计算向量(x1,y1)和(x2,y2)的夹角余弦值*******/ float vectors_cos(float x1,float y1,float x2,float y2){ return ( x1*x2+y1*y2 )/(sqrt(x1*x1+y1*y1)*sqrt(x2*x2+y2*y2)); } /*****用于排序的compare函数 *************/ int comp(point & x1, point &x2){ return x1.angle >x2.angle; } /*****判断p3是否在向量（p1,p2）左边（共线情况不纳入）************/ bool left_rotated(point * p1,point *p2, point *p3 ){ return p1->x*(p2->y-p3->y) + p3->x*(p1->y-p2->y) + p2->x*(p3->y-p1->y)>0; } /*****graham函数将plist中的非顶点去掉，余下的是最小凸多边形的顶点**********************/ void graham( list<point> * plist){ if(plist->size()<3) return; list<point>::iterator it1= plist->begin(),it2=plist->end(),it3=it1; float x0=it1->x , y0=it1->y; /*******找出y值最小的点it3***************/ for ( ; it1!=it2 ;it1++ ){ if ( it1->y < y0 ){ y0=it1->y; x0=it1->x; it3=it1; } else if( it1->y == y0 && it1->x<x0){ x0=it1->x; it3=it1; } } /********计算it3与每个顶点形成的向量与x轴正向的夹角的余弦**/ for ( it1=plist->begin(); it1!=it2 ;it1++ ) it1->angle=vectors_cos(1,0,it1->x-it3->x,it1->y-it3->y); /******特殊化处理，it3点本身对应的夹角是无穷大**/ it3->angle=RAND_max+1; /******将点按照夹角余弦值排序*******/ plist->sort(comp); /********将右旋的点去掉********/ it3=plist->begin();it1=it3; it3++;it2=it3;it3++; plist->push_back(*it1); do{ if(left_rotated(&*it1,&*it2,&*it3)){ it1++;it2++;it3++; } else{ plist->erase(it2);it2=it3;it3++; } }while(it3!=plist->end()); plist->pop_back(); } int main(){ list<point> plist; creatdata(&plist); list<point>::iterator it1= plist.begin(),it2=plist.end(),it3=it1; graham(&plist); for (it1=plist.begin(); it1!=it2;it1++){ printf("%.3f %.3f %.3f/n",it1->x,it1->y,it1->angle); } return 0; }