题目链接:【HDU 3264】
给我们n个圆的圆心坐标(x,y)以及圆的半径r,这n个圆彼此相离(不包含,不相交),求一个圆的最小半径,这个圆满足以下两个条件:
1、它的圆心是这n个圆其中一个圆的圆心
2、它与每一个圆相交部分的面积>=那个圆面积的一半
一开始蒙蔽了,圆相交面积的公式求错了,一直wa
相交面积的公式就是两扇形的面积-四边形的面积
枚举每一个圆的圆心,二分半径
#include <iostream>
#include <algorithm>
#include <cstdio>
#include <cstring>
#include <string>
#include <math.h>
using namespace std;
const double pi=acos(-1.0);
const double eps=1e-8;
int t, n;
struct node
{
double x, y, r;
}a[30];
double dis(int i, int j)
{
return sqrt((a[i].x-a[j].x)*(a[i].x-a[j].x)+(a[i].y-a[j].y)*(a[i].y-a[j].y));
}
bool check(double r1, int id)
{
for(int i=0; i<n; i++)
{
double r2=a[i].r, l=dis(id, i);
if(l+r2<=r1) continue;//已知圆内含于求得的圆
if(r1+r2<=l) return false;//已知圆与求得的圆相离
if(l+r1<=r2)//这是针对n=1的,求得的圆内含于已知圆
{
if(r1*r1*2<r2*r2) return false;
}
double x=(l*l+r2*r2-r1*r1)/2/l;
double s = acos((l-x)/r1)*r1*r1+acos(x/r2)*r2*r2-sqrt(r2*r2-x*x)*l;
if(2*s<pi*r2*r2) return false;
}
return true;
}
double solve(double ri, int id)
{
double li=0, mid;
while(fabs(li-ri)>=eps)
{
mid = (li+ri)/2;
if(check(mid, id)) ri=mid;
else li=mid;
}
return mid;
}
int main()
{
scanf("%d", &t);
while(t--)
{
scanf("%d", &n);
for(int i=0; i<n; i++)
{
scanf("%lf%lf%lf", &a[i].x, &a[i].y, &a[i].r);
}
double rr=1e18;
for(int i=0; i<n; i++)
{
double ri=0;
for(int j=0; j<n; j++)
{
ri = max(ri, dis(i, j)+a[j].r);
}
rr = min(rr, solve(ri, i));
}
printf("%.4lf\n", rr);
}
return 0;
}