计算几何模板.docVIP

#includeiostream#includemath.h#includevector#includealgorithm#includestdio.husing namespace std; const double eps = 1e-10;const double pi = acos( -1.0 );const double inf = 1e100; inline int sig( double x ){??? return (x eps) - (x - eps);} struct point{ //数据结构??? double x, y, dis, angle; }??? point(){}??? point(const double x_,const double y_):x(x_),y(y_){}??? //逻辑运算??? inline bool operator ( const point t)const{??????? return sig ( y - t.y ) 0 || sig ( y - t.y ) == 0 sig( x - t.x ) 0 ;??? }??? inline bool operator == ( const point t )const{??????? return sig ( x - t.x ) == 0 sig ( y - t.y ) == 0;??? }??? //算术运算??? inline point operator - ( const point b)const{??????? return point( x - b.x , y - b.y );??? }??? inline point operator + ( const point b )const{??????? return point( x + b.x, y + b.y );??? }??? inline point operator * ( const double t )const{??????? return point ( x * t , y * t );??? }??? inline point operator / ( const double t )const{??????? return point ( x / t , y / t );??? }??? //求长度??? inline double len(){??????? return sqrt( x * x + y * y );??? } ??? //长度的平方??? inline double len2(){??????? return x * x + y * y ;??? }??? inline double dist( point t ){??????? return sqrt( (x - t.x)*(x - t.x) + (y - t.y)*(y - t.y) );??? }??? //输入输出??? int read(){??????? return scanf(%lf%lf,x,y);??? }??? void out(){??????? printf((%.2lf %.2lf) ,x,y);??? }}; bool _cmp(const point a,const point b ){??? return sig( a.angle - b.angle ) 0 || sig( a.angle - b.angle ) == 0 sig( a.dis - b.dis ) 0;}// 点积inline double dot( point p , point a, point b ){//这种形式更常用??? return ( a.x - p.x ) * (b.x - p.x) + (a.y - p.y) * (b.y - p.y);} //叉积 ,遵循 左手定则（左旋为正）inline double cross( point a, point b){??? return a.y * b.x - a.x * b.y ;}inline double cross( point p , point a, point b){//这种形式更常用??? return (a.y - p.y) *(b.x - p.x) - (a.x - p.x) * (b.y - p.y);} //angle 0 表示逆时针旋转，angle 0 为顺时针//矢量pv以p为顶点逆时针旋转angle并放大sca