程序師世界是廣大編程愛好者互助、分享、學習的平台,程序師世界有你更精彩!
首頁
編程語言
C語言|JAVA編程
Python編程
網頁編程
ASP編程|PHP編程
JSP編程
數據庫知識
MYSQL數據庫|SqlServer數據庫
Oracle數據庫|DB2數據庫
 程式師世界 >> 編程語言 >> C語言 >> C++ >> 關於C++ >> UVA 11817 (球面幾何)

UVA 11817 (球面幾何)

編輯:關於C++

這題主要知道如何把經緯度轉化為三維點坐標即可

代碼:

 

#include 
#include 
#include 
#include 
using namespace std;

const double PI = acos(-1.0);
int t;
double lat[2], lng[2];

//角度轉化為弧度
double torad(double deg) {
    return deg / 180 * PI;
}

struct Point3 {
    double x, y, z;
    Point3() {}
    Point3(double x, double y, double z) {
        this->x = x;
        this->y = y;
        this->z = z;
    }
};

const double R = 6371009;

Point3 get_coord(double R, double lat, double lng) {
    lat = torad(lat);
    lng = torad(lng);
    return Point3(R * cos(lat) * cos(lng), R * cos(lat) * sin(lng), R * sin(lat));
}

double dist(Point3 p, Point3 q) {
    double dx = p.x - q.x;
    double dy = p.y - q.y;
    double dz = p.z - q.z;
    return sqrt(dx * dx + dy * dy + dz * dz);
}

int main() {
    scanf("%d", &t);
    while (t--) {
        scanf("%lf%lf%lf%lf", &lat[0], &lng[0], &lat[1], &lng[1]);
        Point3 p[2];
        p[0] = get_coord(R, lat[0], lng[0]);
        p[1] = get_coord(R, lat[1], lng[1]);
        double d = dist(p[0], p[1]);
        printf("%.0f\n", 2 * asin(d / 2 / R) * R - d);
    }
    return 0;
}


 

  1. 上一頁:
  2. 下一頁:
Copyright © 程式師世界 All Rights Reserved