题目描述
给定平面上 \(n\) 个点,求凸包直径。
输入格式
第一行一个正整数 \(n\)。
接下来 \(n\) 行,每行两个整数 \(x,y\),表示一个点的坐标。
输出格式
输出一行一个整数,表示答案的平方。
#include<iostream>
#include<cstdio>
#include<cmath>
#include<algorithm>
using namespace std;
int n,top;
double ans;
struct data{
double x,y;
}p[50001],s[50001];
inline double dis(data a,data b){return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));}
inline double mul(data p1,data p2,data p0){
return (p1.x-p0.x)*(p2.y-p0.y)-(p2.x-p0.x)*(p1.y-p0.y);
}
inline bool cmp(data a,data b){
if(mul(a,b,p[0])==0)return dis(a,p[0])<dis(b,p[0]);
return mul(a,b,p[0])>0;
}
void graham(){
top=2; int k=0;
for(int i=1;i<n;i++)
if(p[k].y>p[i].y||(p[k].y==p[i].y&&p[k].x>p[i].x))k=i;
swap(p[0],p[k]);
sort(p+1,p+n,cmp);
s[0]=p[0],s[1]=p[1],s[2]=p[2];
for(int i=3;i<n;i++){
while(top&&mul(p[i],s[top],s[top-1])>=0)top--;
s[++top]=p[i];
}
s[++top]=p[0];
for(int i=0,q=1;i<top;i++){
int j=(i+1)%top;
while(abs(mul(s[i],s[j],s[q]))<abs(mul(s[i],s[j],s[(q+1)%top])))q=(q+1)%top;
ans=max(ans,dis(s[i],s[q]));
}
printf("%.0lf\n",ans*ans);
}
signed main(){
scanf("%d",&n);
for(int i=0;i<n;i++)scanf("%lf%lf",&p[i].x,&p[i].y);
graham();
}