#include<stdio.h>
#include<math.h>
#include<algorithm>
using namespace std;
typedef struct
{
double x;
double y;
} POINT;
POINT result[102];//保存凸包上的点,相当于所说的栈S
POINT a[102];
int n,top;
double Distance(POINT p1,POINT p2)//两点间的距离
{
return sqrt((p1.x-p2.x)*(p1.x-p2.x)+(p1.y-p2.y)*(p1.y-p2.y));
}
double Multiply(POINT p1,POINT p2,POINT p3)//叉积
{
return((p2.x-p1.x)*(p3.y-p1.y)-(p2.y-p1.y)*(p3.x-p1.x));
}
bool cmp(POINT aa,POINT bb)//按点乘排序,cos=a*b/|a||b|
{
int ax=aa.x-a[0].x;
int ay=aa.y-a[0].y;
int bx=bb.x-a[0].x;
int by=bb.y-a[0].y;
float cosa=ax/sqrt(ax*ax+ay*ay);
float cosb=bx/sqrt(bx*bx+by*by);
return cosa>cosb;
}
//寻找凸包的过程,p0,p1,p2..的寻找过程在下面main中进行了
void Tubao()
{
int i;
result[0].x=a[0].x;
result[0].y=a[0].y;
result[1].x=a[1].x;
result[1].y=a[1].y;
result[2].x=a[2].x;
result[2].y=a[2].y;
top=2;
for(i=3; i<=n; i++)
{
while(Multiply(result[top-1],result[top],a[i])<=0&&top>=2)
top--;
result[top+1].x=a[i].x;
result[top+1].y=a[i].y;
top++;
}
}
int main()
{
int i,p;
double px,py,len,temp;
while(scanf("%d",&n)!=EOF,n)//
{
for(i=0; i<n; i++)
scanf("%lf%lf",&a[i].x,&a[i].y);
if(n==1)
{
printf("0.00\n");
continue;
}
else if(n==2)
{
printf("%.2lf\n",Distance(a[0],a[1]));
continue;
}
//这里的目的好像是找出y坐标最小的点,之后把他定义为P0
py=-1;
for(i=0; i<n; i++)
{
if(py==-1||a[i].y<py)
{
px=a[i].x;
py=a[i].y;
p=i;
}
else if(a[i].y==py&&a[i].x<px)
{
px=a[i].x;
py=a[i].y;
p=i;
}
}
//swap(a[0],a[p])
temp=a[0].x;
a[0].x=a[p].x;
a[p].x=temp;
temp=a[0].y;
a[0].y=a[p].y;
a[p].y=temp;
sort(a+1,a+n,cmp);
a[n].x=a[0].x;
a[n].y=a[0].y;
//调用tubao()
Tubao();
len=0.0;
for(i=0; i<top; i++)
len=len+Distance(result[i],result[i+1]);
printf("%.2lf\n",len);
}
return 0;
}
hdu-1392(二维凸包模板)
最新推荐文章于 2018-07-27 13:49:48 发布