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

#include<iostream>
#include<cstdio>
using namespace std;
double a[101],b[101];
int pnpoly(int nvert, double *vertx, double *verty, double testx, double testy)///nvert顶点数，vertx顶点x坐标数组，verty顶点y坐标数组,testx,testy判断坐标
{
    int i, j, c = 0;
    for (i = 0, j = nvert-1; i < nvert; j = i++)
    {
        if ( ((verty[i]>testy) != (verty[j]>testy)) &&(testx < (vertx[j]-vertx[i]) * (testy-verty[i]) / (verty[j]-verty[i]) + vertx[i]) )
            c = !c;
    }
    return c;
}

int main()
{
    int n;
    double p1,p2;
    while (scanf("%d", &n) != EOF)
    {
        for (int i=0; i<n; ++i) scanf ("%lf%lf", &a[i], &b[i]);
        int m;
        scanf ("%d", &m);
        while (m--)
        {
            scanf ("%lf%lf", &p1, &p2);

            if (pnpoly(n,a,b,p1,p2)) printf ("Yes\n");
            else printf ("No\n");
        }
    }
    return 0;
}