4219: 跑得比谁都快

3007: 拯救小云公主

三角剖分的解释可以看这里:http://www.cnblogs.com/Enceladus/p/6706444.html

后一道是前一道的弱化版。

如果考虑一下二分答案再判定,判定的过程无非就是把每个点变成一个以答案为半价的圆,然后如果左上和右下被切开了肯定就不合法,那就变成1001里面那种对偶图,跑一下spfa就好。

再看其实就不需要二分答案了,把各个boss连接起来,二分的答案达到边权一半时这两个圆就相交。整条路径会在达到最大边权的时候完全连通,那么最大边权最小的路径就是答案。

看了题解才意识到这一点,我的破题能力是不是吃枣药丸啊?

所以就三角剖分+最小生成树,至于三角剖分我选择拉模版。

貌似得用long double

#include<cmath>
#include<cstdio>
#include<algorithm>
#define eps 1e-12
#define MN 110001
#define ld long double
using namespace std; int read_p,read_ca;
inline int read(){
read_p=;read_ca=getchar();
while(read_ca<''||read_ca>'') read_ca=getchar();
while(read_ca>=''&&read_ca<='') read_p=read_p*+read_ca-,read_ca=getchar();
return read_p;
}
struct po{ld x,y;};
struct na{po p;int n;}p[MN];
struct e{int x,y,ne,pre;}b[MN<<];
struct bi{int x,y;ld z;}B[MN<<];
bool operator < (bi a,bi b){return a.z<b.z;}
bool operator < (na a,na b){return a.p.x==b.p.x?a.p.y<b.p.y:a.p.x<b.p.x;}
int n,m,num=,l[MN],st[MN],top,NUM=,x,y,_l[MN],_num=,fa[MN];ld _A,_B;
inline void _add(int x,int y){b[++num].x=x;b[num].y=y;b[num].ne=l[x];b[l[x]].pre=num;l[x]=num;}
inline void add(int x,int y){_add(x,y);_add(y,x);}
inline void _del(int i){if(l[b[i].x]==i)l[b[i].x]=b[i].ne;else{b[b[i].pre].ne=b[i].ne;if(b[i].ne)b[b[i].ne].pre=b[i].pre;}}
inline void del(int i){_del(i);_del(i^);}
ld cha(po a,po b,po c){return (b.x-a.x)*(c.y-b.y)-(c.x-b.x)*(b.y-a.y);}
bool _cha(po a,po b,po c){return cha(a,b,c)<=;}
inline bool ju(po a,po b,po c,po d){return (_cha(a,b,c)^_cha(a,b,d))&&(_cha(c,d,a)^_cha(c,d,b));}
ld sqr(ld x){return x*x;}
inline ld min(ld a,ld b){return a<b?a:b;}
bool incir(po a,po b,po c,po d){
ld m=(a.x*(b.y-c.y)+b.x*(c.y-a.y)+c.x*(a.y-b.y))*.;
ld x=(sqr(a.x)*(b.y-c.y)+sqr(b.x)*(c.y-a.y)+sqr(c.x)*(a.y-b.y)-(a.y-b.y)*(b.y-c.y)*(c.y-a.y))/m;
ld y=((a.x-b.x)*(b.x-c.x)*(c.x-a.x)-sqr(a.y)*(b.x-c.x)-sqr(b.y)*(c.x-a.x)-sqr(c.y)*(a.x-b.x))/m;
return sqrt(sqr(x-a.x)+sqr(y-a.y))+eps>sqrt(sqr(x-d.x)+sqr(y-d.y));
}
void work(int L,int R){
if (R-L<=){
for (int i=L;i<R;i++)
for (int j=i+;j<=R;j++) add(i,j);
return;
}
int mid=L+R>>;
work(L,mid);work(mid+,R);
top=;
for (int i=L;i<=R;st[++top]=i++)
while (top>&&cha(p[st[top-]].p,p[st[top]].p,p[i].p)<) top--;
int Ld,Rd,id;
for (int i=;i<top;i++)
if (st[i]<=mid&&st[i+]>mid) Ld=st[i],Rd=st[i+];
add(Ld,Rd);
for(;;){
id=-;
for (int i=l[Ld];i;i=b[i].ne)
if (b[i].y!=Rd)
if (cha(p[b[i].y].p,p[Ld].p,p[Rd].p)>&&(id==-||incir(p[Ld].p,p[Rd].p,p[id].p,p[b[i].y].p))) id=b[i].y;
for (int i=l[Rd];i;i=b[i].ne)
if (b[i].y!=Ld)
if (cha(p[Ld].p,p[Rd].p,p[b[i].y].p)>&&(id==-||incir(p[Ld].p,p[Rd].p,p[id].p,p[b[i].y].p))) id=b[i].y;
if (id==-) break; if (id<=mid){
for (int i=l[Ld];i;i=b[i].ne)
if (b[i].y!=Rd)
if (ju(p[id].p,p[Rd].p,p[b[i].y].p,p[Ld].p)) del(i);
Ld=id;
}else{
for (int i=l[Rd];i;i=b[i].ne)
if (b[i].y!=Ld)
if (ju(p[id].p,p[Ld].p,p[b[i].y].p,p[Rd].p)) del(i);
Rd=id;
}
add(Ld,Rd);
}
}
int gf(int x){return fa[x]==x?x:fa[x]=gf(fa[x]);}
int main(){
//freopen("a.in","r",stdin);
//freopen("a.out","w",stdout);
register int i,j;
n=read();scanf("%Lf%Lf",&_A,&_B);
for (i=;i<n;i++) scanf("%Lf%Lf",&p[i].p.x,&p[i].p.y),p[i].n=i+,fa[i]=i;
sort(p,p+n);
work(,n-);
for (i=;i<n;i++)
for (j=l[i];j;j=b[j].ne) if (b[j].y>i) B[NUM].x=i,B[NUM].y=b[j].y,B[NUM++].z=sqrt(sqr(p[i].p.x-p[b[j].y].p.x)+sqr(p[i].p.y-p[b[j].y].p.y))/;
for (i=;i<n;i++) B[NUM].x=i,B[NUM].y=n,B[NUM++].z=min(p[i].p.x,_B-p[i].p.y);
for (i=;i<n;i++) B[NUM].x=i,B[NUM].y=n+,B[NUM++].z=min(_A-p[i].p.x,p[i].p.y);
n+=;fa[n-]=n-;fa[n-]=n-;
sort(B,B+NUM);
for (i=;i<NUM;i++){
int x=gf(B[i].x),y=gf(B[i].y);
if (x==y) continue;
fa[x]=y;
if (gf(n-)==gf(n-)) return printf("%.6Lf\n",B[i].z),;
}
}

upd 2017/12/26

  代码被hack了T_T,详见http://www.cnblogs.com/Enceladus/p/6706444.html

BZOJ:4219: 跑得比谁都快 3007: 拯救小云公主的更多相关文章

  1. 【BZOJ】【3007】拯救小云公主

    思路题 我的naive的做法是二分答案+判定是否有路径可走……但是没有正确理解[走的方向任意]这句话…… 其实就是说想咋走咋走= =360°无死角乱走…… 所以其实是个平面上的问题…… 我们可以换个方 ...

  2. BZOJ 3007 [SDOI2012]拯救小云公主 - 对偶图 + 并查集

    Solution 答案具有单调性, 显然可以二分答案. 有两个注意点 : 英雄是可以随便走的, 也就是不是网格图... 还有坐标不能小于$1$ QAQ 开始时英雄在左下角, 公主在右上角, 我们反过来 ...

  3. BZOJ 3007 解救小云公主 二分答案+对偶图

    题目大意:给定一个矩形和矩形内的一些点.求一条左下角到右上角的路径.使全部点到这条路径的最小距离最大 最小距离最大.果断二分答案 如今问题转化成了给定矩形中的一些圆形障碍物求左下角和右上角是否连通 然 ...

  4. 51Nod 1781 跑的比谁都快

    香港记者跑的比谁都快是众所周知的常识. 现在,香港记者站在一颗有 n 个点的树的根结点上(即1号点),编号为 i 的点拥有权值 a[i] ,数据保证每个点的编号都小于它任意孩子结点的别号. 我们假定这 ...

  5. 跑的比谁都快 51Nod - 1789

    香港记者跑的比谁都快是众所周知的常识.   现在,香港记者站在一颗有  nn 个点的树的根结点上(即1号点),编号为  ii 的点拥有权值  a[i]a[i] ,数据保证每个点的编号都小于它任意孩子结 ...

  6. [51nod1789] 跑得比谁都快

    题面 题解 设\(f[i]\)为根节点到\(i\)的最小耗时 设\(S\)为\(i\)的祖先集合, 可以得到 \[ f[i] = min(f[j] + (i - j)^p),j \in S \] 对于 ...

  7. 新型序列化类库MessagePack,比JSON更快、更小的格式

    MessagePack is an efficient binary serialization format. It lets you exchange data among multiple la ...

  8. 不想业务被中断?快来解锁华为云RDS for MySQL新特性

    摘要:新特性上线!华为云RDS for MySQL又添新技能,实力保障业务连续性. 本文分享自华为云社区<不想业务被中断?快来解锁华为云RDS for MySQL新特性>,作者:Gauss ...

  9. 有趣 GIF 动图集 - 仿佛每张小动图都诉说了一个小笑话或者小故事

    点这里 来自法国南特(Nantes)的 Guillaume Kurkdjian 目前还是个学生.Kurkdjian 擅长创作一些平面动态图像,这些有趣的小动图仿佛每张都诉说了一个小笑话或者小故事,像个 ...

随机推荐

  1. animate.css+wow.js页面滚动即时显示动画

    1.地址引入 <link href="css/animate.min.css" rel="stylesheet" type="text/css& ...

  2. 【python】字符串

    >>> str1="welcom to China">>> str1[2:4]'lc'>>> str1[7]'t'>&g ...

  3. 【JMeter】获取json响应报文中数组长度

    import com.jayway.jsonpath.JsonPath; import com.jayway.jsonpath.Predicate; import net.minidev.json.J ...

  4. Python3 将txt数据转换成列表,进行排序,筛选

    Python 程序员需要知道的 30 个技巧 首先是数据: 将上边的四个数据分别写在新建的txt文件中 1.将txt数据转为列表 with open('james.txt') as jaf: data ...

  5. PXE+Kickstart 全自动安装部署CentOS7.4

    一.简介 1.什么是PXE PXE(preboot execute environment,预启动执行环境)是由Intel公司开发的最新技术,工作于Client/Server的网络模式,支持工作站通过 ...

  6. WPF单位真的与分辨率无关吗?

    转载自http://www.cnblogs.com/helloj2ee/archive/2009/04/21/1440709.htm WPF从发布之日起,一直将"分辨率无关(resoluti ...

  7. mybatis源码分析(一)

    mybatis源码分析(sqlSessionFactory生成过程) 1. mybatis框架在现在各个IT公司的使用不用多说,这几天看了mybatis的一些源码,赶紧做个笔记. 2. 看源码从一个d ...

  8. 用Go校验下载文件之SHA256

    用GO校验下载文件之SHA256 原来对计算机和网络使用安全这块不够重视,用了N多年盗版的操作系统和办公软件,为了破解使用过各种激活软件,也安装使用过很多别人破解过的软件:网络下载的文件从不校验.慢慢 ...

  9. Qt个人研究进展

    1:纯socket通信实现多线程邮件发送,支持多个收件人和附件,通用任何平台,包括ARM.2:纯串口通信AT命令实现多线程短信收发,支持多个收件人和长短信,通用任何平台,包括ARM.3:纯串口通信PO ...

  10. 小白的Python之路 day4 软件目录结构规范

    软件目录结构规范 为什么要设计好目录结构? "设计项目目录结构",就和"代码编码风格"一样,属于个人风格问题.对于这种风格上的规范,一直都存在两种态度: 一类同 ...