到网上查一下,四边形的费马点就:
如果是凸四边形,就是对角线的交点。
如果是凹四边形,就是该凹点。
本来很简单,由于粗心,先是这样想的:枚举四个点的距离,再算出对角线的距离,取最小就行了,这样是错的,必须要判断凸包问题,在分别求。
#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 ;
}