HDU 1756 点在多边形内

题意：

判断点在多边形内

题解：

随机一个区域外的点，与所有线段求交，如果点数为奇数，则在多边形内，反之，在多边形外

需要注意:

1、与多边形的顶点相交时，重新随机一个点

2、点若在多边形上的话，直接返回奇数（我就在这里tle了。。一直随机，不停了。。。）

#include <iostream>
#include <cstdlib>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <ctime>

#define N 222
#define EPS 1e-7

using namespace std;

struct PO
{
    double x,y;
}p[N];

struct LI
{
    PO a,b;
}li[N];

int n,m;

inline int doublecmp(double x)
{
    if(x>EPS) return 1;
    else if(x<-EPS) return -1;
    return 0;
}

inline double cross(PO &a,PO &b,PO &c)
{
    return (c.x-a.x)*(b.y-a.y)-(c.y-a.y)*(b.x-a.x);
}

inline void read()
{
    for(int i=1;i<=n;i++) scanf("%lf%lf",&p[i].x,&p[i].y);
    for(int i=1;i<n;i++) li[i].a=p[i],li[i].b=p[i+1];
    li[n].a=p[n]; li[n].b=p[1];
    scanf("%d",&m);
}

inline int segcross(LI &a,LI &b)
{
    int p1,p2,d1,d2;
    //b跨立a 
    p1=doublecmp(cross(a.a,a.b,b.a));
    p2=doublecmp(cross(a.a,a.b,b.b));
    //a跨立b 
    d1=doublecmp(cross(b.a,b.b,a.a));
    d2=doublecmp(cross(b.a,b.b,a.b));
    if(p1==0||p2==0) return -1;//多边形的边的端点在射线上 
    else if(p1*p2<=0&&d1*d2<=0) return 1;
    else return 0;
}

inline bool inside(PO &a,PO &b,PO &c)//点在直线上 
{
    if(!doublecmp(cross(a,b,c))&&doublecmp((b.x-a.x)*(c.x-a.x))<=0&&doublecmp((b.y-a.y)*(c.y-a.y))<=0) return true;
    return false;
}

inline int getnum()
{
    int ans,i=0,tmp;
    while(i<=n)
    {
        li[i].b.x=rand()+1000.0;
        li[i].b.y=rand()+1000.0;
        ans=0;
        for(i=1;i<=n;i++)
        {
            if(inside(li[0].a,li[i].a,li[i].b)) return 1;//点在多边形上 
            tmp=segcross(li[0],li[i]);
            if(tmp==-1) break;
            else if(tmp==1) ans++;
        }
    }
    return ans;
}

inline void go()
{
    int ans;
    while(m--)
    {
        scanf("%lf%lf",&li[0].a.x,&li[0].a.y);
        if(getnum()&1) puts("Yes");
        else puts("No");
    }
}

int main()
{
    while(scanf("%d",&n)!=EOF) read(),go();
    return 0;
}