Graham算法构造凸包
2016-01-31 16:46
218 查看
double eps = 1e-10; struct P { double x, y; P(double x=0, double y=0):x(x), y(y) {} double add(double a, double b){ if(abs(a+b)<eps*(abs(a)+abs(b))) return 0; return a+b; } P operator + (P p){ return P(add(x, p.x), add(y, p.y)); } P operator - (P p){ return P(add(x, -p.x), add(y, -p.y)); } P operator *(double d){ return P(x*d, y*d); } double dot(P p){ //点积 return add(x*p.x, y*p.y); } double det(P p){ //差积 return add(x*p.y, -y*p.x); } }ps[10000 + 100]; double dist(P a, P b) { return sqrt((b-a).dot(b-a)); } bool cmp_x(const P& p, const P& q) { if(p.x!=q.x) return p.x < q.x; return p.y < q.y; } vector<P> convex_hull(P *ps, int n) { sort(ps, ps+n, cmp_x); int k = 0; //凸包顶点数 vector<P> qs(n*2); //构造凸包的下侧 for(int i=0; i<n; i++) { while(k>1 && (qs[k-1]-qs[k-2]).det(ps[i]-qs[k-1])<=0) k--; qs[k++] = ps[i]; } //构造凸包的上侧 for(int i=n-2,t=k; i>=0; i--) { while(k>t && (qs[k-1]-qs[k-2]).det(ps[i]-qs[k-1])<=0) k--; qs[k++] = ps[i]; } qs.resize(k-1); return qs; }
相关文章推荐
- Material Design 相关资源
- LeetCode算法练习题Java
- 高精度减法
- UML图总结
- java BigInteger源码学习
- 初学最短路算法
- 信息学奥林匹克竞赛-小玉买文具
- linux下常用的几个时间函数:time,gettimeofday,clock_gettime,_ftime
- 关于字符数组与string类的比较(长度以及)
- UML——类图
- Linux下套接字详解(五)----基于fork多进程的TCP套接字(阻塞/同步/并发)
- JQuery.Ajax()的data参数类型
- 举例讲解C#编程中委托的实例化使用
- Android 实现ListView 内部的点击事件
- 查找 2
- Bug 是改不完的……
- 高精度乘法
- UML——对象图
- 自己动手开发一个 Web 服务器(三)
- linux中的阻塞机制及等待队列