【模板】旋转卡壳

link

模板一个。其实感觉说是求凸包的直径倒不如说是求平面内最远点对，毕竟它的输入没保证是个凸包，自己还要再求一遍啊……

旋转卡壳的思想十分优雅易懂，就是先证明平面内最近点对一定是凸包上两点，再根据这个结论，枚举每一条边的同时找出离线段最远的点并更新答案即可。为了尽可能地降低找点的复杂度，对问题进行分析。容易发现距离是单峰的，所以可以用三分查找来做；但事实上没有必要，因为会发现按顺序枚举边的过程中决策点是单调不降的，所以可以直接用高级决策单调性来解决。算法复杂度的瓶颈是凸包部分的排序，故复杂度为 \(O(NlogN)\) 。

有一些小bug。有可能凸包上只有两个点，这就会导致永远都找不到离线段最远的点而死循环，特判即可。另外凸包的边上不应该有点，否则会造成一些错误。

#include<cstdio>
#include<algorithm>
//#define zczc
using namespace std;
const int N=50010;
inline void read(int &wh){
	wh=0;int f=1;char w=getchar();
	while(w<'0'||w>'9'){if(w=='-')f=-1;w=getchar();}
	while(w>='0'&&w<='9'){wh=wh*10+w-'0';w=getchar();}
	wh*=f;return;
}
inline int max(int s1,int s2){
	return s1<s2?s2:s1;
}

struct node{int x,y;}a[N],st[N],s[N];
inline int pow(int wh){return wh*wh;}
inline int dis(node s1,node s2){return pow(s1.x-s2.x)+pow(s1.y-s2.y);}
inline node ad(node s1,node s2){return (node){s2.x-s1.x,s2.y-s1.y}; }
inline int cj(node s1,node s2){return s1.x*s2.y-s2.x*s1.y;} 
inline bool cmp(node s1,node s2){return s1.x==s2.x?s1.y<s2.y:s1.x<s2.x;}
inline bool check(node a,node b,node c,node d){return (b.y-a.y)*(d.x-c.x)<(b.x-a.x)*(d.y-c.y);}
inline int mj(node s1,node s2,node s3){return cj(ad(s1,s2),ad(s1,s3));}
int m,ans,top,cnt;

signed main(){
	
	#ifdef zczc
	freopen("in.txt","r",stdin);
	#endif
	
	read(m);
	for(int i=1;i<=m;i++){read(a[i].x);read(a[i].y);}
	sort(a+1,a+m+1,cmp);
	
	//求凸包
	for(int i=1;i<=m;i++){
		while(top>1&&!check(st[top-1],st[top],st[top],a[i]))top--;st[++top]=a[i];
	}
	for(int i=1;i<top;i++)s[++cnt]=st[i];top=0;
	for(int i=1;i<=m;i++){
		while(top>1&&!check(st[top],a[i],st[top-1],st[top]))top--;st[++top]=a[i];
	}
	for(int i=top;i>1;i--)s[++cnt]=st[i];
	s[++cnt]=s[1];m=cnt;
	
	if(m<=3){
		for(int i=1;i<=m;i++)
			for(int j=1;j<=m;j++)
				ans=max(ans,dis(s[i],s[j]));
		printf("%d\n",ans);
		return 0;
	}
	
	//求直径
	for(int i=1,last=1;i<m;i++){
		while(mj(s[last],s[i],s[i+1])<=mj(s[last%m+1],s[i],s[i+1]))last=last%m+1;
		ans=max(ans,max(dis(s[last],s[i]),dis(s[last],s[i+1])));
	}
	
	printf("%d\n",ans);
	
	return 0;
}