luogu P1452 【模板】旋转卡壳

题目描述

给定平面上 \(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();
    
}