【POJ2587】Airline Hub(经纬度求距离)

最新推荐文章于 2017-08-07 17:39:22 发布

oranges_c

最新推荐文章于 2017-08-07 17:39:22 发布

阅读量437

点赞数

CC 4.0 BY-SA版权

分类专栏： poj 文章标签： poj

本文链接：https://blog.youkuaiyun.com/oranges_c/article/details/54381402

poj 专栏收录该内容

15 篇文章

订阅专栏

本文介绍了一道算法题的解决思路，该题要求找出在一组经纬度坐标中，找到一个点使得它到其它所有点的最大距离最小。文章提供了一个C++实现方案，并使用了球面三角学来计算地球上两点之间的距离。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

记录一个菜逼的成长。。

题目链接
题目大意：
给你n个点的经纬度，求一个点到其他点的距离的最大值里是最小。

具体解释
需要用到经纬度求距离的公式。
当个模板好了。

#include <cstdio>
#include <cmath>
#include <algorithm>
using namespace std;
int INF = 0x3f3f3f3f;
const double pi = acos(-1);
const int maxn = 1000 + 10;
struct node{
    double x,y;
}a[maxn];
double cal(double d)
{
    return d * pi / 180.0;
}
int main()
{
    int n;
    while(~scanf("%d",&n)){
        for( int i = 0; i < n; i++ ){
            scanf("%lf%lf",&a[i].x,&a[i].y);
        }
        double dis = INF * 100;
        double x,y;
        for( int i = 0; i < n; i++ ){
            double tmp = 0;
            for( int j = 0; j < n; j++ ){
                if(i != j){
                    double t1 = cal(a[i].x),t2 = cal(a[j].x);
                    double c = t1 - t2;
                    double b = cal(a[i].y) - cal(a[j].y);
                    double s = 2 * asin(sqrt(pow(sin(c/2),2) + cos(t1) * cos(t2) * pow(sin(b/2),2)));
                    s = s * 6378.137;
                    tmp = max(tmp,s);
                }
            }
            if(tmp < dis){
                dis = tmp;
                x = a[i].x;
                y = a[i].y;
            }
        }
        printf("%.2f %.2f\n",x,y);
    }
    return 0;
}