POJ 3335 半平面相交算法

题意是问在多边形中是不是存在一个点，这个点可以看到多边形内部的任意一个地方。

用的是半平面相交的模板，调试了半天，最后发现竟然是m和n写错了，悲剧~

 

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<iostream>
#include<algorithm>

using namespace std;

const int maxn=105;

struct node
{
   double x,y;
}point[maxn],p[maxn],q[maxn];

int m,n;

//求点x,y形成的直线ax+by+c=0;
void getline(node x,node y,double &a,double &b,double &c)
{
   a=y.y-x.y;
   b=x.x-y.x;
   c=x.y*y.x-x.x*y.y;
}

//求xy这条直线与ax+by+c=0这条直线的交点
node intersect(node x,node y,double a,double b,double c)
{
   double u=fabs(a*x.x+b*x.y+c);
   double v=fabs(a*y.x+b*y.y+c);
   node pt;
   pt.x=(x.x*v+y.x*u)/(u+v);
   pt.y=(x.y*v+y.y*u)/(u+v);
   return pt;
}

void cut(double a,double b,double c)
{
   int cnt=0;
   //这里的m写成了n，害我调试了半天
   for(int i=1;i<=m;i++)
   {
      if(a*p[i].x+b*p[i].y+c>=0) q[++cnt]=p[i];
      else
      {
         if(a*p[i-1].x+b*p[i-1].y+c>0)
            q[++cnt]=intersect(p[i],p[i-1],a,b,c);
         if(a*p[i+1].x+b*p[i+1].y+c>0)
            q[++cnt]=intersect(p[i],p[i+1],a,b,c);
      }
   }
   //cout<<cnt<<endl;
   for(int i=1;i<=cnt;i++) p[i]=q[i];
   p[cnt+1]=q[1],p[0]=p[cnt];
   m=cnt;
}

void solve()
{
   double a,b,c;
   for(int i=1;i<=n;i++)
   {
      getline(point[i],point[i+1],a,b,c);
      cut(a,b,c);
   }
}

int main()
{
   int t;
   scanf("%d",&t);
   while(t--)
   {
      scanf("%d",&n);
      for(int i=1;i<=n;i++)
         scanf("%lf%lf",&point[i].x,&point[i].y);
      m=n;
      point[n+1]=point[1];
      for(int i=1;i<=n+1;i++)
         p[i]=point[i];
      p[0]=p[n];
      solve();
      //cout<<m<<endl;
      if(m<1) puts("NO");
      else puts("YES");
   }
  
   return 0;
}