Grandpa's Estate POJ 1228 （凸包）

于 2019-06-26 19:37:18 发布
阅读量180 收藏
点赞数
分类专栏： POJ 凸包&旋转卡壳
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/xiaoshazheng/article/details/93766097
版权

题目地址：http://poj.org/problem?id=1228
凸包丢失了若干个点，给你剩余的点，求这些剩余的点能不能构成一个唯一的凸包，即剩余点加上任意点（剩余点还在凸包边上），只能构成一个凸包。
在剩余点中，一遍至少三个点才能保持唯一性，两点的话可以在两点之间的外面加一点，三点加一点的话，原先三点中的中间的点会失效。
一种错误的写法
建两个凸包，一个允许多点共线，一个不允许，看这两个凸包是否存在边相等，


#include<stdio.h>
#include<string.h>
#include<math.h>
#include<algorithm>
#include<iostream>
using namespace std;
const int MAXN=50098;
const double PI=acos(-1.0);

struct point
{
    int x,y;
};
point list[MAXN];
int stack1[MAXN],stack2[MAXN],top1,top2;

int cross(point p0,point p1,point p2) //计算叉积  p0p1 X p0p2
{
    return (p1.x-p0.x)*(p2.y-p0.y)-(p1.y-p0.y)*(p2.x-p0.x);
}
double dis(point p1,point p2)  //计算 p1p2的 距离
{
    return sqrt((double)(p2.x-p1.x)*(p2.x-p1.x)+(p2.y-p1.y)*(p2.y-p1.y));
}
bool cmp(point p1,point p2) //极角排序函数 ， 角度相同则距离小的在前面
{
    int tmp=cross(list[0],p1,p2);
    if(tmp>0) return true;
    else if(tmp==0&&dis(list[0],p1)<dis(list[0],p2)) return true;
    else return false;
}
void init(int n) //输入，并把  最左下方的点放在 list[0]  。并且进行极角排序
{
    int i,k;
    point p0;
    scanf("%d%d",&list[0].x,&list[0].y);
    p0.x=list[0].x;
    p0.y=list[0].y;
    k=0;
    for(i=1;i<n;i++)
    {
        scanf("%d%d",&list[i].x,&list[i].y);
        if( (p0.y>list[i].y) || ((p0.y==list[i].y)&&(p0.x>list[i].x)) )
        {
            p0.x=list[i].x;
            p0.y=list[i].y;
            k=i;
        }
    }
    list[k]=list[0];
    list[0]=p0;
    sort(list+1,list+n,cmp);
}


void graham1(int n)
{
    int i;
    if(n==1) {top1=0;stack1[0]=0;}
    if(n==2)
    {
        top1=1;
        stack1[0]=0;
        stack1[1]=1;
    }
    if(n>2)
    {
        for(i=0;i<=1;i++) stack1[i]=i;
        top1=1;

        for(i=2;i<n;i++)
        {
            while(top1>0&&cross(list[stack1[top1-1]],list[stack1[top1]],list[i])<=0) top1--;
            top1++;
            stack1[top1]=i;
        }
    }
}
void graham2(int n)
{
    int i;
    if(n==1) {top2=0;stack2[0]=0;}
    if(n==2)
    {
        top2=1;
        stack2[0]=0;
        stack2[1]=1;
    }
    if(n>2)
    {
        for(i=0;i<=1;i++) stack2[i]=i;
        top2=1;

        for(i=2;i<n;i++)
        {
            while(top2>0&&cross(list[stack2[top2-1]],list[stack2[top2]],list[i])<0) top2--;
            top2++;
            stack2[top2]=i;
        }
    }
}
bool judge()
{
    int j=0;
    top1++;
    top2++;
    for(int i=0;i<top1;i++)
    {
        while(stack1[i]!=stack2[j])
            j=(j+1)%top2;
        if(stack1[i]==stack2[j]&&stack1[(i+1)%top1]==stack2[(j+1)%top2])
            return  false;
    }
    return  true;
}
int main()
{

   int t,n;
   scanf("%d",&t);
   while(t--)
   {
        scanf("%d",&n);
        memset(stack1,0,sizeof(stack1));
           init(n);
            graham1(n);
            graham2(n);
        if(judge())
            printf("YES\n");
        else
            printf("NO\n");
   }


    return 0;
}

一种AC的：
建一个允许多点共线的凸包，求他与相邻两点或其后两点是否共线

#include<stdio.h>
#include<math.h>
#include<algorithm>
#include<iostream>
#include<string.h>
using namespace std;
const int MAXN=1005;
const double PI=acos(-1.0);

struct point
{
    int x,y;
};
point list[MAXN];
int stack[MAXN],top;

int cross(point p0,point p1,point p2) //计算叉积  p0p1 X p0p2
{
    return (p1.x-p0.x)*(p2.y-p0.y)-(p1.y-p0.y)*(p2.x-p0.x);
}
double dis(point p1,point p2)  //计算 p1p2的 距离
{
    return sqrt((double)(p2.x-p1.x)*(p2.x-p1.x)+(p2.y-p1.y)*(p2.y-p1.y));
}
bool cmp(point p1,point p2) //极角排序函数 ， 角度相同则距离小的在前面
{
    int tmp=cross(list[0],p1,p2);
    if(tmp>0) return true;
    else if(tmp==0&&dis(list[0],p1)<dis(list[0],p2)) return true;
    else return false;
}
void init(int n) //输入，并把  最左下方的点放在 list[0]  。并且进行极角排序
{
    int i,k;
    point p0;
    scanf("%d%d",&list[0].x,&list[0].y);
    p0.x=list[0].x;
    p0.y=list[0].y;
    k=0;
    for(i=1;i<n;i++)
    {
        scanf("%d%d",&list[i].x,&list[i].y);
        if( (p0.y>list[i].y) || ((p0.y==list[i].y)&&(p0.x>list[i].x)) )
        {
            p0.x=list[i].x;
            p0.y=list[i].y;
            k=i;
        }
    }
    list[k]=list[0];
    list[0]=p0;
    sort(list+1,list+n,cmp);
}

void graham(int n)
{
    int i;
    if(n==1) {top=0;stack[0]=0;}
    if(n==2)
    {
        top=1;
        stack[0]=0;
        stack[1]=1;
    }
    if(n>2)
    {
        for(i=0;i<=1;i++) stack[i]=i;
        top=1;

        for(i=2;i<n;i++)
        {
            while(top>0&&cross(list[stack[top-1]],list[stack[top]],list[i])<0) top--;
            top++;
            stack[top]=i;
        }
    }
}
bool judge()
{
  //stack[++top]=stack[0];
  //stack[++top]=stack[1];
  //printf("%d\n",top);
 /* for(int i=0;i<=top;i++)
  {
      printf("%d  %d *",list[stack[i]].x,list[stack[i]].y);
  }
  printf("\n");*/
    for(int i=0;i<top;i++)
    {
       if(i==0)
        {
              if(cross(list[stack[i]],list[stack[i+1]],list[stack[i+2]])!=0)
              {
                  //printf("%d**\n",i);
                 // printf("%d %d %d %d %d %d\n",
///list[stack[i]].x,list[stack[i]].y,list[stack[i+1]].x,list[stack[i+1]].y,list[stack[i+2]].x,list[stack[i+2]].y);
                     return false;

              }

        }

        else
        if(cross(list[stack[i]],list[stack[i+1]],list[stack[i+2]])!=0&&cross(list[stack[i-1]],list[stack[i]],list[stack[i+1]])!=0)
                        return false;
    }
    return true;
}
int main()
{

   int t,n;
   scanf("%d",&t);
   while(t--)
   {
        scanf("%d",&n);
        memset(stack,0,sizeof(stack));
           init(n);
            graham(n);
            if(n<6)
            {
                printf("NO\n");
                continue;
            }
            else
            {
                if(judge())
            printf("YES\n");
                  else
            printf("NO\n");
            }

   }
    return 0;
}

发现一点数据bug
1
7
0 0
1 0
2 0
2 1
2 2
1 2
0 2
YES
按照题意是NO AC代码确实YES

真正AC的代码
1将凸包求全了
2判断凸包最后两个点的时候改进了

#include<stdio.h>
#include<math.h>
#include<algorithm>
#include<iostream>
#include<string.h>
using namespace std;
const int MAXN=1005;
const double PI=acos(-1.0);
int flag;
struct point
{
    int x,y;
};
point list[MAXN];
int stack[MAXN],top;

int cross(point p0,point p1,point p2) //计算叉积  p0p1 X p0p2
{
    return (p1.x-p0.x)*(p2.y-p0.y)-(p1.y-p0.y)*(p2.x-p0.x);
}
double dis(point p1,point p2)  //计算 p1p2的 距离
{
    return sqrt((double)(p2.x-p1.x)*(p2.x-p1.x)+(p2.y-p1.y)*(p2.y-p1.y));
}
bool cmp(point p1,point p2) //极角排序函数 ， 角度相同则距离小的在前面
{
    int tmp=cross(list[0],p1,p2);
    if(tmp>0) return true;
    else if(tmp==0&&dis(list[0],p1)<dis(list[0],p2)) return true;
    else return false;
}
void init(int n) //输入，并把  最左下方的点放在 list[0]  。并且进行极角排序
{
    int i,k;
    point p0;
    scanf("%d%d",&list[0].x,&list[0].y);
    p0.x=list[0].x;
    p0.y=list[0].y;
    k=0;
    for(i=1;i<n;i++)
    {
        scanf("%d%d",&list[i].x,&list[i].y);
        if( (p0.y>list[i].y) || ((p0.y==list[i].y)&&(p0.x>list[i].x)) )
        {
            p0.x=list[i].x;
            p0.y=list[i].y;
            k=i;
        }
    }
    list[k]=list[0];
    list[0]=p0;
    sort(list+1,list+n,cmp);
}

void graham(int n)
{
    int i;
    if(n==1) {top=0;stack[0]=0;}
    if(n==2)
    {
        top=1;
        stack[0]=0;
        stack[1]=1;
    }
    if(n>2)
    {
        for(i=0;i<=1;i++) stack[i]=i;
        top=1;

        for(i=2;i<n;i++)
        {
            while(top>0&&cross(list[stack[top-1]],list[stack[top]],list[i])<0) top--;
            top++;
            stack[top]=i;
        }
    }
   // printf("&&\n");
    for(int i=stack[top]-1;i>0;i--)
    {
        if(cross(list[0],list[stack[top]],list[i])==0)
         {

                        stack[++top]=i;
        }
    }
    //printf("%d\n",top);
}
bool judge()
{
  stack[++top]=stack[0];
  stack[++top]=stack[1];
  /*printf("%d\n",top);
  for(int i=0;i<=top;i++)
  {
      printf("%d  %d *",list[stack[i]].x,list[stack[i]].y);
  }
  printf("\n");*/
    for(int i=0;i<top-1;i++)
    {
       if(i==0)
        {
              if(cross(list[stack[i]],list[stack[i+1]],list[stack[i+2]])!=0)
                     return false;
        }
        else
        if(cross(list[stack[i]],list[stack[i+1]],list[stack[i+2]])!=0&&cross(list[stack[i-1]],list[stack[i]],list[stack[i+1]])!=0)
                        return false;
    }
    return true;
}
int main()
{

   int t,n;
   scanf("%d",&t);
   while(t--)
   {
        scanf("%d",&n);
        memset(stack,0,sizeof(stack));
           init(n);
            graham(n);
            if(n<6)
            {
                printf("NO\n");
                continue;
            }
            else
            {
                if(judge())
            printf("YES\n");
                  else
            printf("NO\n");
            }

   }
    return 0;
}
become__better
关注
专栏目录
Grandpa国内开源引擎
06-18
"Grandpa国内开源引擎"是一款专为游戏编程初学者设计的开源游戏开发平台。这款引擎以其易用性和强大的功能，为新手提供了理想的起点，让他们能够快速地进入游戏开发的世界。在"Grandpa"这个亲切的名字背后，蕴含的是...
poj1228——稳定凸包
mumei的博客
09-15 220
题目链接poj.org/problem?id=1228 作为祖父唯一活着的后代，信徒卡姆兰继承了祖父的所有财产。最有价值的是爷爷出生村的一块凸面形农场。农场最初被一根厚厚的绳子与邻近的农场隔开，绳子钩在多边形边界上的一些尖刺（大钉子）。但是，当卡姆兰去参观他的农场时，他注意到绳子和一些尖峰不见了。你的任务是编写一个程序来帮助Kamran决定他的农场的边界是否只能由剩余的峰值决定。‎ ‎输入‎ ...
POJ 1228 （稳定凸包问题)
K键盘里的青春K
04-15 917
转自：http://www.cnblogs.com/xdruid/archive/2012/06/20/2555536.html    这道题算是很好的一道凸包的题吧，做完后会加深对凸包的理解。    题意很关键。。。这英语看了好几遍才差不多看明白了。意思就是给你一堆点，这堆点本来就是某个凸包上的部分点，问你这堆点是否能确定唯一的凸包（大概这意思吧。。。）。后来搜了一下，发现这种凸
poj 1228
weixin_30666943的博客
02-27 90
就是给你一堆点，看这些点能否构成一个 稳定的凸包。 凸包每条边上有3个及以上的点就可以了。 1 #include <cstdio> 2 #include <cstring> 3 #include <cmath> 4 #include <iostream> 5 #include <iomanip> 6 #i...
POJ 1228
weixin_30737433的博客
07-31 76
给题意跪了。。。 题目输入一个凸包上的点（没有凸包内部的点，要么是凸包顶点，要么是凸包边上的点），判断这个凸包是否稳定。所谓稳定就是判断能不能在原有凸包上加点， 得到一个更大的凸包，并且这个凸包包含原有凸包上的所有点。 这样，只需判断每条边是否有大于等于三点即可。注意，一条直线的凸包是NO #include <iostream> #include <cstdi...
POJ 1228 (计算几何 凸包)
morejarphone~
11-12 354
题意是给你n个点，问它们构成的凸包能不能扩大了。 板子题，暴力枚举凸包每条边，只要每条边上都有点（不包含端点），那么这个凸包就不能扩大，否则肯定能扩大。特判1个点和2个点的凸包。 #include #include #include #include #include using namespace std; typedef unsigned long long ll; #defin
Grandpa_visionfall-开源
07-22
这是一款专为万圣节游戏比赛而设计的简短点击式冒险游戏。
Bernie Sanders to Grandpa-crx插件
04-06
【标题】"Bernie Sanders to Grandpa-crx插件"是一个专门为英语用户设计的浏览器扩展程序，其主要功能是将网页上出现的“Bernie Sanders”这个名称替换为“Grandpa”。这款插件以轻松幽默的方式，将美国政治人物伯尼...
五年级英语作文爷爷的愿望Grandpa’sHope
09-09
这篇五年级英语作文“爷爷的愿望Grandpa’s Hope”主要讲述了作者与爷爷之间的亲情故事，体现了家庭关系和代际沟通的重要性。在教育孩子的过程中，这些主题可以帮助孩子们理解和尊重长辈的情感需求，同时培养他们的...
POJ 1228 (稳定凸包问题)
weixin_30478619的博客
08-23 113
<题目链接> <转载于 >>> > 首先来了解什么是稳定的凸包。比如有4个点： 这四个点是某个凸包上的部分点，他们连起来后确实还是一个凸包。但是原始的凸包可能不是这样。 比如： 即这四个点构成的凸包不算做“稳定”的。我们发现，当凸包上存在一条边上的点只有端点两个点的时候，这个凸包不是稳定的，因为它可以在这条边外再引入一个点，构成一个新的...
Grandpa's Estate(poj1228稳定的凸包)
leonharetd的专栏
10-09 1029
题意：给你一个凸包的所有的点，求看其是否能组成唯一的凸包 思路：组成唯一的凸包，那么他的每一条边至少有三个点，因为如果有有两个点，那么可以再加一个点使得凸包扩大，所以每条边至少三个点。 这就和以前组成的凸包不太一样的，以前组凸包有可能没有考虑其边界上有点，因为这对面积没有什么影响，但是这题有影响，所以再求凸包的时候不能要边界上的点，等于说是求得了真正的凸包，每条边只有两个端点，然后枚举每条线段
POJ1228--Grandpa's Estate
SCUT_Pein
02-04 514
Description Being the only living descendant of his grandfather, Kamran the Believer inherited all of the grandpa's belongings. The most valuable one was a piece of convex polygon shaped farm in
poj1228
liuguidongliuguidong的专栏
04-30 600
#include #include #include using namespace std; #define esp 1e-8 struct upoint{double x;double y;}Point[1003]; int stack[1003]; int N; double Xmult(upoint a,upoint b,upoint c) { retu
poj1228-Grandpa's Estate 带边上节点的凸包（稳定凸包）问题
wangheng
07-03 379
Grandpa's Estate Time Limit: 1000MS   Memory Limit: 10000K Total Submissions: 13649   Accepted: 3818 Description Being the only living descendant of his grandfather, Kamr
poj1228Grandpa's Estate
fouzhe的博客
06-07 596
题目描述： DescriptionBeing the only living descendant of his grandfather, Kamran the Believer inherited all of the grandpa’s belongings. The most valuable one was a piece of convex polygon shaped farm in
POJ1228-Grandpa's Estate
ecjtu_acm_菜鸟—鱼尾尾
07-24 668
Grandpa's Estate Time Limit: 1000MS   Memory Limit: 10000K Total Submissions: 10651   Accepted: 2891 Description Being the only living descendant of his grandfather, Kamran
poj1228 (判断凸包的完整性)
无奈滴微笑
03-20 522
Grandpa's Estate Time Limit: 1000MS   Memory Limit: 10000K Total Submissions: 11185   Accepted: 3081 Description Being the only living descendant of his grandfather, Kamr
POJ1228——Grandpa's Estate（计算几何，凸包）
say_c_box的博客
07-29 393
Grandpa's Estate Time Limit: 1000MS   Memory Limit: 10000K Total Submissions: 12618   Accepted: 3540 Description Being the only living descendant of his grandfather, Ka
#include<iostream> using namespace std; class Family { protected: string name1; string name2; int age1; int age2; public: Family(string a,string b) { name1 = a; name2 = b; } virtual void Getage() = 0; virtual void working() = 0; }; class Grand:public Family { public: Grand(string a, string b) : name1(a), name2(b) { cout << "Grandpa's name is " << name1 << "," << "Grandma's name is " << name2 << endl; } void Getage(int x, int y) { age1 = x; age2 = y; cout << "Grandpa's age is " << age1 << "," << "Grandma's age is " << age2 << endl; } void working() { cout << "Grandparents' activity is watering the flowers" << endl; } }; int main() { Grand *g = new Grand("tom", "sara"); g->Getage(65, 34); g->working(); return 0; }
最新发布
06-09
这是一段 C++ 代码，定义了一个家庭类 Family 和一个继承自 Family 的祖父母类 Grand。在 Family 类中，有父母的姓名和年龄属性，并且定义了两个纯虚函数 Getage 和 working，表示获取年龄和工作。...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值