栏目分类:
子分类:
返回
名师互学网用户登录
快速导航关闭
当前搜索
当前分类
子分类
实用工具
热门搜索
名师互学网 > IT > 面试经验 > 面试问答

zoj 1868 Airline Hub

面试问答 更新时间: 发布时间: IT归档 最新发布 模块sitemap 名妆网 法律咨询 聚返吧 英语巴士网 伯小乐 网商动力

zoj 1868 Airline Hub

#include <set>#include <map>#include <queue>#include <stack>#include <math.h>#include <stdio.h>#include <stdlib.h>#include <iostream>#include <limits.h>#include <string.h>#include <string>#include <algorithm>#define MId(x,y) ( ( x + y ) >> 1 )#define L(x) ( x << 1 )#define R(x) ( x << 1 | 1 )#define BUG puts("here!!!")using namespace std;const int MAX = 1010;const double pi = acos(-1.0);const double inf = 1e30;struct point { double lng, lat, x,y;};point p[MAX];const double eps = 1e-6;bool dy(double x,double y){return x > y + eps;}bool xy(double x,double y){return x < y - eps;}bool dyd(double x,double y){ return x > y - eps;}bool xyd(double x,double y){return x < y + eps;} bool dd(double x,double y) {return fabs( x - y ) < eps;}  double torad(double x){return x*pi/180;}double angle_3d(double lng1, double lat1, double lng2, double lat2)  {    return acos(cos(lat1)*cos(lat2)*cos(lng1 - lng2) + sin(lat1)*sin(lat2));  }  double dis(int k,int i){double ang = angle_3d(p[i].lng, p[i].lat, p[k].lng, p[k].lat);return ang;}int solve(int n){double mmin = inf;int t;for(int i=0; i<n; i++){double mmax = 0;for(int k=0; k<n; k++){if( k == i ) continue;if( dy(dis(k,i), mmax) )mmax = dis(k,i);}if( xyd(mmax, mmin) ){mmin = mmax;t = i;}}return t;}int main(){int n;double x, y;while( ~scanf("%d", &n) ){for(int i=0; i<n; i++){scanf("%lf %lf", &x, &y);p[i].lat = torad(x);p[i].lng = torad(y);p[i].x = x; p[i].y = y;}int ans = solve(n);printf("%.2lf %.2lfn", p[ans].x, p[ans].y);}return 0;}
转载请注明:文章转载自 www.mshxw.com
本文地址:https://www.mshxw.com/it/376116.html
我们一直用心在做
关于我们 文章归档 网站地图 联系我们

版权所有 (c)2021-2022 MSHXW.COM

ICP备案号:晋ICP备2021003244-6号