【zznu-2174】

题目链接

 

题目描述

给出一个圆C1的圆心和半径x1,y1,r1和另外一个圆C2的圆心为x2,y2（x，y均为整数 r为正整数  且绝对值不超1e8）。
已知两圆圆心距不超1e8，给出两者相交面积S（大于0），求出圆C2的最小半径r2。

 

输入

每行输出五个整数x1,y1,r1,x2,y2。以及相交面积S（保留2位有效数字）。

 

输出

结果为正整数，输出占一行。

 

 

样例输入

0 0 1 0 1 1.22
0 0 1 1 1 0.57

 

样例输出

1
1

 

 

 

#include <bits/stdc++.h>
using namespace std;
const double pi = acos(-1.0);
struct Circle
{
    double x, y, r;
};

double dis(Circle a, Circle b)
{
    return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}

double solve(Circle a, Circle b)
{
    double d = dis(a, b);
    if (d >= a.r + b.r) return 0;
    if (d <= fabs(a.r - b.r))
    {
        double r = a.r < b.r ? a.r : b.r;
        return pi * r * b.r;
    }

    double ang1 = acos((a.r*a.r+d*d-b.r*b.r)/2.00/a.r/d);
    double ang2 = acos((b.r*b.r+d*d-a.r*a.r)/2.00/b.r/d);
    double ret = ang1*a.r*a.r+ang2*b.r*b.r-d*a.r*sin(ang1);
    return ret;
}

int main()
{
    Circle a, b;
    double s;
    while(~scanf("%lf%lf%lf%lf%lf%lf", &a.x, &a.y, &a.r, &b.x, &b.y, &s))
    {
        double l = 0, r = 1e8;
        while(fabs(r - l) >= 1e-5)
        {
            double mid = (l+r)/2.0;
            b.r = mid;
            double ss = solve(a, b);
            if(ss > s)
                r = mid;
            else
                l = mid;
        }
        printf("%.0f\n", l);
    }
    return 0;
}