POJ 1584 计算几何

思路：

求一遍凸包

用三角形面积（叉积求一下）/边长

求出来高，跟半径比一比

坑点：凸包上三点共线

//By SiriusRen
#include <cmath>
#include <cstdio>
#include <algorithm>
using namespace std;
#define eps 1e-9
const int N=;
int n,top,k,f;double r;
struct Point{double x,y;}point[N],tubao[N],cir;
Point operator-(Point a,Point b){Point c;c.x=a.x-b.x,c.y=a.y-b.y;return c;}
bool operator<(Point a,Point b){return abs(a.x-b.x)<eps?a.y<b.y:a.x<b.x;}
double chaji(Point a,Point b){return a.x*b.y-a.y*b.x;}
double cross(Point a,Point b,Point c){return chaji(a-c,b-c);}
double dis(Point a){return sqrt(a.x*a.x+a.y*a.y);}
bool check(Point a,Point b){return abs(chaji(a-cir,b-cir))/dis(a-b)>r-eps;}
int main(){
    while(scanf("%d",&n)){
        if(n<)break;top=f=;
        scanf("%lf%lf%lf",&r,&cir.x,&cir.y);
        for(int i=;i<=n;i++)scanf("%lf%lf",&point[i].x,&point[i].y);
        sort(point+,point++n);
        for(int i=;i<=n;i++){
            while(top>&&cross(tubao[top],point[i],tubao[top-])<-eps)top--;
            tubao[++top]=point[i];
        }k=top;
        for(int i=n-;i;i--){
            while(top>k&&cross(tubao[top],point[i],tubao[top-])<-eps)top--;
            tubao[++top]=point[i];
        }
        if(top!=n+){puts("HOLE IS ILL-FORMED");continue;}
        for(int i=;i<=n;i++)if(!check(tubao[i],tubao[i+]))f=;
        for(int i=;i<=n;i++)
            if(chaji(tubao[i+]-tubao[i],cir-tubao[i])<)f=;
        printf("PEG WILL %sFIT\n",f?"NOT ":"");
    }
}