POJ 2187美人コンテスト

質問の意味:

タイトルリンクが
与えられた(N- \)\点二つ最も遠い点間の距離を見つけ、最も遠い距離の出力の平方\(N <= 50000 \)

アイデア:

立ち往生回転。

注意事項:

配列名を混同しないでください

コード:

#include<cstdio>
#include<algorithm>
using namespace std;
const int N=50005;
int n,top,per[N],res;
struct point{int x,y;int dist(){return x*x+y*y;}}p[N],q[N];
point operator-(point x,point y){return (point){x.x-y.x,x.y-y.y};}
int operator^(point x,point y){return x.x*y.y-x.y*y.x;}
inline int read()
{
    int s=0,w=1; char ch=getchar();
    for(;'0'>ch||ch>'9';ch=getchar())if(ch=='-')w=-1;
    for(;'0'<=ch&&ch<='9';ch=getchar())s=(s<<1)+(s<<3)+(ch^48);
    return s*w;
}
inline bool cmp(int x,int y)
{
    int dit=(p[x]-p[1])^(p[y]-p[1]);
    if(dit)return dit>0;
    return (p[x]-p[1]).dist()<(p[y]-p[1]).dist();
}
inline void Graham()
{
    for(int i=2;i<=n;++i)
        if(p[1].x>p[i].x||(p[1].x==p[i].x&&p[1].y>p[i].y))
            swap(p[1],p[i]);
    for(int i=1;i<=n;++i)per[i]=i;
    sort(per+2,per+n+1,cmp);
    q[++top]=p[1];
    for(int i=2;i<=n;++i)
    {
        int j=per[i];
        while(top>1&&((p[j]-q[top-1])^(q[top]-q[top-1]))>=0)--top;
        q[++top]=p[j];
    }
    q[top+1]=q[1];
}
inline int nxt(int x){return x==top?1:x+1;}
inline int S(point x,point y){return abs(x^y);}
inline int solve()
{
    if(top==2) return (q[1]-q[2]).dist();
    for(int i=1,j=3;i<=top;++i)
    {
        while(i!=nxt(j)&&S(q[i]-q[j],q[i+1]-q[j])<=S(q[i]-q[j+1],q[i+1]-q[j+1]))
            j=nxt(j);
        res=max(res,(q[i]-q[j]).dist());
        res=max(res,(q[i+1]-q[j]).dist());
    }
    return res;
}
int main()
{
    n=read();
    for(int i=1;i<=n;++i)
        scanf("%d%d",&p[i].x,&p[i].y);
    Graham();
    printf("%d\n",solve());
    return 0;
}

おすすめ

転載: www.cnblogs.com/zmyzmy/p/12228303.html