poj 3990 hdu 3694

来源：互联网发布：淘宝充值号码查询编辑：程序博客网时间：2024/05/22 14:36

题目意思是求四边形的费马点。

这里要严谨一点了，随便一个随机贪心算法会wa。

若四边形为凸边形，则费马点为对角线之交点。

若四边形为凹边形，则费马点为四边形中的一个顶点。

以下是代码：

#include<cstdio>#include<iostream>#include<utility>#include<string>#include<algorithm>#include<cmath>using namespace std;const int M=15; const double inf=1e250;struct point{double x,y;bool operator < (const point a)const{if(y==a.y)return x<a.x;return y<a.y;}}p[M],stac[M];int n,top;double dis(point a,point b){return sqrt((b.x-a.x)*(b.x-a.x)+(b.y-a.y)*(b.y-a.y));}double xmult(double x1,double y1,double x2,double y2){return x1*y2-x2*y1;}double cross(point p0,point p1,point p2){double x1=p1.x-p0.x,y1=p1.y-p0.y;double x2=p2.x-p0.x,y2=p2.y-p0.y;return xmult(x1,y1,x2,y2);}void covex_hull(){sort(p,p+n);stac[0]=p[0];stac[1]=p[1];top=1;for(int i=2;i<n;i++){while(top>0 && cross(stac[top],p[i],stac[top-1])<=0)top--;stac[++top]=p[i];}int temp=top;for(int i=n-2;i>=0;i--){while(top>temp && cross(stac[top],p[i],stac[top-1])<=0)top--;stac[++top]=p[i];}}int main(){while(scanf("%lf%lf%lf%lf%lf%lf%lf%lf",&p[0].x,&p[0].y,&p[1].x,&p[1].y,&p[2].x,&p[2].y,&p[3].x,&p[3].y)==8){if(p[0].x==-1 && p[0].y==-1 && p[1].x==-1 && p[1].y==-1 && p[2].x==-1 && p[2].y==-1 && p[3].x==-1 && p[3].y==-1)break;n=4;covex_hull();double ans;if(top==4){ans=dis(stac[0],stac[2])+dis(stac[1],stac[3]);}else{ans=inf;ans=min(ans,dis(p[0],p[1])+dis(p[0],p[2])+dis(p[0],p[3]));ans=min(ans,dis(p[1],p[0])+dis(p[1],p[2])+dis(p[1],p[3]));ans=min(ans,dis(p[2],p[1])+dis(p[2],p[0])+dis(p[2],p[3]));ans=min(ans,dis(p[3],p[1])+dis(p[3],p[2])+dis(p[3],p[0]));}printf("%.4lf\n",ans);}return 0;}