一.相关定义
在平面内取一个定点O,叫极点,引一条射线Ox,叫做极轴,再选定一个长度单位和角度的正方向(通常取逆时针方向)。对于平面内任何一点M,用ρ表示线段OM的长度(有时也用r表示),θ表示从Ox到OM的角度,ρ叫做点M的极径,θ叫做点M的极角,有序数对 (ρ,θ)就叫点M的极坐标,这样建立的坐标系叫做极坐标系。
在通常情况下,我们需要将这些点按照一个选定的中心顺(逆)时针排序。
二.排序方法
在排序过程中用到结构体和相关函数
struct Point{
double x,y;
Point(){}
Point(double xx,double yy){
x=xx;
y=yy;
}
};
Point arr[2005];
double operator *(Point a,Point b){ //叉乘
return a.x*b.y-a.y*b.x;
}
double compare(Point a,Point b,Point c){ //计算极角
Point A(b.x-a.x,b.y-a.y);
Point B(c.x-a.x,c.y-a.y);
return A*B;
}
方法1:利用atan2()函数极角排序
利用#include<cmath>的atan2(double y,double x)函数排序,此函数返回的是与x轴正方向的
夹角,可以处理四个象限的任意情况,排序结果是3,4,1,2象限
bool cmp(Point a,Point b){
if(atan2(a.y,a.x)!=atan2(b.y,b.x))
return atan2(a.y,a.x)<atan2(b.y,b.x);
return a.x<b.x;
}
方法2:利用叉积按极角排序
bool cmp(Point b,Point c){
Point a(0,0);
if(compare(a,b,c)==0) return b.x<c.x;
return compare(a,b,c)>0;
}
对比:
方法1,利用atan2(double y,double x)排序时间快,但是精度不够
方法2,利用叉积排序,时间慢一点,精度较高
二者,有利有弊,选择合适的即可。