Intersecting Lines POJ - 1269

http://poj.org/problem?id=1269

叉积真的是计算几何利器

首先对于两条直线(p1 p2) (p3 p4) 若两直线共线 则p3-p1与p2-p1的叉积模为0且p4-p1与p2-p1的叉积模为0

否则 用斜率判一下平行

最后 如果相交求交点 设交点为p0 则p1-p0与p2-p0的叉积模为0 p3-p0与p4-p0的叉积模为0

联立后可得a1*x+b1*y+c1=0与a2*x+b2*y+c2=0这两个线性方程组 直接克拉默法则即可求解

 

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
using namespace std;
const double eps=1e-8;

double getval(double x1,double y1,double x2,double y2,double x3,double y3)
{
    return (x2-x1)*(y3-y1)-(x3-x1)*(y2-y1);
}

int main()
{
    double x1,y1,x2,y2,x3,y3,x4,y4;
    double a1,b1,c1,a2,b2,c2,x,y;
    int t;
    scanf("%d",&t);
    printf("INTERSECTING LINES OUTPUT\n");
    while(t--){
        scanf("%lf%lf%lf%lf%lf%lf%lf%lf",&x1,&y1,&x2,&y2,&x3,&y3,&x4,&y4);
        if(fabs(getval(x3,y3,x1,y1,x2,y2))<eps&&fabs(getval(x4,y4,x1,y1,x2,y2))<eps) printf("LINE\n");
        else if(fabs((y2-y1)*(x4-x3)-(y4-y3)*(x2-x1))<eps) printf("NONE\n");
        else{
            a1=y1-y2,b1=x2-x1,c1=x1*y2-x2*y1;
            a2=y3-y4,b2=x4-x3,c2=x3*y4-x4*y3;
            x=(-c1*b2+b1*c2)/(a1*b2-a2*b1);
            y=(-a1*c2+a2*c1)/(a1*b2-a2*b1);
            printf("POINT %.2f %.2f\n",x,y);
        }
    }
    printf("END OF OUTPUT\n");
    return 0;
}

/*
1
5 0 7 6 3 -6 4 -3
*/