计算几何袁伟模板.docVIP

以下函数全是我亲自验证过的 如果有问题 请告诉我 让我也不要犯错 QQ；497330525 谢谢； /*?需要包含的头文件?*/?#include?cmath??/*?常用的常量定义?*/?const?double?INF?=?1E200????const?double?EP?=?1E-10;?const?int?MAXV?=?300;?const?double?PI?=?3?/*?基本几何结构?*/?class?POINT?{?public:double?x;?double?y;?POINT(double?a=0,?double?b=0)?{?x=a;?y=b;}?//constructor?};? //线段:class?LINESEG?{?public:POINT?s;?POINT?e;?LINESEG(POINT?a,?POINT?b)?{?s=a;?e=b;}?LINESEG()?{?}?};? //直线：class?LINE???????????//?直线的解析方程?a*x+b*y+c=0??为统一表示，约定?a?=?0?{????double?a;????double?b;????double?c;?LINE(double?d1=1,?double?d2=-1,?double?d3=0)?{a=d1;?b=d2;?c=d3;}?}; 所有函数： /********************?*????????????????????*?*???点的基本运算?????*?*????????????????????*?********************/? double?dist(POINT?p1,POINT?p2)???????//?返回两点之间欧氏距离?{?return(?sqrt(?(p1.x-p2.x)*(p1.x-p2.x)+(p1.y-p2.y)*(p1.y-p2.y)?)?);?}?bool?equal_point(POINT?p1,POINT?p2)???? //?判断两个点是否重合??{?return?(?(abs(p1.x-p2.x)EP)(fabs(p1.y-p2.y)EP)?);?}?/******************************************************************************?/* 返回值大于0 one-O在two-O顺时针方向endl; 返回值小于0 one-O在two-O逆时针方向endl; 返回值等于0 one-O与two-O在一条直线上endl; */ double multiply(POINT one,POINT two,POINT O)//查乘积 { return((one.x-O.x)*(two.y-O.y)-(two.x-O.x)*(one.y-O.y)); } /*******************************************************************************?r=dotmultiply(p1,p2,op),得到矢量(p1-op)和(p2-op)的点积，如果两个矢量都非零矢量?r0:两矢量夹角为锐角；r=0：两矢量夹角为直角；r0:两矢量夹角为钝角?*******************************************************************************/?double?dotmultiply(POINT?p1,POINT?p2,POINT?p0)?{?return?((p1.x-p0.x)*(p2.x-p0.x)+(p1.y-p0.y)*(p2.y-p0.y));?}?/*?判断点p是否在线段l上，条件：(p在线段l所在的直线上)?(点p在以线段l为对角线的矩形内)?*/?bool?online(LINESEG?l,POINT?p)?{?return((multiply(l.e,p,l.s)==0)?(?(?(p.x-l.s.x)*(p.x-l.e.x)=0?)(?(p.y-l.s.y)*(p.y-l.e.y)=0?)?)?);?}?//?返回点p以点o为圆心逆时针旋转alpha(单位：弧度)后所在的位置?POINT?rotate(POINT?o,double?alpha,POINT?p)?{?POINT?tp;?p.x-=o.x;?p.y-=o.y;?tp.x=p.x*cos(alpha)-p.y*sin(alpha)+o.x;?tp.y=p.