视频1 视频21 视频41 视频61 视频文章1 视频文章21 视频文章41 视频文章61 推荐1 推荐3 推荐5 推荐7 推荐9 推荐11 推荐13 推荐15 推荐17 推荐19 推荐21 推荐23 推荐25 推荐27 推荐29 推荐31 推荐33 推荐35 推荐37 推荐39 推荐41 推荐43 推荐45 推荐47 推荐49 关键词1 关键词101 关键词201 关键词301 关键词401 关键词501 关键词601 关键词701 关键词801 关键词901 关键词1001 关键词1101 关键词1201 关键词1301 关键词1401 关键词1501 关键词1601 关键词1701 关键词1801 关键词1901 视频扩展1 视频扩展6 视频扩展11 视频扩展16 文章1 文章201 文章401 文章601 文章801 文章1001 资讯1 资讯501 资讯1001 资讯1501 标签1 标签501 标签1001 关键词1 关键词501 关键词1001 关键词1501 专题2001
POJ1269判断2条直线的位置关系
2020-11-09 07:19:25 责编:小采
文档

题目大意:给两个点能够确定一条直线,题目给出两条直线(由4个点确定),要求判断出这两条直线的关系:平行,同线,相交。如果相交还要求出交点坐标。 解题思路: 先判断两条直线p1p2, q1q2是否共线, 如果不是,再判断 直线 是否平行, 如果还不是, 则两

题目大意:给两个点能够确定一条直线,题目给出两条直线(由4个点确定),要求判断出这两条直线的关系:平行,同线,相交。如果相交还要求出交点坐标。


解题思路:


先判断两条直线p1p2, q1q2是否共线, 如果不是,再判断 直线 是否平行, 如果还不是, 则两直线相交。


判断共线: p1p2q1 共线 且 p1p2q2 共线 ,共线用叉乘为 0 来判断,


判断 平行: p1p2 与 q1q2 共线


求交点:


直线p1p2上的点 可表示为 p1+t(p2-p1) , 而交点 又在 直线q1q2上, 所以有 (q2-q1)X (p1 + t(p2-p1 ) - q1 ) =0


解得 交点 t = p1 + ( ((q2-q1) X (q1 - p1)) /( (q2-q1) X(p2-p1) ) *(p2-p1) )

-----------------------------------------------------------------------
注意: double 型数据为0 不能直接==0

------------------------------------------------------------------------

叉乘不满足交换律


const double eps = 1e-8 ;

double add(double x , double y){
 if(fabs(x+y) < eps*(fabs(x) + fabs(y))) return 0 ;
 return x + y ;
}

struct Point{
 double x , y ;
 Point(){}
 Point(double _x , double _y):x(_x),y(_y){}
 Point operator + (Point o){
 return Point(add(x , o.x) , add(y , o.y)) ;
 }
 Point operator - (Point o){
 return Point(add(x , -o.x) , add(y , -o.y)) ;
 }
 Point operator * (double o){
 return Point(x*o , y*o) ;
 }
 double operator ^(Point o){
 return add(x*o.y , -y*o.x) ;
 }
 double dist(Point o){
 return sqrt((x-o.x)*(x-o.x) + (y-o.y)*(y-o.y)) ;
 }
 void read(){
 scanf("%lf%lf" ,&x , &y) ;
 }
};

//判断2条直线的位置关系
int twoline(Point p1 , Point p2 , Point q1 , Point q2 , Point &interp){
 if(fabs((p2-p1) ^ (q1-p1)) < eps &&
 fabs((p2-p1) ^ (q2-p1)) < eps ) return 1 ; //共线
 if(fabs((p2-p1) ^ (q2-q1)) < eps) return 2 ; //平行
 double d1 = (q2 - q1) ^ (q1 - p1) ;
 double d2 = (q2 - q1) ^ (p2 - p1) ;
 double t = d1 / d2 ;
 interp = p1 + (p2 - p1) * t ; // 交点 interp
 return 3 ; // 相交
}

int main(){
 int t , k ;
 Point p1 , p2 , q1 , q2 , interp ;
 puts("INTERSECTING LINES OUTPUT") ;
 cin>>t ;
 while(t--){
 p1.read() , p2.read() ;
 q1.read() , q2.read() ;
 k = twoline(p1 , p2 , q1 , q2 , interp) ;
 if(k == 1) puts("LINE") ;
 else if(k == 2) puts("NONE") ;
 else printf("POINT %.2lf %.2lf\n" , interp.x , interp.y) ;
 }
 puts("END OF OUTPUT") ;
 return 0 ;
}

下载本文
显示全文
专题