文章目录
直线的表示
直线有多种表示方法。编程时,可以灵活使用这些方法:
struct Line
{
Point p1,p2; //线上的两个点
Line(){}
Line(Point p1,Point p2):p1(p1),p2(p2){}
//根据一个点和倾斜角angle——>确定直线
//其中:0<=angle<pi (就是0到180度)
Line(Point p,double angle)
{
p1 = p;
if(sgn(angle – pi/2) == 0)
{
p2 = (p1 + Point(0,1));
}
else
{
p2 = (p1 + Point(1,tan(angle)));
}
}
//ax+by+c=0
Line(double a,double b,double c)
{
if(sgn(a) == 0)
{
//y=-c/b
//p1是把x=0带入y=-c/b,得-c/b
//p2是把x=1带入y=-c/b,得-c/b
p1 = Point(0,-c/b);
p2 = Point(1,-c/b);
}
else if(sgn(b) == 0)
{
//x=-c/a
p1 = Point(-c/a,0);
p2 = Point(-c/a,1);
}
else
{
//p2-p1刚好是直线:y=-a/b*x-c/b的斜率
//p1是把x=0带入y=-a/b*x-c/b,得-c/b
//p2是把x=1带入y=-a/b*x-c/b,得(-c-a)/b
p1 = Point(0,-c/b);
p2 = Point(1,(-c-a)/b);
}
}
};
线段的表示
可以用两个点来表示线段,起点是p1,终点是p2。
直接用直线的数据结构来定义线段:typedef Line Segment
;
位置关系
1. 点和直线的位置关系:点在直线左侧、在右侧、在直线上。
2. 点和线段的位置关系:点 p 是否在线段 v 上
点和直线的其它关系
1. 点到直线的距离
2. 点关于直线的对称点
补充:【常规做法】
简单做法:
先求点p在直线上的投影q,再用中点公式求对称点p’
点在直线上的投影
Point Point_line_proj(Point p, Line v)
{
double k=Dot(v.p2-v.p1,p-v.p1)/Len2(v.p2-v.p1);
return v.p1+(v.p2-v.p1)*k;
}
Point Point_line_symmetry(Point p, Line v)
{
Point q = Point_line_proj(p,v);
return Point(2*q.x-p.x,2*q.y-p.y);
}
代码解释
3. 点到线段的距离
我们设现在有一个点p以及线段AB,那么点p到线段AB的距离只要从3个距离中取最小值即可。
4. 两条直线的位置关系
两条直线的位置关系可能有:
重合
平行
相交
int Line_relation(Line v1, Line v2)
{
//判断直线是否相交 为0:不相交
if(sgn(Cross(v1.p2-v1.p1,v2.p2-v2.p1)) == 0)
{
if(Point_line_relation(v1.p1,v2)==0) return 1; //1 重合
else return 0; //0 平行
}
return 2; //2 相交
}
5. 判断两个线段是否相交
关系两线段是否相交,我们可以利用叉积有正负的特点。
如果一条线段的两端在另一条线段的两侧,那么两个端点与另一线段的产生的两个叉积正负相反,也就是说两个叉积相乘为负。如果两条线段互相满足这一点,那么就是相交的。
//Line1:ab, Line2:cd
bool Cross_segment(Point a,Point b,Point c,Point d)
{
double c1=Cross(b-a,c-a),c2=Cross(b-a,d-a);
double d1=Cross(d-c,a-c),d2=Cross(d-c,b-c);
return sgn(c1)*sgn(c2)<=0 && sgn(d1)*sgn(d2)<=0; //1相交;0不相交
}
6. 求两条线段的交点
这个简单:先判断两条线段是否相交。若相交,问题转化成两条直线求交点;若不相交,则无交点