POJ2187

最新推荐文章于 2020-01-19 11:55:15 发布

nkkkkk

最新推荐文章于 2020-01-19 11:55:15 发布

阅读量243

点赞数

分类专栏：题解文章标签：几何

本文链接：https://blog.csdn.net/nkkkkk/article/details/86498674

版权

题解专栏收录该内容

19 篇文章 0 订阅

订阅专栏

求凸包两最远点距离

#include<stdio.h>
#include<math.h>
#include<iostream>
#include<algorithm>
#define pi acos(-1.0)
using namespace std;
const int N=100005;
struct node{int x,y;}P[N],p[N];
int cnt,n;
int X(node a,node b,node c)
{
    return (b.x-a.x)*(c.y-a.y)-(c.x-a.x)*(b.y-a.y);
}
int dis(node a,node b)
{
    return (a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y);
}
bool cmp(node a,node b)
{
    int k=X(p[0],a,b);
    if(k>0) return true;
    if(k<0) return false;
    return dis(p[0],a)<dis(p[0],b);
}
void Graham(int n)
{
    if(n==1)
    {
        cnt=0;
        P[0]=p[0];
    }
    else if(n==2)
    {
        cnt=1;
        P[0]=p[0],P[1]=p[1];
    }
    else
    {
        P[0]=p[0],P[1]=p[1],cnt=1;
        for(int i=2;i<n;i++)
        {
            while(cnt>0&&(X(P[cnt-1],P[cnt],p[i])<=0))
                cnt--;
            cnt++;
            P[cnt]=p[i];
        }
    }
}
void kk(int n)
{
    int q=1;
    int ans=0;
    P[n+1]=P[0];
    for(int i=0;i<cnt;i++)
    {
        while(X(P[i+1],P[q+1],P[i])>X(P[i+1],P[q],P[i]))
            q=(q+1)%n;
        ans=max(ans,max(dis(P[q],P[i]),dis(P[q+1],P[i+1])));
    }
    printf("%d\n",ans);
}
int main()
{
    scanf("%d",&n);
    scanf("%d%d",&p[0].x,&p[0].y);
    for(int i=1;i<n;i++)
    {
        scanf("%d%d",&p[i].x,&p[i].y);
        if(p[i].y<p[0].y)
            swap(p[0],p[i]);
        else if((p[i].y==p[0].y)&&(p[i].x<p[0].x))
        swap(p[0],p[i]);

    }
     sort(p+1,p+n,cmp);
    Graham(n);
    kk(cnt);
    return 0;
}

啊啊啊啊啊啊啊啊啊啊一开始找p[0]的时候第二个比较x写成了y哇了一下午1555555555555555555555555551TAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAT