计算几何 Graham_scan凸包 pku 1113

http://162.105.81.212/JudgeOnline/problem?id=1113

 

 

Graham扫描法寻找凸包

语法：Graham_scan(Point PointSet[],Point ch[],int n,int &len);

参数：

PointSet[]：输入的点集

ch[]：输出的凸包上的点集，按照逆时针方向排列

n：PointSet中的点的数目

len：输出的凸包上的点的个数

 

/*************凸包 1113pku Graham_scan */#include<iostream>//#include<vector>#include<cmath>using namespace std;struct point{ double x,y;}one;//vector<point> PointSet,ch;point PointSet[1005],ch[1005];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 multi(point p1,point p2,point p0){ return (p1.x-p0.x)*(p2.y-p0.y)-(p2.x-p0.x)*(p1.y-p0.y); }int Graham_scan(int n,int &len){ int i,j,k=0,top=2; point tmp; for(i=1;i<n;i++) if ( (PointSet[i].y<PointSet[k].y)|| ((PointSet[i].y==PointSet[k].y)&&(PointSet[i].x<PointSet[k].x)) ) k=i; tmp=PointSet[0]; PointSet[0]=PointSet[k]; PointSet[k]=tmp; for (i=1;i<n-1;i++) { k=i; for (j=i+1;j<n;j++) if ( (multi(PointSet[j],PointSet[k],PointSet[0])>0) || multi(PointSet[j],PointSet[k],PointSet[0])==0 &&dis(PointSet[0],PointSet[j])<dis(PointSet[0],PointSet[k]) ) k=j; tmp=PointSet[i]; PointSet[i]=PointSet[k]; PointSet[k]=tmp; } ch[0]=PointSet[0]; ch[1]=PointSet[1]; top=1; for (i=2;i<n;i++) { while(multi(PointSet[i],ch[top],ch[top-1])>=0&&top>1) top--; ch[++top]=PointSet[i]; } len=top+1; return 0;}int main(){ int m,l,len; double s; cin>>m>>l; for(int i=0;i<m;i++) { cin>>PointSet[i].x>>PointSet[i].y; } Graham_scan(m,len); s=dis(ch[0],ch[len-1]); for(int i=1;i<len;i++) { s+=dis(ch[i],ch[i-1]); } s+=3.1415926*2*l; printf("%.0lf/n",s); system("pause"); return 0;}