poj1269——Intersecting Lines//直线交点

判断两直线的位置：或重合，或平行，或相交。

用叉积搞定。框架，求交点算法自己实现。

#include<iostream>
using namespace std;
const double epx=1e-10;
struct Point
{
    double x;
    double y;
};
double Multi (Point p1,Point p2,Point  p3)
{
	// p1p2 X p1p3
    return (p2.x -p1.x )*(p3.y -p1.y )-(p2.y -p1.y)*(p3.x -p1.x );
}
//判断两直线的关系
/*
在这里有三种关系：1共线 2平行 3相交
1 共线可通过叉积来判断
2 平行通过向量来判断
3 通了上面2种情况的其他情况
  求交点可通过叉积及解二元一次方程来求解
*/
int N;
Point p1,p2,p3,p4;
Point p0;//交点
double a1,b1,c1,a2,b2,c2;
int main()
{
    scanf("%d",&N);
    int i;
    printf("INTERSECTING LINES OUTPUT\n");
    for(i=0;i<N;++i)
    {
        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);
        if(Multi(p3,p4,p1)==0 && Multi(p3,p4,p2)==0)//共线  即判断三点是否共线  做两次叉积 
            printf("LINE\n");
        else
        {
			double b1,b2,D,D1,D2;//求交点
			b1=(p2.y -p1.y )*p1.x +(p1.x -p2.x )*p1.y ;
			b2=(p4.y -p3.y )*p3.x +(p3.x -p4.x )*p3.y ;
			D=(p2.x -p1.x )*(p4.y -p3.y )-(p4.x -p3.x )*(p2.y -p1.y );
            if(D==0 )//平行
                printf("NONE\n");
            else
            {
				Point pt;
				D1=b2*(p2.x -p1.x )-b1*(p4.x -p3.x );
		     	D2=b2*(p2.y -p1.y )-b1*(p4.y -p3.y );
		    	pt.x =D1/D;pt.y =D2/D;//求出交点
                printf("POINT %.2f %.2f\n",pt.x,pt.y);
            }
        }
    }
    printf("END OF OUTPUT\n");
    return 0;
}