【poj2187】最远点对（勉强凑数）

题目简述

输入n个点，及其坐标，n<=50000，所有坐标都是不超过10000的整数组成，没有重点。

问最远点对间的距离的平方是多少

题解

这是一道旋转卡壳的裸题

我们要求这个多边形的直径，这可怎么办呢

首先，最远点对一定在凸包上，我们考虑这样一个凸包

显然的，卡在两个点上，一定可以转化成卡在一个边和一个点上

更显然的，如果卡在一条边上和一个点上，那么这个点一定离这个边最远

那么，这个点和这条边组成的三角形一定是包括这条边的三角形中最大的

再之，假定点i和点i+1卡到了点j

那么随着i增加，j也增加

根据这个单调性，我们可以计算出每条边对应的点了

那么直径也就等于边的端点到对应的点的距离

注意边界即可！

#include<stdio.h>
#include<stdlib.h>
#include<iostream>
#include<string>
#include<string.h>
#include<algorithm>
#include<math.h>
#include<queue>
#include<map>
#include<vector>
#include<set>
#define il inline
#define re register
using namespace std;
typedef double db;
const int N=;
struct P{int x,y;} a[N],s[N];
int n,top=,ans=;
il P operator-(P a,P b){
    return (P){a.x-b.x,a.y-b.y};
}
il int operator*(P a,P b){
    return a.x*b.y-a.y*b.x;
}
il int dis(P a){
    return a.x*a.x+a.y*a.y;
}
il bool cmp(P x,P y){
    if(abs((x-a[])*(y-a[]))>) return (x-a[])*(y-a[])>;
    return dis(x-a[])<dis(y-a[]);
}
il int S(P a,P b,P c){
    return abs((a-b)*(a-c))/;
}
int main(){
    scanf("%d",&n);
    for(int i=;i<=n;i++)
        scanf("%d%d",&a[i].x,&a[i].y);
    for(int i=;i<=n;i++)
        if(a[i].y<a[].y||(a[i].y==a[].y&&a[i].x<a[].x))
            swap(a[],a[i]);
    sort(a+,a+n+,cmp);
    top=;s[]=a[];s[]=a[];
    for(int i=;i<=n;i++){
        while(top>&&(a[i]-s[top-])*(s[top]-s[top-])>=) top--;
        s[++top]=a[i];
    }
    for(int i=;i<top;i++)
        s[i]=s[i+];
    n=top;
    for(int i=,j=;i<n;i++){
        while(S(s[i],s[(i+)%n],s[j])<S(s[i],s[(i+)%n],s[(j+)%n])) j=(j+)%n;
        ans=max(dis(s[i]-s[j]),ans);
    }
    printf("%d",ans);
    return ;
}