poj1228稳定凸包

时间：2022-09-14 10:49:23

就是给一系列点，看这是不是一个稳定凸包

稳定凸包是指一个凸包不能通过加点来使它扩大面积，也就是说每条边最少有三个点

判断的地方写错了，写了两边循环，其实数组s已经排好了序，直接每三个判断就好了

#include<map>

#include<set>

#include<cmath>

#include<queue>

#include<stack>

#include<vector>

#include<cstdio>

#include<cassert>

#include<iomanip>

#include<cstdlib>

#include<cstring>

#include<iostream>

#include<algorithm>

#define pi acos(-1.0)

#define ll long long

#define mod 1000000007

#define ls l,m,rt<<1

#define rs m+1,r,rt<<1|1

#pragma comment(linker, "/STACK:1024000000,1024000000")

using namespace std;

const double g=10.0,eps=1e-;

const int N=+,maxn=+,inf=0x3f3f3f;

struct point{

    double x,y;

};

point p[N],s[N];

int n;

double dir(point p1,point p2,point p3)

{

    return (p3.x-p1.x)*(p2.y-p1.y)-(p3.y-p1.y)*(p2.x-p1.x);

}

double dis(point p1,point p2)

{

    return sqrt((p1.x-p2.x)*(p1.x-p2.x)+(p1.y-p2.y)*(p1.y-p2.y));

}

bool comp(point p1,point p2)

{

    double te=dir(p[],p1,p2);

    if(te<)return ;

    if(te==&&dis(p[],p1)<dis(p[],p2))return ;

    return ;

}

bool graham()

{

    int pos;

    double minx,miny;

    minx=miny=inf;

    for(int i=;i<n;i++)

    {

        if(p[i].y<miny||(p[i].y==miny&&p[i].x<minx))

        {

            minx=p[i].x;

            miny=p[i].y;

            pos=i;

        }

    }

    swap(p[],p[pos]);

    sort(p+,p+n,comp);

    int top=;

    p[n]=p[];

    s[]=p[],s[]=p[],s[]=p[];

    for(int i=;i<=n;i++)

    {

        while(top>=&&dir(s[top-],s[top],p[i])>)top--;

        s[++top]=p[i];

    }

  /*  cout<<endl;

    for(int i=0;i<n;i++)

        cout<<p[i].x<<" "<<p[i].y<<endl;

    cout<<endl;

    for(int i=0;i<top;i++)

        cout<<s[i].x<<" "<<s[i].y<<endl;*/

    for(int i=;i<top-;i++)

    {

        if(dir(s[i-],s[i],s[i+])!=&&dir(s[i],s[i+],s[i+])!=)

            return ;

    }

    return ;

}

int main()

{

    ios::sync_with_stdio(false);

    cin.tie();

    int t;

    cin>>t;

    while(t--){

        cin>>n;

        for(int i=;i<n;i++)

            cin>>p[i].x>>p[i].y;

        if(n<)

        {

            cout<<"NO"<<endl;

            continue;

        }

        if(graham())cout<<"YES"<<endl;

        else cout<<"NO"<<endl;

    }

    return ;

}

/*

1

4

0 1

0 2

0 3

0 4

*/

相关文章