Intersecting Lines（POJ-1269）

---

1 题意

  判断两条直线是否平行，重合或相交，如果相交则输出交点坐标。
  链接：link。

2 思路

  命题1：向量$\vec{a},\vec{b}$平行的充要条件为$\vec{a}\times \vec{b}=0$。
  证明：
  先证必要性，已知$\vec{a}\times \vec{b}=\left|\vec{a}\right|\left|\vec{b}\right|\sin <\vec{a},\vec{b}>=0$。
  1. 若$\left|\vec{a}\right|=0$或$\left|\vec{b}\right|=0$，则说明$\vec{a}$或$\vec{b}$为零向量，由零向量与任何向量平行，可知$\vec{a}//\vec{b}$。
  2. 若$\sin <\vec{a},\vec{b}>=0$，则说明$<\vec{a},\vec{b}>=\pi$或$<\vec{a},\vec{b}>=2\pi$，显然两个向量平行。
  再证必要性，若向量$\vec{a},\vec{b}$平行，则说明$<\vec{a},\vec{b}>=\pi$或$<\vec{a},\vec{b}>=2\pi$，那么$\sin <\vec{a},\vec{b}>=0$，即$\vec{a}\times \vec{b}=0$，得证。

  命题2：两条直线平行或重合的充要条件是其方向向量的向量积为0。
  证明：略。

  命题3：假设两条直线$l_1,l_2$的方向向量分别为$\overrightarrow{P_1{P}_2},\overrightarrow{P_3{P}_4}$，且两条直线平行或重合，即$\overrightarrow{P_1{P}_2}\times \overrightarrow{P_3{P}_4}=0$。如果$\overrightarrow{P_1{P}_4}\times \overrightarrow{P_3{P}_2}=0$，则两条直线重合，否则两条直线平行。
  证明：略。

  命题3：假设两条直线$l_1,l_2$相交，且其方向向量分别为$\overrightarrow{P_1{P}_2},\overrightarrow{P_3{P}_4}$，其中$P_1=(x_1,y_1),P_2=(x_2,y_2),P_3=(x_3,y_3),P_4=(x_4,y_4)$，则两条直线的交点为$(\frac{(x_1{y}_2-x_2{y}_1)(x_4-x_3)-(x_3{y}_4-x_4{y}_3)(x_2-x_1)}{(y_2-y_1)(x_4-x_3)-(y_4-y_3)(x_2-x_1)},\frac{(y_2-y_1)(x_3{y}_4-x_4{y}_3)-(y_4-y_3)(x_1{y}_2-x_2{y}_1)}{(x_2-x_1)(y_4-y_3)-(x_4-x_3)(y_2-y_1)})$。

2.1 时间复杂度分析

  判断直线是否平行或相交的时间复杂度为$\mathcal{O}(1)$，计算交点坐标的时间复杂度为$\mathcal{O}(1)$，所以总的时间复杂度为$\mathcal{O}(1)$。

2.2 实现

#include<iostream>
#include<cstdio>
using namespace std;
int main(){
    int T;scanf("%d",&T);
    puts("INTERSECTING LINES OUTPUT");
    while(T--){
        double x1,y1,x2,y2,x3,y3,x4,y4;
        scanf("%lf %lf %lf %lf %lf %lf %lf %lf",&x1,&y1,&x2,&y2,&x3,&y3,&x4,&y4);
        if((x2-x1)*(y4-y3)==(x4-x3)*(y2-y1)){
            if((x3-x1)*(y2-y1)==(x2-x1)*(y3-y1)) puts("LINE");
            else puts("NONE");
        }else{
            double a1=y2-y1,b1=x2-x1,c1=x1*y2-x2*y1;
            double a2=y4-y3,b2=x4-x3,c2=x3*y4-x4*y3;
            printf("POINT %.2lf %.2lf\n",(c1*b2-c2*b1)/(a1*b2-a2*b1),(a1*c2-a2*c1)/(b1*a2-b2*a1));
        }
    }
    cout<<"END OF OUTPUT"<<endl;
    return 0;
}