poj1269求两条直线的关系（相交（求交点）---叉乘），平行，共线）

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cmath>
#include <cstring>
#include <string.h>
#include <queue>          //poj1269求两条直线的关系（相交（求交点）---叉乘），平行，共线）
using namespace std;
struct line
{
double x1, y1, x2, y2, x3, y3, x4, y4;
}p[15];
int main()
{
int n, t;
double k1, k2, a, b, c, d, x, y, a1, a2, b1, b2, c1, c2;
while(scanf("%d", &n)!=EOF)
{
for(t=0; t<n; ++t)
{
scanf("%lf%lf%lf%lf%lf%lf%lf%lf", &p[t].x1, &p[t].y1, &p[t].x2, &p[t].y2, &p[t].x3, &p[t].y3, &p[t].x4, &p[t].y4);
}
printf("INTERSECTING LINES OUTPUT\n");
for(t=0; t<n; ++t)
{
k1=(p[t].y2-p[t].y1)/(p[t].x2-p[t].x1);
k2=(p[t].y3-p[t].y4)/(p[t].x3-p[t].x4);
a=p[t].x2-p[t].x1;
b=p[t].y2-p[t].y1;
c=p[t].x3-p[t].x4;
d=p[t].y3-p[t].y4;
if(a*d-b*c==0)
{
c=p[t].x3-p[t].x1;
d=p[t].y3-p[t].y1;
if(c*b-a*d==0)
{
printf("LINE\n");
}
else 
{
printf("NONE\n");
}
}
else 
{
a1=p[t].y1-p[t].y2;
b1=p[t].x2-p[t].x1;
c1=p[t].x1*p[t].y2-p[t].x2*p[t].y1;
a2=p[t].y3-p[t].y4;
b2=p[t].x4-p[t].x3;
c2=p[t].x3*p[t].y4-p[t].x4*p[t].y3;
x=(b1*c2-b2*c1)/(a1*b2-a2*b1);
y=(a2*c1-a1*c2)/(a1*b2-a2*b1);
printf("POINT %.2lf %.2lf\n", x, y);
}
}
printf("END OF OUTPUT\n");
}
return 0;
}