【Google Code Jam 2009 round2 problem D】Watering Plants (两圆交点求法详解)

本文主要是介绍【Google Code Jam 2009 round2 problem D】Watering Plants (两圆交点求法详解),希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

【Google Code Jam 2009 round2 problem D】Watering Plants 

你在温室种植了N株植物。为了给这些植物浇水,你购买了2台自动浇水的机器。每株植物i的圆心为(Xi,Yi),半径为Ri。任意两个圆都不相交或相切。每台机器都可以给某个完全包含于半径为r的圆形区域内的植物浇水。求最小的r使得存在一种方案能给所有植物浇水。

限制条件

1<=Xi<=1000

1<=Yi<=1000

1<=Ri<=100

Small

1<=N<=3

Large

1<=N<=40

输入

2
3
20 10 2
20 20 2
40 10 3
3
20 10 3
30 10 3
40 10 3

输出

Case #1: 7.000000
Case #2: 8.000000


分析:

0)该题在《挑战程序竞赛第二版》【1】3.7.4中有所介绍,给出了求解算法,误差分析和临界情况的示例,本文整理了其中的思路并添加了产生临界情况的详细分析和求两圆交点坐标公式的推导

1)若对半径r存在一种浇水方案,则对任意R>r,对半径R肯定也存在一种浇水方案。因此容易想到可以使用二分搜索的方法在解空间搜索满足条件的最小值。

2)现在的问题在于能否找到一种高效的判别算法C(r)判断对半径r是否存在一种浇水方案。当两个浇水器的半径确定时,可以放置的位置有无数种情况,但是我们只需要判断两浇水器处于临界位置时是否存在浇水方案(即2个半径为r的圆是否刚好能覆盖所有植物)。当两个圆半径最小时,由于植物数N>=3,两个圆覆盖的植物数N1,N2>=1,若一个圆覆盖的植物数为1,则它半径最小时必定正好与它覆盖的植物同心。若一个圆覆盖的植物数Ni>=2,则它一定至少与两个植物相切。因此两个圆的圆心存在两种临界位置。求出所有临界位置的圆的圆心(这样的圆共有O(n^2)个)。判断这些临界位置的圆中是否存在两个圆可以覆盖所有的植物。因此存在复杂度为O(n^2)的判断算法

3)求解和两植物相切的浇水器的圆心

设植物圆1的圆心为O1(X1,Y1),半径为r1,植物圆2的圆心为O2(X2,Y2),半径为r2,浇水器半径为R

由于浇水器与植物1,2都相切,浇水器圆心在以O1为圆心,R1=R-r1为半径的圆和以O2为圆心,R2=R-r2为半径的圆的交点处。因此我们需要求解两圆交点的坐标。

设存在两圆X,Y

圆X 圆心A(x,y),半径R1, 圆Y的圆心B(a,b),半径R2,两圆交点为C、D,CD与AB的交点为E,过C作平行于Y轴直线,过E作平行于X轴直线,两者相较于F,过B作平行于Y轴直线,过A作平行X轴直线,两者相交于G,设E坐标为(x0,y0),FE=x4,CF=y4,如下图【2】所示



令dx=a-x dy=b-y  a=AB=sqrt((a-b)^2+(b-y)^2)

CE^2=AC^2-AE^2=BC^2-BE^2=BC^2-(AB-AE)^2

整理得AE=(R1^2-R2^2+a^2)/(2a)

易得x0=dx*AE/AB

令b=AE/AB=(R1^2-R2^2+a^2)/(2a^2)

则x0=dx*b

同理y0=dy*b

易得三角形CFE与三角形AGB相似,因此

FE/CE=BG/AB

CE^2=AC^2-AE^2=R1^2-(b*a)^2

FE=BG/AB*CE=dy/a*sqrt(R1^2-(b*a)^2)=dy*sqrt(R1^2/a^2-b^2)

令d=sqrt(R1^2/a^2-b^2)

则x4=FE=dy*d

同理y4=CF=dy*d

于是

C点坐标为(x0-x4,y0+y4)

D点坐标为(x0+x4,y0-y4)

Q.E.D 

4)代码可以参考【1】中的实现

参考资料

【1】《挑战程序设计竞赛》作者秋叶拓哉 / 岩田阳一/北川宜稔 

【2】图片来自于《两圆相交求交点算法证明》 

http://wenku.baidu.com/link?url=hlnh_qEKfbTg1tEcLLHeLUdhfL_5px9sEj2aLPSmVSjuqdwNIAaM46iSzWX-L3srkdY1LnNJBsUJCwQyZklh3cPasuA1TLLOsmdL2D6GEkm

这篇关于【Google Code Jam 2009 round2 problem D】Watering Plants (两圆交点求法详解)的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/396823

相关文章

Linux换行符的使用方法详解

《Linux换行符的使用方法详解》本文介绍了Linux中常用的换行符LF及其在文件中的表示,展示了如何使用sed命令替换换行符,并列举了与换行符处理相关的Linux命令,通过代码讲解的非常详细,需要的... 目录简介检测文件中的换行符使用 cat -A 查看换行符使用 od -c 检查字符换行符格式转换将

详解C#如何提取PDF文档中的图片

《详解C#如何提取PDF文档中的图片》提取图片可以将这些图像资源进行单独保存,方便后续在不同的项目中使用,下面我们就来看看如何使用C#通过代码从PDF文档中提取图片吧... 当 PDF 文件中包含有价值的图片,如艺术画作、设计素材、报告图表等,提取图片可以将这些图像资源进行单独保存,方便后续在不同的项目中使

Android中Dialog的使用详解

《Android中Dialog的使用详解》Dialog(对话框)是Android中常用的UI组件,用于临时显示重要信息或获取用户输入,本文给大家介绍Android中Dialog的使用,感兴趣的朋友一起... 目录android中Dialog的使用详解1. 基本Dialog类型1.1 AlertDialog(

C#数据结构之字符串(string)详解

《C#数据结构之字符串(string)详解》:本文主要介绍C#数据结构之字符串(string),具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录转义字符序列字符串的创建字符串的声明null字符串与空字符串重复单字符字符串的构造字符串的属性和常用方法属性常用方法总结摘

Java中StopWatch的使用示例详解

《Java中StopWatch的使用示例详解》stopWatch是org.springframework.util包下的一个工具类,使用它可直观的输出代码执行耗时,以及执行时间百分比,这篇文章主要介绍... 目录stopWatch 是org.springframework.util 包下的一个工具类,使用它

Java进行文件格式校验的方案详解

《Java进行文件格式校验的方案详解》这篇文章主要为大家详细介绍了Java中进行文件格式校验的相关方案,文中的示例代码讲解详细,感兴趣的小伙伴可以跟随小编一起学习一下... 目录一、背景异常现象原因排查用户的无心之过二、解决方案Magandroidic Number判断主流检测库对比Tika的使用区分zip

Java实现时间与字符串互相转换详解

《Java实现时间与字符串互相转换详解》这篇文章主要为大家详细介绍了Java中实现时间与字符串互相转换的相关方法,文中的示例代码讲解详细,感兴趣的小伙伴可以跟随小编一起学习一下... 目录一、日期格式化为字符串(一)使用预定义格式(二)自定义格式二、字符串解析为日期(一)解析ISO格式字符串(二)解析自定义

springboot security快速使用示例详解

《springbootsecurity快速使用示例详解》:本文主要介绍springbootsecurity快速使用示例,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝... 目录创www.chinasem.cn建spring boot项目生成脚手架配置依赖接口示例代码项目结构启用s

Python中随机休眠技术原理与应用详解

《Python中随机休眠技术原理与应用详解》在编程中,让程序暂停执行特定时间是常见需求,当需要引入不确定性时,随机休眠就成为关键技巧,下面我们就来看看Python中随机休眠技术的具体实现与应用吧... 目录引言一、实现原理与基础方法1.1 核心函数解析1.2 基础实现模板1.3 整数版实现二、典型应用场景2

一文详解SpringBoot响应压缩功能的配置与优化

《一文详解SpringBoot响应压缩功能的配置与优化》SpringBoot的响应压缩功能基于智能协商机制,需同时满足很多条件,本文主要为大家详细介绍了SpringBoot响应压缩功能的配置与优化,需... 目录一、核心工作机制1.1 自动协商触发条件1.2 压缩处理流程二、配置方案详解2.1 基础YAML