POJ 2187 Beauty Contest（凸包&&平面最远点对）

最新推荐文章于 2020-10-19 17:13:46 发布

原创最新推荐文章于 2020-10-19 17:13:46 发布 · 579 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#计算几何 #凸包 #平面最远点对

计算几何专栏收录该内容

5 篇文章

订阅专栏

本文介绍了一种利用Graham扫描算法寻找二维平面上点集构成的凸包的方法。通过排序和栈操作来高效地确定凸包顶点，最后计算凸包直径。

利用graham算法枚举凸包。

#include<map>
#include<cmath>
#include<queue>
#include<cmath>
#include<cstdio>
#include<stack>
#include<iostream>
#include<cstring>
#include<algorithm>
#define inf 999999999.9
#define eps 1e-8
#define ls l,mid,rt<<1
#define rs mid+1,rt,rt<<1|1
#define LL __int64
using namespace std;
struct P{
    int x,y;
}p[50005],s[50005];
int n,top;
int dis(P a,P b){
    return (a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y);
}

P operator-(P a,P b){
    P t;t.x=a.x-b.x;t.y=a.y-b.y;return t;
}

int operator*(P a,P b){
    return a.x*b.y-a.y*b.x;
}

bool operator<(P a,P b){
    int t=(a-p[1])*(b-p[1]);
    if(t==0)return dis(a,p[1])<dis(b,p[1]);
    return t>0;
}

void Graham(){
    int k = 1;
    for(int i = 2;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[1],p[k]);
    sort(p + 2,p + n + 1);
    s[1]=p[1];
    s[2]=p[2];top = 2;
    for(int i = 3;i <= n; ++ i){
        while(top > 1&&(p[i] - s[top - 1])*(s[top]-s[top-1]) >= 0)
            top--;
        s[++top] = p[i];
    }
}
int RC(){
    int q = 2,ans = 0;
    s[top+1 ]= p[1];
    for(int i = 1;i <= top;++ i){
        while((s[i+1]-s[i])*(s[q+1]-s[i])>(s[i+1]-s[i])*(s[q]-s[i]) )
            q = q%top+1;
        ans = max(ans,max(dis(s[q+1],s[i+1]),dis(s[q],s[i])));
    }
    return ans;
}
int main(){
    int j,i;
    while(~scanf("%d",&n)){
        for(int i = 1;i <= n;++ i){
            scanf("%d%d",&p[i].x,&p[i].y);
        }
        Graham();
        printf("%d",RC());
    }
    return 0;
}