POJ 1410 Intersection 判断线段交和点在矩形内 【计算几何】

最新推荐文章于 2019-07-14 21:16:11 发布
最新推荐文章于 2019-07-14 21:16:11 发布
阅读量701 收藏 2
点赞数
分类专栏： 计算几何
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_40507857/article/details/83473513
版权

ACM博客_kuangbin

POJ 1410 Intersection（判断线段交和点在矩形内）

Intersection

Time Limit: 1000MS Memory Limit: 10000K
Total Submissions: 9996 Accepted: 2632

Description

You are to write a program that has to decide whether a given line segment intersects a given rectangle. 

An example: 
line: start point: (4,9) 
end point: (11,2) 
rectangle: left-top: (1,5) 
right-bottom: (7,1) 

 
Figure 1: Line segment does not intersect rectangle 

The line is said to intersect the rectangle if the line and the rectangle have at least one point in common. The rectangle consists of four straight lines and the area in between. Although all input values are integer numbers, valid intersection points do not have to lay on the integer grid. 

Input

The input consists of n test cases. The first line of the input file contains the number n. Each following line contains one test case of the format: 
xstart ystart xend yend xleft ytop xright ybottom 

where (xstart, ystart) is the start and (xend, yend) the end point of the line and (xleft, ytop) the top left and (xright, ybottom) the bottom right corner of the rectangle. The eight numbers are separated by a blank. The terms top left and bottom right do not imply any ordering of coordinates.

Output

For each test case in the input file, the output file should contain a line consisting either of the letter "T" if the line segment intersects the rectangle or the letter "F" if the line segment does not intersect the rectangle.

Sample Input

1
4 9 11 2 1 5 7 1

Sample Output

F

Source

Southwestern European Regional Contest 1995

题意：

给了一个线段和矩形。如果线段和矩形的边相交，或者线段在矩形内。输出T，否则输出F

思路：

先判断矩形4条边是否与线段相交，然后判断线段两个端点是否在矩形内

AC代码1：

转角法判断点在多边形内：实现简单，效率较低

//判断线段相交和点在矩形内
#include <iostream>
#include <stdio.h>
#include <string.h>
#include <algorithm>
#include <queue>
#include <map>
#include <vector>
#include <set>
#include <string>
#include <math.h>

using namespace std;
const double eps = 1e-8;
int sgn(double x)
{
    if(fabs(x) < eps)return 0;
    if(x < 0)return -1;
    else return 1;
}
struct Point
{
    double x,y;
    Point(){}
    Point(double _x,double _y)
    {
        x = _x;y = _y;
    }
    Point operator -(const Point &b)const
    {
        return Point(x - b.x,y - b.y);
    }
    //叉积
    double operator ^(const Point &b)const
    {
        return x*b.y - y*b.x;
    }
    //点积
    double operator *(const Point &b)const
    {
        return x*b.x + y*b.y;
    }
    //绕原点旋转角度B（弧度值），后x,y的变化
    void transXY(double B)
    {
        double tx = x,ty = y;
        x = tx*cos(B) - ty*sin(B);
        y = tx*sin(B) + ty*cos(B);
    }
};
struct Line
{
    Point s,e;
    double k;
    Line(){}
    Line(Point _s,Point _e)
    {
        s = _s;e = _e;
        k = atan2(e.y - s.y,e.x - s.x);
    }
    //两条直线求交点，
    //第一个值为0表示直线重合，为1表示平行，为2是相交
    //只有第一个值为2时，交点才有意义
    pair<int,Point> operator &(const Line &b)const
    {
        Point res = s;
        if(sgn((s-e)^(b.s-b.e)) == 0)
        {
            if(sgn((s-b.e)^(b.s-b.e)) == 0)
                return make_pair(0,res);//重合
            else return make_pair(1,res);//平行
        }
        double t = ((s-b.s)^(b.s-b.e))/((s-e)^(b.s-b.e));
        res.x += (e.x-s.x)*t;
        res.y += (e.y-s.y)*t;
        return make_pair(2,res);
    }
};
//两点间距离
double dist(Point a,Point b)
{
    return sqrt((a-b)*(a-b));
}
//判断线段相交
bool inter(Line l1,Line l2)
{
    return
        max(l1.s.x,l1.e.x) >= min(l2.s.x,l2.e.x) &&
        max(l2.s.x,l2.e.x) >= min(l1.s.x,l1.e.x) &&
        max(l1.s.y,l1.e.y) >= min(l2.s.y,l2.e.y) &&
        max(l2.s.y,l2.e.y) >= min(l1.s.y,l1.e.y) &&
        sgn((l2.s-l1.s)^(l1.e-l1.s))*sgn((l2.e-l1.s)^(l1.e-l1.s)) <= 0 &&
        sgn((l1.s-l2.s)^(l2.e-l1.s))*sgn((l1.e-l2.s)^(l2.e-l2.s)) <= 0;
}
//判断点在线段上
bool OnSeg(Point p,Line L)
{
    return
    sgn((L.s-p)^(L.e-p))==0&&
    sgn((p.x-L.s.x)*(p.x-L.e.x))<=0&&
    sgn((p.y-L.s.y)*(p.y-L.e.y))<=0;
}

//判断点在凸多边形内，点形成一个凸包，而且按逆时针排序
//如果是顺时针，就把<0改成>0
//点的编号：0~n-1
//返回值：-1点在凸多边形外，0点在凸多边形边上，1点在凸多边形内
int inConvexPoly(Point a,Point p[],int n)
{
    for(int i=0;i<n;i++)
    {
        if(sgn((p[i]-a)^(p[(i+1)%n]-a))<0)
           return -1;
        else if(OnSeg(a,Line(p[i],p[(i+1)%n])))
            return 0;
    }
    return 1;
}

//判断点在任意多边形内
//射线法，poly[]的顶点数要大于等于3，点的编号：0~n-1
//返回值
//-1：点在多边形外
//0：点在多边形边界上
//1：点在多边形内
int inPoly(Point p,Point poly[],int n)
{
    int cnt=0;
    Line ray,side;
    ray.s=p;
    ray.s.y=p.y;
    ray.e.x=-100000000000.0;//-INF,注意取值，防止越界

    for(int i=0;i<n;i++)
    {
        side.s=poly[i];
        side.e=poly[(i+1)%n];
        if(OnSeg(p,side))   return 0;
        //如果平行轴，则不考虑
        if(sgn(side.s.y-side.e.y)==0)
            continue;
        if(OnSeg(side.s,ray)){
            if(sgn(side.s.y-side.e.y)>0)
                cnt++;
        }
        else if(OnSeg(side.e,ray)){
            if(sgn(side.e.y-side.s.y)>0)
                cnt++;
        }
        else if(inter(ray,side))
            cnt++;
    }
    if(cnt%2==1)
        return 1;
    else
        return -1;
}
int main()
{
    int T;
    scanf("%d",&T);
    while(T--){
        Line line;
        scanf("%lf%lf%lf%lf",&line.s.x,&line.s.y,&line.e.x,&line.e.y);
            double x1,y1,x2,y2;
            scanf("%lf%lf%lf%lf",&x1,&y1,&x2,&y2);
            if(x1>x2)   swap(x1,x2);
            if(y1>y2)   swap(y1,y2);
            Point p[5];//0~3
            p[0]=Point(x1,y1);//left_bottom
            p[1]=Point(x2,y1);//right_bottom
            p[2]=Point(x2,y2);//right_top
            p[3]=Point(x1,y2);//left_top

            Line bottom=Line(p[1],p[0]);
            if(inter(line,bottom)){
                printf("T\n");
                continue;
            }

            Line right=Line(p[2],p[1]);
            if(inter(line,right)){
                printf("T\n");
                continue;
            }

            Line top=Line(p[3],p[2]);
            if(inter(line,top)){
                printf("T\n");
                continue;
            }

            Line left=Line(p[3],p[0]);
            if(inter(line,left)){
                printf("T\n");
                continue;
            }

            if(inConvexPoly(line.s,p,4)>=0||inConvexPoly(line.e,p,4)>=0){
                printf("T\n");
                continue;
            }

            printf("F\n");
            continue;
    }
    return 0;
}

AC代码2：

射线法判断点是否在多边形内，貌似要快一点：0 ms

//判断线段相交和点在矩形内
#include <iostream>
#include <stdio.h>
#include <string.h>
#include <algorithm>
#include <queue>
#include <map>
#include <vector>
#include <set>
#include <string>
#include <math.h>

using namespace std;
const double eps = 1e-8;
int sgn(double x)
{
    if(fabs(x) < eps)return 0;
    if(x < 0)return -1;
    else return 1;
}
struct Point
{
    double x,y;
    Point(){}
    Point(double _x,double _y)
    {
        x = _x;y = _y;
    }
    Point operator -(const Point &b)const
    {
        return Point(x - b.x,y - b.y);
    }
    //叉积
    double operator ^(const Point &b)const
    {
        return x*b.y - y*b.x;
    }
    //点积
    double operator *(const Point &b)const
    {
        return x*b.x + y*b.y;
    }
    //绕原点旋转角度B（弧度值），后x,y的变化
    void transXY(double B)
    {
        double tx = x,ty = y;
        x = tx*cos(B) - ty*sin(B);
        y = tx*sin(B) + ty*cos(B);
    }
};
struct Line
{
    Point s,e;
    double k;
    Line(){}
    Line(Point _s,Point _e)
    {
        s = _s;e = _e;
        k = atan2(e.y - s.y,e.x - s.x);
    }
    //两条直线求交点，
    //第一个值为0表示直线重合，为1表示平行，为2是相交
    //只有第一个值为2时，交点才有意义
    pair<int,Point> operator &(const Line &b)const
    {
        Point res = s;
        if(sgn((s-e)^(b.s-b.e)) == 0)
        {
            if(sgn((s-b.e)^(b.s-b.e)) == 0)
                return make_pair(0,res);//重合
            else return make_pair(1,res);//平行
        }
        double t = ((s-b.s)^(b.s-b.e))/((s-e)^(b.s-b.e));
        res.x += (e.x-s.x)*t;
        res.y += (e.y-s.y)*t;
        return make_pair(2,res);
    }
};
//两点间距离
double dist(Point a,Point b)
{
    return sqrt((a-b)*(a-b));
}
//判断线段相交
bool inter(Line l1,Line l2)
{
    return
        max(l1.s.x,l1.e.x) >= min(l2.s.x,l2.e.x) &&
        max(l2.s.x,l2.e.x) >= min(l1.s.x,l1.e.x) &&
        max(l1.s.y,l1.e.y) >= min(l2.s.y,l2.e.y) &&
        max(l2.s.y,l2.e.y) >= min(l1.s.y,l1.e.y) &&
        sgn((l2.s-l1.s)^(l1.e-l1.s))*sgn((l2.e-l1.s)^(l1.e-l1.s)) <= 0 &&
        sgn((l1.s-l2.s)^(l2.e-l1.s))*sgn((l1.e-l2.s)^(l2.e-l2.s)) <= 0;
}
//判断点在线段上
bool OnSeg(Point p,Line L)
{
    return
    sgn((L.s-p)^(L.e-p))==0&&
    sgn((p.x-L.s.x)*(p.x-L.e.x))<=0&&
    sgn((p.y-L.s.y)*(p.y-L.e.y))<=0;
}

//判断点在凸多边形内，点形成一个凸包，而且按逆时针排序
//如果是顺时针，就把<0改成>0
//点的编号：0~n-1
//返回值：-1点在凸多边形外，0点在凸多边形边上，1点在凸多边形内
int inConvexPoly(Point a,Point p[],int n)
{
    for(int i=0;i<n;i++)
    {
        if(sgn((p[i]-a)^(p[(i+1)%n]-a))<0)
           return -1;
        else if(OnSeg(a,Line(p[i],p[(i+1)%n])))
            return 0;
    }
    return 1;
}

//判断点在任意多边形内
//射线法，poly[]的顶点数要大于等于3，点的编号：0~n-1
//返回值
//-1：点在多边形外
//0：点在多边形边界上
//1：点在多边形内
int inPoly(Point p,Point poly[],int n)
{
    int cnt=0;
    Line ray,side;
    ray.s=p;
    ray.s.y=p.y;
    ray.e.x=-100000000000.0;//-INF,注意取值，防止越界

    for(int i=0;i<n;i++)
    {
        side.s=poly[i];
        side.e=poly[(i+1)%n];
        if(OnSeg(p,side))   return 0;
        //如果平行轴，则不考虑
        if(sgn(side.s.y-side.e.y)==0)
            continue;
        if(OnSeg(side.s,ray)){
            if(sgn(side.s.y-side.e.y)>0)
                cnt++;
        }
        else if(OnSeg(side.e,ray)){
            if(sgn(side.e.y-side.s.y)>0)
                cnt++;
        }
        else if(inter(ray,side))
            cnt++;
    }
    if(cnt%2==1)
        return 1;
    else
        return -1;
}
int main()
{
    int T;
    scanf("%d",&T);
    while(T--){
        Line line;
        scanf("%lf%lf%lf%lf",&line.s.x,&line.s.y,&line.e.x,&line.e.y);
            double x1,y1,x2,y2;
            scanf("%lf%lf%lf%lf",&x1,&y1,&x2,&y2);
            if(x1>x2)   swap(x1,x2);
            if(y1>y2)   swap(y1,y2);
            Point p[5];//0~3
            p[0]=Point(x1,y1);//left_bottom
            p[1]=Point(x2,y1);//right_bottom
            p[2]=Point(x2,y2);//right_top
            p[3]=Point(x1,y2);//left_top

            Line bottom=Line(p[1],p[0]);
            if(inter(line,bottom)){
                printf("T\n");
                continue;
            }

            Line right=Line(p[2],p[1]);
            if(inter(line,right)){
                printf("T\n");
                continue;
            }

            Line top=Line(p[3],p[2]);
            if(inter(line,top)){
                printf("T\n");
                continue;
            }

            Line left=Line(p[3],p[0]);
            if(inter(line,left)){
                printf("T\n");
                continue;
            }

            if(inPoly(line.s,p,4)>=0||inPoly(line.e,p,4)>=0){
                printf("T\n");
                continue;
            }

            printf("F\n");
            continue;
    }
    return 0;
}

 

紫芝
关注
专栏目录
计算机图形学中点分割算法法 (线段与矩形关系的判断 向量)
sherlockreal的博客
11-07 6695
线段与矩形的关系我理解的是三种：1、线段完全在矩形内 2、线段完全在矩形外 3、线段有一部分在矩形内 点与矩形的关系比线段与矩形的关系更容易判断，所以以下都是通过端点的情况去判断线段与矩形的关系的。具体也是分为三种： 1、线段两端点都在矩形内 2、线段一个端点在矩形内一个在矩形外 3、线段两个端点都在矩形外1、线段两端点都在矩形内：　　很容易得出这种情况下，线段是完全在矩形内部的。具体关
POJ 1410 Intersection (判断线段是否与矩形相交)
Christine的专栏
07-17 1001
题目大意：给定一条线段的起始和终止点的坐标，和一个矩形的左上和右下的坐标，判断线段是否与矩形相交。注意：矩形是实心       的，而且矩形的两点坐标不一定是按照左上、右下的顺序。 解题思路：首先需要确定矩形左上和右下的坐标，判断线段是否在矩形内部。接着判断线段是否和矩形的某条边相交即可。 代码如下：#include #include #include #include #include
POJ1410
xxx0624
07-25 202
注意：这个矩形是实心！ 其余没什么说的。 判断线段与矩形的位置关系 #include&lt;stdio.h&gt; #include&lt;string.h&gt; #include&lt;stdlib.h&gt; #include&lt;algorithm&gt; using namespace std; const double eps = 1e-8; struct Point{ ...
POJ 1410 Intersection（判断线段交和点在矩形内）
HJ
01-14 540
Description You are to write a program that has to decide whether a given line segment intersects a given rectangle. An example: line: start point: (4,9) end point: (11,2) rectangle: left-top: (
判断线段、折线、多边形是否在矩形中
weixin_34406061的博客
11-04 603
2019独角兽企业重金招聘Python工程师标准>>> ...
POJ 1410 Intersection（判断线段交 点在多边形内）
AI蜗牛之家
05-18 995
Intersection Time Limit:1000MS     Memory Limit:10000KB     64bit IO Format:%I64d & %I64u Submit Status Practice POJ 1410 Description You are to write a program that has to decide wh
POJ 1410 Intersection(判断线段是否在矩形面里)
binwin20
09-06 741
题意：给你线段和矩形，问是否相交，或矩形包含线段； 注意：线段相交时， 和 不相交。。。。郁闷死了。。 #include #include #include #include #include #include using namespace std; const double INF = 1e20; const double EPS = 1e-6; bool zero(doub
计算几何——Intersection(线段与矩形相交)
mumei的博客
07-14 4320
题目链接http://poj.org/problem?id=1410 You are to write a program that has to decide whether a given line segment intersects a given rectangle. An example: line: start point: (4,9) end point: (11,2) rect...
POJ-1410——Intersection（简单计算几何，坑题）
lee371042的博客
01-08 534
Intersection Time Limit: 1000MS   Memory Limit: 10000K Total Submissions: 17761   Accepted: 4556 Description You are to write a program that has to decide whether a given line segmen...
ACM计算几何题目推荐
智障
01-21 5725
把下面的东东都看看，题目刷刷应该就差不多了吧哈。。哈哈。。 其实也谈不上推荐，只是自己做过的题目而已，甚至有的题目尚未AC，让在挣扎中。之所以推荐计算几何题，是因为，本人感觉ACM各种算法中计算几何算是比较实际的算法，在很多领域有着重要的用途（例如本人的专业，GIS）。以后若有机会，我会补充、完善这个列表。 计算几何题的特点与做题要领： 1.大部分不会很难，少部分题目思路
C++编写判断直线是否在矩形内高效算法
11-17
判断直线是否在矩形内，不适用乘除法，高效率算法，仅供参考！
POJ 3449 Geometric Shapes（几何相交）
退役的铜牌狗
10-16 298
DescriptionWhile creating a customer logo, ACM uses graphical utilities to draw a picture that can later be cut into special fluorescent materials. To ensure proper processing, the shapes in the pictur
Intersection POJ - 1410
DeathYmz的博客
08-21 319
Intersection POJ - 1410 分类：计算几何+线段相交 题目： You are to write a program that has to decide whether a given line segment intersects a given rectangle. An example: line: start point: (4,9) end poi...
poj 1410
[CHN]liurui
07-22 487
题目概述给定线段的端点坐标(x1,y1),(x2,y2)，矩形一条对角线上的两个点(x3,y3),(x5,y5)，定义矩形的边及内部均属于矩形，相交为线段和矩形有至少一个公共点，且矩形四条边必然分别平行于两条坐标轴，问线段和矩形是否相交时限1000ms/3000ms输入第一行正整数times，其后times行，每行8个整数x1,y1,x2,y2,x3,y3,x5,y5限制没有限制输出每行一个字符，若
POJ 1410 Intersection(判断线段与矩形位置关系)
不期而遇
07-31 779
题目链接 题目大意：判断一条线段与一个实心矩形（四条边以及中间包含的部分）是否相交。 分析：线段与矩形四条边的不规范相交（交于一条线段的端点或者重合），以及线段在矩形的内部，都算作相交。这题看似简单，其实有很多比较细节的地方，比较坑，调了好久。 #include&lt;map&gt; #include&lt;set&gt; #include&lt;cmath&gt; #include&lt...
poj1410
__Praesidio_
04-02 374
Source Code Problem: 1410   User: Praesidio Memory: 180K   Time: 0MS Language: C++   Result: Accepted Source Code //poj 1410 #include #include using namespace st
POJ - 1410 Intersection
M_Lter的博客
08-25 179
题目链接 判断线段是否和矩形相交或者在矩形内或者在矩形边上，是则输出T，否则输出F。 一开始自己没看清题意，直接判断矩形四条边是否和线段相交，一直wa。看了题解才知道，还有在矩形内和矩形边上，这三种情况都单独讨论。 #include&lt;cstdio&gt; #include&lt;cmath&gt; #include&lt;cstring&gt; #include&lt;algorit...
POJ - 1410 - Intersection（计算几何）
Steve·Curcy
10-27 179
You are to write a program that has to decide whether a given line segment intersects a given rectangle. An example: line: start point: (4,9) end point: (11,2) rectangle: left-top: (1,5) right-bottom...
POJ-1410-Intersection（线段相交）
永远鲜红の幼月的博客
04-10 330
题目链接：http://poj.org/problem?id=1410 题目大意：给出一个线段和一个矩形，判断他们是否相交 思路：因为只有几种情况，全部列出来就好了，注意，线段全部在矩形中不算相交。 ACCode： //#pragma comment(linker, "/STACK:1024000000,1024000000") #include<stdio.h> #i...
计算几何：线段相交与交点判断
- 提供了多个POJ题目，涵盖了线段相交、多边形面积计算、点在多边形内的判断等计算几何问题。 10. 直线的数据结构： - 使用`struct line_t`表示直线，通常包含三个系数`a, b, c`，表示方程`ax + by + c = 0`。 -...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值