计算几何-基本定义（编译通过，性能未测试）

函数列表如下

小于判定 点A（x1，y1） B（x2，y2） 当x1<x2 ||（x1==x2 && y1<y2） 时 return true；
int dcmp(double x)//判断x是正，负，或者0；
double degtorad(double x)//角度转弧度
double radtodeg(double x)//弧度转角度
double arad1(double x)//将一个弧度角转换到 [0,2pi) 内
double arad2(double x)//将一个弧度角转换到 [-pi,pi) 内
double adeg1(double x)//将一个角度制角转换到 [0,360) 内
double adeg2(double x)//将一个角度制角转换到 [-180,180) 内
double polarangle(Vector a)//极角 以弧度制返回
double dot(Vector A,Vector B)//点积（数量积）
double vlen(Vector A)//矢量长度
double angle(Vector A,Vector B)//两矢量夹角
double cross(Vector A,Vector B)//叉积（失量积）
double area3(point A,point B,point C)//三角形面积（三点式）
double arean(point *p,int n)//n边形面积（注意 ***数组 p 编号 从 0 开始***）
Vector rotate(Vector A,double rad)//矢量旋转 rad是弧度
Vector normal(Vector A)//求单位法向量
//直线的参数方程 ： A+(B-A)*t t∈R；射线 t∈ (0,+oo); 线段 t∈(0,1);
point line_line_point(point p,Vector v,point q,Vector w)//求两直线交点 （参数方程传入）
double dis_point_line(point p,point a,point b)//点 p 到直线 ab 的距离
double dis_point_segment(point p,point a,point b)//点 到 线段 最近距离
point point_line_projection(point p,point a,point b)//点 在 直线 的 投影
bool is_segment_intersection(point a1,point a2,point b1,point b2)//判断 是否 两直线（不允许在端点处） 相交；如果相交 返回true；
bool is_on_segment(point p,point a1,point a2)//判断一个点是否在直线上

/*计算几何模板---基本定义*//*0.基本定义*/struct point{    double x,y;    point (double x=0,double y=0):x(x),y(y){};};typedef point Vector;Vector operator +(Vector A,Vector B){return Vector(A.x+B.x,A.y+B.y);}Vector operator -(Vector A,Vector B){return Vector(A.x-B.x,A.y-B.y);}Vector operator *(Vector A,double p){return Vector(A.x*p,A.y*p);}Vector operator /(Vector A,double p){return Vector(A.x/p,A.y/p);}bool operator < (const point &a,const point &b){    return a.x<b.x||(a.x==b.x&&a.y<b.y);}int dcmp(double x){    if(fabs(x)<eps)return 0;    return x<0?-1:1;}bool operator ==(const point &a,const point &b){    return dcmp(a.x-b.x)==0&&dcmp(a.y-b.y)==0;}double const pi=3.14159265358979323846264;/*end-0*//*基本运算*//*1.角度运算*/double degtorad(double x)//角度转弧度{    return x/180.0*pi;}double radtodeg(double x)//弧度转角度{    return x/pi*180;}double arad1(double x)//将一个弧度角转换到 [0,2pi) 内{    return x-floor(x/2.0/pi)*2*pi;}double arad2(double x)//将一个弧度角转换到 [-pi,pi) 内{    double ang=arad1(x);    return ang<pi?ang:ang-2*pi;}double adeg1(double x)//将一个角度制角转换到 [0,360) 内{    return x-floor(x/360.0)*360.0;}double adeg2(double x)//将一个角度制角转换到 [-180,180) 内{    double ang=adeg1(x);    return ang<180?ang:ang-360;}double polarangle(Vector a)//极角 以弧度制返回{    return atan2(a.y,a.x);}/*end-1*//*2.数量积*/double dot(Vector A,Vector B){return A.x*B.x+A.y*B.y;}double vlen(Vector A){return sqrt(dot(A,A));}double angle(Vector A,Vector B){return acos(dot(A,B)/vlen(A)/vlen(B));}/*end-2*//*3.矢量积&面积*/double cross(Vector A,Vector B){return A.x*B.y-A.y*B.x;}double area3(point A,point B,point C){return cross(B-A,C-A)/2.0;}//三角形有向面积double arean(point *p,int n)//n边形面积（注意 ***数组 p 编号 从 0 开始***）{    double s=0;    point p0(0,0);    fr(i,0,n-1)    s+=cross(p0,p[i]);    return s/2;}/*end-3*//*4.矢量旋转*/Vector rotate(Vector A,double rad)//矢量旋转 rad是弧度{    return Vector(A.x*cos(rad)-A.y*sin(rad),A.x*sin(rad)+A.y*cos(rad));}Vector normal(Vector A)//求单位法向量{    double l=vlen(A);    return Vector(-A.y/l,A.x/l);}/*end-4*//*5.点与直线*///直线的参数方程 ： A+(B-A)*t t∈R；射线 t∈ (0,+oo); 线段 t∈(0,1);point line_line_point(point p,Vector v,point q,Vector w)//求两直线交点 （参数方程传入）{    Vector u=p-q;    double t=cross(w,u)/cross(v,w);    return p+v*t;}double dis_point_line(point p,point a,point b)//点 p 到直线 ab 的距离{    Vector v1=b-a,v2=p-a;    return fabs(cross(v1,v2))/vlen(v1);}double dis_point_segment(point p,point a,point b)//点 到 线段 最近距离{    if(a==b)return vlen(p-a);    Vector v1=b-a,v2=p-a,v3=p-b;    if(dcmp(dot(v1,v2))<0)return vlen(v2);    else if(dcmp(dot(v1,v3))>0)return vlen(v3);    else return fabs(cross(v1,v2))/vlen(v1);}point point_line_projection(point p,point a,point b)//点 在 直线 的 投影{    Vector v=b-a;    return a+v*(dot(v,p-a)/dot(v,v));}bool is_segment_intersection(point a1,point a2,point b1,point b2)//判断 是否 两直线（不允许在端点处） 相交；如果相交 返回true；{    double c1=cross(a2-a1,b1-a1),c2=cross(a2-a1,b2-a1),c3=cross(b2-b1,a1-b1),c4=cross(b2-b1,a2-b1);    return dcmp(c1)*dcmp(c2)<0&&dcmp(c3)*dcmp(c4)<0;}bool is_on_segment(point p,point a1,point a2)//判断一个点是否在直线上{    return dcmp(cross(a1-p,a2-p))==0&&dcmp(dot(a1-p,a2-p))<0;}/*end-5*/