计算几何相关（更新中）

1.向量

结构体定义点和直线

struct point{
	double x,y;
};
struct line{
	point a,b;
};

两点相减求向量，重载 - 运算符

struct point{
	int x,y;
	point(){};                               //因为有带参的构造函数，下面无法直接定义，所以加一个无参的构造函数 
	point(int newx,int newy)                 //需要重载-，所以设立一个带参的构造函数 
    	{  
        	x=newx;y=newy;  
    	}  
	point operator -(const point &b)        //两点相减得到向量，故重载一下- 
	{
		return point(x-b.x,y-b.y);
	}
}

向量点积

a•b=|a|*|b|*cosα              P•Q=x1*x2+y1*y2

点积代码

double dot(point a,point b)
{
	return a.x*b.x+a.y*b.y;
}

向量叉积

|PQ|=x1*y2-x2*y1

代码实现

double across(point a,point b)
{
	return a.x*b.y-a.y*b.x;
}

double across(line x,line y)
{
	point temp1=x.a-x.b;
	point temp2=y.a-y.b;
	
	return temp1.x*temp2.y-temp1.y*temp2.x;
}

叉乘可以解决很多线段或直线的关系，若两向量叉乘的积大于0，说明前者在后者的顺时针方向，小于0，说明前者在后者的逆时针方向，等于0属于共线。

2.点定位

判断点是否在线段上

设线段端点为P1，P2，点为Q。

满足（Q-P1)*(P2-P1)=0且Q在以P1，P2为对角顶点的矩阵内

第一条保证Q在直线P1P2上，第二条保证Q不在其延长线或反向延长线上。

代码实现

bool onSegment(point p1,point p2,point q)
{
	if((q.x-p1.x)*(p2.y-p1.y)==(p2.x-p1.x)*(q.y-p1.y)&&min(p1.x,p2.x)<=q.x&&q.x<=max(p1.x,p2.x)&&min(p1.y,p2.y)<=q.y&&q.y<=max(p1.y,p2.y))
	return 1;
	else
	return 0;
}

判断点在三角形的内外

点Q与三角形的其他两点构成三个三角形，若S1+S2+S3==Sabc，则表明Q在三角形内

面积计算公式 S=1/2*|across（a，b）|     向量叉积的一半

代码实现

bool intriangle(triangle p,point q)
{
	double Sabc=1/2*fabs(across(p.a,p.b));
	double Sqab=1/2*fabs(across(q-p.a,q-p.b));
	double Sqac=1/2*fabs(across(q-p.a,q-p.c));
	double Sqbc=1/2*fabs(across(q-p.b,q-p.c));
	
	if(Sqab+Sqac+Sqbc==Sabc)
	return true;
	else
	return false;
}

判断点在多边形的内外

第一种方法   扫描法

懒得解释

代码略

第二种方法   叉乘判别法

类似与凸包算法吧，每个点与第一个点不断做叉乘，累积叉乘，结果若为正则继续，直到遍历结束。

代码可以参考凸包

第三种方法   角度和的判断法

若点在图形内，连接点与每一个顶点，所有角度之和为360°时，代表点在图形内。

代码实现

double angle = 0;
realPointList::iterator iter1 = points.begin();
for (realPointList::iterator iter2 = (iter1 + 1); iter2 < points.end(); ++iter1, ++iter2)
 {
   double x1 = (*iter1).x - p.x;  
   double y1 = (*iter1).y - p.y;  
   double x2 = (*iter2).x - p.x;
   double y2 = (*iter2).y - p.y;  
   angle += angle2D(x1, y1, x2, y2);
 }

if (fabs(angle - span::PI2) < 0.01) return true;
else return false;