已知圆外一点，圆心和半径，求过圆外点的直线与圆的切点算法

CPoint CalcQieDian(CPoint ptCenter, CPoint ptOutside, double dbRadious) {     struct point {double x, y;}; point E,F,G,H; double r=dbRadious; //1. 坐标平移到圆心ptCenter处,求园外点的新坐标E E.x= ptOutside.x-ptCenter.x; E.y= ptOutside.y-ptCenter.y; //平移变换到E  //2. 求园与OE的交点坐标F, 相当于E的缩放变换 double t= r / sqrt (E.x * E.x + E.y * E.y);  //得到缩放比例 F.x= E.x * t;   F.y= E.y * t;   //缩放变换到F //3. 将E旋转变换角度a到切点G，其中cos(a)=r/OF=t, 所以a=arccos(t); double a=acos(t);   //得到旋转角度 G.x=F.x*cos(a) -F.y*sin(a); G.y=F.x*sin(a) +F.y*cos(a);    //旋转变换到G //4. 将G平移到原来的坐标下得到新坐标H H.x=G.x+ptCenter.x; H.y=G.y+ptCenter.y;             //平移变换到H     //5. 返回H return CPoint(int(H.x),int(H.y)); //6. 实际应用过程中，只要一个中间变量E,其他F,G,H可以不用。}

应该说方法有很多种，几何的，代数的，做的过程中也在纸上运算了好长时间，可能也是这些知识很久没用了，生疏很多。解出来的公式都比较复杂，无意间在网上找到这个方法，没经过具体的笔算验证，但是求出来的结果是我预判范围的，暂定为计算结果是正确的。在这边共享下，有需要的童鞋不妨自己验证一下。

对了，这个函数只能返回一个点，理论上有两个切点的，只要把15行改成如下：

 double a=-acos(t);   //得到旋转角度