• <ins id="pjuwb"></ins>
    <blockquote id="pjuwb"><pre id="pjuwb"></pre></blockquote>
    <noscript id="pjuwb"></noscript>
          <sup id="pjuwb"><pre id="pjuwb"></pre></sup>
            <dd id="pjuwb"></dd>
            <abbr id="pjuwb"></abbr>
            簡單的幾何題,先把經緯度換算成球面坐標,再把球面坐標換算成直角坐標,然后求夾角,乘半徑得到球面距離


            /*************************************************************************
            Author: WHU_GCC
            Created Time: 2007-10-2 10:18:27
            File Name: pku3407.cpp
            Description: 
            ***********************************************************************
            */

            #include 
            <iostream>
            #include 
            <cmath>
            using namespace std;
            #define out(x) (cout<<#x<<": "<<x<<endl)
            const int maxint=0x7FFFFFFF;
            typedef 
            long long int64;
            const int64 maxint64 = 0x7FFFFFFFFFFFFFFFLL;
            template
            <class T>void show(T a, int n){for(int i=0; i<n; ++i) cout<<a[i]<<' '; cout<<endl;}
            template
            <class T>void show(T a, int r, int l){for(int i=0; i<r; ++i)show(a[i],l);cout<<endl;}

            const double pi = acos(-1.0);

            int main()
            {
                
            double theta, phi, theta_m, phi_m;
                
            double theta1, phi1, theta2, phi2;
                
            char s1[10], s2[10];
                
            while (scanf("%lf%lf%s%lf%lf%s"&theta, &theta_m, s1, &phi, &phi_m, s2) != EOF)
                
            {
                    theta 
            += theta_m / 60;
                    phi 
            += phi_m / 60;
                    
            if (strcmp(s1, "N"== 0)
                        theta1 
            = 90 - theta;
                    
            else
                        theta1 
            = 90 + theta;
                    
            if (strcmp(s2, "E"== 0)
                        phi1 
            = phi;
                    
            else
                        phi1 
            = 360 - phi;
                    scanf(
            "%lf%lf%s%lf%lf%s"&theta, &theta_m, s1, &phi, &phi_m, s2);
                    theta 
            += theta_m / 60;
                    phi 
            += phi_m / 60;
                    
                    
            if (strcmp(s1, "N"== 0)
                        theta2 
            = 90 - theta;
                    
            else
                        theta2 
            = 90 + theta;
                    
            if (strcmp(s2, "E"== 0)
                        phi2 
            = phi;
                    
            else
                        phi2 
            = 360 - phi;

                    theta1 
            = theta1 * pi / 180;
                    theta2 
            = theta2 * pi / 180;
                    phi1 
            = phi1 * pi / 180;
                    phi2 
            = phi2 * pi / 180;

                    
            double p = 6370.0;
                    
            double x1, y1, z1, x2, y2, z2;
                    x1 
            = p * sin(theta1) * cos(phi1);
                    y1 
            = p * sin(theta1) * sin(phi1);
                    z1 
            = p * cos(theta1);

                    x2 
            = p * sin(theta2) * cos(phi2);
                    y2 
            = p * sin(theta2) * sin(phi2);
                    z2 
            = p * cos(theta2);
                    
                    
            double alpha = acos((x1 * x2 + y1 * y2 + z1 * z2) / ((sqrt(x1 * x1 + y1 * y1 + z1 * z1)) * (sqrt(x2 * x2 + y2 * y2 + z2 * z2))));
                    printf(
            "%.3lf\n", alpha * p);
                }

                
            return 0;
            }
            posted on 2007-10-02 17:55 Felicia 閱讀(633) 評論(1)  編輯 收藏 引用 所屬分類: 計算幾何
            Comments
             
            久久99久久99精品免视看动漫| 人人狠狠综合久久亚洲婷婷| 久久人妻少妇嫩草AV无码蜜桃| 久久国产视频网| 2020久久精品亚洲热综合一本| 色综合久久久久综合体桃花网| 99国产精品久久久久久久成人热| 成人免费网站久久久| 欧美大战日韩91综合一区婷婷久久青草 | 99久久免费国产精品特黄| 色偷偷88888欧美精品久久久| 色综合合久久天天综合绕视看| 性做久久久久久久久| 狠狠色丁香婷婷综合久久来| 一本久久综合亚洲鲁鲁五月天亚洲欧美一区二区 | 99久久精品免费看国产一区二区三区 | 久久久艹| 99精品国产在热久久| 久久久久亚洲AV无码专区首JN | 久久精品国产第一区二区三区| 久久久久亚洲AV成人网人人网站| 国产精品久久久久jk制服| 久久伊人精品一区二区三区| 久久久久国产一区二区三区| 人人狠狠综合久久亚洲婷婷| aaa级精品久久久国产片| 一本久久a久久精品亚洲| 伊人久久大香线蕉综合5g| 99久久国产综合精品五月天喷水| 一本一本久久A久久综合精品| 一日本道伊人久久综合影| 久久人人爽人爽人人爽av| 久久黄视频| 人妻中文久久久久| 久久亚洲中文字幕精品一区四| 国产香蕉97碰碰久久人人| 精品综合久久久久久97超人| 久久亚洲国产中v天仙www| 久久综合中文字幕| 91精品免费久久久久久久久| 久久精品嫩草影院|