判断两条直线 共线或 平行 或相交 求交点 叉乘的应用 poj 1269 Intersecting Lines

#include <cstdlib>
#include <cstring>
#include <algorithm>
#include <cstdio>
#include <cmath>
#include <iostream>
#include <vector>
#include<map>
using namespace std;
typedef long long ll;
const int N =30;
const double PI = 3.1415927;
const double EPS=1e-10;
struct Point{
    double x,y;
    Point(){}
    Point(double x,double y):x(x),y(y){} // 构造函数，方便代码编写
    Point(const Point & p):x(p.x),y(p.y){}
    Point operator +(Point p){
        return Point(x+p.x,y+p.y);
    }
    Point operator-(Point p){
        return Point(x-p.x,y-p.y);
    }
    Point operator*(double d){
        return Point(x*d,y*d);
    }
    double operator*(Point p){  // 内积 点乘
        return x*p.x+ y*p.y;
    }
    double operator^(Point p){//  外积 叉乘
        return x*p.y-y*p.x;
    }
    friend ostream& operator<<(ostream& os,const Point& p ){
        os<<p.x<<" "<<p.y<<endl;
        return os;
    }
    friend istream& operator>>(istream& is, Point&p) {// Point 不能是常量，是变量
        is>>p.x>>p.y;
        return is;
    }
    double dist(Point p){
        return sqrt( (x-p.x)*(x-p.x)+ (y-p.y)*(y-p.y) );
    }
};
Point intersecting(Point p1,Point p2,Point q1,Point q2){ // 计算直线p1p2与q1q2的交点
    double d1=( (q2-q1)^(q1-p1) );
    double d2=( (q2-q1)^(p2-p1) );
    return p1+ (p2-p1)*(d1/d2) ;
}
Point p1,p2,q1,q2;
void solve() // 判断直线p1p2与直线q1q2的关系
{
    double d1=(p2-p1)^(q1-p1);
    double d2=(p2-p1)^(q2-p1);
    if(fabs(d1 ) < EPS && fabs(d2)  < EPS)  // 如果p1p2q1 共线 且 p1p2q2 共线 则直线共线
    {
        puts("LINE");
        return ;
    }
    else
    {
        if(fabs( (p2-p1)^(q2-q1) ) < EPS  )  // 否则，如果 p1p2 与 q1q2 共线， 则直线平行
        {
            puts("NONE");
            return ;
        }
        else                           //  否则，则相交
        {
            Point temp=intersecting(p1,p2,q1,q2);
            printf("POINT %.2lf %.2lf\n",temp.x,temp.y);
        }
    }
}

int main() {
    int n;
    cin>>n;
    printf("INTERSECTING LINES OUTPUT\n");
    for(int i=0;i<n;i++)
    {
        cin>>p1>>p2>>q1>>q2;
        solve();
    }
    printf("END OF OUTPUT\n");
    return 0;
}