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

来源：互联网发布：视频剪切软件中文版编辑：程序博客网时间：2024/04/29 22:33

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

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

本题的结果要转换成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;}