poj2187最远点对(勉强凑数)

Posted ExiledPoet

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了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=100001;
struct P{int x,y;} a[N],s[N];
int n,top=0,ans=0;
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[1])*(y-a[1]))>0) return (x-a[1])*(y-a[1])>0;
    return dis(x-a[1])<dis(y-a[1]);
}
il int S(P a,P b,P c){
    return abs((a-b)*(a-c))/2;
}
int main(){
    scanf("%d",&n);
    for(int i=1;i<=n;i++)
        scanf("%d%d",&a[i].x,&a[i].y);
    for(int i=2;i<=n;i++)
        if(a[i].y<a[1].y||(a[i].y==a[1].y&&a[i].x<a[1].x)) 
            swap(a[1],a[i]);
    sort(a+2,a+n+1,cmp);
    top=2;s[1]=a[1];s[2]=a[2];
    for(int i=3;i<=n;i++){
        while(top>1&&(a[i]-s[top-1])*(s[top]-s[top-1])>=0) top--;
        s[++top]=a[i];
    }
    for(int i=0;i<top;i++)
        s[i]=s[i+1];
    n=top;
    for(int i=0,j=2;i<n;i++){
        while(S(s[i],s[(i+1)%n],s[j])<S(s[i],s[(i+1)%n],s[(j+1)%n])) j=(j+1)%n;
        ans=max(dis(s[i]-s[j]),ans);
    }
    printf("%d",ans);
    return 0;
}

 

以上是关于poj2187最远点对(勉强凑数)的主要内容,如果未能解决你的问题,请参考以下文章

poj 2187旋转卡壳

计算几何-RC-poj2187

POJ 2187 Beauty Contest

POJ2187(凸包+旋转卡壳)

计算几何

POJ2187 Beauty Contest (旋转卡壳算法 求直径)