POJ2187-最远点对

题目：题目链接

 

题意：题目的意思就是说给你N个点，然后让你求出最远的两个点的距离的平方（注意：是平方）；

 

分析：对于50000个点暴力肯定是不行的，那么我们可以联系到凸包，因为最远的点肯定是在凸包上，那么要求最远的

距离，我们就可以先对这些点求一次凸包（复习凸包的写法）。然后对于凸包上面的点可以枚举，因为这时候要枚举的

点的数目已经减少很多了，这时候，就是普通的暴搞了。主要是前面的凸包的写法。

 

代码：

#include <iostream>
#include <cstdio>
#include <string>
#include <string.h>
#include <map>
#include <vector>
#include <cstdlib>
#include <algorithm>
#include <cmath>
#include <queue>
#include <set>
#include <stack>
#include <functional>
#include <fstream>
#include <sstream>
#include <iomanip>
#include <numeric>
#include <cassert>
#include <bitset>
#include <stack>
#include <ctime>
#include <list>
#define INF 0x7fffffff
#define max3(a,b,c) (max(a,b)>c?max(a,b):c)
#define min3(a,b,c) (min(a,b)<c?min(a,b):c)
#define mem(a,b) memset(a,b,sizeof(a))
using namespace std;
#define maxn 50005
struct node
{
    double x, y;
}num[maxn], res[maxn];//;

bool ral(node p1,node p2, node p3)
{
    return (p2.x -p1.x)*(p3.y - p1.y) > (p3.x - p1.x)*(p2.y - p1.y);
}
double dis(int a, int b)
{
    return (res[a].x - res[b].x)*(res[a].x - res[b].x) + (res[a].y - res[b].y)*(res[a].y - res[b].y);
}
int cmp(node p1,node p2)
{
    return p1.y < p2.y || (p1.y == p2.y && p1.x < p2.x);
}

int main()
{
    int n;
    scanf("%d", &n);
    for(int i = 0; i < n; ++i)
    scanf("%lf%lf", &num[i].x, &num[i].y);
    if(n == 2)
    {
        int ans = ((num[0].x - num[1].x)*(num[0].x - num[1].x) + (num[0].y - num[1].y)*(num[0].y - num[1].y));
        printf("%d\n", ans);
    }
    else
    {
        sort(num, num + n, cmp);
        res[0] = num[0];
        res[1] = num[1];
        int top = 1;

        for(int i = 2; i < n; ++i)
        {
            while(top && !ral(res[top], res[top-1], num[i]))
            top--;
            res[++top] = num[i];
        }

        int len = top;
        res[++top] = num[n-2];
        for(int i = n-3; i >= 0; --i)
        {
            while(top!= len && !ral(res[top], res[top-1], num[i]))
            top--;
            res[++top] = num[i];
        }

        int MAX = -INF;

        for(int i = 0; i < top; ++i)
        {
            for(int j = 0; j < top; ++j)
            {
                int ans = dis(i, j);
                if(ans > MAX)
                MAX = ans;
            }
        }
        printf("%d\n", MAX);
    }
    return 0;
}