A Round Peg in a Ground Hole POJ - 1584

A Round Peg in a Ground Hole POJ - 1584
分类：判断凸包+圆在凸包内
题意：给出n，和圆心的圆心点和半径，在给出你n个按顺序的点（注意这里不一定说按照顺时针逆时针的）。判断否是凸包，，是的话圆是否能完全在凸包内，相切是可以的。
思路：向kuangbin大佬模板低头，真好使。整个题套上模板差不多改一下细节就好了。emmm。
1、这里不可以直接用判断点是否在凸包内的哪个，因为：点形成一个凸包，而且按逆时针排序（如果是顺时针把里面的<0改为>0）
2、判断凸包->判断圆心是否在多边形内—>判断点到线段的距离是否大于半径

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cmath>
#include<cstring>
using namespace std;
const int N=200;
const double eps=1e-8; 
int sgn(double x)
{
    if(fabs(x)<eps) return 0;
    if(x<0) return -1;
    else return 1;
}
struct point{
    double x,y;
    point(){}
    point(double a,double b){x=a;y=b;}
    point operator - (const point &b) const
    {
        return point(x-b.x,y-b.y);
    }
    //叉积
    double operator ^ (const point &b) const
    {
        return x*b.y-y*b.x;
    }
    //点积
    double operator * (const point &b) const
    {
        return x*b.x+y*b.y;
    }
    bool operator <(const point &u)const{
        if(y==u.y) return x<u.x;
        return y<u.y;
    }
}p[N];//存储
struct Line
{
    point s,e;
    Line(){}
    Line(point _s,point _e)
    {
        s = _s;e = _e;
    }
};
//判断凸包
bool isconvex(int n)
{
    bool s[3];
    memset(s,false,sizeof(s));
    for(int i=0;i<n;i++)
    {
        s[sgn((p[(i+1)%n]-p[i])^(p[(i+2)%n]-p[i]))+1]=true;
        if(s[0]&&s[2]) return false;
    }
    return true;//是凸包
}
//*判断点在线段上
bool OnSeg(point p0,point p1,point p2)
{
    return
    sgn((p1-p0)^(p2-p0))==0&&
    sgn((p0.x-p1.x)*(p0.x-p2.x))<=0 &&
    sgn((p0.y-p1.y)*(p0.y-p2.y))<=0;
}
//判断点在tu边形内，角度和判断/凸包特有
//*判断点在凸多边形内
//点形成一个凸包，而且按逆时针排序（如果是顺时针把里面的<0改为>0）
//点的编号:0~n-1
//返回值：
//-1:点在凸多边形外
//0:点在凸多边形边界上
//1:点在凸多边形内
int inConvexPoly(point a,int n)
{
    for(int i = 0;i < n;i++)
    {
        if(sgn((p[i]-a)^(p[(i+1)%n]-a)) < 0)return -1;//在外面
        else if(OnSeg(a,p[i],p[(i+1)%n]))return 0;//在边界
    }
    return 1;//在里面
}

//*判断线段相交
bool inter(Line l1,Line l2)
{
    return
    max(l1.s.x,l1.e.x) >= min(l2.s.x,l2.e.x) &&
    max(l2.s.x,l2.e.x) >= min(l1.s.x,l1.e.x) &&
    max(l1.s.y,l1.e.y) >= min(l2.s.y,l2.e.y) &&
    max(l2.s.y,l2.e.y) >= min(l1.s.y,l1.e.y) &&
    sgn((l2.s-l1.e)^(l1.s-l1.e))*sgn((l2.e-l1.e)^(l1.s-l1.e)) <= 0 &&
    sgn((l1.s-l2.e)^(l2.s-l2.e))*sgn((l1.e-l2.e)^(l2.s-l2.e)) <= 0;
}
int inPoly(point p,point poly[],int n)
{
    int cnt;
    Line ray,side;
    cnt = 0;
    ray.s = p;
    ray.e.y = p.y;
    ray.e.x = -100000000000.0;//-INF,注意取值防止越界

    for(int i = 0;i < n;i++)
    {
        side.s = poly[i];
        side.e = poly[(i+1)%n];

        if(OnSeg(p,side.s,side.e))return 0;

        //如果平行轴则不考虑
        if(sgn(side.s.y - side.e.y) == 0)
            continue;

        if(OnSeg(side.s,ray.s,ray.e))
        {
            if(sgn(side.s.y - side.e.y) > 0)cnt++;
        }
        else if(OnSeg(side.e,ray.s,ray.e))
        {
            if(sgn(side.e.y - side.s.y) > 0)cnt++;
        }
        else if(inter(ray,side))
            cnt++;
    }
    if(cnt % 2 == 1)return 1;
    else return -1;
}
//求点到线段距离
//返回点到线段的最近的点
double dis(point p1,point p2)
{
    return sqrt((p1-p2)*(p1-p2));
}
point nearpoint(point p0,point p1,point p2)
{
    point result;
    double t=((p0-p1)*(p2-p1))/((p2-p1)*(p2-p1));
    if(t>=0&&t<=1)
    {
        result.x=p1.x+(p2.x-p1.x)*t;
        result.y=p1.y+(p2.y-p1.y)*t;
    }
    else
    {
        if(dis(p0,p1)<dis(p0,p2))
            result=p1;
        else result=p2;
    }
    return result;
}

//判断是否在凸变形的内部，相切也可以
bool judge(point p1,int n,double r)
{
    int flag=inPoly(p1,p,n);
   if(!flag)//圆心在边上
    {
        if(r==0) return 1;
        else return 0;
    }
    if(flag==-1) return 0;//不在多边形外部
    for(int i=0;i<n;i++) 
    {
         point tmp=nearpoint(p1,p[i],p[(i+1)%n]);
         double distance=dis(p1,tmp);
        if(sgn(distance-r)<0) return 0; //点到直线的距离
    }
    return 1;
}

int main()
{
    int n;
    double r,px,py;
    while(scanf("%d",&n)!=EOF)
    {
        if(n<3) break;
        scanf("%lf%lf%lf",&r,&px,&py);
        for(int i=0;i<n;i++) 
        scanf("%lf%lf",&p[i].x,&p[i].y);
        int flag=isconvex(n);//判断是否为凸边形
        if(!flag) printf("HOLE IS ILL-FORMED\n");
        else
        {
            if(judge(point(px,py),n,r)) printf("PEG WILL FIT\n");
            else printf("PEG WILL NOT FIT\n");
        }
    }
    return 0;
}