URAL 1294 Mars Satellites(圆内接三角形)

URAL 1294 Mars Satellites(圆内接三角形)

题意:

       圆边上按顺序有4个不同的点A,B,C,D,现在已知AC,AD,BC,BD的长度,问你边CD的长度是多少?

分析:

       假设AC,AD,BC,BD长分别为a1,b1,a2,b2.那么根据余弦定理可以得到:  (a1^2+b1^2-c^2)/2a1b1= (a2^2+b2^2-c^2)/2a2b2.

       化简可以得到:

       c^2= {a1b1*(a2^2+b2^2)-a2b2(a1^2+b1^2) }/(a1b1-a2b2)

       所以当分母==0 或c^2<=0时 都是不合法的值.

       最终输出的值需要*1000,因为a1单位是KKm.

AC代码:

#include<cstdio>
#include<cmath>
using namespace std;
const double eps=1e-10;
int dcmp(double x)
{
    if(fabs(x)<eps) return 0;
    return x<0?-1:1;
}

int main()
{
    double a1,b1,a2,b2;
    while(scanf("%lf%lf%lf%lf",&a1,&b1,&a2,&b2)==4)
    {
        if(dcmp(a1*b1-a2*b2)==0)
        {
            printf("Impossible.\n");
            continue;
        }
        double c2= (a1*b1*(a2*a2+b2*b2)-a2*b2*(a1*a1+b1*b1)) / (a1*b1-a2*b2);
        if(dcmp(c2)<=0) printf("Impossible.\n");
        else printf("Distance is %.0lf km.\n",sqrt(c2)*1000.0 );
        //如果上面用 printf("Distance is %d km.\n",int(sqrt(c2)*1000.0) );
        //就是错的,可能int溢出,也可能精度出现问题
    }
    return 0;
}