Beauty Contest

http://acm.hust.edu.cn/vjudge/contest/view.action?cid=28417#problem/F

题目大意:给n个点,求相聚最远距离的平方(输出整形)

集体思路:先求出包围所有点的凸包,然后暴力枚举求解(直接暴力会超时)

#include<iostream>
#include<cmath>
#include<string.h>
#include<string>
#include<stdio.h>
#include<algorithm>
#include<iomanip>

using namespace std;
#define eps 0.0000000001
#define PI acos(-1.0)

//*******************************************************************************
//点和向量
struct Point{
    double x,y;
    Point(double x=0,double y=0):x(x),y(y){}
};
typedef Point Vector;
Vector operator-(Vector a,Vector b){return Vector(a.x-b.x,a.y-b.y);}
bool operator<(const Vector& a,const Vector& b){return a.x<b.x||(a.x==b.x && a.y<b.y);}
int dcmp(double x){
    if(fabs(x)<eps)return 0;
    else return x<0 ? -1:1;
}
double Dot(Vector A,Vector B){return A.x*B.x+A.y*B.y;}//向量点积
double Length(Vector A){return sqrt(Dot(A,A));}//向量模长
double Cross(Vector A,Vector B){return A.x*B.y-A.y*B.x;}
//********************************************************************************
//计算凸包输入点数组p,个数n,输出点数组ch,返回凸包定点数
//输入不能有重复,完成后输入点顺序被破坏
//如果不希望凸包的边上有输入点,把两个<=改成<
//精度要求高时,建议用dcmp比较
//基于水平的Andrew算法-->1、点排序2、删除重复的然后把前两个放进凸包
//3、从第三个往后当新点在凸包前进左边时继续,否则一次删除最近加入的点,直到新点在左边
int ConVexHull(Point* p,int n,Point*ch){
    sort(p,p+n);
    int m=0;
    for(int i=0;i<n;i++){//下凸包
        while(m>1 && Cross(ch[m-1]-ch[m-2],p[i]-ch[m-2])<=0)m--;
        ch[m++]=p[i];
    }
    int k=m;
    for(int i=n-2;i>=0;i--){//上凸包
        while(m>k && Cross(ch[m-1]-ch[m-2],p[i]-ch[m-2])<=0)m--;
        ch[m++]=p[i];
    }
    if(n>1)m--;
    return m;
}
//*******************************************************************
Point x[50005];
Point ch[50005];
int main(){
    for(int n;cin>>n&&n;){
        for(int i=0;i<n;i++)cin>>x[i].x>>x[i].y;
        sort(x,x+n);
        int m=ConVexHull(x,n,ch);
        double maxx=-1,anss;
        for(int i=0;i<m;i++){
            for(int j=0;j<m;j++){
                anss=(ch[i].x-ch[j].x)*(ch[i].x-ch[j].x)+(ch[i].y-ch[j].y)*(ch[i].y-ch[j].y);
                if(anss>maxx)maxx=anss;
            }
        }
        cout<<(int)maxx<<'\n';
    }return 0;
}
View Code

 

 

posted @ 2013-08-07 20:03  beautifulzzzz  阅读(414)  评论(0编辑  收藏  举报