UVa 11817 - Tunnelling the Earth

2023-10-15 02:24:01

将平面上的亮点坐标转化成空间坐标，为球面上的两点。求所对应的弧长与弦长的差值，转化为空间坐标的函数为

double torad(double deg) {return deg / 180 * M_PI;
}
void get_coord(double lat, double lng, double& x, double& y, double& z) {lat = torad(lat);lng = torad(lng);x = R * cos(lat) * cos(lng);y = R * cos(lat) * sin(lng);z = R * sin(lat);
}

/*************************************************************************> File Name: 11817.cpp> Author: AcToy> Mail: ycsgldy@163.com > Created Time: 2013年07月19日 星期五 19时45分23秒************************************************************************/#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include using namespace std;
typedef unsigned int u32;
typedef long long i64;
typedef unsigned long long u64;
typedef vector IV;
typedef vector BV;
typedef pair II;
typedef vector IIV;
#define For(t,v,c) for(t::const_iterator v=c.begin(); v!=c.end(); ++v)
const int INF = 0x7FFFFFFF;
const double eps = 1E-10;
const double PI = acos(-1);
const double R = 6371009;
double torad(double deg) {return deg / 180 * M_PI;
}
void get_coord(double lat, double lng, double& x, double& y, double& z) {lat = torad(lat);lng = torad(lng);x = R * cos(lat) * cos(lng);y = R * cos(lat) * sin(lng);z = R * sin(lat);
}
double sqr(double x) { return x * x; }
double dist(double x1, double y1, double z1, double x2, double y2, double z2) {return sqrt(sqr(x1 - x2) + sqr(y1 - y2) + sqr(z1 - z2));
}
int main() {int Case;double lat1, lng1, lat2, lng2, d1, d2;double x1, y1, z1, x2, y2, z2;scanf("%d", &Case);while(Case--) {scanf("%lf%lf%lf%lf", &lat1, &lng1, &lat2, &lng2);get_coord(lat1, lng1, x1, y1, z1);get_coord(lat2, lng2, x2, y2, z2);d1 = dist(x1, y1, z1, x2, y2, z2);d2 = 2 * asin(d1 / (2 * R)) * R;printf("%.0lf\n", d2 - d1);}return 0;
}

本文来自互联网用户投稿，文章观点仅代表作者本人，不代表本站立场，不承担相关法律责任。如若转载，请注明出处。 如若内容造成侵权/违法违规/事实不符，请点击【内容举报】进行投诉反馈！

标签：技术

上一篇 > hdoj4405-Aeroplane chess
下一篇 > UVa 406-Prime Cuts

Duilib中list控件支持ctrl和shif多行选中的实现

[ICML2015]Batch Normalization:Accelerating Deep Network Training by Reducing Internal Covariate Shif

win10系统微软输入法于eclipse ctrl+shif+f冲突间接处理办法

Codeforces Round #259 (Div. 2) B. Little Pony and Sort by Shif

读LDD3，内存映射与DMA--PAGE_SHIF…

VMware虚拟机安装XP【要先分区，再设置BOOT 启动CD，shif+上移】

更换iBus五笔的左与右Shif

sublime ctrl+shif+f 没用解决办法

idea 对 ctrl + z 的撤销是 ctrl + shif + z

计算机最早的设计师应用于,计算机应用基础选择题doc.doc

win10自带截图神器：Win+Shift+S

Python基础之文件目录操作

python简述目录_Python基础之文件目录操作(示例代码)

tp5 如何做数据采集

任务2-7(服务器字体+阿里巴巴矢量库)

html标签（1)：h1~h6,p,br,pre,hr

TI 电量计介绍与芯片选型指南

几款TI电源芯片简介

TI DSP芯片C2000系列读取FLASH数据

德州仪器(Ti)平台嵌入式开发基础

TI三相电机智能栅极驱动芯片特点分类

省选模拟（12.08） T3 圈圈圈圈圈圈圈圈

Hadoop生态圈技术栈（上）

大数据开发基础入门与项目实战（三）Hadoop核心及生态圈技术栈之6.Impala交互式查询

小猿圈之Linux下Mysql 操作命令

大数据Hadoop生态圈常用面试题

大数据开发基础入门与项目实战（三）Hadoop核心及生态圈技术栈之4.Hive DDL、DQL和数据操作

备战Noip2018模拟赛11（B组）T3 Monogatari 物语

【智能优化算法-圆圈搜索算法】基于圆圈搜索算法Circle Search Algorithm求解单目标优化问题附matlab代码

NYOJ 78 圈水池

递归问题跑道汽车绕圈问题 Python实现

Hadoop生态圈（三）：MapReduce

UVa 11817 - Tunnelling the Earth

相关文章