nyist 78 圈水池

http://acm.nyist.net/JudgeOnline/problem.php?pid=78

圈水池

时间限制：3000 ms  |  内存限制：65535 KB

难度：4

 

描述

有一个牧场，牧场上有很多个供水装置，现在牧场的主人想要用篱笆把这些供水装置圈起来，以防止不是自己的牲畜来喝水，各个水池都标有各自的坐标，现在要你写一个程序利用最短的篱笆将这些供水装置圈起来！（篱笆足够多，并且长度可变）

 

输入

第一行输入的是N,代表用N组测试数据（1<=N<=10）
第二行输入的是m,代表本组测试数据共有m个供水装置（3<=m<=100）
接下来m行代表的是各个供水装置的横纵坐标

输出

输出各个篱笆经过各个供水装置的坐标点，并且按照x轴坐标值从小到大输出，如果x轴坐标值相同，再安照y轴坐标值从小到大输出

样例输入

1
4
0 0
1 1
2 3
3 0

样例输出

0 0
2 3
3 0

 

 

分析：

凸包。

 

 

AC代码：

 

// NYOJ 78 -- 圈水池 Andrew 凸包算法 O(nlogn)
// 将凸包的顶点按 从小到大 先x后y的 顺序输出
//
/*test data
9 100
200 400
300 400
300 300
400 300
400 400
500 400
500 200
350 200
200 200
    == 1628

*/

#include <cstdio>
#include <cmath>
#include <algorithm>
using namespace std;

const int MAXP = 1010;
const double Pzero = 0.00001; 

struct POINT{
    int x,y;
}p[MAXP],CHp[MAXP*2];

double XJ(POINT a, POINT b, POINT c){
    // vector a->b (b.x-a.x, b.y-a.y)
    // vector a->c (c.x-a.x, c.y-a.y)
    return (b.x-a.x)*(c.y-a.y) - (c.x-a.x)*(b.y-a.y);
}

bool comp(POINT a,POINT b){
    return (a.x < b.x) || ((a.x==b.x) && (a.y<b.y));
}

int AndrewConvexHull(POINT *CHp, int n){
    // 计算凸包 CHp[] 逆时针存点 
    sort(p, p + n, comp);
    int nCHp = 0;
    for (int i = 0 ; i < n ; CHp[nCHp++] = p[i++])
        while(nCHp >= 2 && XJ(CHp[nCHp-2], CHp[nCHp-1], p[i]) <= 0)
            nCHp--;
    for (int i = n - 1,m = nCHp + 1 ; i >= 0 ; CHp[nCHp++] = p[i--])
        while(nCHp >= m && XJ(CHp[nCHp-2], CHp[nCHp-1], p[i]) <= 0)
            nCHp--;
    return --nCHp;
}

void init(int n){
    // 输入
    for (int i = 0 ; i < n; i++)
        scanf("%d%d", &p[i].x, &p[i].y);
}

void Output(int nCHp){
    sort(CHp, CHp + nCHp, comp);
    for (int i = 0; i < nCHp ; i++){
        printf("%d %d\n", CHp[i].x, CHp[i].y);
    }
}

int main()
{
//    freopen("in.txt","r",stdin);
    int n;scanf("%d",&n);
    while(~scanf("%d",&n) ){

        init(n);
        int nCHp = AndrewConvexHull(CHp, n);
        Output(nCHp);
    }
    return 0;
}