URAL 1294

原以为是会分分钟搞定的。。结果也是折腾了半个小时

首先，明明列出了余弦公式，还先去求一下角度，这样就不大好了。

一开始看到int类型10000以内， 于是就像用int做，但是想到如果是这样的话中间过程可能会数据过多。

这个判断impossible情况真是比我想象中还要方便，直接对这个公式进行分析就可以了

#include<cmath>
#include<iostream>
#include<algorithm>
#include<cstdio>
#include<string>
using namespace std;
const double eps=1e-10;
double a1,a2,b1,b2;

int dcmp(double x)
{
    if(x>eps)return 1;
    if(x<-eps)return -1;
    return 0;
}

int main()
{
    scanf("%lf%lf%lf%lf",&a1,&a2,&b1,&b2);
    double down=a1*a2-b1*b2,up=a1*a2*(b1*b1+b2*b2)-b1*b2*(a1*a1+a2*a2);
    if(!dcmp(down)||!dcmp(up)){
        cout<<"Impossible.\n";
        return 0;
    }
    double ans=up/down;
    if(dcmp(ans)<0){
        cout<<"Impossible.\n";
        return 0;
    }
    printf("Distance is %.0f km.\n",1000*sqrt(ans));
    return 0;
}