hdu__1756 Cupid's Arrow

http://acm.hdu.edu.cn/showproblem.php?pid=1756

/**
题意:题目是中文的
分析:首先我们输入的多边形有可能是凸多边形，也可能是凹多边形，所以我们这里使用弧长法，
弧长法适用于所有的多边形。接下来我们简单讲解一下弧长法
弧长法要求多边形是有向多边形，因为我们输入多边形的时候是是按照顺序读入的所以符合第一点
有向我们只要自己规定哪边为正向就可以了，一般是顺时针方向，规定沿着正向边的左侧为多边形
的内侧，然后再以被测点作单位园，将全部有向边向单位园作径向投影，并计算其在的单位园上的
弧长的代数和。代数和若为0，则点在多边形外部，若代数和为2π,则点在多边形里面，若代数和为
π，则点在多边形的边上。还有就是被测点和多边形的顶点重合，我们要判断是否重合和。
**/

参考《计算机图形学》

和zoj 的1081算法一样

http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemCode=1081

 

 

 

#include<stdio.h>
#include<math.h>
#define PI acos((double)-1)
const double eps=1e-8;
struct point
{
   double x,y;    
}p[101],d[1001];
double Dis(point a,point b)         //计算两点距离的平方
{
   return  (a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y);    
}
int main()
{
    int n,m,i,j;
 while(scanf("%d",&n)!=EOF){
  //   printf("%.9lf\n",PI);
    for(i=0;i<n;i++)
    scanf("%lf%lf",&p[i].x,&p[i].y);
    p[n]=p[0];
    scanf("%d",&m);
    for(i=1;i<=m;i++)
    {
      scanf("%lf%lf",&d[i].x,&d[i].y);
      double ans=0,t,r;
      for(j=1;j<=n;j++)
      {         //如果目标点在多边形边上，则跳过
      if((d[i].x-p[j-1].x)*(p[j].y-p[j-1].y)==(d[i].y-p[j-1].y)*(p[j].x-p[j-1].x))
      continue;
        t=Dis(d[i],p[j-1])+Dis(d[i],p[j])-Dis(p[j],p[j-1]);
        r=acos(t/(2*sqrt(Dis(d[i],p[j-1]))*sqrt(Dis(d[i],p[j]))));
        //求向量叉乘，以此来判断夹角的方向
        double temp=(p[j-1].x-d[i].x)*(p[j].y-p[j-1].y)-(p[j-1].y-d[i].y)*(p[j].x-p[j-1].x);
        if(temp<0)  r=-r;
        ans+=r;
   }
   if(fabs(ans)==PI)  printf("Yes\n");
   else if(fabs(ans)<eps)  printf("No\n");
   else printf("Yes\n");
       }        
    }
 return 0; 
}