【圆的反演变换】hdu4773

这个blog解释的挺不错:http://blog.csdn.net/acdreamers/article/details/16966369

给定不相交的两个圆以及圆外一点，找一个经过给定点的圆与其他两个圆相切

首先来看反演变换，首先是给定一个圆圆心为O，半径为R

1、圆外一点P与圆内一点P‘会一一对应的反演OP*OP'=R*R

2、经过O的圆，反演后成为不经过O的一条直线

3、不经过O的圆，反演后成为另一个圆，且圆心并不对应

4、不经过O的直线反演后成为一个经过O的圆

5、过O的直线反演后不变

那么这道题就把两个圆反演之后求两圆的公切线，然后反演回去，就会成为一个过O的圆，且与另两圆相切（由反演的过程可看出）

求反演后的圆的圆心，可以将过O的直线所在的那条直径两端点反演，然后再求圆心

这道题貌似不需要求内公切线，且反演中心圆的半径要设大一点，否则会有精度问题

#include <iostream>
#include <cstdio>	
#include <cstdlib>
#include <cstring>
#include <algorithm>
#include <cmath>
#define sqr(x) ((x)*(x))
using namespace std;
const double eps=1e-10;
struct Point{
	double x,y,z;
	Point operator +(const Point &b) const;
	Point operator *(const double &b) const;
	Point operator -() cons