HDU 3264 Open-air shopping malls（计算几何 && 二分优化）

题目链接：【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;
}