POJ - 1269 Intersecting Lines 两直线的位置关系

#include <cmath>

#include <cstdio>

const double T = 1e-8;

struct Node{

    double x, y;

}p1, p2, p3, p4;

double ChaJi(Node a, Node b, Node c){

    return (c.x-a.x)*(b.y-a.y)-(c.y-a.y)*(b.x-a.x);

} 

void Solve(){

    if(ChaJi(p1, p2, p3) == 0 && ChaJi(p1, p2, p4) == 0)        //共线； 

        puts("LINE");

    else if((p1.x-p2.x)*(p3.y-p4.y) == (p1.y-p2.y)*(p3.x-p4.x)) // 平行

        puts("NONE");

    else{                                                       //相交

        double a1 = p1.y-p2.y;

        double b1 = p2.x-p1.x;

        double c1 = p1.x*p2.y-p2.x*p1.y;

        double a2 = p3.y-p4.y;

        double b2 = p4.x-p3.x;

        double c2 =p3.x*p4.y-p4.x*p3.y;

        double x = (c1*b2-c2*b1)/(a2*b1-a1*b2);

        double y = (a2*c1-a1*c2)/(a1*b2-a2*b1);

        printf("POINT %.2f %.2f\n", x, y);

    } 

} 

int main(){

    int n;

    while(scanf("%d", &n) != EOF){

        printf("INTERSECTING LINES OUTPUT\n");

        while(n--){

            scanf("%lf%lf%lf%lf%lf%lf%lf%lf", &p1.x, &p1.y, &p2.x, &p2.y, &p3.x, &p3.y, &p4.x, &p4.y);

            Solve(); 

        }

        printf("END OF OUTPUT\n");

    }

    return 0;

}