HDU 3264 Open-air shopping malls （计算几何-圆相交面积）

传送门：http://acm.hdu.edu.cn/showproblem.php?pid=3264

题意：给你n个圆，坐标和半径，然后要在这n个圆的圆心画一个大圆，大圆与这n个圆相交的面积必须大于等于每个圆面积的一半，问你建在那个圆心半径最小，为多少。

题解：枚举这n个圆，求每个圆的最小半径，通过二分半径来求，然后取这n个的最小值即可，注意点精度就OK了。

AC代码：

 

#include <iostream>

#include <cstdio>

#include <cstring>

#include <string>

#include <cstdlib>

#include <cmath>

#include <vector>

#include <list>

#include <deque>

#include <queue>

#include <iterator>

#include <stack>

#include <map>

#include <set>

#include <algorithm>

#include <cctype>

using namespace std;



#define si1(a) scanf("%d",&a)

#define si2(a,b) scanf("%d%d",&a,&b)

#define sd1(a) scanf("%lf",&a)

#define sd2(a,b) scanf("%lf%lf",&a,&b)

#define ss1(s)  scanf("%s",s)

#define pi1(a)    printf("%d\n",a)

#define pi2(a,b)  printf("%d %d\n",a,b)

#define mset(a,b)   memset(a,b,sizeof(a))

#define forb(i,a,b)   for(int i=a;i<b;i++)

#define ford(i,a,b)   for(int i=a;i<=b;i++)



typedef long long LL;

const int N=33;

const int INF=0x3f3f3f3f;

const double PI=acos(-1.0);

const double eps=1e-8;



int n;



struct xkn

{

    double x,y,r;

    double area;

}p[22],h;



double dis(xkn a,xkn b)

{

    return sqrt( (a.x-b.x)*(a.x-b.x) + (a.y-b.y)*(a.y-b.y) );

}



double fuck(xkn a,xkn b)  //求两圆的相交面积函数。

{

    double d=dis(a,b);

    if(d>=a.r+b.r) return 0;

    double r=(a.r>b.r?b.r:a.r);

    if( d<=fabs(a.r-b.r) ) return PI*r*r;

    double A1=acos( (a.r*a.r+d*d-b.r*b.r)/2/a.r/d );

    double A2=acos( (b.r*b.r+d*d-a.r*a.r)/2/b.r/d );

    double res=A1*a.r*a.r + A2*b.r*b.r;

    res-=sin(A1)*a.r*d;

    return res;

}



bool xiaohao(xkn h)

{

    for(int i=0;i<n;i++)

    {

        double jiao=fuck(p[i],h);

        if(jiao<p[i].area/2)

            return false;

    }

    return true;

}



int main()

{

//    freopen("input.txt","r",stdin);

    int T;

    si1(T);

    while(T--)

    {

        si1(n);

        for(int i=0;i<n;i++)

        {

            sd2(p[i].x,p[i].y);

            sd1(p[i].r);

            p[i].area=PI*p[i].r*p[i].r;

        }

        double Min=55555;

        for(int i=0;i<n;i++)

        {

            double l=0,r=22222,m;

            while((r-l)>eps)

            {

                m=(l+r)/2;

                h=p[i]; h.r=m;

                if(xiaohao(h))

                    r=m;

                else

                    l=m;

            }

            if(Min>m)   Min=m;

        }

        printf("%.4f\n",Min);

    }

    return 0;

}

//9253993	2013-09-30 19:54:36	Accepted	3264	0MS	320K	2260 B	G++	XH_Reventon