poj2187-Beauty Contest(极角排序+Graham Scan)

本题为凸包入门题

一、题意

在输入的一堆散点（坐标为整数）中，找出最远点对距离的平方。即求凸包直径。

二、Something

1.用到的知识：极角排序（叉积）、Graham Scan算法

2.用到的数据结构：栈（用vector实现）

三、说明

1.可以解决共线和点数为2的情况。

2.叉积cross(a, b, c),是以向量(c, a)为开始边，逆时针走到向量(b, c)的角的sin值（正数->角<180度，负数->角>180度，0->点共线）。

四、代码实现

1.极角排序

2.Graham Scan算法构造凸包

3.暴力求凸包上最长点对距离的平方

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cmath>
#include <vector>
using namespace std;
const int MAXN = 50000 + 5;

struct P{ int x, y; };
P p[MAXN];
int n;

int dis(P a, P b)
{
    return (a.x - b.x) * (a.x - b.x) + (a.y - b.y) * (a.y - b.y);
}

int cross(P a, P b, P c)
{
    return (a.x - c.x) * (b.y - c.y) - (b.x - c.x) * (a.y - c.y);
}

bool cmp(P a, P b)
{
    int x = cross(a, b, p[0]);
    if(x > 0 || (x == 0 && dis(a, p[0]) < dis(b, p[0]))) return 1;
    return 0;
}

vector<P> Graham()
{
    vector<P> stk(n);//存凸包
    int k = 0;
    for (int i = 0; i < n; i++) {
        while (k > 1 && (cross(p[i], stk[k-2], stk[k-1])) <= 0)//
            k--;
        stk[k++] = p[i];
    }
    stk.resize(k);//调整正确的凸包顶点数
    return stk;
}

int main()
{
    while (~scanf("%d", &n)) {
        int ans = 0, rp = 0;//rp用于临时保存最左下点的下标
        //输入
        for (int i = 0; i < n; i++) {
            scanf("%d%d", &p[i].x, &p[i].y);
            if(p[rp].y > p[i].y
               || (p[rp].y == p[i].y && p[rp].x > p[i].x))//找最左下的点
               rp = i;
        }

        //极角排序
        swap(p[rp], p[0]);//初始化，方便cmp过程
        sort(p + 1, p + n, cmp);

        //Graham Scan构成凸包
        vector<P> stk = Graham();

        //暴力求解
        for (int i = 0; i < stk.size(); i++) {
            for (int j = 0; j < i; j++)
                ans = max(ans, dis(stk[i], stk[j]));
        }

        cout<<ans<<endl;
    }

    return 0;
}

提供一些测试数据：

8
0 5
1 8
3 4
5 0
6 2
6 6
8 3
8 7

80

5
1 0
0 1
1 3
1 2
1 1

9

2
1 0
0 1

2

五、相关题目

1.poj-wall(求凸包+几何计算) http://poj.org/problem?id=1113

2.poj-Scrambled Polygon(极角排序) http://poj.org/problem?id=2007