poj 3990 Fermat Point in Quadrangle / hdu 3694

// 求四边形费马点，一开始想到的是用模拟退火， 写的不好， WA了两次， 只好去看了一下四边形费马点的规律
//对于三角形：当三个角都小于120度，则费马点在三角形内部与任意两定点的连线构成的角都为120度；
//若存在一个角大于等于120度，则费马点为此角顶点。
//对于平面四边形：若为凸四边形，则费马点为两对角线交点；对于凹四边形，为其凹顶点。
//在不是凸四边形的情况下，只要枚举四个顶点即可。

#include<iostream>

#include<cstdio>
#include<cmath>
#include<algorithm>
using namespace std;
const double EP=1e-8;
const double INF=1e8;
struct Point{double x, y;};
Point p[4], pt[4], c;
double x_mult(Point sp, Point ep, Point op){
    return (sp.x-op.x)*(ep.y-op.y)-(sp.y-op.y)*(ep.x-op.x);
}
double dist(Point p1, Point p2){
    return sqrt((p1.x-p2.x)*(p1.x-p2.x)+(p1.y-p2.y)*(p1.y-p2.y));
}
bool cmp(const Point p1, const Point p2){
    return p1.y<p2.y ||(p1.y==p2.y&&p1.x<p2.x);
}

int graham(Point pp[], Point res[]){
    int i, len1, k=0, top=1;
    sort(pp,pp+4,cmp);
    res[0]=pp[0];
    res[1]=pp[1];
    res[2]=pp[2];
    for(i=2; i<4; i++){
        while(top&&x_mult(pp[i],res[top],res[top-1])> - EP)
            top--;
        res[++top]=pp[i];
    }
    len1=top; res[++top]=pp[2];
    for(i=1; i>=0; i--){
        while(top!=len1&&x_mult(pp[i],res[top],res[top-1])> - EP)
            top--;
            res[++top]=pp[i];
    }
    return top;      
}
Point intersect(Point u1,Point u2,Point v1,Point v2){
      Point res=u1;
      double t=((u1.x-v1.x)*(v1.y-v2.y)-(u1.y-v1.y)*(v1.x-v2.x))/
               ((u1.x-u2.x)*(v1.y-v2.y)-(u1.y-u2.y)*(v1.x-v2.x));
     res.x+=(u2.x-u1.x)*t;
     res.y+=(u2.y-u1.y)*t;
     return res;
}

int main(){
    //freopen("1.txt", "r", stdin);
    while(1){
        double ans=0, s;
        int i, flag, j;
        for(i=0; i<4; i++)
        scanf("%lf%lf", &p[i].x, &p[i].y);
        if(p[0].x<0)break;
        int sum=graham(p, pt);
        if(sum==4){
            c=intersect(pt[0], pt[2], pt[1], pt[3]);
            for(i=0; i<4; i++)
            ans+=dist(p[i], c);
        }
        else {
            ans=INF;
            for(i=0; i<4; i++){
                s=0;
                for(j=0; j<4; j++)
                s+=dist(p[i], p[j]);
                if(s<ans)
                    ans=s;
            }
        }
        printf("%.4f\n", ans);
    }
    return 0;
}