tzc3109 两圆外公切线之交

题目链接：tzc3109

方法：计算几何

思想：当两个圆的半径相同或者两圆关系为内含时不存在外公切线交点。
当存在外公切线交点时，由于知道两半径，而外公切线与两圆心够成的直线
及圆心与切点的连线所够成的图形中，得到交点到两圆圆心的向量之比
等于半径之比，从而(x-o1.x,y-o1.y)=k(x-o2.x,y-o2.y)。可惜这向量的关系
很久没用了，被我给忘记了，用距离+斜率联立方程，既麻烦，又容易出错。

代码：

#include <iostream>
#include <cmath>
using  namespace std;

const double eps=1e-10;
struct point { double x,y; };
struct circle
{ 
	point centre;
	double r;
};
double dist(point a,point b)   
{   
    return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));   
}  

int main()
{
	circle o1,o2;
	int n;
	cin>>n;
	while(n--)
	{
		cin>>o1.centre.x>>o1.centre.y>>o1.r;
		cin>>o2.centre.x>>o2.centre.y>>o2.r;
		if(fabs(o1.r-o2.r)>=dist(o1.centre,o2.centre)//两圆内含 
			||fabs(o1.r-o2.r)<eps) //半径相等 
		{   printf("Error\n");  continue; } 
		double k=o1.r/o2.r;
		//向量公式 
		// (x-o1.x,y-o1.y)=k(x-o2.x,y-o2.y)
		double x=(k*o2.centre.x-o1.centre.x)/(k-1);
		double y=(k*o2.centre.y-o1.centre.y)/(k-1);
		printf("%.2lf %.2lf\n",x,y);
	}
	return 0;
}