zoj 1199 Point of Intersection 求两个圆公切线的交点

来源:互联网 发布:7u分享网络注册 编辑:程序博客网 时间:2024/05/21 22:35

题目来源:

http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemId=199

分析:

两个圆,不存在交点的情况。

1)在 半径相同的情况下, 两条公切线 是平行的, 不存在交点。

2)当一个圆在另一个圆内时, 不存在交点。

当存在交点时, 计算 交点的坐标。 我们用三角形的相似,计算出x,y 值。

x  =  (p2.x * r1 - p1.x * r2 ) / (r1 - r2) ,

y =  (p2.y * r1 - p1.y * r2) / (r1 - r2)) 

 

代码如下:

double add(double a, double b){    return (fabs(a+b) < EPS * (fabs(a) + fabs(b))) ? 0 : (a + b) ;}struct Point{    double x, y ;    Point(){}    Point(double x, double y):x(x),y(y){}} ;double dist(Point a, Point b){    return sqrt(add( (a.x - b.x)*(a.x - b.x) ,(a.y - b.y)*(a.y - b.y) )) ;}double r1 ,r2 ;Point p1, p2;bool is_insection(){    if(r1 == r2)        return 0 ;    if(dist(p1, p2) <= fabs(r1 - r2))        return 0 ;    return 1;}Point tangle_p(){    return Point((p2.x * r1 - p1.x * r2 ) / (r1 - r2) , (p2.y * r1 - p1.y * r2) / (r1 - r2)) ;}int main(){    int t;    cin>>t;    while(t--){        scanf("%lf%lf%lf%lf%lf%lf", &p1.x , &p1.y,  &r1 ,&p2.x , &p2.y , &r2) ;        if(is_insection()){            Point p = tangle_p();            printf("%.2lf %.2lf\n" , p.x , p.y) ;        }        else printf("Impossible.\n");    }}

 

0 0