POJ 2007 叉积排序

题意：

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

题解：

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

View Code

 1 #include <iostream>
 2 #include <cstdio>
 3 #include <cstdlib>
 4 #include <algorithm>
 5 #include <cstring>
 6 #include <cmath>
 7 
 8 #define N 2222222
 9 #define EPS 1e-3 
10 //不能直接极角排序，因为图形可能全部在原点左侧 
11 using namespace std;
12 
13 struct PO
14 {
15     double x,y,angle;
16 }p[N];
17 
18 int n=1;
19 
20 inline double cross(PO &a,PO &b,PO &c)
21 {
22     return (b.x-a.x)*(c.y-a.y)-(b.y-a.y)*(c.x-a.x);
23 }
24 
25 inline int doublecmp(double x)
26 {
27     if(x>EPS) return 1;
28     else if(x<-EPS) return -1;
29     return 0;
30 }
31 
32 inline bool cmp(PO a,PO b)
33 {
34     return doublecmp(cross(p[1],a,b))>=0;
35 }
36 
37 inline void go()
38 {
39     for(int i=2;i<=n;i++)
40         if(doublecmp(p[i].x)==0&&doublecmp(p[i].y)==0)
41         {
42             swap(p[i].x,p[1].x); swap(p[i].y,p[1].y);
43             break;
44         }
45     sort(p+2,p+1+n,cmp);
46     for(int i=1;i<=n;i++) printf("(%.0lf,%.0lf)\n",p[i].x,p[i].y);
47 }
48 
49 int main()
50 {
51     while(scanf("%lf%lf",&p[n].x,&p[n].y)!=EOF) n++;
52     n--; go();
53     return 0;
54 }