POJ1269判断2条直线的位置关系

题目大意：给两个点能够确定一条直线，题目给出两条直线（由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 ;
}