专栏首页饶文津的专栏【POJ 1269】判断两直线相交

【POJ 1269】判断两直线相交

利用叉积解方程

#include <cstdio>
#define MAX 1<<31
#define dd double
int xmult(dd x1,dd y1,dd x2,dd y2,dd x,dd y){
    return (x1-x)*(y2-y)-(x2-x)*(y1-y);
}
int main(){
    int n;
    dd x1,y1,x2,y2,x3,y3,x4,y4;
    scanf("%d",&n);
    puts("INTERSECTING LINES OUTPUT");
    while(n--){
        scanf("%lf%lf%lf%lf%lf%lf%lf%lf",&x1,&y1,&x2,&y2,&x3,&y3,&x4,&y4);
        dd a1=y1-y2;
        dd b1=x2-x1;
        dd c1=xmult(x1,y1,x2,y2,0,0);
        dd a2=y3-y4;
        dd b2=x4-x3;
        dd c2=xmult(x3,y3,x4,y4,0,0);
        if(a1*b2==a2*b1){
            if(xmult(x1,y1,x2,y2,x3,y3)==0)
                puts("LINE");
            else puts("NONE");
        }
        else {
            dd cx,cy;
            cx=(b1*c2-c1*b2)/(a1*b2-a2*b1);
            cy=(a1*c2-c1*a2)/(b1*a2-b2*a1);
            printf("POINT %.2f %.2f\n",cx,cy);
        }
    }
    puts("END OF OUTPUT");
}

利用点斜式解方程

#include <cstdio>
#define MAX 1<<31
#define dd double
struct P{
    dd x,y;
    void input(){
        scanf("%lf%lf",&x,&y);
    }
    void output(){
        printf("POINT %.2f %.2f\n",x,y);
    }
};
struct L{
    P s,e;
    void input(){
        s.input(),e.input();
    }
    dd k(){
        if(s.x==e.x)return MAX;
        return (s.y-e.y)/(s.x-e.x);
    }
}l1,l2;
int xmult(P a,P b,P o){
    return (a.x-o.x)*(b.y-o.y)-(b.x-o.x)*(a.y-o.y);
}
void getCross(L a,L b){
    P c;
    dd ka=a.k(),kb=b.k();
    if(ka==MAX){//a是竖直的
        c.x=a.s.x;
        c.y=(c.x-b.s.x)*kb+b.s.y;
    }
    else{
        if(kb==MAX)
            c.x=b.s.x;
        else
            c.x=(a.s.y-b.s.y-ka*a.s.x+kb*b.s.x)/(kb-ka);
        c.y=(c.x-a.s.x)*ka+a.s.y;
    }
    c.output();
}
int main(){
    int n;
    scanf("%d",&n);
    puts("INTERSECTING LINES OUTPUT");
    while(n--){
        l1.input(),l2.input();
        dd ka=l1.k(),kb=l2.k();
        if(ka==kb){
            if(xmult(l1.s,l1.e,l2.s)==0)
                puts("LINE");
            else puts("NONE");
        }
        else 
            getCross(l1,l2);
    }
    puts("END OF OUTPUT");
}

本文参与腾讯云自媒体分享计划,欢迎正在阅读的你也加入,一起分享。

我来说两句

0 条评论
登录 后参与评论

相关文章

  • 【POJ 3525】Most Distant Point from the Sea(直线平移、半平面交)

    二分这个半径,将所有直线向多边形中心平移r距离,如果半平面交不存在那么r大了,否则r小了。

    饶文津
  • 【HDU 4445】Crazy Tank(暴力)

    高中物理斜抛运动,简单分析一下角度固定下来则可以计算每个cannonball的降落坐标lnd。

    饶文津
  • 【HDU 5105】Math Problem

      b≠0,f为二次函数,最大值点在区间端点或者x0=c/(2*b),当L≤x0≤R时,ans=max{f(L),f(R),f(x0)}。

    饶文津
  • 如何用<dl>标签做表格而不用table标签

      我们都知道很多的内容编辑器(TinyMCE编辑器、fck)都有插入表格功能,快速方便,但是这些表格用到的<table>标签,可以查看html源代码就能发现,...

    ytkah
  • 用re和xpath进行爬虫信息提取

    一般而言,3种提取数据的方法中,re速度最快,但设计正则表达式规则相对复杂;xpath速度其次,其设计规则一定程度上类似有些类似于从sql中查询数据,难度居中;...

    luanhz
  • VS2017中运行MySQL的存储过程

    zls365
  • 分布式跟踪工具-pinpoint安装部署

    博主无意间接触到了这款工具!而且使用的时候觉得特别爽!比其他的检测工具强的不是一丁半点(个人认为)! 所以在这里分享给大家使用!

    互扯程序
  • 百度熊掌号折腾手记(二)chrome控制台报警告性错误。

    自从上次部署完熊掌号之后,就没在管过了,发的几篇文章有收录的有没收录的。前两天小伙伴们在博客112群里说页面加了熊掌号后会在chrome中报警告性错误。网上也都...

    世纪访客
  • 谷歌小姐姐搞出魔法画板:你随便画,补不齐算AI输

    当然,你也可以不让它画猫,改成画狗。只要你设定了一个绘画的目标,之后随便画一笔,AI就能脑补出余下的画面。

    量子位
  • 淘宝镜像错误导致vue项目一系列异常

    新建一个项目的时候,项目提示的这个 indexing paused due to project generation。这个问题当时没有咋当回事,使用Inval...

    用户1956326

扫码关注云+社区

领取腾讯云代金券