HDU 3694 Fermat Point in Quadrangle (计算几何- 四边形的费马点)
来源:互联网 发布:webstorm破解版 mac 编辑:程序博客网 时间:2024/06/06 08:39
Fermat Point in Quadrangle
Problem Description
In geometry the Fermat point of a triangle, also called Torricelli point, is a point such that the total distance from the three vertices of the triangle to the point is the minimum. It is so named because this problem is first raised by Fermat in a private letter. In the following picture, P0 is the Fermat point. You may have already known the property that:
Alice and Bob are learning geometry. Recently they are studying about the Fermat Point.
Alice: I wonder whether there is a similar point for quadrangle.
Bob: I think there must exist one.
Alice: Then how to know where it is? How to prove?
Bob: I don’t know. Wait… the point may hold the similar property as the case in triangle.
Alice: It sounds reasonable. Why not use our computer to solve the problem? Find the Fermat point, and then verify your assumption.
Bob: A good idea.
So they ask you, the best programmer, to solve it. Find the Fermat point for a quadrangle, i.e. find a point such that the total distance from the four vertices of the quadrangle to that point is the minimum.
Alice and Bob are learning geometry. Recently they are studying about the Fermat Point.
Alice: I wonder whether there is a similar point for quadrangle.
Bob: I think there must exist one.
Alice: Then how to know where it is? How to prove?
Bob: I don’t know. Wait… the point may hold the similar property as the case in triangle.
Alice: It sounds reasonable. Why not use our computer to solve the problem? Find the Fermat point, and then verify your assumption.
Bob: A good idea.
So they ask you, the best programmer, to solve it. Find the Fermat point for a quadrangle, i.e. find a point such that the total distance from the four vertices of the quadrangle to that point is the minimum.
Input
The input contains no more than 1000 test cases.
Each test case is a single line which contains eight float numbers, and it is formatted as below:
x1 y1 x2 y2 x3 y3 x4 y4
xi, yi are the x- and y-coordinates of the ith vertices of a quadrangle. They are float numbers and satisfy 0 ≤ xi ≤ 1000 and 0 ≤ yi ≤ 1000 (i = 1, …, 4).
The input is ended by eight -1.
Each test case is a single line which contains eight float numbers, and it is formatted as below:
x1 y1 x2 y2 x3 y3 x4 y4
xi, yi are the x- and y-coordinates of the ith vertices of a quadrangle. They are float numbers and satisfy 0 ≤ xi ≤ 1000 and 0 ≤ yi ≤ 1000 (i = 1, …, 4).
The input is ended by eight -1.
Output
For each test case, find the Fermat point, and output the total distance from the four vertices to that point. The result should be rounded to four digits after the decimal point.
Sample Input
0 0 1 1 1 0 0 11 1 1 1 1 1 1 1-1 -1 -1 -1 -1 -1 -1 -1
Sample Output
2.82840.0000
Source
2010 Asia Fuzhou Regional Contest
Recommend
axubiao
解题思路:对于凸四边形,费马点一定在对角线的连线上,对于凹多边形,费马点在那个凹点上,证明的话,画张图很好证明。
View Code
Problem : 3694 ( Fermat Point in Quadrangle ) Judge Status : Accepted
RunId : 9336690 Language : C++ Author : a1061747415
Code Render Status : Rendered By HDOJ C++ Code Render Version 0.01 Beta
RunId : 9336690 Language : C++ Author : a1061747415
Code Render Status : Rendered By HDOJ C++ Code Render Version 0.01 Beta
#include <iostream>#include <cmath>#include <cstdio>#include <algorithm>using namespace std;struct point{ double x,y; friend bool operator < (point p1,point p2){ if(p1.y!=p2.y) return p1.y<p2.y; else return p1.x<p2.x; } double getdis(point p0){ return sqrt((x-p0.x)*(x-p0.x)+(y-p0.y)*(y-p0.y)); } friend istream& operator >> (istream &input,point &p0){ input>>p0.x>>p0.y; return input; }};point p[10],pp[10];int n=4,top;double getans(int k){ double dis=0; for(int i=0;i<4;i++){ dis+=pp[i].getdis(pp[k]); //cout<<p[i].x<<" "<<p[i].y<<" "<<p[k].x<<" "<<p[k].y<<" "<<dis<<endl; } return dis;}double xchen(point a,point b,point c){ return (b.x-a.x)*(c.y-a.y)-(c.x-a.x)*(b.y-a.y);}double dchen(point a,point b,point c){ return (b.x-a.x)*(c.x-a.x)+(c.y-a.y)*(b.y-a.y);}bool cmp(point a,point b){ if(fabs(xchen(p[0],a,b))<1e-7) return a.getdis(p[0])<b.getdis(p[0]); else return xchen(p[0],a,b)>0;}bool deal(){ top=1; sort(p,p+n); sort(p+1,p+n,cmp); for(int i=2;i<n;i++){ while(top>0 && xchen(p[top-1],p[top],p[i])<=0) top--; p[++top]=p[i]; } if(top>=3) return true; else return false; }void computing(){ //for(int i=0;i<4;i++) cout<<p[i].x<<" "<<p[i].y<<endl; double ans=1e9,tmp; if(deal()){ tmp=p[0].getdis(p[2])+p[1].getdis(p[3]); if(tmp<ans) ans=tmp; }else{ for(int i=0;i<4;i++){ tmp=getans(i); if(tmp<ans) ans=tmp; } } printf("%.4f\n",ans);}int main(){ while(cin>>p[0]>>p[1]>>p[2]>>p[3]){ bool flag=true; for(int i=0;i<4;i++){ pp[i]=p[i]; if(fabs(p[i].x+1)>1e-9) flag=false; if(fabs(p[i].y+1)>1e-9) flag=false; } if(flag) break; computing(); } return 0;}
- HDU 3694 Fermat Point in Quadrangle (计算几何- 四边形的费马点)
- HDU 3694 Fermat Point in Quadrangle(四边形的费马点)
- HDU 3694 Fermat Point in Quadrangle (费马定理求四边形的费马点)
- Fermat Point in Quadrangle POJ 3990 四边形的费马点 数学
- hdu 3694 Fermat Point in Quadrangle
- hdu 3694 Fermat Point in Quadrangle (数学)
- HDU 3694 Fermat Point in Quadrangle (数学-费马点)
- HDU 3694 Fermat Point in Quadrangle (数学-费马点)
- poj 3990 Fermat Point in Quadrangle / hdu 3694
- Hdu 3694 Fermat Point in Quadrangle(三分)
- HDU 3694Fermat Point in Quadrangle(三分 不过网上大多是费马点的题解)
- hdu3694 Fermat Point in Quadrangle 费马点
- HDOJ 3694 Fermat Point in Quadrangle
- poj 3990 Fermat Point in Quadrangle 凸包和费马点
- HDU3694 Fermat Point in Quadrangle 多边形费马点结论
- 2010 Asia Fuzhou Regional Contest HDOJ 3694 Fermat Point in Quadrangle
- UVALive 5102 Fermat Point in Quadrangle 极角排序+找距离二维坐标4个点最近的点
- HDU3694Fermat Point in Quadrangle
- 电子邮件恢复2--PST文件查看与修复
- Cocos2D-X渲染框架之CCSprite
- scanf()和printf()
- gson 使用
- CMFCPropertyGridCtrl 的简单教程
- HDU 3694 Fermat Point in Quadrangle (计算几何- 四边形的费马点)
- HQL语句总结
- yum配置与使用(详解)
- C/c++面试题一
- linux man使用方法
- Centos6.3 安装pptp vpn
- poj3150
- 网上报名系统(浅读)
- 切词分词之判断字符串是否能够拆分为字段中的单词