二分+半平面交——poj1279
/*
二分距离,凸包所有边往左平移这个距离,半平面交后看是否还有核存在
*/ #include<iostream>
#include<cstring>
#include<cstdio>
#include<cmath>
#include<vector>
#include<algorithm>
#include<queue>
using namespace std;
#define N 205 typedef double db;
const db eps=1e-;
const db pi=acos(-1.0);
int sign(db k){
if (k>eps) return ; else if (k<-eps) return -; return ;
}
int cmp(db k1,db k2){return sign(k1-k2);}
int inmid(db k1,db k2,db k3){return sign(k1-k3)*sign(k2-k3)<=;}// k3 在 [k1,k2] 内
struct point{
db x,y;
point(){}
point(db x,db y):x(x),y(y){}
point operator + (const point &k1) const{return point(k1.x+x,k1.y+y);}
point operator - (const point &k1) const{return point(x-k1.x,y-k1.y);}
point operator * (db k1) const{return point(x*k1,y*k1);}
point operator / (db k1) const{return point(x/k1,y/k1);}
int getP() const{return sign(y)==||(sign(y)==&&sign(x)>=);}
point turn90(){return point(-y,x);}
db abs(){return sqrt(x*x+y*y);}
point unit(){db w=abs(); return point(x/w,y/w);}
};
int inmid(point k1,point k2,point k3){return inmid(k1.x,k2.x,k3.x)&&inmid(k1.y,k2.y,k3.y);}
db cross(point k1,point k2){return k1.x*k2.y-k1.y*k2.x;}
db dot(point k1,point k2){return k1.x*k2.x+k1.y*k2.y;}
int comp(point k1,point k2){
if(k1.getP()==k2.getP())return sign(cross(k1,k2))>;
return k1.getP()<k2.getP();
} struct line{
point p[];
line(point k1,point k2){p[]=k1; p[]=k2;}
point& operator [] (int k){return p[k];}
int include(point k){return sign(cross(p[]-p[],k-p[]))>=;}//k在l左端
point dir(){return p[]-p[];}
line push(db eps){//向左边平移eps
point delta = (p[]-p[]).turn90().unit()*eps;
return line(p[]+delta,p[]+delta);
}
}; //输入的点是顺时针:ans<0,逆时针:ans>0
bool judge(vector<point> v){
double ans=;
for(int i=;i<v.size()-;i++)
ans+=cross(v[i]-v[],v[i+]-v[]);
return ans>;
}
point getLL(point k1,point k2,point k3,point k4){
db w1=cross(k1-k3,k4-k3),w2=cross(k4-k3,k2-k3); return (k1*w2+k2*w1)/(w1+w2);
}
point getLL(line k1,line k2){return getLL(k1[],k1[],k2[],k2[]);}
int parallel(line k1,line k2){return sign(cross(k1.dir(),k2.dir()))==;}
int sameDir(line k1,line k2){return parallel(k1,k2)&&sign(dot(k1.dir(),k2.dir()))==;}
int checkpos(line k1,line k2,line k3){return k3.include(getLL(k1,k2));}//k1,k2交点在 k3 左端
int operator<(line k1,line k2){//按极角排序,角度相同的从左到右排
if(sameDir(k1,k2))return k2.include(k1[]);
return comp(k1.dir(),k2.dir());
}
vector<line> getHL(vector<line> L){ // 求半平面交 , 半平面是逆时针方向 , 输出按照逆时针
sort(L.begin(),L.end()); deque<line> q;
for (int i=;i<(int)L.size();i++){
if (i&&sameDir(L[i],L[i-])) continue;
while (q.size()>&&!checkpos(q[q.size()-],q[q.size()-],L[i])) q.pop_back();
while (q.size()>&&!checkpos(q[],q[],L[i])) q.pop_front();
q.push_back(L[i]);
}
while (q.size()>&&!checkpos(q[q.size()-],q[q.size()-],q[])) q.pop_back();
while (q.size()>&&!checkpos(q[],q[],q[q.size()-])) q.pop_front();
vector<line>ans; for (int i=;i<q.size();i++) ans.push_back(q[i]);
return ans;
} vector<point>v;
vector<line>L,LL;
int n; int judge(db eps){
LL.clear();
for(int i=;i<L.size();i++)
LL.push_back(L[i].push(eps));
vector<line> res=getHL(LL);
if(res.size()>)return ;
return ;
} int main(){
while(cin>>n && n){
v.clear();L.clear();
for(int i=;i<=n;i++){
point t;cin>>t.x>>t.y;
v.push_back(t);
}
for(int i=;i<v.size();i++)
L.push_back(line(v[i],v[(i+)%v.size()])); db l=,r=1e9,ans,mid;
while(r-l>eps){
mid=(l+r)/;
if(judge(mid))
l=mid,ans=mid;
else r=mid;
} printf("%.6lf\n",ans);
}
}
二分+半平面交——poj1279的更多相关文章
- UVa 1475 (二分+半平面交) Jungle Outpost
题意: 有n个瞭望塔构成一个凸n边形,敌人会炸毁一些瞭望台,剩下的瞭望台构成新的凸包.在凸多边形内部选择一个点作为总部,使得敌人需要炸毁的瞭望塔最多才能使总部暴露出来.输出敌人需要炸毁的数目. 分析: ...
- [HNOI2012][BZOJ2732] 射箭 [二分+半平面交]
题面 BZOJ题面 思路 半平面交代码讲解戳这里,用的就是这道题 我们射箭的函数形如$y=Ax^2+Bx$ 考虑每一个靶子$(x_0,y_1,y_2)$,实际上是关于$A,B$的不等式限制条件 我们只 ...
- poj 3525Most Distant Point from the Sea【二分+半平面交】
相当于多边形内最大圆,二分半径r,然后把每条边内收r,求是否有半平面交(即是否合法) #include<iostream> #include<cstdio> #include& ...
- POJ3525:Most Distant Point from the Sea(二分+半平面交)
pro:给定凸多边形,求凸多边形内的点到最近边界的最远距离. sol:显然是二分一个圆,使得圆和凸多边形不相交,但是这样很难实现. 由于是凸多边形,我们可以把二分圆转化为二分凸多边形的移动. 如果每一 ...
- POJ 3525/UVA 1396 Most Distant Point from the Sea(二分+半平面交)
Description The main land of Japan called Honshu is an island surrounded by the sea. In such an isla ...
- POJ3525-Most Distant Point from the Sea(二分+半平面交)
Most Distant Point from the Sea Time Limit: 5000MS Memory Limit: 65536K Total Submissions: 3955 ...
- poj3525Most Distant Point from the Sea(半平面交)
链接 求凸多边形内一点距离边最远. 做法:二分+半平面交判定. 二分距离,每次让每条边向内推进d,用半平面交判定一下是否有核. 本想自己写一个向内推进..仔细一看发现自己的平面交模板上自带.. #in ...
- 简单几何(半平面交+二分) LA 3890 Most Distant Point from the Sea
题目传送门 题意:凸多边形的小岛在海里,问岛上的点到海最远的距离. 分析:训练指南P279,二分答案,然后整个多边形往内部收缩,如果半平面交非空,那么这些点构成半平面,存在满足的点. /******* ...
- POJ 3525 Most Distant Point from the Sea [半平面交 二分]
Most Distant Point from the Sea Time Limit: 5000MS Memory Limit: 65536K Total Submissions: 5153 ...
随机推荐
- 机器学习算法中的网格搜索GridSearch实现(以k-近邻算法参数寻最优为例)
机器学习算法参数的网格搜索实现: //2019.08.031.scikitlearn库中调用网格搜索的方法为:Grid search,它的搜索方式比较统一简单,其对于算法批判的标准比较复杂,是一种复合 ...
- 87)PHP,PDO的预编译技术
(1) 比如以下的语句: insert into biao1 values(‘李宁’,‘’): insert into biao1 values(‘安踏’,‘’): insert into biao1 ...
- Emergency
题意:有N个点,M条边,每个点有权值,问从起点到终点最短路的个数以及权值最大的最短路的权值. 分析:修改Dijstra模板. #include<bits/stdc++.h> using n ...
- vscode修改样式
以修改上方滚动条宽度为例 打开开发者工具 help->toggle developer tool 或者快捷键 ctrl+shift+i 选择滚动条,找到css对应文件 鼠标移上去可以看到路径,类 ...
- 51nod 1433:0和5
1433 0和5 题目来源: CodeForces 基准时间限制:1 秒 空间限制:131072 KB 分值: 10 难度:2级算法题 收藏 取消关注 小K手中有n张牌,每张牌上有一个一位数的数, ...
- ETC系列产品非接触式读卡器方案:SI522
随着科技的不断发展,出行上高速这是非常寻常的事.但是在很多节假日高峰时期,在高速路口塞车缴费给很多车主造成很大的烦心.为了解决这一系列的问题,科技发明了ETC这种便捷式缴费技术,让车主们顺畅通过高速路 ...
- 2019年第十届蓝桥杯【C++省赛B组】
试题 A: 组队 本题总分:5 分 作为篮球队教练,你需要从以下名单中选出 1 号位至 5 号位各一名球员,组成球队的首发阵容.每位球员担任 1 号位至 5 号位时的评分如下表所示.请你计算首发阵容 ...
- java不要登录发送邮件noReply
最近在工作中遇到一个问题: 客户需要让我们的系统给他们发邮件,但是不给我们提供用户密码. 这样就导致一个问题,即,一般我们使用的都是smtp方式发邮件,发邮件的时候必须登录,虽然可以修改发件人的邮箱名 ...
- postman 请求get post方法的 区别
1.HTTP的五种请求方法:GET, POST ,HEAD,OPTIONS, PUT, DELETE, TRACE 和 CONNECT 方法. GET请求:请求指定的页面信息,并返回实体主体.(通常用 ...
- 01 DDL(DataDefinitionLanguage)
注: 语句用 ; 或 \g \G 表示结束 . 建库语句 : CREATE DATABASE db_name ; 查询有哪些库 : SHO ...