平面最远点对距离


(求平面集合中最远距离的长度)

//求出的是最远距离的平方;


#define MAXN 100005
struct Point {
    int x, y;
    bool operator < (const Point& _P) const {
        return y<_P.y||(y==_P.y&&x<_P.x);
    };
} pset[MAXN],ch[MAXN];

int cross(Point a,Point b,Point o) {
    return (a.x - o.x) * (b.y - o.y) - (b.x - o.x) * (a.y - o.y);
}
void convex_hull(Point *p,Point *ch,int n,int &len) {
    sort(p, p+n);
    ch[0]=p[0];
    ch[1]=p[1];
    int top=1;
    for(int i=2; i<n; i++) {
        while(top>0&&cross(ch[top],p[i],ch[top-1])<=0)
            top--;
        ch[++top]=p[i];
    }
    int tmp=top;
    for(int i=n-2; i>=0; i--) {
        while(top>tmp&&cross(ch[top],p[i],ch[top-1])<=0)
            top--;
        ch[++top]=p[i];
    }
    len=top;
}

int dist2(Point a,Point b) {
    return (a.x - b.x) * (a.x - b.x) + (a.y - b.y) * (a.y - b.y);
}

int rotating_calipers(Point *ch,int n) {
    int q=1,ans=0;
    ch[n]=ch[0];
    for(int p=0; p<n; p++) {
        while(cross(ch[p+1],ch[q+1],ch[p])>cross(ch[p+1],ch[q],ch[p]))
            q=(q+1)%n;
        ans=max(ans,max(dist2(ch[p],ch[q]),dist2(ch[p+1],ch[q+1])));
    }
    return ans;
}
//pset存放的是点的集合,n个点
int main() {
    int n, len;
    scanf("%d", &n);
     for(int i=0; i < n; i++) {
        scanf("%d%d",&pset[i].x,&pset[i].y);
     }
     convex_hull(pset,ch,n,len);
      printf("%d\n",rotating_calipers(ch,len));
    }
    return 0;
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值