POJ-1269判断线与线的关系

内容来自http://blog.csdn.net/dreamvyps/article/details/6162690

代码是自己写的，WA两次，第一次一小部分的x写成了y，第二次就是编译器的原因，需要用C++提交，G++对于2.lf会出问题。

一、问题描述

http://acm.pku.edu.cn/JudgeOnline/problem?id=1269

题目大意：给两个点能够确定一条直线，题目给出两条直线（由4个点确定），要求判断出这两条直线的关系：平行，同线，相交。如果相交还要求出交点坐标。

二、解题思路

先判断两条直线是不是同线，不是的话再判断是否平行，再不是的话就只能是相交的，求出交点。

如何判断是否同线？由叉积的原理知道如果p1，p2，p3共线的话那么(p2-p1)X(p3-p1)=0。因此如果p1，p2，p3共线，p1，p2，p4共线，那么两条直线共线。direction（）求叉积，叉积为0说明共线。

如何判断是否平行？由向量可以判断出两直线是否平行。如果两直线平行，那么向量p1p2、p3p4也是平等的。即((p1.x-p2.x)*(p3.y-p4.y)-(p1.y-p2.y)*(p3.x-p4.x))==0说明向量平等。

如何求出交点？这里也用到叉积的原理。假设交点为p0(x0,y0)。则有：

(p1-p0)X(p2-p0)=0

(p3-p0)X(p2-p0)=0

展开后即是

(y1-y2)x0+(x2-x1)y0+x1y2-x2y1=0

(y3-y4)x0+(x4-x3)y0+x3y4-x4y3=0

将x0,y0作为变量求解二元一次方程组。

假设有二元一次方程组

a1x+b1y+c1=0;

a2x+b2y+c2=0

那么

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

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

因为此处两直线不会平行，所以分母不会为0。

 

下面贴上自己的代码：

#include<stdio.h>
#include<math.h>
#include<string.h>
struct Point
{
    double x,y;
    Point(double x=0,double y=0):x(x),y(y){}//构造函数
}p[5];
//就算OA和OB的叉积
double Cross(Point A,Point B)
{
    return A.x*B.y-A.y*B.x;
}
//点-点=向量
Point operator|(Point A,Point B)
{
    return Point(A.x-B.x,A.y-B.y);
}
int main()
{
    int n;
    scanf("%d",&n);
    printf("INTERSECTING LINES OUTPUT\n");
    while(n--)
    {scanf("%lf%lf%lf%lf%lf%lf%lf%lf",&p[1].x,&p[1].y,&p[2].x,&p[2].y,&p[3].x,&p[3].y,&p[4].x,&p[4].y);
   // printf("%lf %lf %lf %lf %lf %lf %lf %lf\n",p[1].x,p[1].y,p[2].x,p[2].y,p[3].x,p[3].y,p[4].x,p[4].y);
     if(Cross(p[2]|p[1],p[3]|p[1])==0&&Cross(p[2]|p[1],p[4]|p[1])==0) printf("LINE\n");
else if( Cross(p[1]|p[2],p[3]|p[4]) ==0 ) printf("NONE\n");
else { double a1=p[1].y-p[2].y,a2=p[3].y-p[4].y;
       double b1=p[2].x-p[1].x,b2=p[4].x-p[3].x;
       double c1=Cross(p[1],p[2]),c2=Cross(p[3],p[4]);
       double d=a1*b2-a2*b1;
       double dx=(b1*c2-b2*c1)/d,dy=(a2*c1-a1*c2)/d;
       printf("POINT %.2lf %.2lf\n",dx,dy);
     }

    }
    printf("END OF OUTPUT");
}