POJ 2187 Beauty Contest G++

//#include <iostream>
#include <vector> 
#include <algorithm>
#include <cstdio>
//挑战程序设计竞赛 例题 没理解 
using namespace std;
struct P{
	double x,y;
	P(){}
	P(double x,double y):x(x),y(y){}
	P operator+ (P p){
		return P(x+p.x,y+p.y);
	}
	P operator- (P p){
		return P(x-p.x,y-p.y);
	}
	P operator* (double d){
		return P(x*d,y*d);
	}
	double dot(P p){
		return x*p.x+y*p.y;
	}
	double det(P p){
		return x*p.y-y*p.x;
	}	
};
//字典序比较 
bool cmp_x(const P&p,const P& q)
{
	if(p.x!=q.x)
	{
		return p.x<q.x;
	}
	return p.y<q.y;	
}
//距离的平方
double dist(P p,P q)
{
	return (p-q).dot(p-q);
} 
//求凸包
vector<P> convex_hull(P*ps,int n)
{
	sort(ps,ps+n,cmp_x);
	int k=0;//凸包的顶点数
	vector<P> qs(n*2);//构造中的凸包
	//构造凸包的下侧
	for(int i=0;i<n;i++)
	{
		while(k>1&&(qs[k-1]-qs[k-2]).det(ps[i]-qs[k-1])<=0)
		k--;
		qs[k++]=ps[i];
	}
	//构造凸包的上侧 
	for(int i=n-2,t=k;i>=0;i--)
	{
		while(k>t&&(qs[k-1]-qs[k-2]).det(ps[i]-qs[k-1])<=0)
		{
			k--;
		}
		qs[k++]=ps[i];
	}
	qs.resize(k-1);
	return qs;
}
 
int main()
{
	int N;
	//cin>>N;
	scanf("%d",&N);
	P ps[N];
	for(int i=0;i<N;i++)
	{
		//cin>>ps[i].x>>ps[i].y;
		double x,y;
		scanf("%lf%lf",&x,&y);
		ps[i].x=x;
		ps[i].y=y;
	}
	vector<P> qs=convex_hull(ps,N);
	double res=0;
	for(int i=0;i<qs.size();i++)
	{
		for(int j=0;j<i;j++)
		{
			res=max(res,dist(qs[i],qs[j]));
		}
	}
	printf("%d\n",(int)res);
	//cout<<(int)res;
	return 0;
} 

N个牧场,已知每个牧场的坐标。求两个牧场最远距离的平方。

发布了1118 篇原创文章 · 获赞 9 · 访问量 33万+

猜你喜欢

转载自blog.csdn.net/woniupengpeng/article/details/78173072