求两个圆公切线的模板

最新推荐文章于 2021-09-22 22:14:21 发布

tzteyang

最新推荐文章于 2021-09-22 22:14:21 发布

阅读量415

点赞数

分类专栏：计算几何

原文链接：https://www.cnblogs.com/1625--H/p/12567362.html

版权

计算几何专栏收录该内容

8 篇文章 0 订阅

订阅专栏

// a[i] 和 b[i]存放第 i 条公切线与 圆A 和 园B 的交点
// 返回值为切线的条数 如果没有切线返回值为-1
int getTangents(circle A, circle B, Point*a, Point *b){
    int cnt = 0;
    // 以A为半径更大的那个圆进行计算
    if(A.r < B.r) return getTangents(B, A, b, a);
    db d2 = (A.p-B.p).len2();  // 圆心距平方
    db rdiff = A.r - B.r;		// 半径差
    db rsum = A.r + B.r;		//半径和
    if(d2 < rdiff * rdiff) return 0; 	// 情况1，内含,没有公切线
    Vector AB = B.p - A.p;				// 向量AB，其模对应圆心距
    db base = atan2(AB.y, AB.x);		// 求出向量AB对应的极角
    if(d2 == 0 && A.r == B.r) return -1;// 情况3，两个圆重合，无限多切线
    if(d2 == rdiff * rdiff){ 			// 情况2，内切，有一条公切线
        a[cnt] = A.point(base);			
        b[cnt] = B.point(base);cnt++;
        return 1;
    }
    // 求外公切线
    db ang = acos((A.r - B.r) / sqrt(d2)); //求阿尔法
    // 两条外公切线
    a[cnt] = A.point(base+ang); b[cnt] = B.point(base+ang); cnt++;
    a[cnt] = A.point(base-ang); b[cnt] = B.point(base-ang); cnt++;
    if(d2 == rsum * rsum){  // 情况5，外切，if里面求出内公切线
        a[cnt] = A.point(base); b[cnt] = B.point(pi+base); cnt++;
    }
    else if(d2 > rsum * rsum){	//情况6，相离，再求出内公切线
        db ang = acos((A.r + B.r) / sqrt(d2));
        a[cnt] = A.point(base + ang); b[cnt] = B.point(pi+base+ang);cnt++;
        a[cnt] = A.point(base - ang); b[cnt] = B.point(pi+base-ang);cnt++;
    }
    // 此时，d2 < rsum * rsum 代表情况 4 只有两条外公切线
    return cnt;
}