几何模板

struct Point
 {
     double x,y;
     friend Point operator - (Point a,Point b)
     {
         Point temp;
         temp.x = a.x - b.x;
         temp.y = a.y - b.y;
         return temp;
     }
 };
 
 Point p1,p2,p3,pc,pp;
 double r;
  
Point get_pc1(Point p1, Point p2, Point p3)    //已知三点，求圆心坐标
{
  Point p;
  if(fabs(p1.y-p2.y) > eps)      //因为计算过程中有出现(y2-y1)做分母的情况，所以这里分了两种情况讨论
  {
  //      double t1 = p3.x*p3.x - p1.x*p1.x+p3.y*p3.y-p1.y*p1.y-((p3.y-p1.y)*(p1.x*p1.x-p2.x*p2.x+p1.y*p1.y-p2.y*p2.y)) / (p2.y-p1.y);
    //    double t2 = 2.0*(p3.x-p1.x)-(2*(p3.y-p1.y)*(p2.x-p1.x)) / (p2.y-p1.y);
         double  t1 = (p2.x*p2.x-p1.x*p1.x+p2.y*p2.y-p1.y*p1.y) *(p3.y-p1.y) - (p2.y-p1.y)*(p3.x*p3.x-p1.x*p1.x+p3.y*p3.y-p1.y*p1.y);
         double t2 = 2.0*((p3.y-p1.y)*(p2.x-p1.x) - (p2.y-p1.y)*(p3.x-p1.x));
        p.x = t1 / t2;
        p.y = (p2.x*p2.x-p1.x*p1.x+p2.y*p2.y-p1.y*p1.y-2*(p2.x-p1.x)*p.x) / (2*(p2.y-p1.y));
     }
     else
     {
         p.x = (p2.x*p2.x-p1.x*p1.x+p2.y*p2.y-p1.y*p1.y) / (2*(p2.x-p1.x));
         p.y = (p3.x*p3.x-p1.x*p1.x+p3.y*p3.y-p1.y*p1.y-2*(p3.x-p1.x)*p.x) / (2*(p3.y-p1.y));
     }
     return p;
}