题意:平面上n个点,n<=5e4,求两点间距离的最大值
直接枚举Tle,显然距离最大的两点 肯定为 n个点集所构成的凸包的顶点.(凸包的顶点 相对于n减少很多 但最坏情况还是O(n^2))
旋转卡壳:
凸包上A距离C最远,则A距离边BC或者边CD肯定也是最远的,则枚举凸包上每一条边 找到距离该边最远的点
注意到当我们逆时针枚举边的时候,最远点的变化也是逆时针的,这样就可以不用从头计算最远点,而可以紧接着上一次的最远点继续计算
#include <iostream>
#include <algorithm>
#include <cstdio>
#include <cmath>
#include <cstring>
using namespace std;
typedef long long ll;
const int N=5e5+20;
const int M=1<<20;
struct point{
int x,y;
}p[N];
int n;
int s[N],top;
int cross(point p0,point p1,point p2)//p0p1 * p0p2
{
return (p1.x-p0.x)*(p2.y-p0.y)-(p1.y-p0.y)*(p2.x-p0.x);
}
int dis(point p1,point p2)
{
return (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(p[0],p1,p2);
if(tmp>0) return true;
else if(tmp==0&&dis(p[0],p1)<dis(p[0],p2)) return true;//????,?????
else return false;
}
void Graham()
{
for(int i=0;i<=1;i++)
s[i]=i;
top=1;
for(int i=2;i<n;i++)//求出前i个点集构成的凸包
{
while(top>0&&cross(p[s[top-1]],p[s[top]],p[i])<=0)
top--;
s[++top]=i;
}
top++;//
}
void Rotating_Calipers()//旋转卡壳
{
int q=1,ans=0;
ans=dis(p[s[0]],p[s[1]]);
for(int i=0;i<top;i++)//边[i,i+1] ,转到下一条边->凸包逆时针旋转,对应的最远点也跟着旋转
{
while(abs(cross(p[s[(i+1)%top]],p[s[i%top]],p[s[(q+1)%top]])) > abs(cross(p[s[(i+1)%top]],p[s[i%top]],p[s[q%top]])))
q = (q + 1)%top;
ans = max(ans,max(dis(p[s[(i+1)%top]],p[s[q]]),dis(p[s[i%top]],p[s[q]])));
}
printf("%d\n",ans);
}
int main()
{
while(cin>>n)
{
int low=0;
for(int i=0;i<n;i++)
{
scanf("%d%d",&p[i].x,&p[i].y);
if((p[low].y==p[i].y&&p[low].x>p[i].x)||p[low].y>p[i].y)
low=i;
}
swap(p[0],p[low]);
sort(p+1,p+n,cmp);
Graham();
Rotating_Calipers();
}
return 0;
}