HDU 1756 判断点在多边形内外

简介:

题意:判断点在多边形内外,用来试模版的

#include <iostream>
#include <stdio.h>
#include <math.h>

typedef double DIY;

const DIY EPS = 0;
const DIY N = 1000005;

using namespace std;
struct Point
{
    DIY x,y;
};
struct Segment
{
    Point a;
    Point b;
};
typedef Point point;
Point p[100];
DIY cross(Point A,Point B,Point C)
{
    return (B.x-A.x)*(C.y-A.y)-(B.y-A.y)*(C.x-A.x);
}
bool IsOnline(Point p,Segment s)
{
    if(fabs(cross(s.a,p,s.b))<=EPS&&min(s.a.x,s.b.x)<=p.x&&
            p.x<=max(s.a.x,s.b.x)&&min(s.a.y,s.b.y)<=p.y&&p.y<=max(s.a.y,s.b.y))
        return true;
    return false;
}
bool Segment_crossing(Segment u,Segment v)  /*判断两线段是否相交*/
{
    return((max(u.a.x,u.b.x)>=min(v.a.x,v.b.x))&&
           (max(v.a.x,v.b.x)>=min(u.a.x,u.b.x))&&
           (max(u.a.y,u.b.y)>=min(v.a.y,v.b.y))&&
           (max(v.a.y,v.b.y)>=min(u.a.y,u.b.y))&&
           (cross(v.a,u.b,u.a)*cross(u.b,v.b,u.a)>=0)&&
           (cross(u.a,v.b,v.a)*cross(v.b,u.b,v.a)>=0));
}
bool Inside(Point t,Point p[],int n)     /*判断点p是否在以边点顺序排列的点集s所围成的多边形内*/
{
    Segment seg;
    Segment pp;
    int num=0;
    Point MM;          /*定义无穷点*/
    MM.y=t.y;
    MM.x=N;
    pp.a=t;
    pp.b=MM;           /*构造一条平行于x轴的以p为端点的右射线pp*/
    for(int i=0; i<n; i++)
    {
        seg.a=p[i];
        seg.b=p[(i+1)%n];
        if(IsOnline(t,seg))
            return 1;
        if(fabs(seg.a.y-seg.b.y)<EPS)
            continue;
        if(IsOnline(seg.a,pp))
        {
            if(seg.a.y>seg.b.y)
                num++;
        }
        else if(IsOnline(seg.b,pp))
        {
            if(seg.b.y>seg.a.y)
                num++;
        }
        else if(Segment_crossing(pp,seg))
            num++;
    }
    if(num%2==1)
        return true;
    return false;
}
int main()
{
    Point t;
    int n,m;
    while(cin>>n)
    {
        for(int i=0; i<n; i++)
            scanf("%lf%lf",&p[i].x,&p[i].y);
        cin>>m;
        while(m--)
        {
            scanf("%lf%lf",&t.x,&t.y);
            if(Inside(t,p,n)) puts("Yes");
            else              puts("No");
        }
    }
    return 0;
}


目录
相关文章
|
算法 定位技术 数据安全/隐私保护
简谈百度坐标反转至WGS84的三种思路
文章版权由作者李晓晖和博客园共有,若转载请于明显处标明出处:http://www.cnblogs.com/naaoveGIS/ 1.背景 基于百度地图进行数据展示是目前项目中常见场景,但是因为百度地图是基于BD09坐标系的,GPS坐标(WGS84)或者其他常见的标准坐标是无法准确在地图上进行展示的,但是互联网在线情况下,百度提供了将WGS84经纬度转换成百度经纬度坐标的API,这里不再对其进行研究(离线情况下也有专门方法解决)。
1716 0
|
22天前
|
C++
C++代码来计算一个点围绕另一个点旋转45度后的坐标
C++代码来计算一个点围绕另一个点旋转45度后的坐标
42 0
|
1月前
|
算法 数据建模
平面中判断点在三角形内算法(重心法)
平面中判断点在三角形内算法(重心法)
24 0
|
4月前
|
C++
[C++/PTA] 判断一个点是否在一个圆的内部
[C++/PTA] 判断一个点是否在一个圆的内部
64 0
|
Java
判断顶点凹凸性、判断多边形的凹凸性、填充凹坑将凹多边形处理为凸多边形【java实现+原理讲解】
判断顶点凹凸性、判断多边形的凹凸性、填充凹坑将凹多边形处理为凸多边形【java实现+原理讲解】
207 0
判断顶点凹凸性、判断多边形的凹凸性、填充凹坑将凹多边形处理为凸多边形【java实现+原理讲解】
射线法——判断一个点是否在多边形内部(适用于凸多边形和凹多边形)【关键原理解释+文字伪代码】
射线法——判断一个点是否在多边形内部(适用于凸多边形和凹多边形)【关键原理解释+文字伪代码】
587 0
每日三题-最大正方形 、完全平方数 、目标和
每日三题-最大正方形 、完全平方数 、目标和
78 2
每日三题-最大正方形 、完全平方数 、目标和
PTA 7-1 打印三角形拼图 (15 分)
一个正方形可以用两个等边直角三角形拼出来。给定正方形的边长、两个三角形和对角线所用的符号,请你打印出这两个三角形拼出的正方形。
119 0
|
算法
判断三角形的性质(直角或等腰)简便算法
判断三角形的性质(直角或等腰)简便算法
118 0