极角排序
2016-05-18 21:22
369 查看
这里讲了几种极角排序方法,http://www.csie.ntnu.edu.tw/~u91029/PointLinePlane2.html,我这里只说下atan2还有象限排序。
因为纯粹叉积是不能排360度的,所以需要改一些东西。
atan2是按向量排序的,这个很好啊,可是就是可能会伤精度,那么会死很惨的,虽然我还没遇到用atan2伤精度的题~
atan2的参数是向量的y坐标,x坐标,排序范围是(-180°, +180°]。
象限排序是不损精度的,因为只用到叉积还有一些比较而已。
先判断两个点是在哪两个象限,然后如果不在同一个象限,直接象限比较,如果在同一个象限,那么用叉积比较即可,边界问题考虑下。
大小判定函数。
const double eps = 1e-6;
bool dy(double x,double y) { return x > y + eps;} // x > y
bool xy(double x,double y) { return x < y - eps;} // x < y
bool dyd(double x,double y) { return x > y - eps;} // x >= y
bool xyd(double x,double y) { return x < y + eps;} // x <= y
bool dd(double x,double y) { return fabs( x - y ) < eps;} // x == y
atan2排序
bool cmp(point& a,point& b)
{
double t1 = atan2(a.y - C.y, a.x - C.x);
double t2 = atan2(b.y - C.y, b.x - C.x);
if( dd(t1, t2) ) return xy(fabs(a.x),fabs(b.x));
return xy(a.t, b.t);
}
象限极角排序代码
int quad(point a)// 判断象限的函数,每个象限包括半个坐标轴
{
if( dy(a.x,0) && xyd(a.y,0) ) return 1;
if( xyd(a.x,0) && dy(a.y,0) ) return 2;
if( xy(a.x,0) && xyd(a.y,0) ) return 3;
if( dyd(a.x,0) && xy(a.y,0) ) return 4;
}
bool cmp(point& a,point& b)
{
point p1 = a,p2 = b;
p1.x -= C.x; p1.y -= C.y;
p2.x -= C.x; p2.y -= C.y;
int l1,l2;
l1 = quad(p1); l2 = quad(p2);
if( l1 == l2 )
{
double c = crossProduct(C,a,b);
return xy(c,0) || dd(c,0.0) && xy(fabs(a.x),fabs(b.x));
}
return l1 < l2;
}
因为纯粹叉积是不能排360度的,所以需要改一些东西。
atan2是按向量排序的,这个很好啊,可是就是可能会伤精度,那么会死很惨的,虽然我还没遇到用atan2伤精度的题~
atan2的参数是向量的y坐标,x坐标,排序范围是(-180°, +180°]。
象限排序是不损精度的,因为只用到叉积还有一些比较而已。
先判断两个点是在哪两个象限,然后如果不在同一个象限,直接象限比较,如果在同一个象限,那么用叉积比较即可,边界问题考虑下。
大小判定函数。
const double eps = 1e-6;
bool dy(double x,double y) { return x > y + eps;} // x > y
bool xy(double x,double y) { return x < y - eps;} // x < y
bool dyd(double x,double y) { return x > y - eps;} // x >= y
bool xyd(double x,double y) { return x < y + eps;} // x <= y
bool dd(double x,double y) { return fabs( x - y ) < eps;} // x == y
atan2排序
bool cmp(point& a,point& b)
{
double t1 = atan2(a.y - C.y, a.x - C.x);
double t2 = atan2(b.y - C.y, b.x - C.x);
if( dd(t1, t2) ) return xy(fabs(a.x),fabs(b.x));
return xy(a.t, b.t);
}
象限极角排序代码
int quad(point a)// 判断象限的函数,每个象限包括半个坐标轴
{
if( dy(a.x,0) && xyd(a.y,0) ) return 1;
if( xyd(a.x,0) && dy(a.y,0) ) return 2;
if( xy(a.x,0) && xyd(a.y,0) ) return 3;
if( dyd(a.x,0) && xy(a.y,0) ) return 4;
}
bool cmp(point& a,point& b)
{
point p1 = a,p2 = b;
p1.x -= C.x; p1.y -= C.y;
p2.x -= C.x; p2.y -= C.y;
int l1,l2;
l1 = quad(p1); l2 = quad(p2);
if( l1 == l2 )
{
double c = crossProduct(C,a,b);
return xy(c,0) || dd(c,0.0) && xy(fabs(a.x),fabs(b.x));
}
return l1 < l2;
}
相关文章推荐
- Java文件流 字节流和字符流
- 深入浅出: Java回调机制(异步)
- 我的Unity3D学习日记-01
- JAVA获取本机IP地址
- 从头认识多线程-2.23 静态内部类的同步机制跟普通类相同
- Try Git 译文
- java中重载和重写
- 矩阵快速幂+map 雷神之路 scu dp练习A题
- POJ1426 建议打表做该题/用数组模拟队列(比STL队列速度要快)
- Yii2.0 自定义日志类
- Libgdx之Label Image
- 面向项目(十一)—— 库的使用
- JAVA-007前面所学总结加深
- 第十三周的学习进度表
- IOI2000 邮局
- iOS蓝牙开发(一)蓝牙相关基础知识
- 渣渣的艰难的找实习经历
- 5月17日 AJAX 之 XML
- ab中文手册
- AngularJS之表单验证