hdu1392 Surround the Trees 凸包

第一次做凸包，这道题要特殊考虑下，n=2时的情况，要除以二才行。

我是从最左边的点出发，每次取斜率最大的点，一直到最右边的点。

然后从最左边的点出发，每次取斜率最低的点，一直到最右边的点。

#include<stdio.h>
#include<math.h>
#include<algorithm>
#include<string.h>
using namespace std;
const double eps=1e-9;
struct node
{
	double x,y;
}f[105];
bool flag[105];
int n;
bool cmp(node a,node b)
{
	return fabs(a.x-b.x)<eps?(a.y<b.y):(a.x<b.x);
}
double xielv(int a,int b)//求斜率
{
	return (f[a].y-f[b].y)/(f[a].x-f[b].x);
}
double dis(int a,int b)//求距离
{
	double xx=f[a].x-f[b].x,yy=f[a].y-f[b].y;
	return sqrt(xx*xx+yy*yy);
}
void calans()//计算最短周长
{
	int i,j,now=0,maxi;
	double ans=0;
	while(1)
	{
		double max=-4000000;//斜率最大值
		maxi=now;
		for(i=1;i<n;i++)
		{
			if(flag[i])continue;
			if(-f[i].x+f[now].x>eps)continue;//在当前点的左边，不考虑
			if(fabs(f[now].x-f[i].x)<eps)
			{
				if(f[now].y>f[i].y)//在当前点的正下方，不考虑
				  continue;
		        maxi=i;
		        break;
			}
			else
			{
				double p=xielv(now,i);
				if(p>max)
	            {
	            	maxi=i;
	            	max=p;
	            }
			}
		}
		flag[maxi]=1;
		ans+=dis(now,maxi);
		if(maxi==now)
		{
		  flag[maxi]=0;//最后右边的点要计算两次，去除标记
		  break;
		}
		now=maxi;
	}
	now=0;
	while(1)//每次取斜率最小的点
	{
		double min=4000000;
		int mini=now;
		for(i=1;i<n;i++)
		{
			if(flag[i])continue;
			if(-f[i].x+f[now].x>eps)continue;//在当前点的左边，不考虑

			if(fabs(f[i].x-f[now].x)<eps)
			{
				if(f[i].x+0.1<f[n-1].x)//在当前点的正下方，不考虑

				   continue;
                mini=i;
                break;
			}
				double p=xielv(now,i);
				if(p<min)
	            {
	            	mini=i;
	            	min=p;
	            }
		}
		ans+=dis(now,mini);
		flag[mini]=1;
	    if(mini==now)
	       break;
		now=mini;
	}
	if(n==2)ans/=2;//n=2时，除以二
	printf("%.2f\n",ans);
}
int main()
{
	int i,j;
	while(scanf("%d",&n)!=-1&&n)
	{
		memset(flag,0,sizeof(flag));
		for(i=0;i<n;i++)
		  scanf("%lf%lf",&f[i].x,&f[i].y);
        sort(f,f+n,cmp);
	    calans();
	}
	return 0;
}