POJ_1269_Intersecting Lines_求直线交点
Description
Input
Output
Sample Input
5
0 0 4 4 0 4 4 0
5 0 7 6 1 0 2 3
5 0 7 6 3 -6 4 -3
2 0 2 27 1 5 18 5
0 3 4 0 1 2 2 5
Sample Output
INTERSECTING LINES OUTPUT
POINT 2.00 2.00
NONE
LINE
POINT 2.00 5.00
POINT 1.07 2.20
END OF OUTPUT
我们都知道一双独特的点在平面上定义了一条线,一条线在一个平面相交的三种方法:
1)没有交集,因为它们是平行的,
2)相交于一条线,因为他们是在另一个(即他们是相同的线),
3)相交于一点。在这个问题中，你将使用你的代数知识来创建一个程序来决定两条线的交点。
你的程序将会反复地读入四个点，在xy平面上定义两条直线，并确定直线的交点和位置。这个问题所要求的所有数字都是合理的，比如在-1000和1000之间。
先用两条直线的向量的叉积判断是否平行，然后用一个端点向另外一条直线的两个端点连线求叉积判断是否重合。
都不是就直接求两直线的交点，用平行四边形的面积求。
代码：
#include <stdio.h>
#include <string.h>
#include <algorithm>
using namespace std;
typedef double f2;
#define eps 1e-6
f2 fabs(f2 x) {return x>0?x:-x;}
struct Point {
f2 x,y;
Point() {}
Point(f2 x_,f2 y_):
x(x_),y(y_) {}
Point operator + (const Point &p) const {return Point(x+p.x,y+p.y);}
Point operator - (const Point &p) const {return Point(x-p.x,y-p.y);}
Point operator * (f2 rate) const {return Point(x*rate,y*rate);}
void rd() {scanf("%lf%lf",&x,&y);}
};
typedef Point Vector;
f2 dot(const Point &p1,const Point &p2) {return p1.x*p2.x+p1.y*p2.y;}
f2 cross(const Point &p1,const Point &p2) {return p1.x*p2.y-p1.y*p2.x;}
struct Line {
Point p;
Vector v;
Line() {}
Line(const Point &p_,const Vector &v_):
p(p_),v(v_) {}
};
Point get_point(const Line &l1,const Line &l2) {
Vector u=l1.p-l2.p;
f2 t=cross(l2.v,u)/cross(l1.v,l2.v);
return l1.p+l1.v*t;
}
void solve() {
Point a1,a2,b1,b2;
a1.rd();a2.rd();b1.rd();b2.rd();
Vector A=a1-a2,B=b1-b2;
if(fabs(cross(A,B))<eps) {
if(fabs(cross(a1-b1,a2-b1))<eps&&fabs(cross(a1-b2,a2-b2))<eps) puts("LINE");
else puts("NONE");
}else {
Point ans=get_point(Line(a2,A),Line(b2,B));
printf("POINT %.2lf %.2lf\n",ans.x,ans.y);
}
}
int main() {
puts("INTERSECTING LINES OUTPUT");
int n;
scanf("%d",&n);
while(n--) {
solve();
}
puts("END OF OUTPUT");
}
POJ_1269_Intersecting Lines_求直线交点
原文：https://www.cnblogs.com/suika/p/9017750.html