判断点与多边形的位置关系

hdu1756

Cupid's Arrow

Time Limit: 3000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)

Total Submission(s): 1033    Accepted Submission(s): 382

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

410 1020 1020 510 5215 825 8

 

Sample Output

YesNo
程序：
#include"stdio.h"#include"string.h"#include"iostream"#include"map"#include"string"#include"queue"#include"stack"#include"vector"#include"stdlib.h"#include"algorithm"#include"math.h"#define M 109#define eps 1e-10#define inf 0x3f3f3f3f#define mod 1070000009double const PI=acos(-1.0);struct node{    double x,y;}p[M];struct line{    node s,e;};double max(double x,double y){    return x>y?x:y;}double min(double x,double y){    return x<y?x:y;}double cross(node a,node b,node c){    double px1=b.x-a.x;    double py1=b.y-a.y;    double px2=c.x-a.x;    double py2=c.y-a.y;    return px1*py2-py1*px2;}int paichi(line a,line b){    if(max(a.e.x,a.s.x)>min(b.s.x,b.e.x)        &&max(b.e.x,b.s.x)>min(a.s.x,a.e.x)        &&max(a.s.y,a.e.y)>min(b.s.y,b.e.y)        &&max(b.s.y,b.e.y)>min(a.s.y,b.e.y))        return 1;    return 0;}int kuali(line a,line b){    if(cross(a.s,a.e,b.s)*cross(a.s,a.e,b.e)<=0        &&cross(b.s,b.e,a.s)*cross(b.s,b.e,a.e)<=0)        return 1;    return 0;}int judge(node q,int n){    int flag;    int i=0;    while(i<=n)    {        node po;        po.x=rand()+1000.0;        po.y=rand()+1000.0;        flag=0;        for(i=1;i<=n;i++)        {            line a,b;            a.s=q;            a.e=po;            b.s=p[i-1];            b.e=p[i];            if(q.x>=min(p[i-1].x,p[i].x)               &&q.x<=max(p[i-1].x,p[i].x)               &&q.y>=min(p[i-1].y,p[i].y)               &&q.y<=max(p[i-1].y,p[i].y)&&fabs(cross(p[i-1],q,p[i]))<eps)                return -1;//在边上            else if(fabs(cross(po,q,p[i]))<eps)                break;            else if(kuali(a,b)&&paichi(a,b))                flag++;        }    }    return flag&1;//若是0在外面，否则在内部}int main(){    int n,i,m;    while(scanf("%d",&n)!=-1)    {        for(i=1;i<=n;i++)            scanf("%lf%lf",&p[i].x,&p[i].y);        p[0]=p[n];        scanf("%d",&m);        while(m--)        {            node g;            scanf("%lf%lf",&g.x,&g.y);            if(judge(g,n))                printf("Yes\n");            else                printf("No\n");        }    }    return 0;}