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;