POJ 2187（求凸包内两点之间最长距离）

推荐博客 ：借鉴

/*

分别列举每个凸包的每条边，再列举另一个凸包的点，求得一对对踵点，与当前记录的值相取舍。
遍历到下一条边，从上次便利到的点开始寻找下一对踵点，再取舍......。直至取舍出最长、最短距离。
*/
#include<stdio.h>
#include<math.h>
#include<algorithm>
using namespace std;
const int maxn=50001;
struct node
{
    double  x,y;
}st[maxn],num[maxn];
int n,k,e;
double dis(node p1,node p2)
{
    return (p1.x-p2.x)*(p1.x-p2.x)+(p1.y-p2.y)*(p1.y-p2.y);
}
int multi(node p1,node p2,node p3)//叉积
{
    return (p1.x-p3.x)*(p2.y-p3.y)-(p2.x-p3.x)*(p1.y-p3.y);
}
bool cmp(node p1,node p2)
{
    if(multi(p1,p2,num[0])>0) return true;
    if(multi(p1,p2,num[0])==0&&dis(p1,num[0])<dis(p2,num[0]))return true;
    return false;
}
void Graham()//凸包
{
    e=2;
    sort(num+1,num+n,cmp);
    st[0]=num[0];
    st[1]=num[1];
    st[2]=num[2];
    for(int i=3;i<n;i++)
    {
        while(e>1&&multi(num[i],st[e],st[e-1])>=0)
            e--;
        st[++e]=num[i];
    }
}
double rotating_calipers()
{
    e++;
    int q=1;
    double ans=0;
    st[e]=st[0];
    for(int p=0;p<e;p++)
    {
        while(multi(st[p+1],st[q+1],st[p])>multi(st[p+1],st[q],st[p]))
            q=(q+1)%e;
        ans=max(ans,max(dis(st[p],st[q]),dis(st[p+1],st[q+1])));
    }
    return ans;
}
int main()
{
    scanf("%d",&n);
    k=0;//K表示起始点
    for(int i=0;i<n;i++){
        scanf("%lf %lf",&num[i].x,&num[i].y);
        if(num[i].x<num[k].x||(num[i].x==num[k].x&&num[i].y<num[k].y)) k=i;
    }
    swap(num[0],num[k]);
    //printf("%d\n",k);
    Graham();
    printf("%.0lf\n",rotating_calipers());
    return 0;
}