BZOJ1038 瞭望塔
学习了半平交面。
我这里写的是训练指南中的双端队列,每次判断是否删去更优然后更新。
看hzwer中有一处不太明白就是为何要将两段加入队列
后来对拍出错才知道是因为精度,当两线重合时他们叉积返回值是一个极小值
所以判断一下精度即可。
#include<bits/stdc++.h>
using namespace std;
const int N=;
int n,m,cnt,tot;double ans=1e20,eps=1e-;
struct point{
double x,y;
}a[N],p[N];
point operator -(point a,point b){
point t;t.x=b.x-a.x;t.y=b.y-a.y;
return t;
}
double operator *(point a,point b){
return a.x*b.y-a.y*b.x;
}
struct line{
point a,b;double slop;
bool operator <(const line &B)const{
return slop==B.slop?(b-a)*(B.b-a)>:slop<B.slop;
}
}l[N],s[N];
point inter(line a,line b){
double k1,k2,t;
k1=(b.b-a.a)*(a.b-a.a);
k2=(a.b-a.a)*(b.a-a.a);
t=k1/(k1+k2);
point ans;
ans.x=b.b.x+(b.a.x-b.b.x)*t;
ans.y=b.b.y+(b.a.y-b.b.y)*t;
return ans;
}
bool jud(line a,line b,line t)
{
point pp=inter(a,b);
return (t.b-t.a)*(pp-t.a)<-eps;
}
void hpl()
{
int L=,R=;
for(int i=;i<=cnt;++i)
{
if(l[i].slop!=l[i-].slop)++tot;
l[tot]=l[i];
}
cnt=tot;tot=;
s[++R]=l[];s[++R]=l[];
for(int i=;i<=cnt;++i)
{
while(L<R&&jud(s[R-],s[R],l[i]))R--;
while(L<R&&jud(s[L+],s[L],l[i]))L++;
s[++R]=l[i];
}
while(L<R&&jud(s[R-],s[R],s[L]))R--; while(L<R&&jud(s[L+],s[L],s[R]))L++;
for(int i=L;i<R;++i)
a[++tot]=inter(s[i],s[i+]);
}
void getans()
{
for(int i=;i<=tot;++i)
{
point t;t.x=a[i].x;t.y=-;
for(int j=;j<n;++j)
if(p[j].x<=a[i].x&&p[j+].x>=a[i].x)
ans=min(ans,a[i].y-inter((line){p[j],p[j+]},(line){t,a[i]}).y);
}
for(int i=;i<=n;++i)
{
point t;t.x=p[i].x;t.y=-;
for(int j=;j<tot;++j)
if(a[j].x<=p[i].x&&a[j+].x>=p[i].x)
ans=min(ans,inter((line){a[j],a[j+]},(line){t,p[i]}).y-p[i].y);
}
}
int main()
{
// freopen("1.out","r",stdin);
// freopen("my.out","w",stdout);
scanf("%d",&n);
for(int i=;i<=n;++i)scanf("%lf",&p[i].x);
for(int i=;i<=n;++i)scanf("%lf",&p[i].y);
// p[0].x=p[1].x;p[0].y=1e9;
// p[n+1].x=p[n].x;p[n+1].y=1e9;
for(int i=;i<n;++i)l[++cnt].a=p[i],l[cnt].b=p[i+];
if(cnt<){
puts("0.000");return ;
}
for(int i=;i<=cnt;++i)
l[i].slop=atan2(l[i].b.y-l[i].a.y,l[i].b.x-l[i].a.x);
sort(l+,l++cnt);
hpl();
getans();
printf("%.3lf",ans);
return ;
}
BZOJ1038 瞭望塔的更多相关文章
- 【bzoj1038】瞭望塔
[bzoj1038]瞭望塔 题意 致力于建设全国示范和谐小村庄的H村村长dadzhi,决定在村中建立一个瞭望塔,以此加强村中的治安.我们将H村抽象为一维的轮廓.如下图所示 我们可以用一条山的上方轮廓折 ...
- 【BZOJ1038】[ZJOI2008]瞭望塔 半平面交
[BZOJ1038][ZJOI2008]瞭望塔 Description 致力于建设全国示范和谐小村庄的H村村长dadzhi,决定在村中建立一个瞭望塔,以此加强村中的治安.我们将H村抽象为一维的轮廓.如 ...
- bzoj1038: [ZJOI2008]瞭望塔
Description 致力于建设全国示范和谐小村庄的H村村长dadzhi,决定在村中建立一个瞭望塔,以此加强村中的治安.我们将H村抽象为一维的轮廓.如下图所示 我们可以用一条山的上方轮廓折线(x1, ...
- bzoj千题计划126:bzoj1038: [ZJOI2008]瞭望塔
http://www.lydsy.com/JudgeOnline/problem.php?id=1038 本题可以使用三分法 将点按横坐标排好序后 对于任意相邻两个点连成的线段,瞭望塔的高度 是单峰函 ...
- 【bzoj1038】瞭望塔 半平面交
题目描述 致力于建设全国示范和谐小村庄的H村村长dadzhi,决定在村中建立一个瞭望塔,以此加强村中的治安.我们将H村抽象为一维的轮廓.如下图所示 我们可以用一条山的上方轮廓折线(x1, y1), ( ...
- 【BZOJ1038】【ZJOI2008】瞭望塔 [模拟退火]
瞭望塔 Time Limit: 10 Sec Memory Limit: 162 MB[Submit][Status][Discuss] Description 致力于建设全国示范和谐小村庄的H村村 ...
- 「BZOJ1038」「洛谷P2600」「ZJOI2008」瞭望塔 半平面交+贪心
题目链接 BZOJ/洛谷 题目描述 致力于建设全国示范和谐小村庄的H村村长dadzhi,决定在村中建立一个瞭望塔,以此加强村中的治安. 我们将H村抽象为一维的轮廓.如下图所示: 我们可以用一条山的上方 ...
- [BZOJ1038][ZJOI2008]瞭望塔(半平面交)
1038: [ZJOI2008]瞭望塔 Time Limit: 10 Sec Memory Limit: 162 MBSubmit: 2999 Solved: 1227[Submit][Statu ...
- 刷题总结——瞭望塔(bzoj1038)
题目: Description 致力于建设全国示范和谐小村庄的H村村长dadzhi,决定在村中建立一个瞭望塔,以此加强村中的治安.我们将H村抽象为一维的轮廓.如下图所示 我们可以用一条山的上方轮廓折线 ...
随机推荐
- Javascript的执行过程详细研究
下面我们以更形象的示例来说明JavaScript代码在页面中的执行顺序.如果说,JavaScript引擎的工作机制比较深奥是因为它属于底层行为,那么JavaScript代码执行顺序就比较形象了,因为我 ...
- 开源中国愚人节网页变模糊的js blur代码
<![if !IE]> <script> /* * by moli */ $(document).ready(function(){ if(document.cookie.in ...
- 在外网使用ssh连接内网中的多台Linux服务器
最近因为要对全球工控机网络进行协议扫描,需要在实验室配置几台服务器,因为我们只有一个IP地址,所以是用路由器搭建了一个内网(拓扑结构如下图).但是这样做了之后无法在宿舍通过ssh直接连接服务器,因为那 ...
- 基本控件文档-UITextField属性---iOS-Apple苹果官方文档翻译
本系列所有开发文档翻译链接地址:iOS7开发-Apple苹果iPhone开发Xcode官方文档翻译PDF下载地址 //转载请注明出处--本文永久链接:http://www.cnblogs.com/Ch ...
- HDU 1214 圆桌会议 (找规律)
题目链接 Problem Description HDU ACM集训队的队员在暑假集训时经常要讨论自己在做题中遇到的问题.每当面临自己解决不了的问题时,他们就会围坐在一张圆形的桌子旁进行交流,经过大家 ...
- 9 - Python函数定义-位置参数-返回值
目录 1 函数介绍 1.1 为什么要使用函数 1.2 Python中的函数 2 函数的基本使用 3 函数的参数 3.1 参数的默认值 3.2 可变参数 3.2.1 可变位置传参 3.2.2 可变关键字 ...
- 148.Sort List---链表排序(冒泡、归并)
题目链接 题目大意:对链表进行排序,要求时间复杂度是o(nlgn). 法一:冒泡,不交换结点,而交换结点中的数值.超时了.代码如下: public ListNode sortList(ListNode ...
- 如何使用curl命令指定ip访问url
有时我们需要测试一个url,但域名并没解析,这时为了一个简单的测试而写host或去做域名解析,显然这并不高效,而有些域名甚至是正式的域名,因此我们可有使用curl命令进行测试 方法一 curl url ...
- linux---centos7 安装chromedriver
1.安装浏览器 指定yum 源 wget -O /etc/yum.repos.d/CentOS-Base.repo http://mirrors.aliyun.com/repo/Centos-7.re ...
- IDEA 部署项目的时候出错:Jar not loaded错误
2011-10-18 17:03:52 org.apache.catalina.loader.WebappClassLoader validateJarFile 信息: validateJarFile ...