几何算法模板.docx

几何算法 /*****************************************?COMPUTATIONAL?GEOMETRY?ROUTINES*?WRITTEN?BY?:?LIU?Yu?(C)?2003****************************************///????叉乘//????两个点的距离? //????点到直线距离//????返回直线?Ax?+?By?+?C?=0??的系数?//????线段//????圆//????两个圆的公共面积//????矩形//????根据下标返回多边形的边?//????两个矩形的公共面积?//????多边形??,逆时针或顺时针给出x,y?//????多边形顶点?//????多边形的边?//????多边形的周长?//????判断点是否在线段上?//????判断两条线断是否相交，端点重合算相交?//????判断两条线断是否平行//????判断两条直线断是否相交//????直线相交的交点?//????判断是否简单多边形?//????求多边形面积?//????判断是否在多边形上?//????判断是否在多边形内部//????点阵的凸包，返回一个多边形? //? 最近点对的距离#include?cmath#include?cstdio#include?memory#include?algorithm#include?iostreamusing?namespace?std;typedef?double?TYPE;? //把double 定义为TYPE#define?Abs(x)?(((x)0)?(x):(-(x)))? //用Abs()这个宏定义一个绝对值函数 #define?Sgn(x)?(((x)0)?(-1):(1)) //取相反数#define?Max(a,b)?(((a)(b))?(a):(b)) //取大值#define?Min(a,b)?(((a)(b))?(a):(b)) //取小值#define?Epsilon?1e-10 #define?Infinity?1e+10#define?Pi?3.14159265358979323846 //定义几个常量TYPE?Deg2Rad(TYPE?deg){return?(deg?*?Pi?/?180.0);}? //把角度制转化为弧度制TYPE?Rad2Deg(TYPE?rad){return?(rad?*?180.0?/?Pi);}? //把弧度制转化为角度制TYPE?Sin(TYPE?deg){return?sin(Deg2Rad(deg));} //对弧度制求正弦TYPE?Cos(TYPE?deg){return?cos(Deg2Rad(deg));} //对弧度制求余弦TYPE?ArcSin(TYPE?val){return?Rad2Deg(asin(val));}? //求反正弦TYPE?ArcCos(TYPE?val){?return?Rad2Deg(acos(val));} //求反余弦TYPE?Sqrt(TYPE?val){?return?sqrt(val);}struct?POINT{????TYPE?x;????TYPE?y;????TYPE?z;????POINT()?:?x(0),?y(0),?z(0)?{};????POINT(TYPE?_x_,?TYPE?_y_,?TYPE?_z_?=?0)?:?x(_x_),?y(_y_),?z(_z_)?{};};//?cross?product?of?(o-a)?and?(o-b)//?叉乘?TYPE?Cross(const?POINT??a,?const?POINT??b,?const?POINT??o){return?(a.x?-?o.x)?*?(b.y?-?o.y)?-?(b.x?-?o.x)?*?(a.y?-?o.y);}? //叉积模板//?planar?points?distance//??两个点的距离?TYPE?Distance(const?POINT??a,?const?POINT??b){return?Sqrt((a.x?-?b.x)?*?(a.x?-?b.x)?+?(a.y?-?b.y)?*?(a.y?-?b.y)?+?(a.z?-?b.z)?*?(a.z?-?b.z));}struct?LINE{????POINT?a;????POINT?b;????LINE()?{};????LINE(POINT?_a_,?POINT?_b_)?:?a(_a_