HDU 1392 Surround the Trees

題目大意：

求散點集凸包外圍總長度，結果保留兩位小數

 

代碼：

#include<stdio.h>
#include<stdlib.h>
#include<math.h>
#define MAX 102
#define eps 1e-6

typedef struct point 
{
	int x,y;
}point;

point p[MAX],stack[MAX];

bool isZero(double x)
{
	return (x>0 ? x : -x)<eps;
}

double dis(point a,point b)
{
	return sqrt(1.0*(a.x-b.x)*(a.x-b.x)+1.0*(a.y-b.y)*(a.y-b.y));
}

double crossProd(point a,point b,point c)
{
	return (b.x-a.x)*(c.y-a.y)-(b.y-a.y)*(c.x-a.x); //以最左下的點爲基準點，其他各點（逆時針方向）以極角從小到大的排序規則
}

int cmp(const void *a,const void *b)
{
	double k;
	point *c=(point *)a,*d=(point *)b;
	k=crossProd(p[0],*c,*d);
	if(k<eps || isZero(k) && dis(p[0],*c)>dis(p[0],*d))
		return 1;
	return -1;
}

double Graham(int n)
{
	int x=p[0].x, y=p[0].y, mi=0;
	int i,top;
	double len;
	point tmp;
	for(i=1;i<n;i++)
		if(p[i].x<x || p[i].x==x && p[i].y<y)
		{
			x=p[i].x;
			y=p[i].y;
			mi=i;
		} //找出最左下的一個點
	tmp=p[mi]; p[mi]=p[0]; p[0]=tmp;
	qsort(p+1,n-1,sizeof(point),cmp); //從p中第二個元素開始快排
	p[n]=p[0];
	stack[0]=p[0];
	stack[1]=p[1];
	stack[2]=p[2];
	top=2;
	for(i=3;i<=n;i++)
	{
		while(crossProd(stack[top-1],stack[top],p[i])<=eps && top>=2)
			top--; //加入一個點後,向右偏拐或共線，則上一個點不在凸包內，則--top，該過程直到不向右偏拐或沒有三點共線的點 
		stack[++top]=p[i]; //在當前情況下符合凸包的點，入棧
	}
	len=0;
	for(i=0;i<top;i++)
		len+=dis(stack[i],stack[i+1]);
	return len;
}

int main()
{
	int T,i;
	double len;

	while(scanf("%d",&T)!=EOF && T)
	{
		for(i=0;i<T;i++)
			scanf("%d%d",&p[i].x,&p[i].y);
		if(T==1) printf("0.00\n");
		else if(T==2) printf("%.2lf\n",dis(p[0],p[1]));
		else 
		{
			len=Graham(T);
			printf("%.2lf\n",len);
		}
	}
	return 0;
}

注意比較函數中精度的處理