二维几何模板

最新推荐文章于 2021-06-16 21:01:39 发布
最新推荐文章于 2021-06-16 21:01:39 发布
阅读量820 收藏 1
点赞数
分类专栏： 模板 文章标签： 二维计算几何 模板 点与向量 二维图形
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/sinianluoye/article/details/44203215
版权

终于完成了，先挂函数列表，然后是代码

所在行函数（结构）声明功能
23int sign(double d)返回浮点数的符号
24bool isint(double d)判断一个浮点数是不是整数
25bool Equal(double a,double b)判断两个浮点数是否相等
28struct Point点定义
13NAP不是一个点
34bool isnap(const Point& P)判断一个点是不是NAP
40typedef Point Vector;向量定义
42Vector operator + (Vector A,Vector B)向量相加
43Vector operator - (Vector A,Vector B)向量相减
44Vector operator * (Vector A,double t)数乘向量
45Vector operator * (double t,Vector A)数乘向量
46Vector operator / (Vector A,double t)向量除数
49bool operator < (const Point& a,const Point& b)坐标从左到右从下到上的比较
50bool operator == (const Point& a,const Point& b)坐标从左到右从下到上的比较
51bool operator > (const Point& a,const Point& b)坐标从左到右从下到上的比较
52bool operator <= (const Point& a,const Point& b)坐标从左到右从下到上的比较
53bool operator >= (const Point& a,const Point& b)坐标从左到右从下到上的比较
55double Dot(const Point& A,const Point& B)点积
56double Cross(const Point& A,const Point& B)叉积
57double Length(const Point& A)模
58double Length(const Point& A,const Point& B)两点间距
59double Distance(const Point& A,const Point& B)两点间距
60double Angle(const Point& A,const Point& B)两向量夹角
61double Area2(const Point& A,const Point& B,const Point& C)三角形有向面积
63Vector Rotate(const Vector& A,double rad)向量绕起点旋转rad弧度
64Vector Normal(const Vector& A)单位法向量
65Vector Unit(const Vector& A)单位向量
67bool Collinear(const Point& A,const Point& B,const Point& C)判断三点是否共线
72struct Segment线段定义
77struct Line直线定义
14NAS不是一个线段（直线）
84bool isnas(const Segment& S)判断是否为NAS
91bool Online(const Point& P,const Line& L)点在直线上
92bool Online(const Point& P,const Segment& L)点在线段上 不包括端点
93bool Online_ex(const Point& P,const Segment& L)点在线段上 包括端点
95int Side(const Point& P,Line L)点对于直线的位置
96int Side(const Point& P,Segment L)点对于直线的位置
97bool SameSide(const Point& P1,const Point& P2,const Line& L)两点在直线同侧判定 不包含在直线上
98bool SameSide(const Point& P1,const Point& P2,const Segment& L)两点在线段同侧判定 不包含在线段所在直线上
99bool SameSide_ex(const Point& P1,const Point& P2,const Line& L)两点在直线同侧判定 包含在直线上
100bool SameSide_ex(const Point& P1,const Point& P2,const Segment& L)两点在线段同侧判定 包含在线段所在直线上
103bool Parallel(const Line& L1,const Line& L2)平行或重合判定
104Equal(const Line& L1,const Line& L2)重合判定
105bool operator == (const Line& L1,const Line& L2)重合判定
106bool Vertical(const Line& L1,const Line& L2)垂直判定
109bool Parallel(const Segment& L1,const Segment& L2)线段平行或共线
110bool Collinear(const Segment& L1,const Segment& L2)判定共线
111bool Equal(const Segment& L1,const Segment& L2)线段相等判定
112bool operator == (const Segment& L1,const Segment& L2)线段相等判定
113bool Intersec(const Segment& L1,const Segment& L2)判断两线段是否相交 不包括端点
114bool Intersec_ex(const Segment& L1,const Segment& L2)判断两线段是否相交 包括端点
122double Distance(const Line& L,const Point& P)点到直线距离
123double Distance(const Segment& S,const Point& P)点到线段距离
131double Length(const Segment& S)线段长度
133Point Intersection(const Line& L1,const Line& L2)两直线交点
140Point Intersection(const Segment& L1,const Segment& L2)两线段交点 不包含端点
146Point Intersection_ex(const Segment& L1,const Segment& L2)两线段交点 包含端点
153Point LineSymmetry(const Line&symmetryline,const Point& P)点关于直线的对称点
160Point LineSymmetry(const Segment&symmetryline,const Point& P)点关于线段的对称点
167Line LineSymmetry(const Line&symmetryline,const Line& L)直线关于直线的对称线
170Segment LineSymmetry(const Line&symmetryline,const Segment& L)线段关于直线的对称线
173Line LineSymmetry(const Segment&symmetryline,const Line& L)直线关于线段的对称线
176Segment LineSymmetry(const Segment&symmetryline,const Segment& L)线段关于线段的对称线
179Point PointSymmetry(const Point&symmetrypoint,const Point& P)中心对称
181Segment PointSymmetry(const Point&Symmetrypoint,const Segment& L)线段关于点的中心对称
184Line PointSymmetry(const Point&Symmetrypoint,const Line& L)直线关于点的中心对称
189Point Projection(const Line& L,const Point& P)点在直线上的投影
196Point Projection(const Segment& L,const Point& P)点对线段的投影
203Segment Projection(const Line& L,const Segment& S)线段在直线上的投影
210Point Midpoint(const Segment& S)中点
212Line Perpendicular(const Line& S,const Point& P)过定点的垂线
213Line PerpendicularBisector(const Segment& S)垂直平分线
215Line NormalLine(const Line& S,const Point& P)过定点垂线
216Line Midnormal(const Segment& S)垂直平分线
218Line AngleBisector(const Line& L1,const Line& L2)角平分线
232double Area(const Point& A,const Point& B,const Point& C)三角形面积
234Point Circumcenter(const Point& A,const Point& B,const Point& C)三角形外心
237Point Incenter(const Point& A,const Point& B,const Point& C)三角形内心
241Point Perpencenter(const Point& A,const Point& B,const Point& C)三角形垂心
244Point Barycentry(const Point& A,const Point& B,const Point& C)三角形重心
247struct Triangle三角形定义
254double Area(const Triangle& T)三角形面积
255Point Circumcenter(const Triangle& T)三角形外心
256Point Incenter(const Triangle& T)三角形内心
257Point Perpencenter(const Triangle& T)三角形垂心
258Point Barycentry(const Triangle& T)三角形重心
260int Type(const Triangle& T)返回三角形的类型
263bool InTriangle_ex(const Point& P,const Triangle& T)返回点是否在三角形内部（含边）
266bool InTriangleLine(const Point& P,const Triangle& T)返回点是否在三角形边上（含端点）
269bool InTriangle(const Point& P,const Triangle& T)返回点是否在三角形内部(不含边)
276struct Polygon多边形定义
312double Area(const Polygon& P)多边形的面积
322Polygon ConvexHull_ex(Polygon P)凸包扫描算法(含平角)
343int ConvexHull_ex(int n,Point *P,Point *ans)凸包扫描算法（含平角）
351Polygon ConvexHull(Polygon P)凸包扫描算法(不含平角)
372int ConvexHull(int n,Point *P,Point *ans)凸包扫描算法（不含平角）
380bool IsConvex(const Polygon& P)判断是否为凸多边形 （凸多边形定义包含平角）
389bool IsConvex(int n,Point *P)判断是否为凸多边形
394int InPolygon(const Point& P,const Polygon& POLY)点是否在多边形内部
412struct Circle圆的定义
419double Area(const Circle& C)圆的面积
421int Position(const Point& P,const Circle& C)返回点与圆的位置关系
426int Position(const Line& L,const Circle& C)返回直线与圆的位置关系
431Segment Intersection(const Line& L,const Circle& C)返回直线与圆的割线
453int Position(const Circle& C1,const Circle& C2)圆与圆的位置关系
471Segment Intersection(const Circle& C1,const Circle& C2)圆与圆的公共弦

欢迎纠错帮忙改进，转载请注明出处 

/*************************************二维几何模板***************************************************/


/*************************** lc思念落叶（sinianluoye、JLU_Lichuang） 2015-3-10***********************/


#include <cmath>
#include <cstdio>
#include <iostream>
#include <cstdlib>
#include <algorithm>


#define eps 1e-8
#define sqr(x) ((x)*(x))
#define NAP Point(NAN,NAN)
#define NAS Segment(NAP,NAP)




using namespace std;


/***1.点与向量的基本定义***/


const double pi=acos(-1);


int sign(double d){if(fabs(d)<eps)return 0;return (d<0)?-1:1;} ///返回浮点数的符号
bool isint(double d){return fabs(floor(d+0.5)-d)<eps;} ///判断一个浮点数是不是整数
bool Equal(double a,double b){return fabs(a-b)<eps;}///判断两个浮点数是否相等




struct Point
{
    double x,y;
    Point(double x=0,double y=0):x(x),y(y){}
};


bool isnap(const Point& P)
{
    return isnan(P.x)||isnan(P.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 t){return Vector(A.x*t,A.y*t);}///数乘向量
Vector operator * (double t,Vector A){return Vector(A.x*t,A.y*t);}///数乘向量
Vector operator / (Vector A,double t){return Vector(A.x/t,A.y/t);}///向量除数


/**以下比较基于点的坐标从左到右从下到上**/
bool operator < (const Point& a,const Point& b){if(sign(a.x-b.x))return a.x<b.x;return a.y<b.y;}
bool operator == (const Point& a,const Point& b){return !(sign(a.x-b.x)||sign(a.y-b.y));}
bool operator > (const Point& a,const Point& b){return !((a<b)||(a==b));}
bool operator <= (const Point& a,const Point& b){return !(a>b);}
bool operator >= (const Point& a,const Point& b){return !(a<b);}


double Dot(const Point& A,const Point& B){return A.x*B.x+A.y*B.y;}///点积
double Cross(const Point& A,const Point& B){return A.x*B.y-A.y*B.x;}///叉积
double Length(const Point& A){return sqrt(Dot(A,A));}///模
double Length(const Point& A,const Point& B){return Length(A-B);}///两点间距
double Distance(const Point& A,const Point& B){return Length(A-B);}///两点间距
double Angle(const Point& A,const Point& B){return acos(Dot(A,B)/Length(A)/Length(B));}///两向量夹角
double Area2(const Point& A,const Point& B,const Point& C){return 0.5*Cross(B-A,C-A);}///三角形的有向面积


Vector Rotate(const Vector& A,double rad){return Vector(A.x*cos(rad)-A.y*sin(rad),A.x*sin(rad)+A.y*cos(rad));}///向量绕起点旋转rad弧度
Vector Normal(const Vector& A){double L=Length(A);return Vector(-A.y/L,A.x/L);}///单位法向量
Vector Unit(const Vector& A){double L=Length(A);return Vector(A.x/L,A.y/L);}///单位向量


bool Collinear(const Point& A,const Point& B,const Point& C){return !sign(Cross(A-B,A-C));}///判断三点是否共线


/***2.直线射线线段相关***/
///a,b为线段的两端点（直线上的两点） v为方向向量


struct Segment
{
    Point a,b,v;
    Segment(Point a=Point(),Point b=Point()):a(a),b(b),v(b-a){}
};///线段
struct Line
{
    Point a,b,v;
    Line(Point a=Point(),Point b=Point()):a(a),b(b),v(b-a){}
    Line(Segment S):a(S.a),b(S.b),v(S.v){}
};///直线


bool isnas(const Segment& S)
{
    return isnap(S.a)||isnap(S.b);
}


///点与直线线段的位置关系


bool Online(const Point& P,const Line& L){return !sign(Cross(P-L.a,L.v));}///点在直线上
bool Online(const Point& P,const Segment& L){return (Online(P,Line(L))&&(sign(Dot(P-L.a,L.v))>0)&&(sign(Dot(L.b-P,L.v))>0));}///点在线段上 不包括端点
bool Online_ex(const Point& P,const Segment& L){return (Online(P,Line(L))&&(sign(Dot(P-L.a,L.v))>=0)&&(sign(Dot(L.b-P,L.v))>=0));}///点在线段上 包括端点


int Side(const Point& P,Line L){if(L.a>L.b)swap(L.a,L.b);return sign(Area2(L.a,L.b,P));}///返回点对于直线的位置 1（左（上）边） 0（在直线上） -1在（右（下）边）
int Side(const Point& P,Segment L){if(L.a>L.b)swap(L.a,L.b);return sign(Area2(L.a,L.b,P));}///返回点对于直线的位置 1（左（上）边） 0（在线段所在直线上） -1在（右（下）边）
bool SameSide(const Point& P1,const Point& P2,const Line& L){return Side(P1,L)*Side(P2,L)==1;}///两点在直线同侧判定 不包含在直线上
bool SameSide(const Point& P1,const Point& P2,const Segment& L){return Side(P1,L)*Side(P2,L)==1;}///两点在线段同侧判定 不包含在线段所在直线上
bool SameSide_ex(const Point& P1,const Point& P2,const Line& L){return Side(P1,L)*Side(P2,L)>=0;}///两点在直线同侧判定 包含在直线上
bool SameSide_ex(const Point& P1,const Point& P2,const Segment& L){return Side(P1,L)*Side(P2,L)>=0;}///两点在线段同侧判定 包含在线段所在直线上


///直线位置关系
bool Parallel(const Line& L1,const Line& L2){return !sign(Cross(L1.v,L2.v));}///平行和重合都返回1 注意重合！！！
bool Equal(const Line& L1,const Line& L2){Line T(L1.a,L2.a);return Parallel(L1,L2)&&Parallel(L1,T);}///重合判定
bool operator == (const Line& L1,const Line& L2){return Equal(L1,L2);}
bool Vertical(const Line& L1,const Line& L2){return !sign(Dot(L1.v,L2.v));}///垂直判定


///线段位置关系
bool Parallel(const Segment& L1,const Segment& L2){return !sign(Cross(L1.v,L2.v));}///线段平行或共线都返回1 注意共线！！！
bool Collinear(const Segment& L1,const Segment& L2){return Parallel(L1,L2)&&Collinear(L1.a,L1.b,L2.a);}///判定共线
bool Equal(const Segment& L1,const Segment& L2){return L1.a==L2.a&&L1.b==L2.b;}///线段相等判定
bool operator == (const Segment& L1,const Segment& L2){return Equal(L1,L2);}
bool Intersec(const Segment& L1,const Segment& L2){return ((Side(L1.a,L2)*Side(L1.b,L2))<0)&&(Side(L2.a,L1)*Side(L2.b,L1)<0);}///判断两线段是否相交不包括端点
bool Intersec_ex(const Segment& L1,const Segment& L2)
{
    int s1=Side(L1.a,L2),s2=Side(L1.b,L2),s3=Side(L2.a,L1),s4=Side(L2.b,L1);
    return (s1||s2)&&(s3||s4);
}///判断两线段是否相交包括端点


///交点、对称点、距离


double Distance(const Line& L,const Point& P){return fabs(Cross(P-L.a,L.v)/Length(L.v));}///点到直线距离
double Distance(const Segment& S,const Point& P)///点到线段距离
{
    if(sign(Dot(P-S.a,S.v))<=0)
        return Distance(P,S.a);
    if(sign(Dot(P-S.b,S.v))<=0)
        return Distance(P,S.b);
    return fabs(Cross(P-S.a,S.v)/Length(S.v));
}
double Length(const Segment& S){return Length(S.a,S.b);}///线段长度


Point Intersection(const Line& L1,const Line& L2)///两直线交点 请先判断是否平行
{
    Vector u=L1.a-L2.a;
    double t=Cross(L2.v,u)/Cross(L1.v,L2.v);
    return L1.a+L1.v*t;
}


Point Intersection(const Segment& L1,const Segment& L2)///两线段交点 不包含端点 不相交的线段返回（NAN,NAN）
{
    if(Intersec(L1,L2))
        return Intersection(Line(L1),Line(L2));
    return NAP;
}
Point Intersection_ex(const Segment& L1,const Segment& L2)///两线段交点 包含端点 不相交的线段返回（NAN,NAN）
{
    if(Intersec_ex(L1,L2))
        return Intersection(Line(L1),Line(L2));
    return NAP;
}


Point LineSymmetry(const Line&symmetryline,const Point& P)///点关于直线的对称点
{
    Vector u=P-symmetryline.a,v=Normal(symmetryline.v),w=symmetryline.v;
    double t=Cross(w,u)/Cross(v,w);
    return P+2*t*v;
}


Point LineSymmetry(const Segment&symmetryline,const Point& P)///点关于线段的对称点
{
    Vector u=P-symmetryline.a,v=Normal(symmetryline.v),w=symmetryline.v;
    double t=Cross(w,u)/Cross(v,w);
    return P+2*t*v;
}


Line LineSymmetry(const Line&symmetryline,const Line& L)///直线关于直线的对称线
{return Line(LineSymmetry(symmetryline,L.a),LineSymmetry(symmetryline,L.b));}


Segment LineSymmetry(const Line&symmetryline,const Segment& L)///线段关于直线的对称线
{return Segment(LineSymmetry(symmetryline,L.a),LineSymmetry(symmetryline,L.b));}


Line LineSymmetry(const Segment&symmetryline,const Line& L)///直线关于线段的对称线
{return Line(LineSymmetry(symmetryline,L.a),LineSymmetry(symmetryline,L.b));}


Segment LineSymmetry(const Segment&symmetryline,const Segment& L)///线段关于线段的对称线
{return Segment(LineSymmetry(symmetryline,L.a),LineSymmetry(symmetryline,L.b));}


Point PointSymmetry(const Point&symmetrypoint,const Point& P){return 2*symmetrypoint-P;}///中心对称


Segment PointSymmetry(const Point&Symmetrypoint,const Segment& L)///线段关于点的中心对称
{return Segment(PointSymmetry(Symmetrypoint,L.a),PointSymmetry(Symmetrypoint,L.b));}


Line PointSymmetry(const Point&Symmetrypoint,const Line& L)///直线关于点的中心对称
{return Segment(PointSymmetry(Symmetrypoint,L.a),PointSymmetry(Symmetrypoint,L.b));}






Point Projection(const Line& L,const Point& P)///点在直线上的投影
{
    Vector u=P-L.a,v=Normal(L.v),w=L.v;
    double t=Cross(w,u)/Cross(v,w);
    return P+t*v;
}


Point Projection(const Segment& L,const Point& P)///点对线段的投影（不一定在线段上）
{
    Vector u=P-L.a,v=Normal(L.v),w=L.v;
    double t=Cross(w,u)/Cross(v,w);
    return P+t*v;
}


Segment Projection(const Line& L,const Segment& S)///线段在直线上的投影
{
    return Segment(Projection(L,S.a),Projection(L,S.b));
}




///线段相关
Point Midpoint(const Segment& S){return 0.5*(S.a+S.b);}///中点


Line Perpendicular(const Line& S,const Point& P){return Line(P,P+Normal(S.v));}///过定点的垂线
Line PerpendicularBisector(const Segment& S){return Perpendicular(S,Midpoint(S));}///垂直平分线


Line NormalLine(const Line& S,const Point& P){return Line(P,P+Normal(S.v));}///过定点垂线(上边的名太长了 这个好记些 下同)
Line Midnormal(const Segment& S){return Perpendicular(S,Midpoint(S));}///垂直平分线


Line AngleBisector(const Line& L1,const Line& L2)///角平分线
///注意这里L1 L2为射线 方向分别为 L1.v即(L1.b-L1.a),L2.v即(L2.b-L2.a) 函数输入只需保证方向正确且L1，L2不平行即可
{
    Point P=Intersection(L1,L2);
    Vector V=Unit(L1.v)+Unit(L2.v);
    if(Length(V)>eps)
    return Line(P,P+V);
    return NormalLine(L1,L1.a);
}


/********************************二维图形相关**********************************************/


///三角形


double Area(const Point& A,const Point& B,const Point& C){return fabs(Area2(A,B,C));}///三角形面积


Point Circumcenter(const Point& A,const Point& B,const Point& C)
{return Intersection(Midnormal(Segment(A,B)),Midnormal(Segment(A,C)));}///三角形外心


Point Incenter(const Point& A,const Point& B,const Point& C)
{return Intersection(AngleBisector(Line(A,B),Line(A,C)),AngleBisector(Line(B,C),Line(B,A)));}///三角形内心




Point Perpencenter(const Point& A,const Point& B,const Point& C)
{return Intersection(Perpendicular(Line(A,B),C),Perpendicular(Line(A,C),B));}/// 三角形垂心


Point Barycentry(const Point& A,const Point& B,const Point& C)
{return (A+B+C)/3.0;}///三角形重心


struct Triangle
{
    Point A,B,C;
    Triangle(Point A=Point(),Point B=Point(),Point C=Point())
    :A(A),B(B),C(C){}
};


double Area(const Triangle& T){return Area(T.A,T.B,T.C);}///面积
Point Circumcenter(const Triangle& T){return Circumcenter(T.A,T.B,T.C);}///外心
Point Incenter(const Triangle& T){return Incenter(T.A,T.B,T.C);}///内心
Point Perpencenter(const Triangle& T){return Perpencenter(T.A,T.B,T.C);}///垂心
Point Barycentry(const Triangle& T){return Barycentry(T.A,T.B,T.C);}///重心


int Type(const Triangle& T)///返回三角形的类型 -1 钝角 0直角 1 锐角
{return sign(Dot(T.A-T.B,T.A-T.C)*Dot(T.B-T.A,T.B-T.C)*Dot(T.C-T.A,T.C-T.B));}


bool InTriangle_ex(const Point& P,const Triangle& T)///返回点是否在三角形内部（含边）
{return !sign(Area(P,T.A,T.B)+Area(P,T.A,T.C)+Area(P,T.B,T.C)-Area(T));}


bool InTriangleLine(const Point& P,const Triangle& T)///返回点是否在三角形边上（含端点）
{return Online_ex(P,Segment(T.A,T.B))||Online_ex(P,Segment(T.A,T.C))||Online_ex(P,Segment(T.C,T.B));}


bool InTriangle(const Point& P,const Triangle& T)///返回点是否在三角形内部(不含边)
{return InTriangle_ex(P,T)&&!InTriangleLine(P,T);}






///多边形


struct Polygon
{
    int n;
    Point *P;
    Polygon(int nums=0,Point *point=NULL)
    {
        if(nums>2)
        {
            P=new Point[n=nums];
            for(int i=0;i<n;i++)
                P[i]=point[i];
        }
        else
        {
            n=0,P=NULL;
        }
    }


    Polygon(const Polygon& T)
    {
        n=T.n;
        P=new Point[n];
        for(int i=0;i<n;i++)
            P[i]=T[i];
    }


    Point& operator [](int key)
    {
        return P[key];
    }
    const Point& operator [](int key)const
    {
        return P[key];
    }
};


double Area(const Polygon& P)///多边形的面积
{
    double ret=0;
    for(int i=1;i<P.n-1;i++)
    {
        ret+=Area2(P[0],P[i],P[i+1]);
    }
    return fabs(ret);
}


Polygon ConvexHull_ex(Polygon P)///凸包扫描算法(含平角)
{
    sort(P.P,P.P+P.n);
    Polygon ret(P.n,P.P);
    int n=P.n,&m=P.n;
    m=0;
    for(int i=0;i<n;i++)
    {
        while(m>1&&sign(Cross(ret[m-1]-ret[m-2],P[i]-ret[m-2]))<0)m--;
        ret[m++]=P[i];
    }
    int k=m;
    for(int i=n-2;i>=0;i--)
    {
        while(m>k&&sign(Cross(ret[m-1]-ret[m-2],P[i]-ret[m-2]))<0)m--;
        ret[m++]=P[i];
    }
    if(m>1)m--;
    return ret;
}


int ConvexHull_ex(int n,Point *P,Point *ans)///凸包扫描算法（含平角）
{
    Polygon ret=ConvexHull_ex(Polygon(n,P));
    for(int i=0;i<ret.n;i++)
        ans[i]=ret[i];
    return ret.n;
}


Polygon ConvexHull(Polygon P)///凸包扫描算法(不含平角)
{
    sort(P.P,P.P+P.n);
    Polygon ret(P.n,P.P);
    int n=P.n,&m=P.n;
    m=0;
    for(int i=0;i<n;i++)
    {
        while(m>1&&sign(Cross(ret[m-1]-ret[m-2],P[i]-ret[m-2]))<=0)m--;
        ret[m++]=P[i];
    }
    int k=m;
    for(int i=n-2;i>=0;i--)
    {
        while(m>k&&sign(Cross(ret[m-1]-ret[m-2],P[i]-ret[m-2]))<=0)m--;
        ret[m++]=P[i];
    }
    if(m>1)m--;
    return ret;
}


int ConvexHull(int n,Point *P,Point *ans)///凸包扫描算法（不含平角）
{
    Polygon ret=ConvexHull(Polygon(n,P));
    for(int i=0;i<ret.n;i++)
        ans[i]=ret[i];
    return ret.n;
}


bool IsConvex(const Polygon& P)///判断是否为凸多边形 （凸多边形定义包含平角）
{
    if(P.n<=3)return 1;
    int t=sign(Area2(P[0],P[1],P[2]));
    for(int i=1;i<P.n-2;i++)
        if(t*sign(Area2(P[i],P[i+1],P[i+2]))<0)return 0;
    return 1;
}


bool IsConvex(int n,Point *P)///判断是否为凸多边形
{
    return IsConvex(Polygon(n,P));
}


int InPolygon(const Point& P,const Polygon& POLY)///返回点是（1）否（0）在多边形内部 负数为在PLOY[-i]和[(-i+1)%n]的线段上
{
    int ret=0,n=POLY.n;
    for(int i=0;i<n;i++)
    {
        if(Online_ex(P,Segment(POLY[i],POLY[(i+1)%n])))
            return -i;
        int t=sign(Cross(POLY[(i+1)%n]-POLY[i],P-POLY[i]));
        int d1=sign(POLY[i].y-P.y);
        int d2=sign(POLY[(i+1)%n].y-P.y);
        if(t>0&&d1<=0&&d2>0)ret++;
        if(t<0&&d2<=0&&d1>0)ret--;
    }
    return ret!=0;
}


///圆


struct Circle
{
    Point C;
    double r;
    Circle(Point C=Point(),double r=0.0):C(C),r(r){}
};


double Area(const Circle& C){return pi*sqr(C.r);}


int Position(const Point& P,const Circle& C)///返回点与圆的位置关系 -1 内部 0 圆上 1 外部
{
    return sign(Length(P,C.C)-C.r);
}


int Position(const Line& L,const Circle& C)///返回直线与圆的位置关系 -1 相离 0 相切 1 相交
{
    return -sign(Distance(L,C.C)-C.r);
}


Segment Intersection(const Line& L,const Circle& C)
///返回直线与圆的交点 相交时返回割线（两端点即为交点） 相切时返回切点（返回线段的两个端点都相同）
///相离时返回 NAS
{
    double d=sqr(C.r)-sqr(Distance(L,C.C));
    if(sign(d)>0)
    {
        d/=Dot(L.v,L.v);
        d=sqrt(d);
        double dx=L.v.x*d,dy=L.v.y*d;
        Point t=Projection(L,C.C),D(dx,dy);
        return Segment(t+D,t-D);
    }
    else if(sign(d)==0)
    {
        Point ret=Projection(L,C.C);
        return Segment(ret,ret);
    }
    else
        return NAS;
}


int Position(const Circle& C1,const Circle& C2)
/// 圆与圆的位置关系
///0 相交 1 外切 -1内切 2 相离 -2 内含
{
    double d=Distance(C1.C,C2.C);
    switch(sign(d-C1.r-C2.r))
    {
    case -1:return 2;
    case 0:return 1;
    }
    switch(sign(fabs(C1.r-C2.r)-d))
    {
    case 0:return -1;
    case 1:return -2;
    }
    return 0;
}


Segment Intersection(const Circle& C1,const Circle& C2)
/// 圆与圆的公共弦 不相交返回 NAS
{
    int key=Position(C1,C2);
    if(!key)
    {
        double d=Distance(C1.C,C2.C);
        double s=C1.r+C2.r+d;
        s/=2;
        double h2=4.0*s*(s-C1.r)*(s-C2.r)*(s-d)/d/d;
        double l=sqrt(sqr(C1.r)-h2);
        Point D=(C2.C-C1.C)*(l/d)+C1.C;
        double h=sqrt(h2);
        Line L=NormalLine(Line(C1.C,C2.C),D);
        return Intersection(L,C1);
    }
    else if(abs(key)==1)
    {
        Point ret=C1.r/(Distance(C1.C,C2.C))*(C2.C-C1.C)+C1.C;
        return Segment(ret,ret);
    }
    else return NAS;
}


/*************copyright by sinianluoye (JLU_LiChuang)***********/


u013613980
关注
  • 0
    点赞
  • 踩
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
图形学实验三《二维图形几何变换算法》代码+报告
12-23
（1）掌握二维平移、比例、旋转几何变换矩阵 （2）掌握相对于任一参考点的比例变换和旋转变换 （3）了解定时器的使用方法和边界碰撞检测方法 （4）了解静态切分视图框架设计方法 实验概括 　对于一个二维图形作平移...
【C++代码 简单好用的 计算几何学 模板】一些点 线 圆的基本操作
溢流眼泪的博客
04-11 424
之前学了一下下C++ 的OOP，然后跑过来跟着竞赛书把计算几何学的模板补好了！不得不说代码量是真的大！ 利用了类的封装，方法重载，操作符重载等。 这里的范数应该是范2（模长概念？） ，书上的范数是模长的平方 assert()那里，要保证有相交才能保证有交点。 ///存储：点 向量 线段 直线 圆 ///点操作：內积，外积，equals与== 相等，+操作 -操作 ///直线与直线的正交判断dot...
几何模板1 ( C++ )
3000bookworm
06-16 108
#include<iostream> #include<cstdio> #include<cmath> using namespace std; typedef Point Vector; const double pi=acos(-1.0); const double eps=1E-8; class Point { public: double x,y; Point() {} Point(double x, double y):x(x),y(y
几何模板
jkchen's Haven
07-20 404
#include&amp;amp;amp;amp;amp;amp;lt;iostream&amp;amp;amp;amp;amp;amp;gt; #include&amp;amp;amp;amp;amp;amp;lt;cstdio&amp;amp;amp;amp;amp;amp;gt; #include&amp;amp;amp;amp;amp;amp;lt;cmath&amp;amp;amp;amp;amp;amp;gt; #include&amp;a
算法系列之九：计算几何与图形学有关的几种常用算法（二）
weixin_33957648的博客
12-27 133
算法系列之九：计算几何与图形学有关的几种常用算法（二） 分类： 算法系列 2011-12-25 23:04 1149人阅读 评论(5) 收藏 举报 3.6 用矢量的叉积判断直线段是否有交 矢量叉积计算的另一个常用用途是直线段求交。求交算法是计算机图形学的核心算法，也是体现速度和稳定性的...
计算几何
weixin_30505485的博客
02-15 50
点积、叉积 \(u×v<0\) ： \(u\) 在 \(v\) 的逆时针方向 判断线段相交，求交点 首先判断两条线段是否相交 图中，设两线段交点为 \(E\) \(\because \frac{|CE|}{|DE|}=\frac{△ACD在CD边上的高}{△BCD在CD边上的高}\;\;=\frac{S_{△ACD}}{S_{△BCD}}=\frac{CD×CA}{CB×CD}=\frac...
ACM专用模板，包括数据结构，图论，字符串，数论，几何计算，KD树，ST表，二叉搜索树，二维ST表（任意矩形，正方形），二维树状数组模板，分块ACM模板.rar
07-26
ACM专用模板，包括数据结构，图论，字符串，数论，几何计算，KD树，ST表，二叉搜索树，二维ST表（任意矩形，正方形），二维树状数组模板，分块，权值线段树 区间最早出现问题，树状数组模板，线段树模板，线段树扫描...
kuangbin acm模板超级好用
09-05
3.2.2 二维 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 57 3.3 树链剖分 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 59 3.3.1 点权 . . ....
win8风格二维码解说PPT模板
03-20
这是一个win8风格二维码解说PPT，第... 二维码 ，又称二维条码，二维条形码最早发明于日本，它是用某种特定的几何图形按一定规律在平面（二维方向上）分布的黑白相间的图形记录数据符号信息的，通过图象输入设备或光电
浙江大学ACM模板 计算几何，图论，数据结构，经典题的模板
09-18
1.9三维几何 22 1.10 凸包 30 1.11 网格 32 1.12 圆 33 1.13 矢量运算求几何模板 35 1.14结构体表示几何图形 47 1.15四城部分几何模板 52 1.16 一些代码 54 1.16.1 最小圆覆盖_zju1450 54 1.16.2 直线旋转_两凸包的...
计算几何模板(仿照刘汝佳大白书风格)
PSU_VJD的专栏
11-08 1892
不完全版，有时间再补充
【几何】-几何基础大模板（更新中）
邂逅黛玛姑娘
05-29 645
四省赛的选拔赛出了一堆几何题
计算几何题目总结
hechenghai
11-08 574
一。基础题目  1.1 有固定算法的题目  A， 最近点对问题 最近点对问题的算法基于扫描线算法。 ZOJ     2107    Quoit Design    典型最近点对问题 POJ    3714    Raid    变种最近点对问题 B，最小包围圆 最小包围圆的算法是一种增量算法，期望是O(n)。 ZOJ    1450    Minimal Circle
node-v12.20.1-sunos-x64.tar.xz
最新发布
04-27
Node.js，简称Node，是一个开源且跨平台的JavaScript运行时环境，它允许在浏览器外运行JavaScript代码。Node.js于2009年由Ryan Dahl创立，旨在创建高性能的Web服务器和网络应用程序。它基于Google Chrome的V8 JavaScript引擎，可以在Windows、Linux、Unix、Mac OS X等操作系统上运行。 Node.js的特点之一是事件驱动和非阻塞I/O模型，这使得它非常适合处理大量并发连接，从而在构建实时应用程序如在线游戏、聊天应用以及实时通讯服务时表现卓越。此外，Node.js使用了模块化的架构，通过npm（Node package manager，Node包管理器）,社区成员可以共享和复用代码，极大地促进了Node.js生态系统的发展和扩张。 Node.js不仅用于服务器端开发。随着技术的发展，它也被用于构建工具链、开发桌面应用程序、物联网设备等。Node.js能够处理文件系统、操作数据库、处理网络请求等，因此，开发者可以用JavaScript编写全栈应用程序，这一点大大提高了开发效率和便捷性。 在实践中，许多大型企业和组织已经采用Node.js作为其Web应用程序的开发平台，如Netflix、PayPal和Walmart等。它们利用Node.js提高了应用性能，简化了开发流程，并且能更快地响应市场需求。
基于Springboot+Vue的乡政府管理系统-毕业源码案例设计.zip
04-27
网络技术和计算机技术发展至今，已经拥有了深厚的理论基础，并在现实中进行了充分运用，尤其是基于计算机运行的软件更是受到各界的关注。加上现在人们已经步入信息时代，所以对于信息的宣传和管理就很关键。系统化是必要的，设计网上系统不仅会节约人力和管理成本，还会安全保存庞大的数据量，对于信息的维护和检索也不需要花费很多时间，非常的便利。 网上系统是在MySQL中建立数据表保存信息，运用SpringBoot框架和Java语言编写。并按照软件设计开发流程进行设计实现。系统具备友好性且功能完善。 网上系统在让售信息规范化的同时，也能及时通过数据输入的有效性规则检测出错误数据，让数据的录入达到准确性的目的，进而提升数据的可靠性，让系统数据的错误率降至最低。 关键词：vue；MySQL；SpringBoot框架 【引流】 Java、Python、Node.js、Spring Boot、Django、Express、MySQL、PostgreSQL、MongoDB、React、Angular、Vue、Bootstrap、Material-UI、Redis、Docker、Kubernetes
一名合格的程序猿修炼手册.md
04-27
一名合格的程序猿修炼手册.
5MHz 函数发生器使用说明书
04-27
5MHz 函数发生器使用说明书
99- 矿山工业互联网平台解决方案.pptx
04-27
99- 矿山工业互联网平台解决方案.pptx
基于Python大学生社会实践申报系统的设计与实现带vue前后端分离毕业源码案例设计.zip
04-27
网络技术和计算机技术发展至今，已经拥有了深厚的理论基础，并在现实中进行了充分运用，尤其是基于计算机运行的软件更是受到各界的关注。加上现在人们已经步入信息时代，所以对于信息的宣传和管理就很关键。系统化是必要的，设计网上系统不仅会节约人力和管理成本，还会安全保存庞大的数据量，对于信息的维护和检索也不需要花费很多时间，非常的便利。 网上系统是在MySQL中建立数据表保存信息，运用SpringBoot框架和Java语言编写。并按照软件设计开发流程进行设计实现。系统具备友好性且功能完善。 网上系统在让售信息规范化的同时，也能及时通过数据输入的有效性规则检测出错误数据，让数据的录入达到准确性的目的，进而提升数据的可靠性，让系统数据的错误率降至最低。 关键词：vue；MySQL；SpringBoot框架 【引流】 Java、Python、Node.js、Spring Boot、Django、Express、MySQL、PostgreSQL、MongoDB、React、Angular、Vue、Bootstrap、Material-UI、Redis、Docker、Kubernetes
输入一个表示几何形状名称的字符串，再在一行内输入这种图形的数据，根据表示名称的字符串选择合适的公式计算几何形状的(表)面积和体积，若为二维图形，只计算面积，若为三维图形，计算其表面积与体积，结果严格保留2位小数。 模板程序给出了长方形和长方体的代码，参考这些代码，完成圆形、球、圆柱体、圆锥和正三棱柱这5种形状的计算程序。
05-30
好的，这是一个计算几何形状面积和体积的程序，支持圆形、球、圆柱体、圆锥和正三棱柱这5种形状的计算。以下是完整的代码： ```python import math # 计算长方形面积 def area_rectangle(width, height): return ...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值