欢迎您访问程序员文章站本站旨在为大家提供分享程序员计算机编程知识!
您现在的位置是: 首页

POJ2187

程序员文章站 2022-03-02 10:20:36
...

求凸包两最远点距离

#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

相关标签: 几何