HDU 1756 Cupid's Arrow
Problem Description
传说世上有一支丘比特的箭,凡是被这支箭射到的人,就会深深的爱上射箭的人。
世上无数人都曾经梦想得到这支箭。Lele当然也不例外。不过他想,在得到这支箭前,他总得先学会射箭。
日子一天天地过,Lele的箭术也越来越强,渐渐得,他不再满足于去射那圆形的靶子,他开始设计各种各样多边形的靶子。
不过,这样又出现了新的问题,由于长时间地练习射箭,Lele的视力已经高度近视,他现在甚至无法判断他的箭射到了靶子没有。所以他现在只能求助于聪明的Acmers,你能帮帮他嘛?
Input
本题目包含多组测试,请处理到文件结束。
在每组测试的第一行,包含一个正整数N(2<N<100),表示靶子的顶点数。
接着N行按顺时针方向给出这N个顶点的x和y坐标(0<x,y<1000)。
然后有一个正整数M,表示Lele射的箭的数目。
接下来M行分别给出Lele射的这些箭的X,Y坐标(0<X,Y<1000)。
Output
对于每枝箭,如果Lele射中了靶子,就在一行里面输出"Yes",否则输出"No"。
Sample Input
4
10 10
20 10
20 5
10 5
2
15 8
25 8
Sample Output
Yes
No
【思路分析】
该题也是判断点是否在多边形上,这里介绍射线法,虽然相对于二分法的时间效率较慢,但是射线法适用于任何多边形的判断:
射线法就是选取多边形外一点和所需判断的点连成一条线,若该直线和多边形的交点个数为奇数,则该点在多边形内,否则,该点在多边形之外。但是要考虑以下几种情况:
即射线过多边形的顶点和边(1、2),射线过多边形的顶点(3),射线过多边形的顶点且和多边形的边共线(4)。
当满足以下的条件时射线与多边形的交点的个数num才自加一:(设多边形的边为bc,要判断的点为a)
1、c.y != d.y,即针对上图中的4。当射线与边共线时,不计入交点个数(若射线平行边时则该射线不可能和该边有交点);
2、a.y > min(c.y,d.y) && a.y <= max(c.y,d.y),即针对上图中的1、2、3。当射线和多边形的顶点相交时,若顶点为该条边上y值最小的点时,则不计入交点个数。例如1中的边x和边y的y值最小的点和射线相交,则不计入交点个数;2中边x的顶点不为x边上y值最小的点,则该点计入交点个数,边y同理,边z同理;3中的边x计入,y不计入。
3、射线与该边有交点。
代码如下:
#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <cstring>
#include <string>
#include <algorithm>
using namespace std;
#define maxn 100005
#define EPS 1e-8
struct Point
{
double x;
double y;
}p[maxn];
double Max(double a,double b)
{
return a > b ? a : b;
}
double Min(double a,double b)
{
return a < b ? a : b;
}
double Abs(double a)
{
return a > 0 ? a : (- a);
}
double cross(double x1,double y1,double x2,double y2)
{
return x1 * y2 - y1 * x2;
}
bool isOnline(Point a,Point b,Point c)
{
if(Abs(cross(a.x - b.x,a.y - b.y,c.x - a.x,c.y - a.y)) < EPS)
{ //判断a是否在直线bc上,即判断向量ab,ca是否共线
if(a.x >= Min(b.x,c.x) && a.x <= Max(b.x,c.x) && a.y >= Min(b.y,c.y) && a.y <= Max(b.y,c.y))
{//判断a点是否在线段bc上
return true;
}
}
return false;
}
int main()
{
int n,m,num,flag;
Point a,b,c,d;
while(scanf("%d",&n) != EOF)
{
for(int i = 0;i < n;i++)
{
scanf("%lf %lf",&p[i].x,&p[i].y);
}
scanf("%d",&m);
while(m--)
{
flag = num = 0;
scanf("%lf%lf",&a.x,&a.y);
b.x = -100;
b.y = a.y;//ab即为穿过多边形的射线
for(int i = 0;i < n;i++)
{
c.x = p[i].x;
c.y = p[i].y;
d.x = p[(i + 1) % n].x;
d.y = p[(i + 1) % n].y;
if(isOnline(a,c,d) == true)//判断a是否在cd边上
{
flag = 1;
break;
}
if(c.y != d.y)//要求射线不与cd共线(若平行也不符合条件)
{
if(a.y > Min(c.y,d.y) && a.y <= Max(c.y,d.y))//射线不过多边形的“凹顶点”
{
if(cross(d.x - a.x,d.y - a.y,c.x - a.x,c.y - a.y)
* cross(d.x - b.x,d.y - b.y,c.x - b.x,c.y - b.y) < EPS)//判断ab和cd是否相交
{
num++;
}
}
}
}
if(flag == 1)
printf("Yes\n");
else
{
if(num & 1 == 1)//num为奇数个时
printf("Yes\n");
else
printf("No\n");
}
}
}
return 0;
}