P1452 Beauty Contes（旋转卡壳版）

题目背景

此处省略1W字^ ^

题目描述

贝茜在牛的选美比赛中赢得了冠军”牛世界小姐”。因此,贝西会参观N(2 < = N < = 50000)个农场来传播善意。世界将被表示成一个二维平面,每个农场位于一对整数坐标(x,y),各有一个值范围在-10000…10000。没有两个农场共享相同的一对坐标。

尽管贝西沿直线前往下一个农场，但牧场之间的距离可能很大,所以她需要一个手提箱保证在每一段旅程中她有足够吃的食物。她想确定她可能需要旅行的最大可能距离,她要知道她必须带的手提箱的大小。帮助贝西计算农场的最大距离。

输入输出格式

输入格式：

第一行：一个整数n

第2~n+1行：xi yi 表示n个农场中第i个的坐标

输出格式：

一行：最远距离的[b]平方[/b]

输入输出样例

输入样例#1：

4
0 0
0 1
1 1
1 0

输出样例#1：

2

说明

NONE

想了一下还是写写旋转卡壳吧，

毕竟做题功利性不能太强，

但是。。。。

旋转卡壳居然和暴力一样快，，，，，，，，，，，，，，，，，，，，

 #include<iostream>
 #include<cstdio>
 #include<cstring>
 #include<cmath>
 #include<algorithm>
 #define vec point
 using namespace std;
 const double eps=1e-;
 const int MAXN=;
 int n;
 void read(int &n)
 {
     char c='+';int x=;bool flag=;
     while(c<''||c>''){c=getchar();if(c=='-')flag=;}
     while(c>=''&&c<=''){x=x*+(c-);c=getchar();}
     flag==?n=-x:n=x;
 }
 inline int dcmp(double x)
 {
     if(fabs(x)<eps)    return ;
     else return x>?:-;
 }
 struct point
 {
     double x,y;
     inline point(double x=,double y=):x(x),y(y){};
 }p[MAXN],ch[MAXN];
 vec operator - (vec a,vec b) {return vec(a.x-b.x,a.y-b.y);}
 vec operator + (vec a,vec b) {return vec(a.x+b.x,a.y+b.y);}
 vec operator == (vec a,vec b){return (dcmp(a.x-b.x)==&&dcmp(a.y-b.y)==);}
 int comp(const point &a,const point &b)
 {
     if(a.x==b.x)    return a.y<b.y;
     else    return a.x<b.x;
 }
 int stack[MAXN];
 double cross(vec a,vec b){return a.x*b.y-a.y*b.x;}
 int convex_hull()
 {
     sort(p+,p+n+,comp);
     int top=;
     for(int i=;i<=n;i++)
     {
         while(top>&& dcmp(cross(ch[top-]-ch[top-], p[i]-ch[top-]))<=)    top--;
         ch[top++]=p[i];
     }
     int tmp=top+;
     for(int i=n-;i>=;i--)
     {
         while(top+>tmp&& dcmp(cross(ch[top-]-ch[top-], p[i]-ch[top-]))<=)    top--;
         ch[top++]=p[i];
     }
     if(n>) top--;
     return top;
 }
 int dis(point a,point b)
 {
     return    ((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
 }
 int num=;
 int ans=;
 int Cross(point a,point b,point c)
 {
     return (b.x-a.x)*(c.y-a.y)-(c.x-a.x)*(b.y-a.y);
 }
 inline void xzqk()
 {
     if(num==)    ans=dis(ch[],ch[]);
     int j=;
     for(int i=;i<=num;i++)
     {
         while(Cross(ch[i],ch[i+],ch[j])<Cross(ch[i],ch[i+],ch[j+]))
             j=(j+)%num;// 防止溢出
         ans=max(ans,max(dis(ch[i],ch[j]),dis(ch[i+],ch[j])));
     }
 }
 int main()
 {
     //freopen("fc.in","r",stdin);
     //freopen("fc.out","w",stdout);
     read(n);
     //ios::sync_with_stdio(0);
     for(int i=;i<=n;i++)
     {
         double a,b;
         cin>>a>>b;
         p[i]=point(a,b);
     }
     num=convex_hull();
     xzqk();

     printf("%d",ans);
     return ;
 }