poj 3990 四边形费马点

到网上查一下，四边形的费马点就：

    如果是凸四边形，就是对角线的交点。

    如果是凹四边形，就是该凹点。

本来很简单，由于粗心，先是这样想的：枚举四个点的距离，再算出对角线的距离，取最小就行了，这样是错的，必须要判断凸包问题，在分别求。

#include<algorithm>
#include<cstdio>
#include<cmath>
#include<cstring>
using namespace std ;
struct node{
	double x,y; 
}p[4],rec[4];

double Min;
double dis(node a,node b){
	return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}
bool cross(node a,node b,node c){
	return (a.x-c.x)*(b.y-c.y)>=(a.y-c.y)*(b.x-c.x);
}
bool cmp(node a,node b){
	if(a.y!=b.y) return a.y<b.y;
	return a.x<b.x;
}

int graham(int n)
{
	int i,pos=1;
	sort(p,p+n,cmp);
	if(n==0) return 0;
	rec[0]=p[0];
	if(n==1) return 1;
	rec[1]=p[1];
	if(n==2) return 2;
	rec[2]=p[2];
	for(i=2;i<n;i++){
		while(pos&&cross(p[i],rec[pos],rec[pos-1])) pos--;
		rec[++pos]=p[i];
	}
	int len=pos;
	rec[++pos]=p[n-2];
	for(i=n-3;i>=0;i--){
		while(pos!=len&&cross(p[i],rec[pos],rec[pos-1])) pos--;
		rec[++pos]=p[i];
	}
	return pos;
}

void make()
{
	double d;
	Min=9999999;
	for(int i=0;i<4;i++){
		d=0;
		for(int j=0;j<4;j++){
			d+=dis(p[i],p[j]);
		}
		if(d<Min) Min=d;
	}
}
int main()
{
   double t1,t2;
   int op;
   while(1){
		op=0;
		for(int i=0;i<4;i++){
			scanf("%lf %lf",&p[i].x,&p[i].y);
			if(p[i].x==-1&&p[i].y==-1) op++;
		}
		if(op==4) break;
		int cnt=graham(4);
	    if(cnt==4){
		 Min=0;
	     Min+=dis(rec[0],rec[2]);
	     Min+=dis(rec[1],rec[3]);
		}	
		else make();   
	    printf("%.4f\n",Min);	//又忘记了poj不能用lf，贡献了一次WA！
   }
 	return 0 ;
}