凸多边形定点逆时针排序

#include
#include
#include
#include
#include


using namespace std;
#define maxn 20000
int top = 2;
struct Point
{
    double x,y,len;
}Pt[maxn],Stack[maxn],Point_A;


double Cross( Point a, Point b, Point c)
{
    return (b.x - a.x)*(c.y - a.y) - (b.y - a.y)*(c.x - a.x);
}


double Dis(Point a, Point b)
{
    return sqrt((a.x - b.x) * (a.x - b.x) + (a.y - b.y)*(a.y - b.y));
}


void FindPoint( int n)
{
    int i,tempNumber = 0;
    Point tempPoint ;


    Point_A = Pt[1];


    for(  i = 1; i<=n; i++)
    {
        if(Pt[i].y < Point_A.y || Pt[i].y == Point_A.y && Pt[i].x < Point_A.x)
        {
            tempNumber = i;
            Point_A = Pt[i];
        }
    }


    tempPoint = Pt[1];
    Pt[1] = Point_A;
    Pt[tempNumber] = tempPoint;
}


bool cmp( Point a, Point b)
{
    double k = Cross(Point_A, a, b);
    if( k > 0) return true;//顺时针方向
    if( k < 0) return false;
    a.len = Dis(Point_A,a);
    b.len = Dis(Point_A,b);
    if(a.len < b.len)//共线
       return true;
}


void Graham( int n)
{
    int i;


    Stack[1] = Pt[1];
    Stack[2] = Pt[2];
    cout<     for( i = 3; i<=n; i++)
    {
        while(Cross(Stack[top-1],Stack[top],Pt[i]) < 0 && top > 1) top--;
        Stack[++top] = Pt[i];
    }


}




int main()
{
    int i,Num;
    while(scanf("%d",&Num)!= EOF)
    {
        for( int i = 1; i <= Num; i++)
        {
            scanf("%lf %lf",&Pt[i].x , &Pt[i].y);
        }


        FindPoint(Num);
        sort(Pt + 1 , Pt + Num + 1,cmp);
        Graham(Num);


//        for( int i = 1 ;i <= top; i++)
//         cout<     }
    return 0;
}

你可能感兴趣的:(平面几何)