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

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 ;
}