您的位置:首页 > 其它

凸包入门题 HDU 1392 凸包果题,求周长

2016-07-26 20:45 274 查看
最近开始学习凸包,先来两道凸包果题。

#include<cstdio>
#include<cmath>
#include<algorithm>
#include<iostream>
using namespace std;
const int MAXN=110;

struct point
{
double x,y;
};
point list[MAXN];
int stack[MAXN],top;

double cross(point p0,point p1,point p2)
{
return ((p1.x-p0.x)*(p2.y-p0.y)-(p1.y-p0.y)*(p2.x-p0.x));
}
double dis(point p1,point p2)
{
return sqrt((double)(p2.x-p1.x)*(p2.x-p1.x)+(double)(p2.y-p1.y)*(p2.y-p1.y));
}
bool cmp(point p1,point p2) //极角排序函数 , 角度相同则距离小的在前面
{
double tmp=cross(list[0],p1,p2);
if(tmp>0) return true;
else if(tmp==0&&dis(list[0],p1)<dis(list[0],p2)) return true;
else return false;
}
void init(int n) //输入,并把  最左下方的点放在 list[0]  。并且进行极角排序
{
int i,k;
point p0;
scanf("%lf%lf",&list[0].x,&list[0].y);
p0.x=list[0].x;
p0.y=list[0].y;
k=0;
for(i=1;i<n;i++)
{
scanf("%lf%lf",&list[i].x,&list[i].y);
if( (p0.y>list[i].y) || ((p0.y==list[i].y)&&(p0.x>list[i].x)) )
{
p0.x=list[i].x;
p0.y=list[i].y;
k=i;
}
}
list[k]=list[0];
list[0]=p0;

sort(list+1,list+n,cmp);
}

void graham(int n)
{
int i;
if(n==1) {top=0;stack[0]=0;}
if(n==2)
{
top=1;
stack[0]=0;
stack[1]=1;
}
if(n>2)
{
for(i=0;i<=1;i++) stack[i]=i;
top=1;

for(i=2;i<n;i++)
{
while(top>0&&cross(list[stack[top-1]],list[stack[top]],list[i])<=0) top--;
top++;
stack[top]=i;
}
}
}

int main()
{
int n;
double ans;
while(scanf("%d",&n)&&n)
{
ans=0.0;
init(n);
graham(n);
if(n==1)
{
printf("0.00\n");
continue;
}
else if(n==2)
{
printf("%.2f\n",dis(list[stack[0]],list[stack[1]]));
continue;
}
for(int i=1;i<=top;i++)
ans+=dis(list[stack[i]],list[stack[i-1]]);
ans+=dis(list[stack[0]],list[stack[top]]);
printf("%.2f\n",ans);
}
return 0;
}
内容来自用户分享和网络整理,不保证内容的准确性,如有侵权内容,可联系管理员处理 点击这里给我发消息
标签: