求凸包

求凸包代码模板

#include <cstdio>
#include <algorithm>
using namespace std;
int const MAXN = 1005;
struct point
{
    int x, y;
};
point list[MAXN], list2[MAXN];
int stack[MAXN], top;

int cross(point p0, point p1, point p2)//计算叉积 p0p1 X p0p2
{
    return (p1.x - p0.x) * (p2.y - p0.y) - (p2.x - p0.x) * (p1.y - p0.y);
}

double dis(point p1, point p2)//计算p1 p2 的距离
{
    return sqrt((double)(p1.x - p2.x) * (p1.x - p2.x) + (double)(p1.y - p2.y) * (p1.y - p2.y);
}

bool cmp(point p1, point p2)//极角排序，角度相同，则距离小的在前面
{
    int tmp = cross(list[0], p1, p2);
    if(tmp > 0) return true;
    else if(tmp == 0 && dis(list[0], p1) < dis(list[0], p2))   return true;
    else return false;
}
void init(int n)//输入，并把左下方的点放在list[0]。并且进行极角排序
{
    int i, k;
    pount p0;
    scanf("%d %d", &list[0].x, list[0].y);
    p0.x = list[0].x;
    p0.y = list[0].y;
    k = 0;
    for(int i = 1; i < n; i++){
        scanf("%d %d", &list[i].x, &list[i].y);
        if((p0.y > list[0].y) || ((p0.y == list[i].y) && (p0.x > list[i].x)) ){//寻找一个y值最小的点作为起点，如果y一样，选x小的
            p0.x = list[i].x;
            p0.y = list[i].y;
            k = i;
        }
    }
    //list[0]和list[k]进行交换
    list[k] = list[0];
    list[0] = p0;
    sort(list + 1, list + n, cmp);
}

void graham(int n)
{
    int i;
    if(n == 1){top = 0; stack[0] = 0;}
    if(n == 2)
    {
        top = 1; 
        stack[0] = 0;
        stack[1] = 1;
    }
    if(n > 2){
        for(int i = 0; i <= 1; i++) stack[i] = i;
        top = 1;
        for(int i = 2; i < n; i++){
            while(top > 0 && cross(list[stack[top-1]], list[stack[top]], list[i]) <= 0)//满足条件的点不断出栈
                top--;
            top++;
            stack[top] = i;//扫描到的点入栈
        }
    }
}

更多用的Andrew法可以看这篇http://blog.csdn.net/x_y_q_/article/details/52040978

圈水池

时间限制： 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

Andrew法模板题。。

#include <iostream>
#include <algorithm>
using namespace std;
struct point{
    int x, y;
    point(int x = 0, int y = 0):x(x), y(y){}
}p[105], ch[105];
typedef point Vector;
Vector operator -(Vector const &A, Vector const &B){
    return Vector(A.x - B.x, A.y - B.y);
}
int cross(point a, point b){
    return a.x * b.y - a.y * b.x;
}
bool cmp(point a, point b){
    return a.x < b.x || (a.x == b.x && a.y < b.y);
}
int n;

int hull()
{
    int  m = 0;
    sort(p, p + n, cmp);
    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;
}

int main()
{
    int N;
    cin >> N;
    while(N--){
       cin >> n;
       for(int i = 0; i < n; i++)
            cin >> p[i].x >> p[i].y;
       int m = hull();
       sort(ch, ch + m, cmp);
       for(int i = 0; i < m; i++)
            cout << ch[i].x << ' ' << ch[i].y << '\n';
    }
    return 0;
}