计算几何ACM——一些代码整理

实数的gcd

double gcd(double x,double y)
{
    while(fabs(x)>eps && fabs(y)>eps)//此处的eps根据需要取
    {
        if(x>y)
            x-=floor(x/y)*y;
        else
            y-=floor(y/x)*x;
    }
    return x+y;
}

求向量叉积

double cross(const Point& p1,const Point& p2,const Point& q1,const Point& q2)
{
    return (p2.x-p1.x)*(q2.y-q1.y)-(p2.y-p1.y)*(q2.x-q1.x);
}//vector cross product

极角排序的一种实现方式

bool cmp(const Point& a,const Point& b)
{
    if(a.y==0 && b.y==0 && a.x*b.x<=0)
        return a.x>b.x;
    if(a.y==0 && a.x>=0 && b.y!=0)
        return true;
    if(b.y==0 && b.x>=0 && a.y!=0)
        return false;
    if(b.y*a.y<=0)
        return a.y>b.y;
    return cross(O,a,O,b)>0 || (cross(O,a,O,b)==0 && a.x<b.x);

}

得到垂心

Point get_circumcenter(Point a,Point b,Point c)
{
    double a1=b.x-a.x,b1=b.y-a.y,c1=(a1*a1+b1*b1)/2;
    double a2=c.x-a.x,b2=c.y-a.y,c2=(a2*a2+b2*b2)/2;
    double d=a1*b2-a2*b1;
    return Point(a.x+(c1*b2-c2*b1)/d,a.y+(a1*c2-a2*c1)/d);
}