计算二维凸包的Graham_Scan算法.pdfVIP

计算二维凸包的Graham_Scan 算法 凸包的概念： 点集Q 的凸包(convex hull)是指一个最小凸多边形，满足Q 中的点或者在多 边形边上或者在其内。下图中由红色线段表示的多边形就是点集Q={p0,p1,...p12} 的凸包。 凸包的求法： 现在已经证明了凸包算法的时间复杂度下界是O(n*logn)，但是当凸包的顶 点数h 也被考虑进去的话，Krikpatrick 和Seidel 的剪枝搜索算法可以达到 O(n*logh)，在渐进意义下达到最优。最常用的凸包算法是Graham 扫描法和Jarvis 步进法。本文只简单介绍一下Graham 扫描法，其正确性的证明和Jarvis 步进法 的过程大家可以参考《算法导论》。 对于一个有三个或以上点的点集Q，Graham 扫描法的过程如下： 令p0 为Q 中Y-X 坐标排序下最小的点； 设p1,p2,...pm为对其余点按以p0 为中心的极角逆时针排序所得的点集（如 果有多个点有相同的极角，除了距p0 最远的点外全部移除； 压p0 进栈S； 压p1 进栈S； 压p2 进栈S； for i ← 3 to m { do while{ 由S 的栈顶元素的下一个元素、S 的栈顶元素以及pi 构成的折 线段不拐向左侧}对S 弹栈; 压pi 进栈S; } return S; 此过程执行后，栈S 由底至顶的元素就是Q 的凸包顶点按逆时针排列的点 序列。需要注意的是，我们对点按极角逆时针排序时，并不需要真正求出极角， 只需要求出任意两点的次序就可以了。而这个步骤可以用前述的矢量叉积性质实 现。 程序示例 #include iostream #include cmath using namespace std; /* PointSet[]：输入的点集 ch[]：输出的凸包上的点集，按照逆时针方向排列 n：PointSet中的点的数目 len：输出的凸包上的点的个数 */ struct Point { float x,y; }; //小于,说明向量p0p1的极角大于p0p2的极角 float multiply(Point p1,Point p2,Point p0) { return ((p1.x-p0.x)*(p2.y-p0.y)- (p2.x-p0.x)*(p1.y-p0.y)); } float dis(Point p1,Point p2) { return (sqrt((p1.x-p2.x)*(p1.x-p2.x)+(p1.y-p2.y)*(p1.y-p2.y))); } void Graham_scan(Point PointSet[],Point ch[],int n,int len) { int i,j,k=0,top=2; Point tmp; //找到最下且偏左的那个点 for(i=1;in;i++) if ((PointSet[i].yPointSet[k].y)||((PointSet[i].y==PointSet[k].y)(PointSet[i].xPoi ntSet[k].x))) k=i; //将这个点指定为PointSet[0] tmp=PointSet[0]; PointSet[0]=PointSet[k]; PointSet[k]=tmp; //按极角从小到大,距离偏短进行排序 for (i=1;in-1;i++) { k=i; for (j=i+1;jn;j++) if( (multiply(PointSet[j],PointSet[k],PointSet[0])0) ||((multiply(PointSet[j],PointSet[k],PointSet[0])==0) (dis(PointSet[0],PointSet[j])dis(PointSet[0],PointSet[k]))) ) k=j;//k保存极角最小的那个点,或者相同距离原点最近 tmp=PointSet[i]; PointSet[i]=Point