计算几何 || 半平面交 (模板收集)
来源:互联网 发布:淘宝优惠券链接微信 编辑:程序博客网 时间:2024/05/16 17:50
#include<cstdio>#include<algorithm>#include<queue>#include<cmath>using namespace std;const double eps = 1e-8;struct Point{ double x,y; Point(double xx=0.0,double yy=0.0):x(xx),y(yy){} Point operator - (const Point &b)const{ return Point(x-b.x,y-b.y); } Point operator +(const Point &b)const{ return Point(x+b.x,y+b.y); } Point operator /(const double &b)const{ return Point(x/b,y/b); } Point operator *(const double &b)const{ return Point(x*b,y*b); } double operator ^(const Point &b)const{ return x*b.y-y*b.x; }};typedef Point myvec;double cross(myvec a,myvec b){ return a^b;}struct Line{ Point p; myvec v; double ang; Line(){} Line( Point pp,myvec vv):p(pp),v(vv){} bool operator < (const Line &l)const{ return ang < l.ang; }};//点p在有向直线L的左边(线上不算)bool on_left( Line l,Point p){ return cross(l.v,p-l.p)>0;}//直线交点 假设交点唯一存在Point get_inter_section(Line a,Line b){ myvec u = a.p - b. p; double t = cross(b.v,u)/cross(a.v,b.v); return a.p+a.v*t;}int half_plane_inter_section(Line *L,int n,Point *poly){ sort(L,L+n);//级角排序 int fir,lst;//双向队列的第一个元素和最后一个元素的下标 Point *p = new Point[n];//p[i] 为q[i]和q[i+1]的交点 Line *q = new Line[n];//双端队列 q[ fir = lst = 0 ] = L[0];//双端队列初始化为只有一个半平面的L[0] for( int i =1; i <n ; ++i) { while( fir < lst && !on_left(L[i],p[lst-1]) ) lst--; while( fir<lst && !on_left(L[i],p[fir]) ) fir++; q[++lst] = L[i]; if( fabs( cross(q[lst].v,q[lst-1].v) ) < eps ){//两向量平行且同向 取内侧一个 lst--; if( on_left(q[lst],L[i].p) ) q[lst] = L[i]; } if( fir < lst ) p[lst-1] = get_inter_section(q[lst-1],q[lst]); } while( fir< lst && !on_left(q[fir],p[lst-1])) lst--;//删除无用的平面 if(lst - fir <=1 ) return 0;//空集 p[lst] = get_inter_section(q[lst],q[fir]);//计算首尾两个半平面的交点 //从 deque 复制到输出中 int m = 0 ; for( int i = fir;i<=lst;++i) poly[m++] = p[i]; return m;}int main(){ return 0;}
0 0
- 计算几何 || 半平面交 (模板收集)
- 【计算几何】半平面交
- 计算几何 || 半平面交
- <模板><计算几何>半平面求交学习小记
- 计算几何之半平面交算法模板及应用
- POJ 3335-Rotating Scoreboard(计算几何-半平面交顺时针模板)
- POJ 1279 Art Gallery (计算几何+半平面交)
- 计算几何_半平面交
- 计算几何学习之半平面交
- 计算几何学习之半平面交
- 计算几何之二分半平面交
- POJ 3335 (计算几何+半平面交)
- bzoj 2618(计算几何+半平面交)
- POJ 3130-How I Mathematician Wonder What You Are!(计算几何-星形-半平面交逆时针模板)
- 平面计算几何模板
- 【半平面交】【计算几何】[BZOJ1038][ZJOI2008]瞭望塔
- BSOJ1156:CQOI2006 凸多边形 计算几何--半平面交
- BSOJ1850:POJ1474摄像头 计算几何 半平面交
- Android Handler机制
- 安装ggplot2
- 学习alert()和confirm()的区别
- 约瑟夫环
- 我的毕业设计-在线五子棋设计
- 计算几何 || 半平面交 (模板收集)
- svm理解
- 【浅谈ASP.NET】——WebForms UnobtrusiveValidationMode 需要“jquery”ScriptResourceMapping
- POJ2082---Terrible Sets(单调栈)
- 4--黑马程序员--技术总结之数组
- 【leetcode】Search for a Range
- 【java】贪吃蛇小游戏
- Mini2440使用USB无线网卡(wifi)
- VS2013如何自己编写C++类