CF1059D Nature Reserve (精度处理,计算几何,二分)
Posted paulshi
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了CF1059D Nature Reserve (精度处理,计算几何,二分)相关的知识,希望对你有一定的参考价值。
CF1059D Nature Reserve (精度处理,计算几何,二分)
题目链接:CF1059D
首先处理无解情况,如果在 $x$ 轴两侧都有点,则无解。
我们在将所有 $y$ 值都变为正数方便处理
如果圆与 $x$ 轴相切,则该圆的一条半径垂直于 $x$ 轴。
于是我们可以二分半径 $R$
那么圆心的纵坐标是确定的,那么我们如何判断该半径能否覆盖所有圆呢?
如图我们根据勾股定理算出圆心的 $横坐标$ 的的集合
我们可以算出 $X_圆$ $le$ $X+sqrt{R^{2}-(Y-R)^{2}} $ $ps$:上图有点 $bug$
$X_圆$ $ge$ $X-sqrt{R^{2}-(Y-R)^{2}} $
则所有点的并集为空则不成立,反之成立
以上写在 $check$ 函数中
然后就是二分答案
精度问题:
因为题目要求是与正确答案相对相差1e-6即可,所以我们将eps设为期望答案大小的1e-7被即可,防止超时,卡精度
1 #include<bits/stdc++.h> 2 #define MAXN 100010 3 using namespace std; 4 struct Point{ 5 int x,y; 6 }a[MAXN]; 7 int n; 8 bool f1,f2; 9 bool check (double k) 10 { 11 double L=-1e15,R=1e15; 12 for (int i=1;i<=n;i++) 13 { 14 double tmp=a[i].y*(2*k-a[i].y); 15 if (tmp<0) return 0; 16 tmp=sqrt (tmp); 17 L=max (L,a[i].x-tmp);R=min (R,a[i].x+tmp); 18 } 19 return L<=R; 20 } 21 int main() 22 { 23 scanf ("%d",&n); 24 int Min_B=1e9,Max_B=-1e9,High=-1e9; 25 for (int i=1;i<=n;i++) 26 { 27 scanf ("%d%d",&a[i].x,&a[i].y); 28 Min_B=min (Min_B,a[i].x); 29 Max_B=max (Max_B,a[i].x); 30 High=max (High,abs (a[i].y)); 31 if (a[i].y>0) f1=1; 32 if (a[i].y<0) f2=1; 33 } 34 if (f1&&f2) 35 { 36 printf ("-1"); 37 return 0; 38 } 39 if (f2) for (int i=1;i<=n;i++) a[i].y=-a[i].y; 40 double l=0,r=1e14; 41 double eps=max (High,(Max_B-Min_B))*1e-8; 42 while (r-l>eps) 43 { 44 double mid=(l+r)/2; 45 if (check (mid)) r=mid; 46 else l=mid; 47 } 48 printf ("%.7lf",r); 49 return 0; 50 }
以上是关于CF1059D Nature Reserve (精度处理,计算几何,二分)的主要内容,如果未能解决你的问题,请参考以下文章
Codeforces-1059D:Nature Reserve问最大的圆包含全部点
Codeforces Round #514 (Div. 2) D. Nature Reserve
为啥调用vector.reserve(required + 1) 比vector.reserve(required) 快?