[HDU1756]Cupid's Arrow（计算几何-点在多边形内部）

题目：

我是超链接

题解：

HDU的中文题？！
一看这个题不就是点在多边形内的裸题吗，为什么要发呢？
因为这个题目的图形不一定是凸多边形啊！
有什么关系呢？就是我们判断的时候不能只判断P点到AB的距离是不是0了，即不能仅仅判断是否共线，而是需要进一步用点积判断，就是下面这个图啦

那正确的判断姿势应该是。。

bool is_online(po p,po a,po b)
{
    return dcmp(cj(p-a,p-b))==0 && dot(a-p,b-p)<=0;
}

至于那个

        if (a>0 && b<=0 && c>0) cnt++; 
        if (a<0 && b>0 && c<=0) cnt--;

这样比较好，因为

就都不用记录进去了

代码：

#include <cstdio>
#include <cmath>
using namespace std;
const double eps=1e-9;
int dcmp(double x)
{
    if (x<=eps && x>=-eps) return 0;
    return (x>0)?1:-1;
}
struct po
{
    double x,y;
    po (double X=0,double Y=0){x=X;y=Y;}
}dian[10005];int n;
po operator -(po x,po y){return po(x.x-y.x,x.y-y.y);}
double cj(po x,po y){return x.x*y.y-x.y*y.x;}
double dot(po x,po y){return x.x*y.x+x.y*y.y;}
bool is_online(po p,po a,po b)
{
    return dcmp(cj(p-a,p-b))==0 && dot(a-p,b-p)<=0;
}
bool Pinp(po p)
{
    int cnt=0,a,b,c;
    for (int i=1;i<=n;i++)
    {
        if (is_online(p,dian[i],dian[i%n+1])) return 0;
        a=dcmp(cj(dian[i%n+1]-dian[i],p-dian[i]));
        b=dcmp(dian[i].y-p.y); 
        c=dcmp(dian[i%n+1].y-p.y);
        if (a>0 && b<=0 && c>0) cnt++;
        if (a<0 && b>0 && c<=0) cnt--;
    }
    if (cnt) return 1;
    return 0;
}
int main()
{
    while (~scanf("%d",&n))
    {
        int m;double x,y;
        for (int i=1;i<=n;i++) scanf("%lf%lf",&x,&y),dian[i]=po(x,y);
        scanf("%d",&m);
        while (m--)
        {
            scanf("%lf%lf",&x,&y);
            if (Pinp(po(x,y))) printf("Yes\n");
            else printf("No\n");
        }
    }
}