Grandpa's Estate - POJ 1228(稳定凸包)

刚开始看这个题目不知道是什么东东，后面看了大神的题解才知道是稳定凸包问题，什么是稳定凸包呢？所谓稳定就是判断能不能在原有凸包上加点，得到一个更大的凸包，并且这个凸包包含原有凸包上的所有点。知道了这个东西就简单了，直接求出来凸包后，然后判断每条边上的点是否超过三点就行了。

代码如下：

============================================================================================================================

#include<stdio.h>
#include<string.h>
#include<algorithm>
#include<math.h>
using namespace std;

const double EPS = 1e-;
const int MAXN = ;
const int oo = 1e9+;

int sta[MAXN], top;

struct point
{
    double x, y;

    point(double x=, double y=):x(x), y(y){}
    point operator - (const point &t)const{
        return point(x-t.x, y-t.y);
    }
    double operator *(const point &t)const{
        return x*t.x + y*t.y;
    }
    double operator ^(const point &t)const{
        return x*t.y - y*t.x;
    }
}p[MAXN];
int Sign(double t)
{
    if(t > EPS)return ;
    if(fabs(t) < EPS)return ;
    return -;
}
double Dist(point a, point b)
{
    return sqrt((a-b)*(a-b));
}
bool cmp(point a, point b)
{
    int t = Sign((a-p[])^(b-p[]));

    if(t == )
        return Dist(a, p[]) < Dist(b, p[]);
    return t > ;
}
void Graham(int N)
{
    int k=;

    for(int i=; i<N; i++)
    {
        if(p[k].y>p[i].y || (Sign(p[k].y-p[i].y)== && p[k].x > p[i].x))
            k = i;
    }
    swap(p[], p[k]);
    sort(p+, p+N, cmp);

    sta[]=, sta[]=, top=;

    if(N < )
    {
        top = N-;
        return ;
    }

    for(int i=; i<N; i++)
    {
        while(top> && Sign((p[i]-p[sta[top]])^(p[sta[top-]]-p[sta[top]])) <= )
            top--;
        sta[++top] = i;
    }
}

int main()
{
    int T, N, i, j;

    scanf("%d", &T);

    while(T--)
    {
        scanf("%d", &N);

        for(i=; i<N; i++)
            scanf("%lf%lf", &p[i].x, &p[i].y);

        Graham(N);
        sta[++top] = sta[];

        for(i=; i<top; i++)
        {
            int s=sta[i], e=sta[i+];
            for(j=; j<N; j++)
            {
                if(j==s || j==e)
                    continue;
                if(Sign( (p[s]-p[e])^(p[j]-p[e]) ) == )
                    break;
            }

            if(j == N)
                break;
        }

        if(i < top || top < )
            printf("NO\n");
        else
            printf("YES\n");
    }

    return ;
}