转于 http://www.cnblogs.com/devtang/archive/2012/02/01/2334977.html
POJ 2007(计算几何初步——极角排序)
先介绍几种极角排序:
1.利用叉积的正负来作cmp.(即是按逆时针排序).此题就是用这种方法
1 bool cmp(const point &a, const point &b)//逆时针排序
2 {
3 point origin;
4 origin.x = origin.y = 0;
5 return cross(origin,b,origin,a) < 0;
6 }
2.利用complex的内建函数。
1 #include<complex>
2 #define x real()
3 #define y imag()
4 #include<algorithm>
5 using namespace std;
6
7 bool cmp(const Point& p1, const Point& p2)
8 {
9 return arg(p1) < arg(p2);
10 }
3.利用arctan计算极角大小。(范围『-180,180』)
1 bool cmp(const Point& p1, const Point& p2)
2 {
3 return atan2(p1.y, p1.x) < atan2(p2.y, p2.x);
4 }
4.利用象限加上极角,叉积。
1 bool cmp(const point &a, const point &b)//先按象限排序,再按极角排序,再按远近排序
2 {
3 if (a.y == 0 && b.y == 0 && a.x*b.x <= 0)return a.x>b.x;
4 if (a.y == 0 && a.x >= 0 && b.y != 0)return true;
5 if (b.y == 0 && b.x >= 0 && a.y != 0)return false;
6 if (b.y*a.y <= 0)return a.y>b.y;
7 point one;
8 one.y = one.x = 0;
9 return cross(one,a,one,b) > 0 || (cross(one,a,one,b) == 0 && a.x < b.x);
10 }
好了,差不多了。
顺便推荐个网站:http://www.csie.ntnu.edu.tw/~u91029/PointLinePlane2.html(上面讲的很详细)
这题的数据应该只有一组:
1 #include<iostream>
2 #include<cmath>
3 #include<complex>
4 #include<algorithm>
5 #define max(a,b) (a)>(b)?(a):(b)
6 #define min(a,b) (a)<(b)?(a):(b)
7 #define EPS 1e-8
8 using namespace std;
9 struct point {
10 double x,y;
11 };
12 point convex[50];
13
14 double cross(const point &p1, const point &p2, const point &q1, const point &q2)
15 {
16 return (q2.y - q1.y)*(p2.x - p1.x) - (q2.x - q1.x)*(p2.y - p1.y);
17 }
18
19 bool cmp(const point &a, const point &b)
20 {
21 point origin;
22 origin.x = origin.y = 0;
23 return cross(origin,b,origin,a) < 0;
24 }
25
26
27 int main()
28 {
29 int cnt = 0;
30 while (scanf("%lf%lf",&convex[cnt].x,&convex[cnt].y) != EOF) {
31 ++cnt;
32 }
33 sort(convex+1,convex+cnt,cmp);
34 for (int i(0); i<cnt; ++i) {
35 cout<<"("<<convex[i].x<<","<<convex[i].y<<")"<<endl;
36 }
37 return 0;
38 }