计蒜客 羊圈复原 (凸包模板)

来源:互联网 发布:消费观念 知乎 编辑:程序博客网 时间:2024/05/15 06:51

题目链接:https://www.jisuanke.com/course/615/28221

题目大意:判断羊圈是否能够复原。

题目思路:首先先求出凸包,如果凸包每条边上都有点存在的话,那么这凸包和原来是一样的,否则有可能和原来不一样。也就是说你需要来判断两点之间是否是紧挨着的,如果是则不满足。

坑点:用graham_scan的算法,会忽略一开始时三点共线的情况,所以要考虑这种情况,否则会出现各种莫名其妙的错误。三点共线!!!三点共线!!!三点共线!!!

代码:

#include <bits/stdc++.h>using namespace std;const int N=10001;const double eps=1e-10;struct node{    double x,y;}point[N],ch[N];int n,len;double multi(node a,node b,node c) //a在cb左边则小于0 ,a在cb的右边则大于0{    return (a.x-c.x)*(b.y-c.y)-(a.y-c.y)*(b.x-c.x);}double dis(node a,node b)  //距离相同去最短{    return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));}void graham_scan(node point[N],node ch[N],int n){    int i,j,k,top;    struct node t;    k=0; //找到最下且偏左的那个点    for(i=1;i<n;i++)        if(point[i].y<point[k].y||(point[i].y==point[k].y&&(point[i].x<point[k].x)))            k=i;    t=point[0];//将这个点指定为point[0];    point[0]=point[k];    point[k]=t;     //按极角从小到大,距离偏短进行排序    for(i=1;i<n-1;i++)    {        k=i;//选择排序,每次至少找到一个极角最小的点        for(j=i+1;j<n;j++)           if(multi(point[j],point[k],point[0])>0||(fabs(multi(point[j],point[k],point[0]))<=eps&&  //multi为toleft检测,如果大于0说明被测的在右边极角小               (dis(point[0],point[j])<dis(point[0],point[k]))))               k=j; //k保存极角最小的那个点,或者相同距离原点最近       t=point[i];       point[i]=point[k];       point[k]=t;    }      //第三个点先入栈    ch[0]=point[0];    ch[1]=point[1];    ch[2]=point[2];    top=2; //判断与其余所有点的关系    for(i=3;i<n;i++)    {     //不满足向左转的关系,栈顶元素出栈      while(multi(point[i],ch[top],ch[top-1])>0||fabs(multi(point[i],ch[top],ch[top-1]))<=eps)          top--;      //当前点与栈内所有点满足向左关系,因此入栈      ch[++top]=point[i];    }    len=top+1;}bool judge(node a,node b){    if(a.x==b.x&&a.y==b.y) return true;    else return false;}int main(){     int i;     scanf("%d",&n);     for(i=0;i<n;i++)         scanf("%lf%lf",&point[i].x,&point[i].y);     graham_scan(point,ch,n);     //for(int i=0;i<len;i++) cout<<ch[i].x<<" "<<ch[i].y<<endl;     int flag=0;     int num;     if(fabs(multi(ch[2],ch[1],ch[0]))<=eps) i=3,num=3;     else i=0,num=0;     for(;i<len-1;i++){        if(flag) break;        while(!judge(point[num],ch[i])) num++;        if(judge(point[num+1],ch[i+1])) flag=1;        num++;     }     if(flag) printf("NO\n");     else printf("YES\n"); return 0;}


原创粉丝点击