UVa 10075 - Airlines（球面距离+Floyd）

这题要求的是球面的航线距离，不是直线距离，而是球面的弧线距离

我先把球面坐标转换为三维坐标，然后计算出两点间距离，再用余弦定理求出扇形角度，从而求出弧线距离

本题的结果要转换成int，否则会有误差

#include <stdio.h>
#include <string.h>
#include <math.h>
const double pi=3.141592653589793;
const double R=6378;
typedef struct
{
    char s[25];
    double x;
    double y;
    double z;
}air;
air a[110];
int d[110][110],inf;
double dis(air a,air b)
{
    return (a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y)+(a.z-b.z)*(a.z-b.z);
}
int ld(air a,air b)
{
    double dist=dis(a,b);
    double ang=acos((R*R*2.0-dist)/(2.0*R*R));
    return (int)(R*ang+0.5);
}
int ff(char *s,int n)
{
    for(int i=0;i<n;i++)
    {
        if(strcmp(a[i].s,s)==0)
            return i;
    }
}
int main()
{
    int n,m,q,cas=0;
    double l,r;
    char st[25],ed[25];

    while(scanf("%d%d%d",&n,&m,&q)==3)
    {
        if(!n&&!m&&!q)return 0;
        if(cas)printf("\n");
        memset(d,0x3f,sizeof(d));
        inf=d[0][0];
        for(int i=0;i<n;i++)
        {
           scanf("%s %lf %lf",a[i].s,&l,&r);
           l=l*pi/180.0;
           r=r*pi/180.0;
           a[i].x=R*cos(l)*cos(r);
           a[i].y=R*cos(l)*sin(r);
           a[i].z=R*sin(l);
        }
        for(int i=0;i<m;i++)
        {
            scanf("%s%s",st,ed);
            int A=ff(st,n);
            int B=ff(ed,n);
            d[A][B]=ld(a[A],a[B]);
        }
        for(int k=0;k<n;k++)
        {
            for(int i=0;i<n;i++)
            {
                for(int j=0;j<n;j++)
                {
                    if(d[i][k]+d[k][j]<d[i][j])
                     d[i][j]=d[i][k]+d[k][j];
                }
            }
        }
        printf("Case #%d\n",++cas);
        for(int i=0;i<q;i++)
        {
            scanf("%s%s",st,ed);
            if(d[ff(st,n)][ff(ed,n)]==inf)puts("no route exists");
            else printf("%d km\n",d[ff(st,n)][ff(ed,n)]);
        }
    }
    return 0;
}