POJ 1269 Interesting Lines

题目大意：

给出四个点坐标，确定两条直线，判断两条直线，是否平行，是否共线，如果有交点，求出焦点。

解题思路：

从特殊到一般。

如何确定是否共线：

假设此时有A B C D四个点，如果A B C共线，且A B D共线，则说明两条直线共线。

如何确定是否平行：

两条直线的叉积是否相等，如果相等则平行。

一定有交点怎么算交点：

假设交点为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<cstdio>
#include<cstring>
#include<string>
#include<cmath>
#include<iostream>
#define LL long long
#define db double
#define maxn 10000000
#define EPS 1e-15
#define max(a,b) ((a>b)?(a):(b))
#define min(a,b) ((a<b)?(a):(b))
using namespace std;

struct point{
	db x,y;
}a,b,c,d;
db cross(point a,point b,point c){
	return((b.x-a.x)*(c.y-a.y)-(c.x-a.x)*(b.y-a.y));
}
void solve(point a,point b,point c,point d){
	if (fabs(cross(a,b,c))<=EPS && fabs(cross(a,b,d))<=EPS) // 叉积=0 共线
		printf("LINE\n");
	else if((b.x-a.x)*(d.y-c.y)==(d.x-c.x)*(b.y-a.y)) //平行
		printf("NONE\n");
	else   //相交  
    {  
        double a1=a.y-b.y;  
        double b1=b.x-a.x;  
        double c1=a.x*b.y-b.x*a.y;  
  
        double a2=c.y-d.y;  
        double b2=d.x-c.x;  
        double c2=c.x*d.y-d.x*c.y;  
  
        double x=(b1*c2-b2*c1)/(a1*b2-a2*b1);  
        double y=(a2*c1-a1*c2)/(a1*b2-a2*b1);  
  
        printf("POINT %.2lf %.2lf\n",x,y);  
    }  
}
int main(){
	int n;
	while(~scanf("%d",&n)){
	printf("INTERSECTING LINES OUTPUT\n");
	for (int i=1;i<=n;i++){
		scanf("%lf%lf%lf%lf%lf%lf%lf%lf",&a.x,&a.y,&b.x,&b.y,&c.x,&c.y,&d.x,&d.y);
		solve(a,b,c,d);
		}
	printf("END OF OUTPUT\n");
	}
	return 0;
}