极角系Graham-Scan算法.docVIP

#includeiostream#includecmathusing namespace std;/*PointSet[]：输入的点集ch[]：输出的凸包上的点集，按照逆时针方向排列n：PointSet中的点的数目len：输出的凸包上的点的个数*/typedef struct { float x,y;}Point;//小于,说明向量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].xPointSet[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保存极角最小的那个点,或者相同距离原点最近 if(i!=k) { tmp=PointSet[i]; PointSet[i]=PointSet[k]; PointSet[k]=tmp; } } //前三个点先入栈 ch[0]=PointSet[0]; ch[1]=PointSet[1]; ch[2]=PointSet[2]; //判断与其余所有点的关系 for (i=3;in;i++) { //不满足向左转的关系,栈顶元素出栈 while(multiply(PointSet[i],ch[top],ch[top-1])=0) top--; //当前点与栈内所有点满足向左关系,因此入栈. ch[++top]=PointSet[i]; } len=top+1;}const int maxN=1000;Point PointSet[maxN];Point ch[maxN];int n;int len;int main() { int n=5; float x[]={0,3,4,2,1}; float y[]={0,0,0,3,1}; for(int i=0;in;i++) { PointSet[i].x=x[i]; PointSet[i].y=y[i]; } Graham_scan(PointSet,ch,n,len); for(int i=0;ilen;i++) coutch[i].x ch[i].yendl; system(pause); return 0;}