河南省ACM第一届第一题

灾区已经非常困难，灾民需要帐篷、衣物、食品和血浆。可通往灾区的道路到处都是塌方，70%以上的路面损坏，桥梁全部被毁。中国空军立即启动应急预案，展开史上最大强度非作战空运行动，准备向灾区空投急需物资。由于余震不断，天气恶劣，怎样知道空投的物资是否落在某灾区的区域内呢？
经过空中观测，某灾区为一凸多边形，空投的物资落在P(X,Y)点。你能否给出一个正确判断？

第1行： N X Y (N为凸形的边数，X，Y为空投物资的坐标 )
第2行： X1 Y1 X2 Y2…….Xn Yn （逆序给出的N个顶点坐标 ）

YES （物资落在灾区的区域内 ）
或
NO （物资落在灾区的区域外 ）

（1） 3 ≤N ≤ 20
（2） 所有的坐标Xi和Yi为整数 -10000≤ Xi ，Yi ≤10000
（3） X,Y为整数 -10000≤ X ，Y ≤10000
（4） 时间限制： 1000MS

标准输入
4 2 2
-1 0 3 -2 10 9 -1 110

标准输出

YES

题目分析：求一个给定的点A(x0,y0)是否在凸多边型p0p1...pn-1（逆序给出）内部，我们首先想到的是判断一个三角形三个点是顺时针还是逆时针的，我们可以利用这个结论。

看下图：

如果点A在多边形内部，则三角形APiP(i+1%n)必为逆时针，而如果在外部，则有些点和A组成的三角形是顺序的，多以，我们只需要依次判断APiP(i+1%n)是否为逆时针即可，一旦发现不是逆时针，则说明在凸多边型外部！

那么我们来看怎么判断三角形的正序还是逆序呢？我们有以下结论：定义向量AB=（x1,y1),AC=（x2,y2），则：
AB X AC>0  逆时针

AB X AC<0  顺时针，其中运算X（向量叉乘）定义为x1*y2-x2*y1。

#include<fstream>
#include<iostream>
using namespace std;
const int N=25;
struct node{
    int x;
    int y;
}p[N];
int fun(struct node a,struct node b,struct node c)//向量的叉乘表示逆顺时针
{
    int abx,aby,acx,acy,ans;
    abx=b.x-a.x;
    aby=b.y-a.y;
    acx=c.x-a.x;
    acy=c.y-a.y;
    ans=abx*acy-aby*acx;
    if(ans>0)
        return 1;
    else
        return 0;
}
int main()
{
ifstream cin("aaa.txt");
int n;
struct node a,b,c;
while(cin>>n>>a.x>>a.y)
{
   for(int i=0;i<n;i++)
            cin>>p[i].x>>p[i].y;
        int flag=1;
        for(int i=0;i<n;i++)
         {
                b=p[i];
                c=p[(i+1)%n];
                flag=fun(a,b,c);
                if(flag==0)
                    break;
         }
         if(flag==1)
            cout<<"YES"<<endl;
         else
            cout<<"NO"<<endl;
}
return 0;
}