poj1228

本文介绍了一个基于凸包算法的实现,该算法用于处理一系列二维坐标点,并确定这些点是否能够构成一个凸多边形。通过定义关键函数如向量叉乘、距离计算等,文章详细解释了如何筛选出构成凸包的点集并检查是否存在共线点。此外,还提供了一个完整的C++程序示例。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

//赚点积分下资源
#include <iostream>
#include <algorithm>
#include <cmath>
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)
{
	return (a.x-c.x)*(b.y-c.y)-(b.x-c.x)*(a.y-c.y);
}

double dis(upoint a,upoint b)
{
	return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}


bool dayu(double a,double b)
{
	return a-esp>b;
}
bool cmp(upoint a,upoint b)
{
	return a.x<b.x||(a.x==b.x&&a.y<b.y);
}
int tubao()
{
	sort(Point,Point+N,cmp);
	int top=1;
	for(int i=0;i<3;i++)stack[i]=i;

	for(int i=2;i<N;i++)
	{
		while(top)
		{
			double x=Xmult(Point[i],Point[stack[top]],Point[stack[top-1]]);
			double d1=dis(Point[i],Point[stack[top-1]]);
			double d2=dis(Point[stack[top]],Point[stack[top-1]]);
			if(x>0)top--;
			else if(x==0&&dayu(d1,d2))top--;
			else break;
		}
		stack[++top]=i;
	}

	int t=top;
	stack[++top]=N-2;
	for(int i=N-3;i>=0;i--)
	{
		while(top>t)
		{
			double x=Xmult(Point[i],Point[stack[top]],Point[stack[top-1]]);
			double d1=dis(Point[i],Point[stack[top-1]]);
			double d2=dis(Point[stack[top]],Point[stack[top-1]]);
			if(x>0)top--;
			else if(x==0&&dayu(d1,d2))top--;
			else break;
		}
		stack[++top]=i;
	}
	return top;
}
bool manzu(int i,int j,int k)
{
	return 0==(Point[i].x-Point[k].x)*(Point[j].y-Point[k].y)-(Point[i].y-Point[k].y)*(Point[j].x-Point[k].x);
}
bool nothave(int i,int j)
{
	for(int k=0;k<N;k++)
		if(k!=i&&k!=j)
		{
			if(manzu(i,j,k))return false;
		}
	return true;
}
int main()
{

	int kase;
	cin>>kase;
	while(kase--)
	
	{
		cin>>N;
		for(int i=0;i<N;i++)
		{
			int a,b;
			cin>>a>>b;
			Point[i].x=1.0*a;
			Point[i].y=1.0*b;
		}
		if(N<=5)
		{
			cout<<"NO\n";
			continue;
		}
		int top=tubao();
		if(top<3)
		{
			cout<<"NO\n";
			continue;
		}
		int yes=1;
		for(int i=1;i<top;i++)
		{
			int j=i-1;
			if(nothave(stack[i],stack[j])){yes=0;break;}
		}
		if(nothave(stack[top-1],stack[0]))yes=0;
		if(yes)cout<<"YES\n";
		else cout<<"NO\n";
	}

return 0;
}
//0 1 0 2 0 3 0 4 0 5 NO
//N<=5 NO
//TOP<3 NO~
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值