POJ 2187 /// 凸包入门 旋转卡壳

题目大意：

求最远点对距离

求凸包上的最远点对

挑战263页

#include <cstdio>
#include <string.h>
#include <algorithm>
#include <vector>
#include <cmath>
using namespace std;
const int N=5e5+;
const double eps=1e-;

int n;
double add(double a,double b) {
    if(abs(a+b)<eps*(abs(a)+abs(b))) return ;
    return a+b;
}
struct P {
    double x,y;
    P(){};
    P(double _x,double _y):x(_x),y(_y){};
    P operator - (P p) {
        return P(add(x,-p.x),add(y,-p.y)); };
    P operator + (P p) {
        return P(add(x,p.x),add(y,p.y)); };
    P operator * (double d) {
        return P(x*d,y*d); };
    double dot(P p) {
        return add(x*p.x,y*p.y); };
    double det(P p) {
        return add(x*p.y,-y*p.x); };
}p[N];

bool cmp(const P &a,const P &b){
    if(a.x==b.x) return a.y<b.y;
    return a.x<b.x;
}
vector <P> andrew()
{
    sort(p,p+n,cmp);
    int k=;
    vector <P> q(n*);
    for(int i=;i<n;i++) {
        while(k> && (q[k-]-q[k-]).det(p[i]-q[k-])<=) k--;
        // >0时q[k-1]才位于线段q[k-2]p[i]的左侧 形成下凸
        q[k++]=p[i];
    } /// 求下凸
    for(int i=n-,t=k;i>=;i--) {
        while(k>t && (q[k-]-q[k-]).det(p[i]-q[k-])<=) k--;
        q[k++]=p[i];
    } /// 求上凸
    q.resize(k-);
    return q;
}

double dist(P a,P b) {
    return (a-b).dot(a-b);
} // 点ab距离的平方

void solve()
{
    vector <P> q=andrew();
    int n=q.size();
    if(n==) {
        printf("%.0f\n",dist(q[],q[]));
        return ;
    }

    int i=, j=;
    for(int k=;k<n;k++) {
        if(!cmp(q[i],q[k])) i=k; // 凸包的最左下
        if(cmp(q[j],q[k])) j=k; // 凸包的最右上
    }

    double ans=;
    int si=i, sj=j;
    while(i!=sj || j!=si) { // 转半圈就能判断到所有对踵点对
        ans=max(ans,dist(q[i],q[j]));
        /*判断旋转
        <0两者未达到平行 那么先转i更易达到平行
        >0两者超过平行 那么转j更易达到平行
        */
        if((q[(i+)%n]-q[i]).det(q[(j+)%n]-q[j])<)
            i=(i+)%n;
        else j=(j+)%n;
    }
    printf("%.0f\n",ans);
}

int main()
{
    scanf("%d",&n);
    for(int i=;i<n;i++)
        scanf("%lf%lf",&p[i].x,&p[i].y);
    solve();

    return ;
}