凸多边形

题目描述

灾区已经非常困难，灾民需要帐篷、衣物、食品和血浆。可通往灾区的道路到处都是塌方，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

当物资降落在区域边缘时，输出YES。

样例输入

4 2 2
-1 0 3 -2 10 9 -1 110

样例输出

YES

如果投放物资的坐标在多边形的内部，则形成的各个三角形的面积之和与多边形的面积之和肯定相等，如果投放的物资的坐标在多边形的外部，则形成的各个三角形的面积之和一定大于多边形的面积之和。

#include<stdio.h>
#include<math.h>
double xx,yy;
double S(int x[],int y[],int n)//逆时针输入n个顶点的坐标求面积
{
	int i;
    double s=0;
    x[n]=x[0];
    y[n]=y[0];
    for(i=0;i<n;i++)
    {
        s+=(x[i]*y[i+1]-x[i+1]*y[i]);
    }
    return s/2;
}
double hai(double a,double b,double c,double d)//海伦公式
{
	double s1,s2,s3,p,sum;
	s1=sqrt((a-c)*(a-c)+(b-d)*(b-d));
	s2=sqrt((a-xx)*(a-xx)+(b-yy)*(b-yy));
	s3=sqrt((c-xx)*(c-xx)+(d-yy)*(d-yy));
	p=(s1+s2+s3)/2;
	sum=sqrt(p*(p-s1)*(p-s2)*(p-s3));
	return sum;
}
int main()
{
    int n,i,x[100],y[100];
    double s,sum;
    while(scanf("%d%lf%lf",&n,&xx,&yy)!=EOF)
    {
    	sum=0;
        for(i=0;i<n;i++)
            scanf("%d%d",&x[i],&y[i]);
        s=S(x,y,n);
        x[n]=x[0];y[n]=y[0];
        for(i=0;i<n;i++)
        {
        	sum+=hai(x[i],y[i],x[i+1],y[i+1]);
        	//printf("%lf\n",sum);
		}
		if(sum-s<1e-10)
			printf("YES\n");
		else
			printf("NO\n");
//        printf("%lf\n",s);
//        printf("%lf\n",sum);
    }    
    return 0;
 }