题目背景

此处省略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 ;
}

最新文章

  1. 使用Monit监控本地进程
  2. 如何在制作jar包时引用第三方jar包
  3. HTML a标签、4个伪类、常用属性(下载)、锚链接(待扩展:邮件、电话、短信、GPS)
  4. Unity3D 第一人称控制器 C#脚本
  5. IIS7下ajax报未定义错误
  6. Python科学计算包模块的安装(ubuntu)
  7. 制作HP MicroServer Gen8可用的ESXi 5.x SD/TF卡启动盘
  8. 浅析SQL Server 2005中的主动式通知机制
  9. Ossec常用命令
  10. 【转】angular通过$http与服务器通信
  11. 将Cell中的视图取出传递到根视图
  12. javascript 判断整数
  13. jquery 中ajax的参数
  14. Linux cat和EOF的使用
  15. C语言可变参数函数的编写
  16. 使用diff和patch指令生成文件差异和还原文件
  17. 06jQuery-01-基本选择器
  18. git 对 Microsoft Word 进行版本控制
  19. Android WebView 缓存机制和模式详解
  20. Tomcat 报错的解决方法:The APR based Apache Tomcat Native library which allows optimal

热门文章

  1. m_Orchestrate learning system---十八、mo项目的启示是什么
  2. 洛谷P2118 比例简化(暴力)
  3. SQL SERVER 提取字符串中汉字
  4. PC比价软件
  5. CDR X6设计师的福利,3折特惠!
  6. 创建一个 Django 项目
  7. iptables 简单介绍及应用 Linux防火墙
  8. MySQL常用技能篇
  9. [洛谷P1835]素数密度
  10. C语言Huffman压缩和解压