POJ 1269 计算几何 求直线交点

题意：

此题给出N组直线，每组2条直线。

如果共线则输出LINE，相交则输入点坐标，否则输出NONE（平行）。

 

题解：

点积判断平行和相交，有向面积（分点公式）求交点（不知道写的对不对，至少能过这个题，貌似数据很弱的样子）。

 

#include <iostream>
#include <cstring>
#include <cstdio>
#include <cstdlib>
#include <algorithm>
#include <cmath>

#define EPS 1e-7

using namespace std;

struct PO
{
    double x,y;
}li[3][2];

inline double cross(const PO &o,const PO &a,const PO &b)
{
    return (a.x-o.x)*(b.y-o.y)-(b.x-o.x)*(a.y-o.y);
}

inline void read()
{
    for(int i=1;i<=2;i++)
        scanf("%lf%lf%lf%lf",&li[i][0].x,&li[i][0].y,&li[i][1].x,&li[i][1].y);
}

inline bool judge_parallel(PO s[],PO t[])
{
    return fabs((s[1].y-s[0].y)*(t[1].x-t[0].x)-(s[1].x-s[0].x)*(t[1].y-t[0].y))<EPS;
}

inline bool judge_same(PO s[],PO t[])
{
    return fabs((s[1].y-s[0].y)*(t[1].x-s[0].x)-(s[1].x-s[0].x)*(t[1].y-s[0].y))<EPS;
}

inline PO get_intersection(PO s[],PO t[])
{
    PO ans;
    double k2=1.0;
    double k1=cross(s[0],t[1],t[0])/cross(s[1],t[0],t[1]);
    ans.x=(k2*s[0].x+k1*s[1].x)/(k1+k2);
    ans.y=(k2*s[0].y+k1*s[1].y)/(k1+k2);
    return ans;
}

inline void go()
{
    if(judge_parallel(li[1],li[2]))
    {
        if(judge_same(li[1],li[2])) puts("LINE");
        else puts("NONE");
        return;
    }
    PO ans=get_intersection(li[1],li[2]);
    printf("POINT %.2lf %.2lf\n",ans.x,ans.y);
}

int main()
{
    int cas; scanf("%d",&cas);
    puts("INTERSECTING LINES OUTPUT");
    while(cas--) read(),go();
    puts("END OF OUTPUT");
    return 0;
}