码迷,mamicode.com
首页 > 其他好文 > 详细

poj 1269 Intersecting Lines——叉积求直线交点坐标

时间:2018-09-19 23:28:16      阅读:254      评论:0      收藏:0      [点我收藏+]

标签:题目   output   参数   end   lse   col   fc7   href   style   

题目:http://poj.org/problem?id=1269

相关知识:

叉积
求面积:https://www.cnblogs.com/xiexinxinlove/p/3708147.html
什么是叉积:https://blog.csdn.net/sunbobosun56801/article/details/78980467
        其二维:https://blog.csdn.net/qq_38182397/article/details/80508303
计算交点:
    方法1:面积的比例:https://blog.csdn.net/dgq8211/article/details/7952825
    方法2:叉积:https://blog.csdn.net/hktkfly6/article/details/69218777
            (其原因:https://www.jianshu.com/p/3468c9967fc7
例题:POJ1269:https://www.xuebuyuan.com/1979631.html

此题的链接使用的就是那个面积比例的方法。但自己不是很懂,所以用了那个另一种方法。

利用面积算交点坐标应该还可以:https://www.cnblogs.com/Narh/p/9663099.html

  设直线有两个参数:p->起点坐标;t->终点坐标减起点坐标。设有线段 a 和 b 。设交点为 r 。

  double d=( (a.p-b.p)×a.t ) / ( b.t×a.t ) ,r = a.p + d * a.t ;

  那个 d 就是面积比,即交点与一条线段的相对该线段起点的坐标值之比。

    画一画,分子就是三角形面积的两倍,分母是那样的平行四边形的面积,但连一下发现它就是两线交叉的那个四边形的面积的两倍。面积比即高的比即那两段线段之比。

关于这种方法的 a、b、c 怎么算,可以这样想:y=k*x+b 又 a*x+b*y+c=0 ,所以 a : b : c = -k : 1 : b;

  已知一条直线的两个点的坐标,联立可得 k=( y1-y2 ) / ( x1-x2 ) ; b=( x1*y2 - y1*x2 ) / ( x1-x2 ) ;

  所以 a : b : c = ( y2-y1) : ( x1-x2 ) : ( x1*y2 - y1*x2 ) 。可以说它们就等于这个,因为另一侧是0,可以随便乘。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define db double
using namespace std;
int T,tx,ty;
struct Pt{
    db x,y;
    Pt(db x=0,db y=0):x(x),y(y) {}
}ta,tb;
struct Line{
    Pt p,t;
}a,b;
int operator* (const Pt &a,const Pt &b)
{return a.x*b.y-a.y*b.x;}
Pt operator- (const Pt &a,const Pt &b)
{return Pt(a.x-b.x,a.y-b.y);}
int check(Line a,Line b)
{
    int d1=(b.p-a.p)*(a.t-a.p);
    int d2=(b.t-a.p)*(a.t-a.p);
    if(!d1&&!d2) return -1;//line
    if(d1==d2) return 1;//none
    if(d1*d2<=0) return 0;//point ==0 for 端点相交
}
Pt find(Line a,Line b)
{
    db x1=a.p.x,y1=a.p.y,x2=a.t.x,y2=a.t.y;
    db a0=y2-y1,b0=x1-x2,c0=y1*x2-x1*y2;
    x1=b.p.x,y1=b.p.y,x2=b.t.x,y2=b.t.y;
    db a1=y2-y1,b1=x1-x2,c1=y1*x2-x1*y2;
    int tmp=a0*b1-a1*b0;
    return Pt((b0*c1-b1*c0)/tmp,(-a0*c1+a1*c0)/tmp);
}
int main()
{
    scanf("%d",&T);
    puts("INTERSECTING LINES OUTPUT");
    while(T--)
    {
        scanf("%d%d",&tx,&ty); ta=Pt(tx,ty);
        scanf("%d%d",&tx,&ty); tb=Pt(tx,ty);
        a.p=ta; a.t=tb;
        scanf("%d%d",&tx,&ty); ta=Pt(tx,ty);
        scanf("%d%d",&tx,&ty); tb=Pt(tx,ty);
        b.p=ta; b.t=tb;
        int d=check(a,b);
        if(d==-1) puts("LINE");
        else if(d==1) puts("NONE");
        else
        {
            printf("POINT ");
            Pt c=find(a,b);
            printf("%.2lf %.2lf\n",c.x,c.y);
        }
    }
    puts("END OF OUTPUT");
    return 0;
}

 

poj 1269 Intersecting Lines——叉积求直线交点坐标

标签:题目   output   参数   end   lse   col   fc7   href   style   

原文地址:https://www.cnblogs.com/Narh/p/9678068.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!