http://poj.org/problem?id=2187
求任意两点的最远距离平方,枚举所以的点肯定超时。
最远两点肯定是凸包上的两点,所以先求凸包,再对凸包上的点两两枚举,
求最大平方值
#include<iostream>
#include<cmath>
#include<string>
#include<algorithm>
#include<queue>
#include<cstring>
#include<cstdio>
using namespace std;
const double PI=acos(-1);
const int N=50005;
struct node
{
int x,y;
}mem[N];
int used_point[N];
int dist(int i,int j)//求两点距离平方
{
return (mem[i].x-mem[j].x)*(mem[i].x-mem[j].x)
+(mem[i].y-mem[j].y)*(mem[i].y-mem[j].y);
}
bool cmp(node l1,node l2)//叉积排序 如果相等 近的优先
{
if((l1.x-mem[0].x)*(l2.y-mem[0].y)==(l2.x-mem[0].x)*(l1.y-mem[0].y))
return ((l1.x-mem[0].x)*(l1.x-mem[0].x)+(l1.y-mem[0].y)*(l1.y-mem[0].y))<
((mem[0].x-l2.x)*(mem[0].x-l2.x)+(mem[0].y-l2.y)*(mem[0].y-l2.y));
return (l1.x-mem[0].x)*(l2.y-mem[0].y)>(l2.x-mem[0].x)*(l1.y-mem[0].y);
}
bool Left_turn(int i,int j,int l)//判断是否左转
{
int x1=mem[j].x-mem[i].x;
int y1=mem[j].y-mem[i].y;
int x2=mem[l].x-mem[j].x;
int y2=mem[l].y-mem[j].y;
if(x1*y2>x2*y1)//如果左转 true (这种情况为不要直线上多余的点 如果想要应用 >=)
return true;
return false;
}
int main()
{
int n;
while(scanf("%d",&n)!=EOF)
{
int k=0;
for(int i=0;i<n;++i)
{
scanf("%d %d",&mem[i].x,&mem[i].y);
if(mem[i].y<mem[k].y||(mem[i].y==mem[k].y&&mem[i].x<mem[k].x))
k=i;
}
if(k!=0)
{
swap(mem[0].x,mem[k].x);
swap(mem[0].y,mem[k].y);
}
sort(mem+1,mem+n,cmp);
int I=0;
used_point[I++]=0;
used_point[I++]=1;
used_point[I++]=2;
for(int i=3;i<n;)
{
if(I<2||Left_turn(used_point[I-2],used_point[I-1],i))//是左转则入栈加一 (I<2)的情况是对于开始就有直线又不要直线上多余点的情况,
//防止越界。如果直线上多余的点也要则可不写,但会有多余的点使得效率低
used_point[I++]=i++;
else//否则出栈
--I;
}
used_point[I]=used_point[0];
int ans=0;
for(int i=0;i<I;++i)
{
for(int j=i+1;j<I;++j)
{
ans=max(ans,dist(used_point[i],used_point[j]));//枚举最大距离平方
}
}
printf("%d\n",ans);
}
return 0;
}