POJ 2007 叉积排序

题意：

给出凸包上的点，其中一个点是(0,0)，要求从(0,0)按照逆时针输出所有点。

题解：

不能有极角排序，因为多边形可能在原点的左侧。

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <algorithm>
#include <cstring>
#include <cmath>

#define N 2222222
#define EPS 1e-3 
//不能直接极角排序，因为图形可能全部在原点左侧 
using namespace std;

struct PO
{
    double x,y,angle;
}p[N];

int n=1;

inline double cross(PO &a,PO &b,PO &c)
{
    return (b.x-a.x)*(c.y-a.y)-(b.y-a.y)*(c.x-a.x);
}

inline int doublecmp(double x)
{
    if(x>EPS) return 1;
    else if(x<-EPS) return -1;
    return 0;
}

inline bool cmp(PO a,PO b)
{
    return doublecmp(cross(p[1],a,b))>=0;
}

inline void go()
{
    for(int i=2;i<=n;i++)
        if(doublecmp(p[i].x)==0&&doublecmp(p[i].y)==0)
        {
            swap(p[i].x,p[1].x); swap(p[i].y,p[1].y);
            break;
        }
    sort(p+2,p+1+n,cmp);
    for(int i=1;i<=n;i++) printf("(%.0lf,%.0lf)\n",p[i].x,p[i].y);
}

int main()
{
    while(scanf("%lf%lf",&p[n].x,&p[n].y)!=EOF) n++;
    n--; go();
    return 0;
}