凸包详解.docVIP

ACM计算几何之凸包问题 整理编辑----Blue Seven 何谓凸包问题 所谓凸包，是指平面上一点集Q，将其中某些点连起来得到一凸边形，使得满足Q中的所有的点都在该凸边形内（或其边上），则称该凸边形为点集Q的凸包(convex hull)。如上图中灰色线段组成的凸边形即为点集Q={p0，p1，p2，...，p11，p12}的凸包。 一组平面上的点，求一个包含所有点的最小凸边形，既是凸包问题。 形象地说：在一平木板（平面）上钉若干钉子（点），将一橡皮筋套上去后，会把钉子圈起来，形成一个凸边形，即为该点集的凸包。 解决凸包问题的常用算法及分析 常用来解决凸包问题的算法有Graham扫描法和Jarvis步进法。接下来简单介绍Graham扫描法(极角序)，至于Jarvis步进法可参考《算法导论》 由此看出，Graham扫描法很好得解决了凸包问题，接下来看看其时间复杂度： 显然，在快速排序时要O(nlogn)，而在扫描入栈出栈时为O(n)，所以总的时间复杂度是O(nlogn)。可以得知，Graham扫描法还算是比较高效的。实际上，Graham扫描法并非解决凸包问题最高效的算法，Melkman算法（一种Graham扫描法的改进）则更高效，为O(n)，但是最常用的还是Graham扫描法。 凸包问题应用举例 求点集中最大直径。 给定一个点集，求其面积最小的外接矩形。 ...等等 1113 1228 2007 2079 2187 Graham扫描法模板 typedef struct { double x,y; }Point; Point point[1001],stack[1001]; int n,top; double l; int Multi(Point p0,Point p1,Point p2) { return ((p1.x-p0.x)*(p2.y-p0.y)-(p2.x-p0.x)*(p1.y-p0.y)); } double Dis(Point p1,Point p2) { return sqrt((p1.x-p2.x)*(p1.x-p2.x)+(p1.y-p2.y)*(p1.y-p2.y)); } int cmp(Point a,Point b) { if(Multi(point[0],a,b)0)return 1; if(Multi(point[0],a,b)0)return 0; if(Dis(point[0],a)Dis(point[0],b)) return 1; return 0; } void Graham() { int i,k=0; Point tmp; for(i=1;in;i++) if((point[i].ypoint[k].y)||((point[i].y==point[k].y) (point[i].x point[k].x))) k=i; tmp=point[0]; point[0]=point[k]; point[k]=tmp; sort(point+1,point+n,cmp); stack[0]=point[0]; stack[1]=point[1]; i=top=1; while(in) if(Multi(stack[top],stack[top-1],point[i])0) //=0时共线点也入栈 stack[++top]=point[i++]; else top--; }