凸包,有利于理解凸包极角排序poj1228

题意时判断是不是稳定的凸包,什么叫稳定的凸包,就是求凸包后每条边上要有不少于三个点。如果一条边上只有两个点,则可以加入一个点,扩大这个凸包。

如图,每条边上有至少三个点,要特别考虑基点和最后一个点。

#include<cstdio>
#include<iostream>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;
const int maxn=10005;
const double PI=acos(-1);
struct point
{
	double x;
	double y;
	int id;
};
point p[maxn];
point s[maxn];
double cross(point p0,point p1,point p2)  //叉积.
{
	return (p1.x-p0.x)*(p2.y-p0.y)-(p2.x-p0.x)*(p1.y-p0.y);
}
double dis(point a,point b)
{
	return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}
bool cmp(point p1,point p2)  //按照斜率排序。
{
	if(cross(p1,p2,p[0])>0||(cross(p1,p2,p[0])==0&&(dis(p1,p[0])<dis(p2,p[0]))))
		return true;
    return false;
}
double max(double a,double b)
{
	if(a>b) return a;
    return b;
}
int graham(int n)
{
	int top,i,tmp=0;
	for(i=1;i<n;i++)  //找基点。
	{
		if((p[i].y<p[tmp].y)||((p[i].y==p[tmp].y)&&(p[i].x<p[tmp].x)))
			tmp=i;
	}
	swap(p[0],p[tmp]);
	sort(p+1,p+n,cmp); //按斜率排序。
	for(i=0;i<n;i++) //对排序后的点编号。
	 p[i].id=i;
	top=1;
	s[0]=p[0];
	s[1]=p[1];
	for(i=2;i<n;i++)
	{
		while(top&&cross(s[top-1],s[top],p[i])<=0)
			top--;
		s[++top]=p[i];
	}
	return top;
}
int main()
{
	//freopen("Input.txt","r",stdin);
	int n,i,tt,ncase;
	scanf("%d",&ncase);
	while(ncase--)
	{
	    scanf("%d",&n);
		for(i=0;i<n;i++)
			scanf("%lf%lf",&p[i].x,&p[i].y);
		tt=graham(n);
		int flag=0;
		if(tt<2) flag=1;
		for(i=1;i<=tt;i++)
        {
            if(s[i].id-s[i-1].id>=2) continue;  //看编号。
             flag=1;
            if(flag) break;
        }
        if(cross(p[0],p[n-1],p[n-2])!=0) flag=1; //特别判断部分
        if(flag) puts("NO");
        else puts("YES");
	}
	return 0;
}


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值