POJ 1584[A Round Peg in a Ground Hole]题解

题目梗概

题目给出了一个多边形，和一个已知半径和圆心的圆，判断：
1.该多边形是否是凸多边形。
2.这个圆能否被这个多边形完全覆盖。
（传送门）

那么需要解决的，有3个问题。
1.该多边形是否是凸多边形。
2.圆心（某个点）是否在凸多边形内
3.圆心（某个点）对所有线段的距离是否都小于圆的半径。

关于问题1：在学习笔记2中有提过，打个广告；
关于问题2，本来是有通解的，但是既然验证时凸多边形通过滴，只需利用叉积判断是否在所有边的同一侧即可。
关于问题3：在学习笔记1中有提过，打个广告；

搞定
哦还有，打个广告，不，膜拜—— Orz ZH大佬 Orz Lynstery

代码

#include<cstdio>
#include<cmath>
using namespace std;
int n;
double R,eps=1e-10;
struct data
{
  double x,y;
  data(double x=0,double y=0):x(x),y(y) { }
}o,a[10005];
double absf(double x){return (x>0)?x:-x;}
int fcmp(double x,double y) {if (absf(x-y)<eps) return 0; return (x<y)?-1:1;}
data operator + (const data a,const data b){return data(a.x+b.x,a.y+b.y);}
data operator - (const data a,const data b){return data(a.x-b.x,a.y-b.y);}
data operator * (const data a,double b){return data(a.x*b,a.y*b);}
double dot(const data a,const data b) {return a.x*b.x+a.y*b.y;}
double leng(const data a){return sqrt(dot(a,a));}
double cross(const data a,const data b) {return a.x*b.y-a.y*b.x;}
double dislne(data p,data a,data b){
  data Va=b-a,Vb=p-a;
  return absf(cross(Va,Vb))/leng(Va);
}
bool check_Tu()
{
    bool f=true;
    for (int i=1;i<=n;i++)
      if (fcmp(cross(a[i%n+1]-a[i],a[(i-2+n)%n+1]-a[i]),0)<0) {f=false; break;}
    if (f) return true;
    for (int i=1;i<=n;i++)
      if (fcmp(cross(a[(i+n-2)%n+1]-a[i],a[i%n+1]-a[i]),0)<0) return false;
    return true;
}
bool _check()
{
    bool ga,gb; ga=gb=true;
    for (int i=1;i<=n;i++) if (fcmp(cross(o-a[i],o-a[i%n+1]),0)<0) {ga=false; break;}
    for (int i=1;i<=n;i++) if (fcmp(cross(o-a[i%n+1],o-a[i]),0)<0) {gb=false; break;}
    if (!ga&&!gb) return false;
    for (int i=1;i<=n;i++) if (fcmp(dislne(o,a[i],a[i%n+1]),R)<0) return false;
    return true;
}
int main()
{
  double x,y;
  for (scanf("%d",&n);n>2;scanf("%d",&n))
  {
    scanf("%lf%lf%lf",&R,&x,&y); o=data(x,y);
    for (int i=1;i<=n;i++){scanf("%lf%lf",&x,&y); a[i]=data(x,y);}
    if (!check_Tu()) printf("HOLE IS ILL-FORMED\n");
    else if (!_check()) printf("PEG WILL NOT FIT\n");
    else printf("PEG WILL FIT\n");
  }
  return 0;
}