POJ 1269 Intersecting Lines（几何）

程序员宏春技术 2022年11月15日

0 收藏 363 点赞 3,668 浏览 1443 个字

题意：给你两条线段的起点和终点，一共四个点，让你求交点坐标，如果这四个点是共线的，输出“LINE”，如果是平行的就输出“NONE”。

思路：照着ZN留下的模板果然好用，直接套上模板了事儿，不过在判断是否共线的时候，其实还有另一种方法，直接将平行和共线一起判断了，我是判断三个点三个点的判断是否是共线。

 //
 #include <stdio.h>
 #include <string.h>
 #include <iostream> using namespace std ; #define eps 1e-8
 #define zero(x) (((x)>0?(x):-(x))<eps)
 typedef struct
 {
     double x,y ;
 } point ;
 typedef struct
 {
     point a,b ;
 } line ; //计算cross product (P1-P0)x(P2-P0)
 double xmult(double x1,double y1,double x2,double y2,double x0,double y0)
 {
     return (x1-x0)*(y2-y0)-(x2-x0)*(y1-y0);
 } //判三点共线
 int dots_inline(double x1,double y1,double x2,double y2,double x3,double y3)
 {
     return zero(xmult(x1,y1,x2,y2,x3,y3));
 } //线段交点请另外判线段相交(同时还是要判断是否平行!)
 point intersection(point u1,point u2,point v1,point v2)
 {
     point ret=u1;
     double t=((u1.x-v1.x)*(v1.y-v2.y)-(u1.y-v1.y)*(v1.x-v2.x))/((u1.x-u2.x)*(v1.y-v2.y)-(u1.y-u2.y)*(v1.x-v2.x));
     ret.x+=(u2.x-u1.x)*t;
     ret.y+=(u2.y-u1.y)*t;
     return ret;
 } int parallel(point u1,point u2,point v1,point v2)
 {
     return zero((u1.x-u2.x)*(v1.y-v2.y)-(v1.x-v2.x)*(u1.y-u2.y));
 } int main()
 {
     int n ;
     scanf("%d",&n) ;
     point p1,q1,p2,q2 ;
     printf("INTERSECTING LINES OUTPUT\n") ;
     for(int i =  ; i < n ; i++)
     {
         scanf("%lf %lf %lf %lf",&p1.x,&p1.y,&p2.x,&p2.y) ;
         scanf("%lf %lf %lf %lf",&q1.x,&q1.y,&q2.x,&q2.y) ;
         if(dots_inline(p1.x,p1.y,p2.x,p2.y,q1.x,q1.y)&&dots_inline(p1.x,p1.y,p2.x,p2.y,q2.x,q2.y))
             printf("LINE\n") ;
         else if(parallel(p1,p2,q1,q2))
             printf("NONE\n") ;
         else
         {
             point ret = intersection(p1,p2,q1,q2) ;
             printf("POINT %.2lf %.2lf\n",ret.x,ret.y) ;
         }
     }
     printf("END OF OUTPUT\n") ;
     return  ;
 }